#HT100  Sigy=885  Siga=400
#SHY685 Sigy=685  Siga=330
#SM570  Sigy=450  Siga=240
#SM490  Sigy=315  Siga=175
#SM400  Sigy=235  Siga=130
# ./f90_ams_d  DD  t0  Sigy  Siga
# write(6,'(5f10.3)') DD,t0,Sigy,Siga,Pk

gfortran -o f90_ams_d f90_ams_d.f90
./f90_ams_d 4600 21 315 175
./f90_ams_d 4600 22 315 175
./f90_ams_d 4600 23 315 175
./f90_ams_d 4600 24 315 175
./f90_ams_d 4600 25 315 175
./f90_ams_d 4600 26 315 175
./f90_ams_d 4600 27 315 175
./f90_ams_d 4600 28 315 175
./f90_ams_d 4600 29 315 175
./f90_ams_d 4600 30 315 175
./f90_ams_d 4600 31 315 175
./f90_ams_d 4600 32 315 175
./f90_ams_d 4600 33 315 175
./f90_ams_d 4600 34 315 175
./f90_ams_d 4600 35 315 175
./f90_ams_d 4600 36 315 175
./f90_ams_d 4600 37 315 175
./f90_ams_d 4600 38 315 175
./f90_ams_d 4600 39 315 175
./f90_ams_d 4600 40 315 175