#!/bin/bash
#Usage Example:   ./runKt -2.3 

#Fetch the history  and make rain and range counts
# This only needs to be done once  for all Kts.  Then comment out.
#Examples from Yin-Conle paper:
#./rangeCount 3 < ../heulerSeegerStress.his  | grep -v \#2D >hs_StressHist.rg
#./pdrain 3 < heulerSeegerStress.his  > hs_StressHist.rain

# Using non-overload A36  file
./saefcalc2  merged_a36_fitted.html  $1 <hs_StressHist.rain  > hs_A36_Kt=$1.init
./makeRepSaefcalc2 hs_A36_Kt=$1.init

# Using the improvised  A36 Periodic overload file :  
./saefcalc2  merged_a36_w_POL_fitted.html  $1 <hs_StressHist.rain  > hs_A36+POL_Kt=$1.init
./makeRepSaefcalc2 hs_A36+POL_Kt=$1.init



# Using  Heuler Seeger weldmetal
./saefcalc2  heulerSeeger_WeldMetal_fitted.html  $1 <hs_StressHist.rain  > hs_WeldMetal_Kt=$1.init
./makeRepSaefcalc2 hs_WeldMetal_Kt=$1.init

# Using improvised POL  file for Heuler Seeger weld metal
./saefcalc2  heulerSeeger_WeldMetal_POL_fitted.html  $1 <hs_StressHist.rain  > hs_WeldMetal_POL_Kt=$1.init
./makeRepSaefcalc2 hs_WeldMetal_POL_Kt=$1.init


# Added June 3
# Using POL  with non-damaging cycle compensation a la DuQuesnay, Pompetski, Topper

./saefcalc2 merged_a36_w_POL_fitted.html $1 < hs_StressHist-Decay.rain > hs_A36+POL-Decay_kt=$1.ini
./makeRepSaefcalc2 hs_A36+POL-Decay_kt=$1.ini

./saefcalc2 heulerSeeger_WeldMetal_POL_fitted.html $1 < hs_StressHist-Decay.rain > hs_WeldMetal+POL-Decay_kt=$1.ini
./makeRepSaefcalc2  hs_WeldMetal+POL-Decay_kt=$1.ini