Mon Feb 28 10:43:34 US/Mountain 2000 Diffsrc executing from /home/foster/tiegcm/tgcm13mt/modsrc.snoe Source directory = /home/tgcm/tgcm13mt ======================================================================== Diff of /home/tgcm/tgcm13mt/cons.h and cons.h: 7c7 < | difk,tbound,dift,decayf,decayh,xmue --- > | difk,tbound,dift,decayf,decayh,xmue,alfalp,efluxlp,sfeps 15c15 < 6TBOUND,DIFT(ZKMXP),XMUE(ZKMXP), --- > 6ALFALP,EFLUXLP,TBOUND,DIFT(ZKMXP),XMUE(ZKMXP), 17c17 < 8ZBA(ZJMX),TBA(ZJMX),UBA(ZJMX),VBA(ZJMX),BNDA(ZIMXP) --- > 8ZBA(ZJMX),TBA(ZJMX),UBA(ZJMX),VBA(ZJMX),BNDA(ZIMXP),SFEPS ======================================================================== Diff of /home/tgcm/tgcm13mt/crates.h and crates.h: 36c36,39 < !DIR$ TASKCOMMON rates_priv,metas,fsaray --- > real :: tvib,disn2p > common/rates_priv1/ tvib(zimxp,zkmxp),disn2p(zimxp,zkmxp) > ! > !DIR$ TASKCOMMON rates_priv,metas,fsaray,rates_priv1 ======================================================================== WARNING: cannot find file /home/tgcm/tgcm13mt/othend.h ======================================================================== Diff of /home/tgcm/tgcm13mt/advnce.f and advnce.f: 43c43 < if (iadvda == 1) call advnceday(c) --- > if (iadvda == 1) call advnceday 53a54 > C call taila (sec) 122c123 < !MIC$+ PRIVATE(j,jx,jj,i,racs,ifrst) --- > !MIC$+ PRIVATE(j,jx,jj,i,racs) 146d146 < if(jx.eq.jmax) ifrst=1 168a169 > ifrst=1 204c205 < subroutine advnceday(c) --- > subroutine advnceday 205a207,208 > include "params.h" > include "cons.h" 208,209d210 < ! Args: < real,intent(inout) :: c(120) 213c214 < real :: pi,delta --- > real :: pi,delta,thet0 234a236,246 > C **** SFEPS IS 6% VARIATION IN SOLAR OUTPUT OVER A YEAR > C **** CAUSED BY THE ORBITAL ECCENTRICITY > ! SFEPS is in common in cons.h. It is initialized to 1 in con.f. > ! If calendar day is being advanced, SFEPS is initially set in > ! input.f, then updated here in advnce.f at day boundaries. It is > ! referenced by qrj and xray. > ! > PI = 3.14159265358979 > THET0 = 2.*PI*FLOAT(IIDAY)/365. > SFEPS = 1.000110+0.034221*COS(THET0)+0.001280*SIN(THET0) > 1 +0.000719*COS(2.*THET0)+0.000077*SIN(2.*THET0) 238d249 < PI = 3.14159265358979 239a251,255 > C **** CCM3 CALCULATION OF DELTA > C DELTA = 0.006918-0.399912*COS(THET0)+0.070257*SIN(THET0) > C 1 -0.006758*COS(2.*THET0)+0.000907*SIN(2.*THET0) > C 2 -0.002697*COS(3.*THET0)+0.001480*SIN(3.*THET0) > C **** 243a260 > write(6,"(' Recalculated SFEPS = ',e12.4)") sfeps ======================================================================== WARNING: cannot find file /home/tgcm/tgcm13mt/altv.f ======================================================================== Diff of /home/tgcm/tgcm13mt/aurht.f and aurht.f: 7a8,10 > > ! include "phys.h" ! for debug only > 17a21,24 > c > c Low energy auroral proton input (see ALFALP and EFLUXLP from input): > real :: alfa_lp,eflux_lp,qteaur > common/aurlp/ alfa_lp(zimxp),eflux_lp(zimxp),qteaur(zimxp) 24c31 < real :: pi --- > real :: pi,alfap1,alfap2,pe1,pe2,alfap0,e0p 86a94,96 > > ! write(6,"('aurht: j=',i3,' e0=',e12.4)") j,e0 > 95c105,118 < --- > c > c Low-energy proton inputs: > ALFAP1 = 5. > ALFAP2 = 15. > C PE1 = 0.15 > C PE2 = 0.4 > PE1 = 1.E-20 > PE2 = 1.E-20 > ALFAP0 = 0.5*(ALFAP1+ALFAP2) > alfa_lp(i+2) = ALFAP0*(1-RALFA2*WK1(I)) > E0P = 0.5*(PE1+PE2) > eflux_lp(i+2) = E0P*(1.-REE*WK1(I))*exp(-(wk3(i)/wk2(i))**2) > c > qteaur(i+2) = -7.e+8*exp(-(wk3(i)/wk2(i))**2) 96a120,133 > c > c Periodic points for alfa_lp and eflux_lp: > alfa_lp(1) = alfa_lp(im+1) > alfa_lp(2) = alfa_lp(im+2) > alfa_lp(im+3) = alfa_lp(3) > alfa_lp(im+4) = alfa_lp(4) > eflux_lp(1) = eflux_lp(im+1) > eflux_lp(2) = eflux_lp(im+2) > eflux_lp(im+3) = eflux_lp(3) > eflux_lp(im+4) = eflux_lp(4) > qteaur(1) = qteaur(im+1) > qteaur(2) = qteaur(im+2) > qteaur(im+3) = qteaur(3) > qteaur(im+4) = qteaur(4) ======================================================================== Diff of /home/tgcm/tgcm13mt/cmpn2d.f and cmpn2d.f: 46,47c46,47 < 1(RA1(I,1)*F(I,NNOPK)*0.8+RA3(I,1)*F(I,NN2PK)*0.9)*SQRT(F(I,NEK)* < 2F(I,NEK+1)) --- > 1(RA1(I,1)*F(I,NNOPK)*0.85+RA3(I,1)*F(I,NN2PK)*0.9)*SQRT(F(I,NEK) > 2*F(I,NEK+1)) ======================================================================== Diff of /home/tgcm/tgcm13mt/cmpn4s.f and cmpn4s.f: 80c80 < 6*F(I,NNOPK)*0.2+RA3(I,1)*F(I,NN2PK)*1.1) --- > 6*F(I,NNOPK)*0.15+RA3(I,1)*F(I,NN2PK)*1.1) ======================================================================== Diff of /home/tgcm/tgcm13mt/cmpo2o.f and cmpo2o.f: 78c78 < S4(I,1) = RK3(I,1)*F(NN2PK)+RK8(I,1)*F(I,NNPK) --- > S4(I,1) = RK3(I,1)*F(I,NN2PK)+RK8(I,1)*F(I,NNPK) ======================================================================== Diff of /home/tgcm/tgcm13mt/comp.f and comp.f: 24c24 < DATA SMALL/1.E-6/ --- > DATA SMALL/1.E-20/ ======================================================================== Diff of /home/tgcm/tgcm13mt/con.f and con.f: 297a298,304 > ! > ! Intitialize SFEPS: Default solar flux variation due to orbital > ! eccentricity assumes no calendar day advance, thus no variation. > ! If calendar day *is* being advanced, SFEPS will be calculated > ! for each new day in advnce.f. SFEPS is in common in cons.h. > ! > SFEPS = 1. ======================================================================== Diff of /home/tgcm/tgcm13mt/dfact.f and dfact.f: /home/tgcm/tgcm13mt/dfact.f and dfact.f are identical ======================================================================== Diff of /home/tgcm/tgcm13mt/dt.f and dt.f: 372c372,374 < S5(I,1) = S5(I,1)+S4(I,1)*3. --- > S5(I,1) = S5(I,1)+S4(I,1)*1.5 > C S5(I,1) = S5(I,1)+S4(I,1)*2. > C S5(I,1) = S5(I,1)+S4(I,1)*3. 655a658,664 > C **** > C **** SET MINIMUM VALUE OF TN AT 100K > C **** > NTK = NJNP+NT > DO I=1,LEN3 > F(I,NTK) = merge(F(I,NTK),100.,F(I,NTK)-100.>=0.) > ENDDO ======================================================================== Diff of /home/tgcm/tgcm13mt/dynamics.f and dynamics.f: 46a47 > C call altv ======================================================================== Diff of /home/tgcm/tgcm13mt/elden.f and elden.f: 145c145,148 < S5(I,1) = merge(S5(I,1),1.E+2,S5(I,1)-1.E+2>=0.) --- > C S5(I,1) = merge(S5(I,1),3.3E+3,S5(I,1)-3.3E+3>=0.) > S5(I,1) = merge(S5(I,1),3.1E+3,S5(I,1)-3.1E+3>=0.) > C S5(I,1) = merge(S5(I,1),2.8E+3,S5(I,1)-2.8E+3>=0.) > C S5(I,1) = merge(S5(I,1),2.5E+3,S5(I,1)-2.5E+3>=0.) 163a167,168 > C F(I,NEK)= 1.5*S5(I,1)-0.5*S5(I,2) > C F(I,NEK+KMAX)=1.5*S5(I,KMAX)-0.5*S5(I,KMAXM1) ======================================================================== WARNING: cannot find file /home/tgcm/tgcm13mt/euvac.f ======================================================================== Diff of /home/tgcm/tgcm13mt/getgpi.f and getgpi.f: 260,267c260,267 < if (ifrst.eq.0) then < write(6,"(' ')") < write(6,"(' iter iyd hr:min:sec f107 f107a', < + ' Kp ctpoten hpower')") < write(6,"(72('-'))") < endif < write(6,"(i6,i7,i5,':',i2,':',i2,2x,5f8.3)") < + iter,iyd_cur,ih,im,is,f107,f107a,fkp,ctpoten,hpower --- > ! if (ifrst.eq.0) then > ! write(6,"(' ')") > ! write(6,"(' iter iyd hr:min:sec f107 f107a', > ! + ' Kp ctpoten hpower')") > ! write(6,"(72('-'))") > ! endif > ! write(6,"(i6,i7,i5,':',i2,':',i2,2x,5f8.3)") > ! + iter,iyd_cur,ih,im,is,f107,f107a,fkp,ctpoten,hpower ======================================================================== Diff of /home/tgcm/tgcm13mt/heelis.f and heelis.f: /home/tgcm/tgcm13mt/heelis.f and heelis.f are identical ======================================================================== Diff of /home/tgcm/tgcm13mt/input.f and input.f: 290a291,299 > c > c 6/4/98 btf: Low energy proton inputs (single floats, in cons.h): > c ALFALP Average alpha energy of low proton aurora (Kev) > c (2 < ALFALP < 40) (ALFALP=2 for Maxwellian) > c (default = 10.) > c EFLUXLP Incident energy flux for protons (erg/cm-2/s-1) > c (default = 0.4) > c (see sub low_proton in proton.f, called from orora, and > c common/aurlp/ in aurht and orora) 359c368 < 8, ALFA1(4,NLEX),ALFA2(4,NLEX) --- > 8, ALFA1(4,NLEX),ALFA2(4,NLEX),THET0,PI 1086a1096,1101 > C **** SFEPS IS 6% VARIATION IN SOLAR OUTPUT OVER A YEAR > C **** CAUSED BY THE ORBITAL ECCENTRICITY > PI = 3.14159265358979 > THET0 = 2.*PI*FLOAT(yrda(2))/365. > SFEPS = 1.000110+0.034221*COS(THET0)+0.001280*SIN(THET0) > 1 +0.000719*COS(2.*THET0)+0.000077*SIN(2.*THET0) 1344a1360,1363 > c > c Low energy proton: > CALL LEXCON(ALFALP ,8HALFALP ,ICNT(45)) > CALL LEXCON(EFLUXLP ,8HEFLUXLP ,ICNT(46)) 1574a1594,1602 > c > c Low energy protons in aurora: > c > c if (icnt(45) == 0) alfalp = 10. > c if (icnt(46) == 0) efluxlp = 0.4 > c alfalp = 10. > c efluxlp = 0.4 > alfalp = 10. > efluxlp = 1.e-20 ======================================================================== Diff of /home/tgcm/tgcm13mt/lamdas.f and lamdas.f: 42c42,43 < DATA DIPMIN/.20/ --- > DATA DIPMIN/.005/ > C DATA DIPMIN/0.17/ 53c54,55 < T4(I)=SIN(DIPMAG(I,J)) --- > T4(I)=SIN(merge(DIPMAG(I,J),DIPMIN,DIPMAG(I,J)-DIPMIN>=0.)) > C T4(I)=SIN(DIPMAG(I,J)) 67a70 > V3(I,1,1)=merge(V3(I,1,1),1.E-20,V3(I,1,1)-1.E-20>=0.) 135,138c138,147 < V4(I,1,1)=V6(I,1,1)*(2.33E-11*(1.-F(I,NPS1K)-F(I,NPS2K))/RMASS(3)* < 1F(I,NTEK)*(1.-1.21E-4*F(I,NTEK))+1.82E-10*F(I,NPS1K)/RMASS(1)* < 2V5(I,1,1)*(1.+3.6E-2*V5(I,1,1))+8.9E-11*F(I,NPS2K)/RMASS(2)* < 3V5(I,1,1)*(1.+5.7E-4*F(I,NTEK))) --- > ! V4(I,1,1)=V6(I,1,1)*(2.33E-11*(1.-F(I,NPS1K)-F(I,NPS2K))/RMASS(3)* > ! 1F(I,NTEK)*(1.-1.21E-4*F(I,NTEK))+1.82E-10*F(I,NPS1K)/RMASS(1)* > ! 2V5(I,1,1)*(1.+3.6E-2*V5(I,1,1))+8.9E-11*F(I,NPS2K)/RMASS(2)* > ! 3V5(I,1,1)*(1.+5.7E-4*F(I,NTEK))) > V4(I,1,1)=V6(I,1,1)*(2.33E-11* > 1 merge((1.-F(I,NPS1K)-F(I,NPS2K))/RMASS(3),1.E-20, > 2 (1.-F(I,NPS1K)-F(I,NPS2K))/RMASS(3)-1.E-20>=0.)* > 3F(I,NTEK)*(1.-1.21E-4*F(I,NTEK))+1.82E-10*F(I,NPS1K)/RMASS(1)* > 4V5(I,1,1)*(1.+3.6E-2*V5(I,1,1))+8.9E-11*F(I,NPS2K)/RMASS(2)* > 5V5(I,1,1)*(1.+5.7E-4*F(I,NTEK))) 199a209,219 > C **** > C **** > C **** LEVELS 1 AND KMAXP1 > C F(I,NLXXK)= 1.5*V3(I,1,1)-0.5*V3(I,2,1) > C F(I,NLYYK)= 1.5*V4(I,1,1)-0.5*V4(I,2,1) > C F(I,NLXYK)= 1.5*V5(I,1,1)-0.5*V5(I,2,1) > C F(I,NLXXK+KMAX)=1.5*V3(I,KMAX,1)-0.5*V3(I,KMAXM1,1) > C F(I,NLYYK+KMAX)=1.5*V4(I,KMAX,1)-0.5*V4(I,KMAXM1,1) > C F(I,NLXYK+KMAX)=1.5*V5(I,KMAX,1)-0.5*V5(I,KMAXM1,1) > C **** > C **** ======================================================================== Diff of /home/tgcm/tgcm13mt/magdyn.f and magdyn.f: 26a27 > C SIN10=0.17 ======================================================================== Diff of /home/tgcm/tgcm13mt/new.f and new.f: 38a39,40 > C RNU(I,1,2,2)=3.0E-11*SQRT(TR(I,1))*(1.-0.135*ALOG10(TR(I,1)* > C 1 0.001))**2 ======================================================================== WARNING: cannot find file /home/tgcm/tgcm13mt/nextlu.f ======================================================================== Diff of /home/tgcm/tgcm13mt/opflux.f and opflux.f: 24c24,31 < DATA PHID,PHIN/1.5E+8,-1.5E+8/ --- > C DATA PHID,PHIN/1.5E+8,-1.5E+8/ > C DATA PHID,PHIN/0.,-3.0E+8/ > C DATA PHID,PHIN/0.,-1.5E+8/ > DATA PHID,PHIN/-1.5E+8,-3.0E+8/ > C DATA PHID,PHIN/-2.5E+8,-3.0E+8/ > C DATA PHID,PHIN/-5.0E+8,-6.0E+8/ > C DATA PHID,PHIN/-7.0E+8,-3.0E+8/ > C DATA PHID,PHIN/-9.0E+8,-3.0E+8/ 28c35,37 < DATA PPOLAR/+1.E+8/ --- > C DATA PPOLAR/+1.E+8/ > C DATA PPOLAR/+5.E+7/ > DATA PPOLAR/0./ ======================================================================== Diff of /home/tgcm/tgcm13mt/oplus.f and oplus.f: /home/tgcm/tgcm13mt/oplus.f and oplus.f are identical ======================================================================== Diff of /home/tgcm/tgcm13mt/orora.f and orora.f: 17a18 > include "crates.h" 23a25,28 > c > c Low energy auroral proton input (see ALFALP and EFLUXLP from input): > real :: alfa_lp,eflux_lp,qteaur > common/aurlp/ alfa_lp(zimxp),eflux_lp(zimxp),qteaur(zimxp) 34a40,43 > ! Ion production from low energy proton source > ! (see call low_proton below. Sub low_proton is in proton.f) > real :: qia(zimxp,zkmxp,5) > ! 55c64,73 < IF (ISHUNK.EQ.0) return --- > IF(ISHUNK.EQ.0) then > qia = 1.e-20 > qteaur = 1.e-20 > ! call addfsech('QIA1_N2+',qia(:,:,1),zimxp,zkmxp,zkmx,j) > ! call addfsech('QIA2_O2+',qia(:,:,2),zimxp,zkmxp,zkmx,j) > ! call addfsech('QIA3_O+ ',qia(:,:,3),zimxp,zkmxp,zkmx,j) > ! call addfsech('QIA4_NO+',qia(:,:,4),zimxp,zkmxp,zkmx,j) > ! call addfsech('QIA5_N+ ',qia(:,:,5),zimxp,zkmxp,zkmx,j) > RETURN > endif 79,81d96 < ! 9/99: aion is a cpu consumer because of **real < ! and exphf. It is located in inline.f. < ! 85a101,116 > ! > ! Contribution of low energy protons: > ! qia(zimxp,zkmxp,5) is output ion production for N2+, O2+, O+, NO+, N+ > ! (top level zkmxp not defined) > ! > ! call low_proton(alfa_lp,eflux_lp,qia,j) > qia = 0. > ! > ! Put qia on secondary histories: > ! > ! call addfsech('QIA1_N2+',qia(:,:,1),zimxp,zkmxp,zkmx,j) > ! call addfsech('QIA2_O2+',qia(:,:,2),zimxp,zkmxp,zkmx,j) > ! call addfsech('QIA3_O+ ',qia(:,:,3),zimxp,zkmxp,zkmx,j) > ! call addfsech('QIA4_NO+',qia(:,:,4),zimxp,zkmxp,zkmx,j) > ! call addfsech('QIA5_N+ ',qia(:,:,5),zimxp,zkmxp,zkmx,j) > ! 121a153,158 > c > c 6/98: Include production due to low energy protons. > c qia(zimxp,zkmxp,5) is output ion production from low_proton > c (called above) for N2+, O2+, O+, NO+, N+ respectively > c (currently, NO+ is not defined) > c 125,129c162,167 < S13(I,1)=S1(I,1)*F(I,NPS1K)/(RMASS(1)*S3(I,1)) < S14(I,1)=S1(I,1)*(0.5*F(I,NPS1K)/RMASS(1)+ < 10.56*F(I,NPS2K)/RMASS(2))/S3(I,1) < S15(I,1)=S1(I,1)*0.7*merge(1.-F(I,NPS1K)-F(I,NPS2K),0.,1.- < 1F(I,NPS1K)-F(I,NPS2K)>=0.)/(RMASS(3)*S3(I,1)) --- > S13(I,1)=S1(I,1)*F(I,NPS1K)/(RMASS(1)*S3(I,1))+qia(i,1,2) > S14(I,1)=S1(I,1)*(0.5*F(I,NPS1K)/RMASS(1)+ > 1 0.56*F(I,NPS2K)/RMASS(2))/S3(I,1)+qia(i,1,3) > S15(I,1)=S1(I,1)*0.7*merge(1.-F(I,NPS1K)-F(I,NPS2K),1.E-8, > 1 (1.-F(I,NPS1K)-F(I,NPS2K))-1.E-8>0.)/(RMASS(3)*S3(I,1)) > 2 +qia(i,1,1) 157a196,204 > DO 41 K=2,KMAX > DO 41 I=1,LEN1 > DISN2P(I,K)=DISN2P(I,K)+SQRT(S15(I,K-1)*S15(I,K)) > 41 CONTINUE > DO 42 I=1,LEN1 > DISN2P(I,1) = DISN2P(I,1)+1.5*S15(I,1)-0.5*S15(I,2) > DISN2P(I,KMAXP1) = DISN2P(I,KMAXP1)+1.5*S15(I,KMAX)-0.5* > 1 S15(I,KMAX-1) > 42 CONTINUE ======================================================================== WARNING: cannot find file /home/tgcm/tgcm13mt/othend.f ======================================================================== WARNING: cannot find file /home/tgcm/tgcm13mt/proton.f ======================================================================== Diff of /home/tgcm/tgcm13mt/qinite.f and qinite.f: 29a30,36 > C 5AL/5.E8,5.E7,5.E7/ > C 5AL/5.E7,5.E6,5.E6/ > 5AL/2.0E7,2.0E6,2.0E6/ > C 5AL/1.0E7,1.0E6,1.0E6/ > C 5AL/1.5E7,1.5E6,1.5E6/ > C 5AL/8.0E6,8.0E5,8.0E5/ > C 5AL/5.5E6,5.5E5,5.5E5/ 31c38,40 < 5AL/5.E4,5.E3,5.E3/ --- > C 5AL/5.E5,5.E4,5.E4/ > C 5AL/1.E5,1.E4,1.E4/ > C 5AL/5.E4,5.E3,5.E3/ 122,125c131,132 < F(I,NQO2PK)=F(I,NQO2PK)+0.67*S6(I,1)*S5(I,1) < F(I,NQOPK)=F(I,NQOPK)+(0.33*S6(I,1)+S7(I,1))*S5(I,1) < F(I,NQN2PK)=F(I,NQN2PK)+0.86*S8(I,1)*S5(I,1) < F(I,NQNPK)=F(I,NQNPK)+0.14*S8(I,1)*S5(I,1) --- > C F(I,NQO2PK)=F(I,NQO2PK)+0.67*S6(I,1)*S5(I,1) > F(I,NQO2PK)=F(I,NQO2PK)+0.67*S6(I,1)*S5(I,1)+1.E+2 ======================================================================== Diff of /home/tgcm/tgcm13mt/qrj.f and qrj.f: 213c213 < 1 S3(I,1)/RMASS(3) --- > 1 S3(I,1)/RMASS(3)*2. 300a301 > DISN2P(I,1) = S10(I,1)*S14(I,1)*S5(I,1) 376c377 < + F(I,NNO2K)-1.E18>=0.0)*(1.+0.11*(C(61)-65.)/165.) --- > + F(I,NNO2K)-1.E18>=0.0)*(1.+0.11*(C(61)-65.)/165.)*SFEPS 383a385 > implicit none 390a393,394 > ! > ! Local: 393a398 > real :: EUVFLX(37),wave1(lmax),wave2(lmax) ! output from euvac 396,397c401,402 < real :: flya < integer :: iscale --- > real :: flya,f107,f107a,hlybr,fexvir,hlya,heiew,xuvfac > integer :: iscale,n,nn 402,407c407,434 < C **** < C **** THE FLAG ISCALE CHOOSES THE FLUX MODELING METHOD < C **** FROM THE SUBROUTINE SSFLUX. ISCALE=0 DOES SOLOMON'S, < C **** ISCALE=1 DOES A LIN. INTERP.,AND ISCALE=2 DOES TOBISKA' < C **** < CALL SSFLUX(ISCALE,F107,F107A,FLYA,SFLUX) --- > HLYBR = 0. > FEXVIR = 0. > HLYA = 3.E+11+0.4E+10*(C(61)-70.) > HEIEW = 0. > C XUVFAC =0. > C XUVFAC = 2.0 - (C(61)-68.0) / (243.0-68.0) > XUVFAC = 4.0 - (C(61)-68.0) / (243.0-68.0) > IF (XUVFAC .LT. 1.0) XUVFAC = 1.0 > C > C ISCALE =0 for Hinteregger contrast ratio method > C =1 for Hinteregger linear interpolation > C =2 for Tobiska EUV91 model > C =3 for Woods & Rottman 10 Nov. 1988 measurement > C =4 for Woods & Rottman 20 Jun. 1989 measurement > C F107 daily 10.7 cm flux (1.E-22 W m-2 Hz-1) > C F107A 81-day centered average 10.7 cm flux > C HLYBR ratio of H Ly-b 1026A flux to solar minimum value (optional) > C FEXVIR ratio of Fe XVI 335A flux to solar minimum value (optional) > C HLYA H Lyman-alpha flux (photons cm-2 s-1) (optional) > C HEIEW He I 10830A equivalent width, (milliAngstroms) (optional) > C XUVFAC factor for scaling flux 16-250A (optional) > C WAVE1 longwave bound of spectral intervals (Angstroms) > C WAVE2 shortwave bound of intervals (= WAVE1 for indiv. lines) > C SFLUX scaled solar flux returned by subroutine (photons cm-2 s-1) > C > CALL SSFLUX(ISCALE,F107,F107A,HLYBR,FEXVIR,HLYA, > | HEIEW, XUVFAC, WAVE1, WAVE2, SFLUX) > CALL EUVAC(F107,F107A,EUVFLX) 412c439 < FEUV(N) = SFLUX(N) --- > FEUV(N) = SFLUX(N)*SFEPS 415c442 < FSRC(N) = SFLUX(N) --- > FSRC(N) = SFLUX(N)*SFEPS 416a444,447 > DO 55 N=1,37 > NN = N+15 > FEUV(NN) = EUVFLX(N)*SFEPS > 55 CONTINUE 551,552c582,583 < C DATA EUVEFF/29*0.05/ < DATA EUVEFF/29*0.10/ --- > DATA EUVEFF/29*0.05/ > C DATA EUVEFF/29*0.10/ 556a588 > implicit none 565c597 < parameter (lmax=59) --- > integer,parameter :: lmax=59 570a603,610 > real :: wleuv1,wleuv2,sigao,sigao2,sigan2,sigio,sigio2, > | sigin2,brop4s,brop2d,brop2p,sigop2p,sigop2d,sigop4s, > | sigin,brn2np,bro2op > COMMON/EUV/ WLEUV1(37),WLEUV2(37),SIGAO(37),SIGAO2(37), > 1 SIGAN2(37),SIGIO(37),SIGIO2(37),SIGIN2(37), > 2 BROP4S(37),BROP2D(37),BROP2P(37), > 3 SIGOP2P(37),SIGOP2D(37),SIGOP4S(37), > 4 SIGIN(37),BRN2NP(37),BRO2OP(37) 571a612,614 > ! Local: > integer :: m,n,nn > ! 585a629,651 > ! > ! DO 51 loop moved here from original euvac. > DO 51 N=1,37 > BROP2P(N) = 0. > IF(N.GT.14) BROP2P(N) = 1.-BROP2D(N)-BROP4S(N) > SIGOP2P(N)=SIGIO(N)*BROP2P(N) > SIGOP2D(N)=SIGIO(N)*BROP2D(N) > SIGOP4S(N)=SIGIO(N)*BROP4S(N) > 51 CONTINUE > DO 56 N = 1,37 > NN = N+15 ! 16:52 > SIGEUV(1,NN) = SIGAO2(N) > SIGEUV(2,NN) = SIGAO(N) > SIGEUV(3,NN) = SIGAN2(N) > SIGMAS(1,NN) = SIGIO2(N) > SIGMAS(2,NN) = SIGOP4S(N) > SIGMAS(3,NN) = SIGIN2(N) > SIGMAS(4,NN) = SIGIN(N) > SIGMAS(5,NN) = SIGOP2D(N) > SIGMAS(6,NN) = SIGOP2P(N) > BRN2(NN) = BRN2NP(N) > BRO2(NN) = BRO2OP(N) > 56 CONTINUE ======================================================================== Diff of /home/tgcm/tgcm13mt/qtieff.f and qtieff.f: 11a12 > include "crates.h" 21a23 > C F(I,NQTEFK) = F(I,NQTEFK)+1.57*DISN2P(I,1) ======================================================================== Diff of /home/tgcm/tgcm13mt/rates.f and rates.f: 38,39c38,39 < RK1(I,1)=(((9.65E-16*S15(I,1)-5.17E-14)*S15(I,1)+1.073E-12)* < 1S15(I,1)-7.74E-12)*S15(I,1)+2.82E-11 --- > ! RK1(I,1)=(((9.65E-16*S15(I,1)-5.17E-14)*S15(I,1)+1.073E-12)* > ! 1S15(I,1)-7.74E-12)*S15(I,1)+2.82E-11 40a41,48 > C **** NEW RATE COEFFICIENT FOR O+ + O2 -> O + O2+ > C **** > C RK1(I,1)=1.7E-11*(300./F(I,NTK))**0.77+8.54E-11*exphf(-3467./ > C 1 F(I,NTK)) > C **** NEW RATE COEFFICIENT FOR O+ + O2 -> O + O2+ (HIERL) > C **** > RK1(I,1)=1.6E-11*S15(I,1)**(-0.52)+5.5E-11*exphf(-22.85/S15(I,1)) > C **** 44a53,61 > C **** NEW RATE COEFFICIENT FOR O+ + N2 -> N + NO+ > C **** > C **** > TVIB(I,1) = F(I,NTK) > S5(I,1) = exphf(-3353./TVIB(I,1)) > RK2(I,1) = (((((270.*S5(I,1)+220.)*S5(I,1)+85.)*S5(I,1)+38.) > 1 *S5(I,1)+1.)*RK2(I,1)*S5(I,1)+RK2(I,1)) > 2 *(1.-S5(I,1)) > C **** 70c87 < | exphf(-1.E-8*F(I,NNO2K)**0.38) --- > | exphf(-1.E-8*F(I,NNO2K)**0.38)*SFEPS 72c89 < | exphf(-8.E-21*F(I,NNO2K)) --- > | exphf(-8.E-21*F(I,NNO2K))*SFEPS 74c91 < | exphf(-8.E-21*F(I,NNVO2K)) --- > | exphf(-8.E-21*F(I,NNVO2K))*SFEPS 120c137,138 < BETA4(I,1) = 6.9E-13 --- > C BETA4(I,1) = 6.9E-13 > BETA4(I,1) = 5.0E-13 ======================================================================== Diff of /home/tgcm/tgcm13mt/runmds.f and runmds.f: 32c32,33 < DATA AVTO/3.0E-12/ --- > ! DATA AVTO/3.0E-12/ > DATA AVTO/4.0E-12/ 38c39,40 < DATA HOR/36*.5/ --- > C DATA HOR/36*.5/ > DATA HOR/36*.25/ ======================================================================== Diff of /home/tgcm/tgcm13mt/settei.f and settei.f: 20a21,22 > real :: alfa_lp,eflux_lp,qteaur > common/aurlp/ alfa_lp(zimxp),eflux_lp(zimxp),qteaur(zimxp) 30c32 < DATA FPOLAR/-1.0E+9/ --- > DATA FPOLAR/-3.0E+9/ 87,88c89,94 < T4(I) =( -5.0E+7*F107TE*T3(I)-4.0E+7*F107TE)/4. < T5(I) = T4(I)/3. --- > C T4(I) =( -5.0E+7*F107TE*T3(I)-4.0E+7*F107TE)/2.5 > T4(I) =( -5.0E+7*F107TE*T3(I)-4.0E+7*F107TE)*1.2 > C T5(I) = T4(I)/3. > T5(I) = T4(I)/2. > T4(I) = T4(I)+QTEAUR(I) > T5(I) = T5(I)+QTEAUR(I) 199c205,206 < T2(I) = merge(T2(I),.01,T2(I)-.01>=0.) --- > C T2(I) = merge(T2(I),.01,T2(I)-.01>=0.) > T2(I) = merge(T2(I),.10,T2(I)-.10>=0.) ======================================================================== Diff of /home/tgcm/tgcm13mt/ssflux.f and ssflux.f: 1c1,3 < SUBROUTINE SSFLUX(ISCALE,F107,F107A,FLYA,SFLUX) --- > ! > SUBROUTINE SSFLUX (ISCALE, F107, F107A, HLYBR, FEXVIR, HLYA, > > HEIEW, XUVFAC, WAVE1, WAVE2, SFLUX) 3,4c5 < C < integer,parameter :: LM=59 --- > integer,parameter :: lmax=59 8,10c9,10 < real,intent(in) :: f107,f107a < real,intent(out) :: sflux(lm) < real,intent(inout) :: flya --- > real,intent(in) :: f107,f107a,hlybr,fexvir,hlya,heiew,xuvfac > real,intent(out) :: wave1(lmax),wave2(lmax),sflux(lmax) 13,19c13,21 < integer :: l,j < real :: WAVE1(LM), WAVE2(LM), frat,f107md,sfnorm, < + WAVEL(LM), WAVES(LM), RFLUX(LM), XFLUX(LM), < + SCALE1(LM), SCALE2(LM), SRA(8), SRB(8), B1(3), B2(3), < + CHINT(LM),CHSLP(LM),COINT(LM),COSLP(LM),W1(LM),W2(LM), < + WCHR(LM),FLUXCH(LM),FLUXCO(LM),WAVEAV(59),r1,r2 < --- > real :: > > WAVEL(LMAX), WAVES(LMAX), RFLUX(LMAX), XFLUX(LMAX), > > SCALE1(LMAX), SCALE2(LMAX), > > TCHR0(LMAX), TCHR1(LMAX), TCHR2(LMAX), > > TCOR0(LMAX), TCOR1(LMAX), TCOR2(LMAX), > > WAR1(LMAX), WAR2(LMAX), > > B1(3), B2(3) > real :: frat,r1,r2,hlymod,heimod,xuvf > integer :: l 20a23 > C regression coefficients which reduce to solar min. spectrum: 22a26,29 > C 'best fit' regression coefficients, commented out, for reference: > C DATA B1/1.31, 0.01106, 0.00492/, B2/-6.618, 0.66159, 0.38319/ > C > C 24,32c31,39 < + 1450.00, 1400.00, 1350.00, 1300.00, 1250.00, 1215.67, < + 1200.00, 1150.00, 1100.00, 1050.00, 1031.91, 1025.72, < + 1000.00, 977.02, 950.00, 900.00, 850.00, 800.00, < + 789.36, 770.41, 765.15, 750.00, 703.31, 700.00, < + 650.00, 629.73, 609.76, 600.00, 584.33, 554.37, < + 550.00, 500.00, 465.22, 450.00, 400.00, 368.07, < + 350.00, 303.78, 303.31, 300.00, 284.15, 256.30, < + 250.00, 200.00, 150.00, 100.00, 50.00, 32.00, < + 23.00, 16.00, 8.00, 4.00, 2.00/ --- > > 1450.00, 1400.00, 1350.00, 1300.00, 1250.00, 1215.67, > > 1200.00, 1150.00, 1100.00, 1050.00, 1031.91, 1025.72, > > 1000.00, 977.02, 950.00, 900.00, 850.00, 800.00, > > 789.36, 770.41, 765.15, 750.00, 703.31, 700.00, > > 650.00, 629.73, 609.76, 600.00, 584.33, 554.37, > > 550.00, 500.00, 465.22, 450.00, 400.00, 368.07, > > 350.00, 303.78, 303.31, 300.00, 284.15, 256.30, > > 250.00, 200.00, 150.00, 100.00, 50.00, 32.00, > > 23.00, 16.00, 8.00, 4.00, 2.00/ 34,72c41,79 < + 1400.00, 1350.00, 1300.00, 1250.00, 1200.00, 1215.67, < + 1150.00, 1100.00, 1050.00, 1000.00, 1031.91, 1025.72, < + 950.00, 977.02, 900.00, 850.00, 800.00, 750.00, < + 789.36, 770.41, 765.15, 700.00, 703.31, 650.00, < + 600.00, 629.73, 609.76, 550.00, 584.33, 554.37, < + 500.00, 450.00, 465.22, 400.00, 350.00, 368.07, < + 300.00, 303.78, 303.31, 250.00, 284.15, 256.30, < + 200.00, 150.00, 100.00, 50.00, 32.00, 23.00, < + 16.00, 8.00, 4.00, 2.00, 1.00/ < DATA RFLUX/ 370.45, 203.69, 96.00, 69.71, 50.70, 26.67, < + 17.21, 8.26, 12.86, 4.10, 5.20, 333.80, < + 2.78, 0.70, 3.07, 3.64, 3.18, 4.38, < + 1.78, 5.96, 4.22, 4.43, 1.93, 0.87, < + 0.79, 0.24, 0.20, 0.17, 0.39, 0.22, < + 0.17, 1.50, 0.45, 0.48, 1.58, 0.80, < + 0.51, 0.31, 0.18, 0.39, 0.21, 0.74, < + 0.87, 6.00, 0.24, 0.84, 0.10, 0.27, < + 0.92, 1.84, 0.13, 0.38, 0.0215, 0.0067, < + 0.0009, 0.0003, 1.E-6, 3.E-9, 1.E-11/ < DATA XFLUX/ 464.20, 241.50, 131.50, 101.90, 81.32, 48.71, < + 37.16, 21.14, 30.70, 11.20, 12.00, 438.80, < + 6.50, 1.60, 6.40, 8.66, 9.04, 13.12, < + 4.42, 13.18, 12.03, 13.29, 5.01, 2.18, < + 1.59, 0.67, 0.43, 0.43, 0.72, 0.46, < + 0.48, 3.02, 1.46, 1.02, 4.86, 1.59, < + 1.57, 1.67, 0.36, 0.99, 2.20, 1.39, < + 5.63, 11.28, 2.50, 4.14, 3.16, 0.59, < + 3.70, 4.85, 0.34, 1.15, 0.18, 0.08, < + 0.025, 0.03, 1.E-3, 3.E-5, 1.E-6/ < DATA SCALE1/35347.5, 33095.6, 18040.6, 13733.0, 12564.2, 7121.38, < + 6608.74, 5779.89, 8009.80, 3186.34, 3033.78, 47555., < + 1692.09, 405.95, 1516.20, 2731.70, 3314.57, 4375.00, < + 1316.91, 3621.91, 3908.56, 4432.54, 1541.21, 531.73, < + 364.83, 0.00, 116.00, 129.41, 162.48, 94.07, < + 41.29, 709.50, 0.00, 268.47, 1561.05, 367.64, < + 290.06, 184.36, 0.00, 86.15, 7.50, 0.00, < + 0.00, 2220.00, 0.00, 61.00, 0.00, 86.95, < + 206.00, 135.89, 60.35, 157.12, 7.06, 0.75, < + 0.00, 0.00, 0.00, 0.00, 0.00/ --- > > 1400.00, 1350.00, 1300.00, 1250.00, 1200.00, 1215.67, > > 1150.00, 1100.00, 1050.00, 1000.00, 1031.91, 1025.72, > > 950.00, 977.02, 900.00, 850.00, 800.00, 750.00, > > 789.36, 770.41, 765.15, 700.00, 703.31, 650.00, > > 600.00, 629.73, 609.76, 550.00, 584.33, 554.37, > > 500.00, 450.00, 465.22, 400.00, 350.00, 368.07, > > 300.00, 303.78, 303.31, 250.00, 284.15, 256.30, > > 200.00, 150.00, 100.00, 50.00, 32.00, 23.00, > > 16.00, 8.00, 4.00, 2.00, 1.00/ > DATA RFLUX/ 322.00, 168.00, 95.00, 62.00, 44.00, 25.00, > > 16.90, 11.80, 19.50, 4.10, 11.10, 249.00, > > 2.78, 0.70, 3.07, 3.64, 3.18, 4.38, > > 1.78, 5.96, 4.22, 4.43, 1.93, 0.87, > > 0.79, 0.24, 0.20, 0.17, 0.39, 0.22, > > 0.17, 1.50, 0.45, 0.48, 1.58, 0.80, > > 0.51, 0.31, 0.18, 0.39, 0.21, 0.74, > > 0.87, 6.00, 0.24, 0.84, 0.10, 0.27, > > 0.92, 1.84, 0.13, 0.38, 0.0215, 0.0067, > > 1.E-3, 2.E-3, 1.E-5, 5.E-8, 1.E-10/ > DATA XFLUX/ 354.00, 191.00, 110.00, 76.00, 55.00, 28.00, > > 19.60, 14.30, 25.30, 5.00, 17.20, 401.00, > > 6.26, 1.51, 6.11, 8.66, 9.04, 13.12, > > 4.42, 13.18, 12.03, 13.29, 5.01, 2.18, > > 1.59, 0.67, 0.43, 0.43, 0.72, 0.46, > > 0.48, 3.02, 1.46, 1.02, 4.86, 1.59, > > 1.57, 1.67, 0.36, 0.99, 2.20, 1.39, > > 5.63, 11.28, 2.50, 4.14, 3.16, 0.59, > > 3.70, 4.85, 0.34, 1.15, 0.18, 0.08, > > 2.5E-2, 5.E-2, 8.E-4, 3.E-5, 5.E-7/ > DATA SCALE1/ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 1692.09, 405.95, 1516.20, 2731.70, 3314.57, 4375.00, > > 1316.91, 3621.91, 3908.56, 4432.54, 1541.21, 531.73, > > 364.83, 0.00, 116.00, 129.41, 162.48, 94.07, > > 41.29, 709.50, 0.00, 268.47, 1561.05, 367.64, > > 290.06, 184.36, 0.00, 86.15, 7.50, 0.00, > > 0.00, 2220.00, 0.00, 61.00, 0.00, 86.95, > > 206.00, 135.89, 60.35, 157.12, 7.06, 0.75, > > 0.00, 0.00, 0.00, 0.00, 0.00/ 74,191c81,175 < + 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, < + 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, < + 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, < + 0.00, 5.34, 0.00, 0.00, 0.00, 0.54, < + 3.30, 0.00, 12.60, 0.00, 0.00, 0.00, < + 5.34, 11.63, 2.28, 5.56, 24.93, 8.16, < + 60.69, 0.00, 28.20, 45.90, 40.80, 1.27, < + 35.47, 42.80, 1.12, 6.19, 1.26, 0.69, < + 0.23, 0.30, 0.01, 3.E-4, 1.E-5/ < DATA SRA/ 0.536, 0.216, 0.203, 0.184, 0.175, 0.126, < + 0.114, 0.073/ < DATA SRB/ 334.0, 189.0, 82.2, 57.2, 38.8, 18.1, < + 9.46, 3.30/ < DATA CHINT/0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.18223E-01, < + -0.14763E-02, 0.13765E-02, 0.10044E-01, 0.48687E-01, < + 0.79848E-02, 0.16334E-02, 0.99890E-02, 0.12362E-01, < + 0.38317E-02, 0.00000E+00, 0.22214E-02, 0.87883E-03, < + 0.65437E-02, 0.20197E-02, 0.11897E-02, 0.25303E-01, < + 0.00000E+00, 0.69859E-02, 0.15090E-02, 0.15679E-01, < + 0.71813E-02, 0.18873E-02, 0.00000E+00, 0.44750E-02, < + 0.43206E-03, 0.00000E+00, 0.00000E+00, 0.24958E+00, < + 0.00000E+00, 0.17434E-03, 0.00000E+00, 0.11585E-01, < + 0.11532E-01, 0.98362E-02, 0.24461E-02, 0.91696E-02, < + 0.86457E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00/ < DATA CHSLP/0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.13097E-01, < + 0.15670E-01, 0.20808E-01, 0.66142E-02, 0.18084E-01, < + 0.20603E-01, 0.24699E-01, 0.91077E-02, 0.49109E-02, < + 0.79057E-03, 0.00000E+00, 0.73979E-03, 0.83584E-03, < + 0.11269E-02, 0.68429E-03, 0.34588E-03, 0.54965E-02, < + 0.00000E+00, 0.23091E-02, 0.13033E-01, 0.32354E-02, < + 0.27878E-02, 0.17867E-02, 0.00000E+00, 0.10032E-02, < + 0.96861E-04, 0.00000E+00, 0.00000E+00, 0.35650E-01, < + 0.00000E+00, 0.96698E-03, 0.00000E+00, 0.16549E-02, < + 0.44742E-02, 0.38025E-02, 0.23655E-02, 0.10185E-01, < + 0.96072E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00/ < DATA COINT/0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.74083E-03, < + 0.00000E+00, 0.13995E-02, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.13215E-02,-0.52657E-03, 0.00000E+00, < + 0.17231E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + -0.57828E-02,-0.14201E-01, 0.44921E-02, 0.81510E-03, < + -0.38050E-01, 0.25868E-01,-0.74031E-01, 0.00000E+00, < + -0.49781E-01,-0.65154E-01,-0.93662E-01, 0.00000E+00, < + -0.59300E-01, 0.12726E-01, 0.32644E-02,-0.64290E-02, < + -0.17720E-01,-0.17972E-01, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00/ < DATA COSLP/0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.35613E-04, < + 0.00000E+00, 0.67265E-04, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.71233E-05, 0.45645E-04, 0.00000E+00, < + 0.20078E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, < + 0.98801E-04, 0.23753E-03, 0.45868E-04, 0.12901E-03, < + 0.64519E-03, 0.20540E-03, 0.17408E-02, 0.00000E+00, < + 0.90311E-03, 0.17093E-02, 0.13948E-02, 0.00000E+00, < + 0.15455E-02, 0.23847E-02, 0.87499E-04, 0.79916E-03, < + 0.34005E-03, 0.31881E-03, 0.00000E+00, 0.00000E+00, < + 0.00000E+00, 0.00000E+00, 0.00000E+00/ < DATA W1/ 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, < + 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, < + 0.00000, 0.00000, 0.00000, 0.50000, 0.50000, 0.50000, < + 0.50000, 0.50000, 0.50000, 0.50000, 0.50000, 0.00000, < + 0.50000, 0.00000, 0.50000, 0.50000, 0.50000, 0.90000, < + 4.50000, 0.50000, 0.00000, 0.50000, 0.50000, 0.50000, < + 0.55000, 0.40000, 0.00000, 0.07000, 0.02000, 0.00000, < + 0.00000, 1.00000, 0.00000, 0.18000, 0.00000, 1.00000, < + 0.33750, 1.12500, 1.12500, 0.80000, 0.24000, 0.00000, < + 0.00000, 0.00000, 0.00000, 0.00000, 0.00000/ < DATA W2/ 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, < + 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, < + 0.00000, 0.00000, 0.00000, 0.22750, 0.22750, 0.22750, < + 0.22750, 0.22750, 0.22750, 0.22750, 0.22750, 0.00000, < + 0.22750, 1.50000, 0.22750, 0.22750, 0.22750, 0.30000, < + 3.00000, 0.22750, 0.80000, 0.22750, 0.22750, 0.22750, < + 0.25000, 0.90000, 1.30000, 0.13000, 1.18000, 0.50000, < + 0.80000, 0.22750, 2.00000, 2.07000, 1.05000, 0.22750, < + 1.91250, 1.12500, 1.12500, 1.20000, 1.76000, 2.00000, < + 0.00000, 0.00000, 0.00000, 0.00000, 0.00000/ < DATA WCHR/ 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, < + 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, < + 0.00000, 0.00000, 0.00000, 1.00000, 1.00000, 1.10000, < + 2.00000, 0.30000, 1.50000, 1.80000, 2.00000, 0.00000, < + 8.00000, 1.00000, 3.50000, 1.30000, 0.10000, 1.00000, < + 1.00000, 0.10000, 1.00000, 2.50000, 0.30000, 0.10000, < + 1.00000, 1.00000, 1.00000, 1.00000, 1.00000, 1.00000, < + 1.00000, 0.10000, 1.00000, 1.00000, 1.00000, 1.00000, < + 1.00000, 1.00000, 1.00000, 1.00000, 1.00000, 1.00000, < + 0.00000, 0.00000, 0.00000, 0.00000, 0.00000/ < DATA WAVEAV/ 0.000, 0.000, 0.000, 0.000, < + 0.000, 0.000, 0.000, 0.000, < + 0.000, 0.000, 0.000, 0.000, < + 0.000, 0.000, 0.000, 1025.500, < + 1031.900, 1025.700, 975.500, 977.020, < + 925.500, 875.500, 825.500, 775.005, < + 787.710, 770.410, 765.150, 725.500, < + 703.360, 675.150, 626.930, 629.730, < + 609.760, 576.985, 584.330, 554.370, < + 525.000, 476.185, 465.220, 418.920, < + 377.915, 368.070, 326.580, 303.780, < + 303.310, 275.300, 284.150, 256.320, < + 224.600, 174.340, 124.470, 75.255, < + 39.620, 24.070, 0.000, 0.000, < + 0.000, 0.000, 0.000/ --- > > 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, > > 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, > > 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, > > 0.00, 5.34, 0.00, 0.00, 0.00, 0.54, > > 3.30, 0.00, 12.60, 0.00, 0.00, 0.00, > > 5.34, 11.63, 2.28, 5.56, 24.93, 8.16, > > 60.69, 0.00, 28.20, 45.90, 40.80, 1.27, > > 35.47, 42.80, 1.12, 6.19, 1.26, 0.69, > > 0.23, 0.46, 7.6E-3, 2.9E-4, 4.8E-6/ > DATA TCHR0/ > > 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 0.0, 0.0, 0.0,-4.290E+00,-5.709E+00,-8.493E+00, > >-1.161E+00,-3.429E+00,-5.464E+00,-6.502E+00,-1.912E+00,-4.034E-01, > >-1.448E-01, 0.000E+00,-9.702E-02,-6.591E-02,-2.338E-02,-1.273E-01, > >-2.406E-01,-3.351E-01, 0.000E+00,-1.465E+00,-2.405E+00,-7.975E-02, > >-4.197E-01,-1.971E-01, 0.000E+00,-5.895E-02,-5.815E-03, 0.000E+00, > > 0.000E+00, 2.138E-01, 0.000E+00,-7.713E-02, 0.000E+00,-3.035E-02, > >-2.039E-01,-1.749E-01,-1.041E-01,-2.638E-01,-1.094E-02, 0.000E+00, > > 0.0, 0.0, 0.0, 0.0, 0.0/ > DATA TCHR1/ > > 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 0.0, 0.0, 0.0,-3.023E-13,-3.745E-13,-5.385E-13, > >-1.211E-13,-3.868E-13,-3.646E-13,-4.125E-13,-1.527E-13,-4.753E-14, > >-3.411E-14, 0.000E+00,-1.190E-14,-1.034E-14,-1.343E-14,-1.539E-14, > >-5.174E-14,-6.934E-14, 0.000E+00,-1.215E-13,-1.537E-13,-2.024E-14, > >-4.596E-14,-1.562E-14, 0.000E+00,-1.221E-14,-1.123E-15, 0.000E+00, > > 0.000E+00,-2.263E-13, 0.000E+00,-1.508E-14, 0.000E+00,-1.744E-14, > >-2.100E-14,-1.805E-14,-8.224E-15,-1.919E-14,-7.944E-16, 0.000E+00, > > 0.0, 0.0, 0.0, 0.0, 0.0/ > DATA TCHR2/ > > 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 0.0, 0.0, 0.0, 3.275E-11, 4.057E-11, 6.160E-11, > > 1.312E-11, 4.189E-11, 4.167E-11, 4.716E-11, 1.654E-11, 5.150E-12, > > 3.901E-12, 0.000E+00, 1.289E-12, 1.120E-12, 1.455E-12, 1.667E-12, > > 5.604E-12, 7.931E-12, 0.000E+00, 1.317E-11, 1.757E-11, 2.194E-12, > > 4.978E-12, 1.693E-12, 0.000E+00, 1.324E-12, 1.285E-13, 0.000E+00, > > 0.000E+00, 2.586E-11, 0.000E+00, 1.724E-12, 0.000E+00, 1.889E-12, > > 2.400E-12, 2.063E-12, 8.911E-13, 2.193E-12, 9.090E-14, 0.000E+00, > > 0.0, 0.0, 0.0, 0.0, 0.0/ > DATA TCOR0/ > > 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 0.0, 0.0, 0.0, 0.000E+00, 0.000E+00, 0.000E+00, > > 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00,-6.060E-02, > > 0.000E+00,-3.399E-02, 0.000E+00, 0.000E+00, 0.000E+00, 4.866E-02, > >-1.762E-01, 0.000E+00,-2.412E-01, 0.000E+00, 0.000E+00, 0.000E+00, > >-4.743E-01,-9.713E-01, 5.891E-02,-1.263E-01,-1.246E+00, 2.870E-01, > >-4.659E+00, 0.000E+00,-1.058E+00,-3.821E+00,-1.874E+00, 0.000E+00, > >-1.896E+00,-8.505E-01,-2.101E-04,-2.012E-01,-6.097E-02,-2.925E-02, > >-4.875E-03, 0.0, 0.0, 0.0, 0.0/ > DATA TCOR1/ > > 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 0.0, 0.0, 0.0, 0.000E+00, 0.000E+00, 0.000E+00, > > 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00, 2.877E-03, > > 0.000E+00, 1.760E-03, 0.000E+00, 0.000E+00, 0.000E+00, 3.313E-04, > > 3.643E-03, 0.000E+00, 5.225E-03, 0.000E+00, 0.000E+00, 0.000E+00, > > 4.085E-03, 1.088E-02, 8.447E-04, 3.237E-03, 1.907E-02, 2.796E-03, > > 4.460E-02, 0.000E+00, 1.007E-02, 3.481E-02, 1.604E-02, 0.000E+00, > > 2.029E-02, 2.160E-02, 6.342E-04, 3.594E-03, 5.503E-04, 2.687E-04, > > 4.479E-05, 0.0, 0.0, 0.0, 0.0/ > DATA TCOR2/ > > 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, > > 0.0, 0.0, 0.0, 0.000E+00, 0.000E+00, 0.000E+00, > > 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00, 1.846E-03, > > 0.000E+00, 1.127E-03, 0.000E+00, 0.000E+00, 0.000E+00, 1.891E-04, > > 2.326E-03, 0.000E+00, 2.801E-03, 0.000E+00, 0.000E+00, 0.000E+00, > > 2.446E-03, 7.121E-03, 5.204E-04, 1.983E-03, 1.204E-02, 1.721E-03, > > 2.911E-02, 0.000E+00, 7.177E-03, 2.272E-02, 9.436E-03, 0.000E+00, > > 1.316E-02, 1.398E-02, 4.098E-04, 2.328E-03, 3.574E-04, 1.745E-04, > > 2.909E-05, 0.0, 0.0, 0.0, 0.0/ > DATA WAR1/ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, > > 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, > > 0.00, 0.00, 3.80, 6.25, 4.93, 6.06, > > 2.70, 7.07, 8.62, 9.60, 4.54, 2.37, > > 0.82, 0.33, 0.24, 0.67, 0.28, 0.55, > > 1.56, 1.11, 0.77, 1.32, 1.71, 0.44, > > 1.11, 0.95, 0.39, 0.81, 2.00, 1.49, > > 6.81, 5.07, 1.63, 5.62, 2.08, 0.59, > > 3.89, 5.19, 0.35, 1.18, 0.099, 0.04, > > 0.007, 0.00, 0.00, 0.00, 0.00/ > DATA WAR2/ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, > > 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, > > 0.00, 0.00, 20.80, 17.90, 9.30, 14.30, > > 6.90, 12.00, 15.60, 18.60, 10.10, 4.30, > > 12.40, 8.00, 3.60, 1.80, 0.50, 1.40, > > 3.90, 2.60, 1.60, 3.40, 4.10, 0.70, > > 4.30, 4.30, 3.80, 2.60, 6.08, 1.35, > > 12.60, 9.78, 2.96, 10.20, 4.11, 6.68, > > 6.62, 8.07, 0.47, 1.73, 0.17, 0.075, > > 0.012, 0.00, 0.00, 0.00, 0.00/ 193c177 < C ************** HINTEREGGER METHOD ************************** --- > C Linear Interpolation between SC#21REFW and F79050: 194a179,185 > FRAT = (F107-68.) / (243.-68.) > DO 200 L=1,LMAX > SFLUX(L) = RFLUX(L) + (XFLUX(L)-RFLUX(L)) * FRAT > 200 CONTINUE > C > C Hinteregger contrast ratio method: > C 196,200c187,188 < R1 = B1(1) + B1(2)*(F107A-71.5) + B1(3)*(F107-F107A+3.9) < R2 = B2(1) + B2(2)*(F107A-71.5) + B2(3)*(F107-F107A+3.9) < DO 100 L=1,LM < IF (L .LT. 9) THEN < SFLUX(L) = SRA(L) * F107 + SRB(L) --- > IF (HLYBR .GT. 0.001) THEN > R1 = HLYBR 202,204c190,210 < IF (L .EQ. 12) THEN < C SFLUX(L) = CHINT(L) + (CHSLP(L)*WCHR(L)*FLYA) < SFLUX(L) = 332. + 0.6 * (F107 - 65.) --- > R1 = B1(1) + B1(2)*(F107A-71.5) + B1(3)*(F107-F107A+3.9) > ENDIF > IF (FEXVIR .GT. 0.001) THEN > R2 = FEXVIR > ELSE > R2 = B2(1) + B2(2)*(F107A-71.5) + B2(3)*(F107-F107A+3.9) > ENDIF > DO 100 L=13,LMAX > SFLUX(L) = (RFLUX(L) + ((R1-1.)*SCALE1(L) > > + (R2-1.)*SCALE2(L)) / 1000.) > 100 CONTINUE > ENDIF > C > C Tobiska EUV91 Method: > C > IF (ISCALE .EQ. 2) THEN > IF (HLYA .GT. 0.001) THEN > HLYMOD = HLYA > ELSE > IF (HEIEW .GT. 0.001) THEN > HLYMOD = HEIEW * 3.77847E9 + 8.40317E10 206,208c212,213 < SFLUX(L) = (RFLUX(L) + ((R1-1.)*SCALE1(L) < + + (R2-1.)*SCALE2(L)) / 1000.) < IF (SFLUX(L) .LT. 0.0) SFLUX(L) = 0.0 --- > HLYMOD = 8.70E8 * F107 + 1.90E11 > HLYMOD = 8.70E8 * F107 + 1.90E11 211,212c216,225 < 100 CONTINUE < ELSE --- > IF (HEIEW .GT. 0.001) THEN > HEIMOD = HEIEW * 3.77847E9 + 8.40317E10 > ELSE > HEIMOD = HLYMOD > ENDIF > DO 500 L=16,55 > SFLUX(L) = TCHR0(L) + TCHR1(L)*HLYMOD + TCHR2(L)*HEIMOD > > + TCOR0(L) + TCOR1(L)*F107 + TCOR2(L)*F107A > 500 CONTINUE > ENDIF 214c227 < C **************** LINEAR INTERPOLATION METHOD ******************** --- > C Woods and Rottman (10 Nov. 1988) spectrum: 216,219c229,232 < FRAT = (F107-68.) / (243.-68.) < DO 200 L=1,LM < SFLUX(L) = RFLUX(L) + (XFLUX(L)-RFLUX(L)) * FRAT < 200 CONTINUE --- > IF (ISCALE .EQ. 3) THEN > DO 550 L=15,55 > SFLUX(L) = WAR1(L) > 550 CONTINUE 222c235 < C ********************* TOBISKA AND BARTH METHOD ******************** --- > C Woods and Rottman (20 June 1989) spectrum: 224,240c237,240 < F107MD = 0 < IF (ISCALE.EQ.2) THEN < IF (FLYA.LT.1) THEN < FLYA = (F107 + 218.88)/1.05453E-9 < F107MD = F107 < ELSE < F107MD = (1.05453E-9*FLYA) - 218.88 < ENDIF < FLYA = FLYA*((12400*1.602E-12)/(1215.67)) !TO ERGS < DO 500 J = 16,54 < FLUXCH(J)=CHINT(J)+(CHSLP(J)*WCHR(J)*FLYA) < FLUXCO(J)=COINT(J)+COSLP(J)*((W1(J)*F107MD)+(W2(J)*F107)) < SFLUX(J) = FLUXCH(J) + FLUXCO(J) < SFLUX(J) = SFLUX(J)*((WAVEAV(J))/(12400*1.602E-12)) !TO PHOTONS < SFLUX(J) = SFLUX(J)*1.E-9 !TO COMPENSATE FOR *1.E9 THAT FOLLOWS < 500 CONTINUE < FLYA = FLYA*((1215.67)/(12400*1.602E-12)) !BACK TO PHOTONS --- > IF (ISCALE .EQ. 4) THEN > DO 560 L=15,55 > SFLUX(L) = WAR2(L) > 560 CONTINUE 242,244d241 < C **************************************************************** < SFNORM = 2. - (F107-68.) / (243.-68.) < IF (SFNORM .LT. 1.0) SFNORM = 1.0 246,252c243,262 < DO 300 L=1,LM < WAVE1(L) = WAVEL(L) < WAVE2(L) = WAVES(L) < IF (WAVE1(L) .LT. 251. .AND. WAVE2(L) .GT. 0.1) < + SFLUX(L) = SFLUX(L) * SFNORM < SFLUX(L) = SFLUX(L) * 1.E9 < 300 CONTINUE --- > C Substitute in H Lyman-alpha and XUVFAC if provided: > C > IF (HLYA .GT. 0.001) SFLUX(12) = HLYA / 1.E9 > IF (XUVFAC .GT. 0.001) THEN > XUVF = XUVFAC > ELSE > XUVF = 1.0 > ENDIF > C > C Convert from gigaphotons to photons, etc.: > C > DO 600 L=1,LMAX > WAVE1(L) = WAVEL(L) > WAVE2(L) = WAVES(L) > IF (SFLUX(L) .LT. 0.0) SFLUX(L) = 0.0 > IF (WAVEL(L).LT.251.0 .AND. WAVES(L).GT.15.0) > > SFLUX(L)=SFLUX(L)*XUVF > SFLUX(L) = SFLUX(L) * 1.E9 > 600 CONTINUE > C 255d264 < C ======================================================================== Diff of /home/tgcm/tgcm13mt/start.f and start.f: 368a369,374 > ! > ! Get geographic coords at ends of magnetic field lines: > ! (subroutine othend defines /othend_com/ in othend.h) > ! > call othend > ! ======================================================================== Diff of /home/tgcm/tgcm13mt/tail.f and tail.f: 79c79 < IF (IFRST .GT. 1 .AND. IPR .EQ. 0) RETURN --- > C IF (IFRST .GT. 1 .AND. IPR .EQ. 0) RETURN 141,142c141,144 < H1 = AMIN1 (2.35, 0.83 + 0.33 * PLEVEL ) < H2 = 2.87 + 0.15 * PLEVEL --- > C H1 = AMIN1 (2.35, 0.83 + 0.33 * PLEVEL ) > C H2 = 2.87 + 0.15 * PLEVEL > H1 = 3.+0.1*HP > H2 = 3.0+0.1*HP 150,151c152 < ! C(61) = F107 < ! C(62) = F107A --- > ! (only if they were not defined by getgpi) 177,178c178,181 < EC = 0.01+0.09*HP/100. < ALFAC = 0.5 --- > C EC = 0.01+0.09*HP/100. > C ALFAC = 0.5 > EC = 0.5 > ALFAC = 1.0 180,181c183,186 < ED = 0.01+0.2*HP/100. < ALFAD = 0.75 --- > C ED = 0.01+0.2*HP/100. > C ALFAD = 0.75 > ED = 0.5 > ALFAD = 2.0 204a210,212 > C ALFA1 = 2.0 > C ALFA2 = 3.0 > C **** SNOE AURORA PARTICLE ALPHA 206c214 < ALFA2 = 3.0 --- > ALFA2 = 2.0 281a290,292 > C ALFA1 = 2.0 > C ALFA2 = 3.0 > C **** SNOE AURORA PARTICLE ALPHA 283c294 < ALFA2 = 3.0 --- > ALFA2 = 2.0 332,333c343,362 < E1 = AMAX1( 0.50, -2.15 + 0.62 * PLEVEL) < E2 = 0.95 + 0.117 * HP --- > C E1 = AMAX1( 0.50, -2.15 + 0.62 * PLEVEL) > C E2 = 0.95 + 0.117 * HP > C **** SNOE PARTICLE ENERGY FLUX > C E1 = AMAX1( 0.50, -2.15 + 0.62 * PLEVEL)*1.5 > C E1 = (0.95 + 0.117 * HP) > C E2 = (0.95 + 0.117 * HP) > C E1 = (1.0 + 0.15* HP) > C E2 = (1.0 + 0.15* HP) > E1 = (1.0 + 0.25* HP) > E2 = (1.0 + 0.25* HP) > C H1 = AMIN1( 2.35, 0.83 + 0.33 * PLEVEL ) > C H2 = 2.87 + 0.15 * PLEVEL > H1 = 3.+0.1*HP > H2 = 3.0+0.1*HP > C > write(6,"('tail: ut=',f8.3,' hp=',e12.4,' plevel=',e12.4, > | ' cp=',e12.4)") ut,hp,plevel,cp > write(6,"(' e1=',e12.4,' e2=',e12.4,' h1=',e12.4,' h2=',e12.4)") > | e1,e2,h1,h2 > C 335,336d363 < H1 = AMIN1( 2.35, 0.83 + 0.33 * PLEVEL ) < H2 = 2.87 + 0.15 * PLEVEL 382,393c409,420 < IF (ABS(UT-IUT) .LT. 0.001) WRITE (6,603) UT, HP, CP, BYIMF, < | RHP, RCP, RAUR, (THETA0(I),I=1,51), DISPC,ARAD(1),ARAD(2), < | TWA6,TWA21,JSWOLDA,JDIDK,JDITH < 603 FORMAT(1X,'TAIL: ',4X,'UT HP CP BYIMF RHP RCP RAUR =', 7F6.1/ < |1X,'PARAMD=THETA0(2),OFFA(2),OFFC(2),DSKOFA(2),DSKOFC(2),PHID(2)'/ < |1X,' PHIN(2),PHIDP0(2),PHIDM0(2),PHINP0(2),PHINM0(2),PSIM(2)'/ < |1X,' PSIE(2),PCEN(2),ARAD(2),H1,H2,ROTH,E1,E2,ROTE,EC,ED,ALFA1'/ < |1X,' ALFA2,TWAK,TWA6,TWA21,ROT6,ROT21,D6,D21,H6,H21,T6,T21'/ < | 5X, 'PARAMD(1-10) =', 10E10.3/ 5X, 'PARAMD(11-20)=',10E10.3/ < | 5X, 'PARAMD(21-30)=', 10E10.3/ 5X, 'PARAMD(31-40)=',10E10.3/ < | 5X, 'PARAMD(41-50)=', 10E10.3/ 5X, 'PARAMD(51)=',E10.3/ < | 5X, 'DISPC ARAD(1,2) TWA6 TWA21 JSWOLDA,DIDK,DITH =', 5F6.3,3I3) --- > C IF (ABS(UT-IUT) .LT. 0.001) WRITE (6,603) UT, HP, CP, BYIMF, > C | RHP, RCP, RAUR, (THETA0(I),I=1,51), DISPC,ARAD(1),ARAD(2), > C | TWA6,TWA21,JSWOLDA,JDIDK,JDITH > C 603 FORMAT(1X,'TAIL: ',4X,'UT HP CP BYIMF RHP RCP RAUR =', 7F6.1/ > C |1X,'PARAMD=THETA0(2),OFFA(2),OFFC(2),DSKOFA(2),DSKOFC(2),PHID(2)'/ > C |1X,' PHIN(2),PHIDP0(2),PHIDM0(2),PHINP0(2),PHINM0(2),PSIM(2)'/ > C |1X,' PSIE(2),PCEN(2),ARAD(2),H1,H2,ROTH,E1,E2,ROTE,EC,ED,ALFA1'/ > C |1X,' ALFA2,TWAK,TWA6,TWA21,ROT6,ROT21,D6,D21,H6,H21,T6,T21'/ > C | 5X, 'PARAMD(1-10) =', 10E10.3/ 5X, 'PARAMD(11-20)=',10E10.3/ > C | 5X, 'PARAMD(21-30)=', 10E10.3/ 5X, 'PARAMD(31-40)=',10E10.3/ > C | 5X, 'PARAMD(41-50)=', 10E10.3/ 5X, 'PARAMD(51)=',E10.3/ > C | 5X, 'DISPC ARAD(1,2) TWA6 TWA21 JSWOLDA,DIDK,DITH =', 5F6.3,3I3) ======================================================================== Diff of /home/tgcm/tgcm13mt/vdrift2.f and vdrift2.f: 51c51 < EEZ(I,K) = EZ(I,J,K)/(F(I+2,NZK+1)-F(I+2,NZK-1)) --- > EEZ(I,K) = EZ(I,J,K)/(F(I+2,NZK+1)-F(I+2,NZK-1)) 58,59c58,59 < EEZ(I,1) = 2.*EEZ(I,2)-EEZ(I,3) < EEZ(I,KMAXP1) = 2.*EEZ(I,KMAX)-EEZ(I,KMAX-1) --- > EEZ(I,1) = 2.*EEZ(I,2)-EEZ(I,3) > EEZ(I,KMAXP1) = 2.*EEZ(I,KMAX)-EEZ(I,KMAX-1) ======================================================================== Diff of /home/tgcm/tgcm13mt/xray.f and xray.f: 12a13 > include "crates.h" 18a20,23 > ! 11/30/99: add SFEPS (cons.h): > ! SFEPS is initialized to 1 in con.f, and recalculated in > ! advnce.f only if calendar day is being advanced and at > ! new day boundaries. 20c25,27 < EX=0.3+0.5*(C(61)-67.)/(176.*SF) --- > C **** SNOE X-RAYS > EX=(0.3+0.5*(C(61)-70.)/240.)*SFEPS > C EX=0.6+0.5*(C(61)-67.)/(176.*SF) 87a95 > DISN2P(I,1)=DISN2P(I,1)+0.64*S3(I,1)