c irirtam-test.for c --------------------------------------------------------------- c Test program for the IRI subroutines READIRTAMCOF, FOUT1, c and GAMMA2 that read the IRTAM coefficients and calculate c NmF2, hmF2, and B0 for the Real-Time IRI. c --------------------------------------------------------------- C 2016.01 06/17/16 First release C 2017.01 11/03/17 Added B1 input c --------------------------------------------------------------- c Real ff2(1064),fh2(1064),fb0(1064),fb1(1064) COMMON/CONST/UMR,PI pi=ATAN(1.0)*4. UMR=pi/180. idate=20160523 alati=-11.95 along=283.13 ryear=2016.6 c call igrf_dip(alati,along,ryear,300.0,dec,dip,xmagbr,xmodip) xmodip=-0.37 hourut=0.0 tov=12 itovhm=1200 call READIRTAMCOF(0,idate,itovhm,1064,ff2) call READIRTAMCOF(1,idate,itovhm,1064,fh2) call READIRTAMCOF(2,idate,itovhm,1064,fb0) c call READIRTAMCOF(3,idate,itovhm,1064,fb1) 1122 fof2rt=FOUT1(XMODIP,ALATI,ALONG,HOURUT,TOV,ff2) hmf2rt=FOUT1(XMODIP,ALATI,ALONG,HOURUT,TOV,fh2) B0rt=FOUT1(XMODIP,ALATI,ALONG,HOURUT,TOV,fb0) c B1rt=FOUT1(XMODIP,ALATI,ALONG,HOURUT,TOV,fb1) print*,idate,hourut,itovhm,fof2rt,hmf2rt,B0rt c print*,idate,hourut,itovhm,fof2rt,hmf2rt,B0rt,B1rt hourut=hourut+1.0 if(hourut.lt.18.0) goto 1122 stop end c