NRLMSISE-00 distribution package There are three parts: 1) FORTRAN code for NRLMSISE-00 atmospheric model 2) FORTRAN code for test program 3) test program output $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ C----------------------------------------------------------------------- SUBROUTINE GTD7(IYD,SEC,ALT,GLAT,GLONG,STL,F107A,F107,AP,MASS,D,T) C C NRLMSISE-00 C ----------- C Neutral Atmosphere Empirical Model from the surface to lower C exosphere C C NEW FEATURES: C *Extensive satellite drag database used in model generation C *Revised O2 (and O) in lower thermosphere C *Additional nonlinear solar activity term C *"ANOMALOUS OXYGEN" NUMBER DENSITY, OUTPUT D(9) C At high altitudes (> 500 km), hot atomic oxygen or ionized C oxygen can become appreciable for some ranges of subroutine C inputs, thereby affecting drag on satellites and debris. We C group these species under the term "anomalous oxygen," since C their individual variations are not presently separable with C the drag data used to define this model component. C C SUBROUTINES FOR SPECIAL OUTPUTS: C C HIGH ALTITUDE DRAG: EFFECTIVE TOTAL MASS DENSITY C (SUBROUTINE GTD7D, OUTPUT D(6)) C For atmospheric drag calculations at altitudes above 500 km, C call SUBROUTINE GTD7D to compute the "effective total mass C density" by including contributions from "anomalous oxygen." C See "NOTES ON OUTPUT VARIABLES" below on D(6). C C PRESSURE GRID (SUBROUTINE GHP7) C See subroutine GHP7 to specify outputs at a pressure level C rather than at an altitude. C C OUTPUT IN M-3 and KG/M3: CALL METERS(.TRUE.) C C INPUT VARIABLES: C IYD - YEAR AND DAY AS YYDDD (day of year from 1 to 365 (or 366)) C (Year ignored in current model) C SEC - UT(SEC) C ALT - ALTITUDE(KM) C GLAT - GEODETIC LATITUDE(DEG) C GLONG - GEODETIC LONGITUDE(DEG) C STL - LOCAL APPARENT SOLAR TIME(HRS; see Note below) C F107A - 81 day AVERAGE OF F10.7 FLUX (centered on day DDD) C F107 - DAILY F10.7 FLUX FOR PREVIOUS DAY C AP - MAGNETIC INDEX(DAILY) OR WHEN SW(9)=-1. : C - ARRAY CONTAINING: C (1) DAILY AP C (2) 3 HR AP INDEX FOR CURRENT TIME C (3) 3 HR AP INDEX FOR 3 HRS BEFORE CURRENT TIME C (4) 3 HR AP INDEX FOR 6 HRS BEFORE CURRENT TIME C (5) 3 HR AP INDEX FOR 9 HRS BEFORE CURRENT TIME C (6) AVERAGE OF EIGHT 3 HR AP INDICIES FROM 12 TO 33 HRS PRIOR C TO CURRENT TIME C (7) AVERAGE OF EIGHT 3 HR AP INDICIES FROM 36 TO 57 HRS PRIOR C TO CURRENT TIME C MASS - MASS NUMBER (ONLY DENSITY FOR SELECTED GAS IS C CALCULATED. MASS 0 IS TEMPERATURE. MASS 48 FOR ALL. C MASS 17 IS Anomalous O ONLY.) C C NOTES ON INPUT VARIABLES: C UT, Local Time, and Longitude are used independently in the C model and are not of equal importance for every situation. C For the most physically realistic calculation these three C variables should be consistent (STL=SEC/3600+GLONG/15). C The Equation of Time departures from the above formula C for apparent local time can be included if available but C are of minor importance. c C F107 and F107A values used to generate the model correspond C to the 10.7 cm radio flux at the actual distance of the Earth C from the Sun rather than the radio flux at 1 AU. The following C site provides both classes of values: C ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SOLAR_RADIO/FLUX/ C C F107, F107A, and AP effects are neither large nor well C established below 80 km and these parameters should be set to C 150., 150., and 4. respectively. C C OUTPUT VARIABLES: C D(1) - HE NUMBER DENSITY(CM-3) C D(2) - O NUMBER DENSITY(CM-3) C D(3) - N2 NUMBER DENSITY(CM-3) C D(4) - O2 NUMBER DENSITY(CM-3) C D(5) - AR NUMBER DENSITY(CM-3) C D(6) - TOTAL MASS DENSITY(GM/CM3) C D(7) - H NUMBER DENSITY(CM-3) C D(8) - N NUMBER DENSITY(CM-3) C D(9) - Anomalous oxygen NUMBER DENSITY(CM-3) C T(1) - EXOSPHERIC TEMPERATURE C T(2) - TEMPERATURE AT ALT C C NOTES ON OUTPUT VARIABLES: C TO GET OUTPUT IN M-3 and KG/M3: CALL METERS(.TRUE.) C C O, H, and N are set to zero below 72.5 km C C T(1), Exospheric temperature, is set to global average for C altitudes below 120 km. The 120 km gradient is left at global C average value for altitudes below 72 km. C C D(6), TOTAL MASS DENSITY, is NOT the same for subroutines GTD7 C and GTD7D C C SUBROUTINE GTD7 -- D(6) is the sum of the mass densities of the C species labeled by indices 1-5 and 7-8 in output variable D. C This includes He, O, N2, O2, Ar, H, and N but does NOT include C anomalous oxygen (species index 9). C C SUBROUTINE GTD7D -- D(6) is the "effective total mass density C for drag" and is the sum of the mass densities of all species C in this model, INCLUDING anomalous oxygen. C C SWITCHES: The following is for test and special purposes: C C TO TURN ON AND OFF PARTICULAR VARIATIONS CALL TSELEC(SW), C WHERE SW IS A 25 ELEMENT ARRAY CONTAINING 0. FOR OFF, 1. C FOR ON, OR 2. FOR MAIN EFFECTS OFF BUT CROSS TERMS ON C FOR THE FOLLOWING VARIATIONS C 1 - F10.7 EFFECT ON MEAN 2 - TIME INDEPENDENT C 3 - SYMMETRICAL ANNUAL 4 - SYMMETRICAL SEMIANNUAL C 5 - ASYMMETRICAL ANNUAL 6 - ASYMMETRICAL SEMIANNUAL C 7 - DIURNAL 8 - SEMIDIURNAL C 9 - DAILY AP 10 - ALL UT/LONG EFFECTS C 11 - LONGITUDINAL 12 - UT AND MIXED UT/LONG C 13 - MIXED AP/UT/LONG 14 - TERDIURNAL C 15 - DEPARTURES FROM DIFFUSIVE EQUILIBRIUM C 16 - ALL TINF VAR 17 - ALL TLB VAR C 18 - ALL TN1 VAR 19 - ALL S VAR C 20 - ALL TN2 VAR 21 - ALL NLB VAR C 22 - ALL TN3 VAR 23 - TURBO SCALE HEIGHT VAR C C To get current values of SW: CALL TRETRV(SW) C DIMENSION D(9),T(2),AP(7),DS(9),TS(2) DIMENSION ZN3(5),ZN2(4),SV(25) COMMON/GTS3C/TLB,S,DB04,DB16,DB28,DB32,DB40,DB48,DB01,ZA,T0,Z0 & ,G0,RL,DD,DB14,TR12 COMMON/MESO7/TN1(5),TN2(4),TN3(5),TGN1(2),TGN2(2),TGN3(2) COMMON/LOWER7/PTM(10),PDM(10,8) COMMON/PARM7/PT(150),PD(150,9),PS(150),PDL(25,2),PTL(100,4), $ PMA(100,10),SAM(100) COMMON/DATIM7/ISD(3),IST(2),NAM(2) COMMON/DATIME/ISDATE(3),ISTIME(2),NAME(2) COMMON/CSW/SW(25),ISW,SWC(25) COMMON/MAVG7/PAVGM(10) COMMON/DMIX/DM04,DM16,DM28,DM32,DM40,DM01,DM14 COMMON/PARMB/GSURF,RE COMMON/METSEL/IMR SAVE EXTERNAL GTD7BK DATA MN3/5/,ZN3/32.5,20.,15.,10.,0./ DATA MN2/4/,ZN2/72.5,55.,45.,32.5/ DATA ZMIX/62.5/,ALAST/99999./,MSSL/-999/ DATA SV/25*1./ IF(ISW.NE.64999) CALL TSELEC(SV) C Put identification data into common/datime/ DO 1 I=1,3 ISDATE(I)=ISD(I) 1 CONTINUE DO 2 I=1,2 ISTIME(I)=IST(I) NAME(I)=NAM(I) 2 CONTINUE C C Test for changed input V1=VTST7(IYD,SEC,GLAT,GLONG,STL,F107A,F107,AP,1) C Latitude variation of gravity (none for SW(2)=0) XLAT=GLAT IF(SW(2).EQ.0) XLAT=45. CALL GLATF(XLAT,GSURF,RE) C XMM=PDM(5,3) C C THERMOSPHERE/MESOSPHERE (above ZN2(1)) ALTT=AMAX1(ALT,ZN2(1)) MSS=MASS C Only calculate N2 in thermosphere if alt in mixed region IF(ALT.LT.ZMIX.AND.MASS.GT.0) MSS=28 C Only calculate thermosphere if input parameters changed C or altitude above ZN2(1) in mesosphere IF(V1.EQ.1..OR.ALT.GT.ZN2(1).OR.ALAST.GT.ZN2(1).OR.MSS.NE.MSSL) $ THEN CALL GTS7(IYD,SEC,ALTT,GLAT,GLONG,STL,F107A,F107,AP,MSS,DS,TS) DM28M=DM28 C metric adjustment IF(IMR.EQ.1) DM28M=DM28*1.E6 MSSL=MSS ENDIF T(1)=TS(1) T(2)=TS(2) IF(ALT.GE.ZN2(1)) THEN DO 5 J=1,9 D(J)=DS(J) 5 CONTINUE GOTO 10 ENDIF C C LOWER MESOSPHERE/UPPER STRATOSPHERE [between ZN3(1) and ZN2(1)] C Temperature at nodes and gradients at end nodes C Inverse temperature a linear function of spherical harmonics C Only calculate nodes if input changed IF(V1.EQ.1..OR.ALAST.GE.ZN2(1)) THEN TGN2(1)=TGN1(2) TN2(1)=TN1(5) TN2(2)=PMA(1,1)*PAVGM(1)/(1.-SW(20)*GLOB7S(PMA(1,1))) TN2(3)=PMA(1,2)*PAVGM(2)/(1.-SW(20)*GLOB7S(PMA(1,2))) TN2(4)=PMA(1,3)*PAVGM(3)/(1.-SW(20)*SW(22)*GLOB7S(PMA(1,3))) TGN2(2)=PAVGM(9)*PMA(1,10)*(1.+SW(20)*SW(22)*GLOB7S(PMA(1,10))) $ *TN2(4)*TN2(4)/(PMA(1,3)*PAVGM(3))**2 TN3(1)=TN2(4) ENDIF IF(ALT.GE.ZN3(1)) GOTO 6 C C LOWER STRATOSPHERE AND TROPOSPHERE [below ZN3(1)] C Temperature at nodes and gradients at end nodes C Inverse temperature a linear function of spherical harmonics C Only calculate nodes if input changed IF(V1.EQ.1..OR.ALAST.GE.ZN3(1)) THEN TGN3(1)=TGN2(2) TN3(2)=PMA(1,4)*PAVGM(4)/(1.-SW(22)*GLOB7S(PMA(1,4))) TN3(3)=PMA(1,5)*PAVGM(5)/(1.-SW(22)*GLOB7S(PMA(1,5))) TN3(4)=PMA(1,6)*PAVGM(6)/(1.-SW(22)*GLOB7S(PMA(1,6))) TN3(5)=PMA(1,7)*PAVGM(7)/(1.-SW(22)*GLOB7S(PMA(1,7))) TGN3(2)=PMA(1,8)*PAVGM(8)*(1.+SW(22)*GLOB7S(PMA(1,8))) $ *TN3(5)*TN3(5)/(PMA(1,7)*PAVGM(7))**2 ENDIF 6 CONTINUE IF(MASS.EQ.0) GOTO 50 C LINEAR TRANSITION TO FULL MIXING BELOW ZN2(1) DMC=0 IF(ALT.GT.ZMIX) DMC=1.-(ZN2(1)-ALT)/(ZN2(1)-ZMIX) DZ28=DS(3) C ***** N2 DENSITY **** DMR=DS(3)/DM28M-1. D(3)=DENSM(ALT,DM28M,XMM,TZ,MN3,ZN3,TN3,TGN3,MN2,ZN2,TN2,TGN2) D(3)=D(3)*(1.+DMR*DMC) C ***** HE DENSITY **** D(1)=0 IF(MASS.NE.4.AND.MASS.NE.48) GOTO 204 DMR=DS(1)/(DZ28*PDM(2,1))-1. D(1)=D(3)*PDM(2,1)*(1.+DMR*DMC) 204 CONTINUE C **** O DENSITY **** D(2)=0 D(9)=0 216 CONTINUE C ***** O2 DENSITY **** D(4)=0 IF(MASS.NE.32.AND.MASS.NE.48) GOTO 232 DMR=DS(4)/(DZ28*PDM(2,4))-1. D(4)=D(3)*PDM(2,4)*(1.+DMR*DMC) 232 CONTINUE C ***** AR DENSITY **** D(5)=0 IF(MASS.NE.40.AND.MASS.NE.48) GOTO 240 DMR=DS(5)/(DZ28*PDM(2,5))-1. D(5)=D(3)*PDM(2,5)*(1.+DMR*DMC) 240 CONTINUE C ***** HYDROGEN DENSITY **** D(7)=0 C ***** ATOMIC NITROGEN DENSITY **** D(8)=0 C C TOTAL MASS DENSITY C IF(MASS.EQ.48) THEN D(6) = 1.66E-24*(4.*D(1)+16.*D(2)+28.*D(3)+32.*D(4)+40.*D(5)+ & D(7)+14.*D(8)) IF(IMR.EQ.1) D(6)=D(6)/1000. ENDIF T(2)=TZ 10 CONTINUE GOTO 90 50 CONTINUE DD=DENSM(ALT,1.,0,TZ,MN3,ZN3,TN3,TGN3,MN2,ZN2,TN2,TGN2) T(2)=TZ 90 CONTINUE ALAST=ALT RETURN END C----------------------------------------------------------------------- SUBROUTINE GTD7D(IYD,SEC,ALT,GLAT,GLONG,STL,F107A,F107,AP,MASS, $ D,T) C C NRLMSISE-00 C ----------- C This subroutine provides Effective Total Mass Density for C output D(6) which includes contributions from "anomalous C oxygen" which can affect satellite drag above 500 km. This C subroutine is part of the distribution package for the C Neutral Atmosphere Empirical Model from the surface to lower C exosphere. See subroutine GTD7 for more extensive comments. C C INPUT VARIABLES: C IYD - YEAR AND DAY AS YYDDD (day of year from 1 to 365 (or 366)) C (Year ignored in current model) C SEC - UT(SEC) C ALT - ALTITUDE(KM) C GLAT - GEODETIC LATITUDE(DEG) C GLONG - GEODETIC LONGITUDE(DEG) C STL - LOCAL APPARENT SOLAR TIME(HRS; see Note below) C F107A - 81 day AVERAGE OF F10.7 FLUX (centered on day DDD) C F107 - DAILY F10.7 FLUX FOR PREVIOUS DAY C AP - MAGNETIC INDEX(DAILY) OR WHEN SW(9)=-1. : C - ARRAY CONTAINING: C (1) DAILY AP C (2) 3 HR AP INDEX FOR CURRENT TIME C (3) 3 HR AP INDEX FOR 3 HRS BEFORE CURRENT TIME C (4) 3 HR AP INDEX FOR 6 HRS BEFORE CURRENT TIME C (5) 3 HR AP INDEX FOR 9 HRS BEFORE CURRENT TIME C (6) AVERAGE OF EIGHT 3 HR AP INDICIES FROM 12 TO 33 HRS PRIOR C TO CURRENT TIME C (7) AVERAGE OF EIGHT 3 HR AP INDICIES FROM 36 TO 57 HRS PRIOR C TO CURRENT TIME C MASS - MASS NUMBER (ONLY DENSITY FOR SELECTED GAS IS C CALCULATED. MASS 0 IS TEMPERATURE. MASS 48 FOR ALL. C MASS 17 IS Anomalous O ONLY.) C C NOTES ON INPUT VARIABLES: C UT, Local Time, and Longitude are used independently in the C model and are not of equal importance for every situation. C For the most physically realistic calculation these three C variables should be consistent (STL=SEC/3600+GLONG/15). C The Equation of Time departures from the above formula C for apparent local time can be included if available but C are of minor importance. c C F107 and F107A values used to generate the model correspond C to the 10.7 cm radio flux at the actual distance of the Earth C from the Sun rather than the radio flux at 1 AU. C C OUTPUT VARIABLES: C D(1) - HE NUMBER DENSITY(CM-3) C D(2) - O NUMBER DENSITY(CM-3) C D(3) - N2 NUMBER DENSITY(CM-3) C D(4) - O2 NUMBER DENSITY(CM-3) C D(5) - AR NUMBER DENSITY(CM-3) C D(6) - TOTAL MASS DENSITY(GM/CM3) [includes anomalous oxygen] C D(7) - H NUMBER DENSITY(CM-3) C D(8) - N NUMBER DENSITY(CM-3) C D(9) - Anomalous oxygen NUMBER DENSITY(CM-3) C T(1) - EXOSPHERIC TEMPERATURE C T(2) - TEMPERATURE AT ALT C DIMENSION D(9),T(2),AP(7),DS(9),TS(2) COMMON/METSEL/IMR CALL GTD7(IYD,SEC,ALT,GLAT,GLONG,STL,F107A,F107,AP,MASS,D,T) C TOTAL MASS DENSITY C IF(MASS.EQ.48) THEN D(6) = 1.66E-24*(4.*D(1)+16.*D(2)+28.*D(3)+32.*D(4)+40.*D(5)+ & D(7)+14.*D(8)+16.*D(9)) IF(IMR.EQ.1) D(6)=D(6)/1000. ENDIF RETURN END C----------------------------------------------------------------------- SUBROUTINE GHP7(IYD,SEC,ALT,GLAT,GLONG,STL,F107A,F107,AP, $ D,T,PRESS) C FIND ALTITUDE OF PRESSURE SURFACE (PRESS) FROM GTD7 C INPUT: C IYD - YEAR AND DAY AS YYDDD C SEC - UT(SEC) C GLAT - GEODETIC LATITUDE(DEG) C GLONG - GEODETIC LONGITUDE(DEG) C STL - LOCAL APPARENT SOLAR TIME(HRS) C F107A - 3 MONTH AVERAGE OF F10.7 FLUX C F107 - DAILY F10.7 FLUX FOR PREVIOUS DAY C AP - MAGNETIC INDEX(DAILY) OR WHEN SW(9)=-1. : C - ARRAY CONTAINING: C (1) DAILY AP C (2) 3 HR AP INDEX FOR CURRENT TIME C (3) 3 HR AP INDEX FOR 3 HRS BEFORE CURRENT TIME C (4) 3 HR AP INDEX FOR 6 HRS BEFORE CURRENT TIME C (5) 3 HR AP INDEX FOR 9 HRS BEFORE CURRENT TIME C (6) AVERAGE OF EIGHT 3 HR AP INDICIES FROM 12 TO 33 HRS PRIOR C TO CURRENT TIME C (7) AVERAGE OF EIGHT 3 HR AP INDICIES FROM 36 TO 59 HRS PRIOR C TO CURRENT TIME C PRESS - PRESSURE LEVEL(MB) C OUTPUT: C ALT - ALTITUDE(KM) C D(1) - HE NUMBER DENSITY(CM-3) C D(2) - O NUMBER DENSITY(CM-3) C D(3) - N2 NUMBER DENSITY(CM-3) C D(4) - O2 NUMBER DENSITY(CM-3) C D(5) - AR NUMBER DENSITY(CM-3) C D(6) - TOTAL MASS DENSITY(GM/CM3) C D(7) - H NUMBER DENSITY(CM-3) C D(8) - N NUMBER DENSITY(CM-3) C D(9) - HOT O NUMBER DENSITY(CM-3) C T(1) - EXOSPHERIC TEMPERATURE C T(2) - TEMPERATURE AT ALT C COMMON/PARMB/GSURF,RE COMMON/METSEL/IMR DIMENSION D(9),T(2),AP(7) SAVE DATA BM/1.3806E-19/,RGAS/831.4/ DATA TEST/.00043/,LTEST/12/ PL=ALOG10(PRESS) C Initial altitude estimate IF(PL.GE.-5.) THEN IF(PL.GT.2.5) ZI=18.06*(3.00-PL) IF(PL.GT..75.AND.PL.LE.2.5) ZI=14.98*(3.08-PL) IF(PL.GT.-1..AND.PL.LE..75) ZI=17.8*(2.72-PL) IF(PL.GT.-2..AND.PL.LE.-1.) ZI=14.28*(3.64-PL) IF(PL.GT.-4..AND.PL.LE.-2.) ZI=12.72*(4.32-PL) IF(PL.LE.-4.) ZI=25.3*(.11-PL) IDAY=MOD(IYD,1000) CL=GLAT/90. CL2=CL*CL IF(IDAY.LT.182) CD=1.-IDAY/91.25 IF(IDAY.GE.182) CD=IDAY/91.25-3. CA=0 IF(PL.GT.-1.11.AND.PL.LE.-.23) CA=1.0 IF(PL.GT.-.23) CA=(2.79-PL)/(2.79+.23) IF(PL.LE.-1.11.AND.PL.GT.-3.) CA=(-2.93-PL)/(-2.93+1.11) Z=ZI-4.87*CL*CD*CA-1.64*CL2*CA+.31*CA*CL ENDIF IF(PL.LT.-5.) Z=22.*(PL+4.)**2+110 C ITERATION LOOP L=0 10 CONTINUE L=L+1 CALL GTD7(IYD,SEC,Z,GLAT,GLONG,STL,F107A,F107,AP,48,D,T) XN=D(1)+D(2)+D(3)+D(4)+D(5)+D(7)+D(8) P=BM*XN*T(2) IF(IMR.EQ.1) P=P*1.E-6 DIFF=PL-ALOG10(P) IF(ABS(DIFF).LT.TEST .OR. L.EQ.LTEST) GOTO 20 XM=D(6)/XN/1.66E-24 IF(IMR.EQ.1) XM = XM*1.E3 G=GSURF/(1.+Z/RE)**2 SH=RGAS*T(2)/(XM*G) C New altitude estimate using scale height IF(L.LT.6) THEN Z=Z-SH*DIFF*2.302 ELSE Z=Z-SH*DIFF ENDIF GOTO 10 20 CONTINUE IF(L.EQ.LTEST) WRITE(6,100) PRESS,DIFF 100 FORMAT(1X,29HGHP7 NOT CONVERGING FOR PRESS, 1PE12.2,E12.2) ALT=Z RETURN END C----------------------------------------------------------------------- SUBROUTINE GLATF(LAT,GV,REFF) C CALCULATE LATITUDE VARIABLE GRAVITY (GV) AND EFFECTIVE C RADIUS (REFF) REAL LAT SAVE DATA DGTR/1.74533E-2/ C2 = COS(2.*DGTR*LAT) GV = 980.616*(1.-.0026373*C2) REFF = 2.*GV/(3.085462E-6 + 2.27E-9*C2)*1.E-5 RETURN END C----------------------------------------------------------------------- FUNCTION VTST7(IYD,SEC,GLAT,GLONG,STL,F107A,F107,AP,IC) C Test if geophysical variables or switches changed and save C Return 0 if unchanged and 1 if changed DIMENSION AP(7),IYDL(2),SECL(2),GLATL(2),GLL(2),STLL(2) DIMENSION FAL(2),FL(2),APL(7,2),SWL(25,2),SWCL(25,2) COMMON/CSW/SW(25),ISW,SWC(25) SAVE DATA IYDL/2*-999/,SECL/2*-999./,GLATL/2*-999./,GLL/2*-999./ DATA STLL/2*-999./,FAL/2*-999./,FL/2*-999./,APL/14*-999./ DATA SWL/50*-999./,SWCL/50*-999./ VTST7=0 IF(IYD.NE.IYDL(IC)) GOTO 10 IF(SEC.NE.SECL(IC)) GOTO 10 IF(GLAT.NE.GLATL(IC)) GOTO 10 IF(GLONG.NE.GLL(IC)) GOTO 10 IF(STL.NE.STLL(IC)) GOTO 10 IF(F107A.NE.FAL(IC)) GOTO 10 IF(F107.NE.FL(IC)) GOTO 10 DO 5 I=1,7 IF(AP(I).NE.APL(I,IC)) GOTO 10 5 CONTINUE DO 7 I=1,25 IF(SW(I).NE.SWL(I,IC)) GOTO 10 IF(SWC(I).NE.SWCL(I,IC)) GOTO 10 7 CONTINUE GOTO 20 10 CONTINUE VTST7=1 IYDL(IC)=IYD SECL(IC)=SEC GLATL(IC)=GLAT GLL(IC)=GLONG STLL(IC)=STL FAL(IC)=F107A FL(IC)=F107 DO 15 I=1,7 APL(I,IC)=AP(I) 15 CONTINUE DO 16 I=1,25 SWL(I,IC)=SW(I) SWCL(I,IC)=SWC(I) 16 CONTINUE 20 CONTINUE RETURN END C----------------------------------------------------------------------- SUBROUTINE GTS7(IYD,SEC,ALT,GLAT,GLONG,STL,F107A,F107,AP,MASS,D,T) C C Thermospheric portion of NRLMSISE-00 C See GTD7 for more extensive comments C C OUTPUT IN M-3 and KG/M3: CALL METERS(.TRUE.) C C INPUT VARIABLES: C IYD - YEAR AND DAY AS YYDDD (day of year from 1 to 365 (or 366)) C (Year ignored in current model) C SEC - UT(SEC) C ALT - ALTITUDE(KM) (>72.5 km) C GLAT - GEODETIC LATITUDE(DEG) C GLONG - GEODETIC LONGITUDE(DEG) C STL - LOCAL APPARENT SOLAR TIME(HRS; see Note below) C F107A - 81 day AVERAGE OF F10.7 FLUX (centered on day DDD) C F107 - DAILY F10.7 FLUX FOR PREVIOUS DAY C AP - MAGNETIC INDEX(DAILY) OR WHEN SW(9)=-1. : C - ARRAY CONTAINING: C (1) DAILY AP C (2) 3 HR AP INDEX FOR CURRENT TIME C (3) 3 HR AP INDEX FOR 3 HRS BEFORE CURRENT TIME C (4) 3 HR AP INDEX FOR 6 HRS BEFORE CURRENT TIME C (5) 3 HR AP INDEX FOR 9 HRS BEFORE CURRENT TIME C (6) AVERAGE OF EIGHT 3 HR AP INDICIES FROM 12 TO 33 HRS PRIOR C TO CURRENT TIME C (7) AVERAGE OF EIGHT 3 HR AP INDICIES FROM 36 TO 57 HRS PRIOR C TO CURRENT TIME C MASS - MASS NUMBER (ONLY DENSITY FOR SELECTED GAS IS C CALCULATED. MASS 0 IS TEMPERATURE. MASS 48 FOR ALL. C MASS 17 IS Anomalous O ONLY.) C C NOTES ON INPUT VARIABLES: C UT, Local Time, and Longitude are used independently in the C model and are not of equal importance for every situation. C For the most physically realistic calculation these three C variables should be consistent (STL=SEC/3600+GLONG/15). C The Equation of Time departures from the above formula C for apparent local time can be included if available but C are of minor importance. c C F107 and F107A values used to generate the model correspond C to the 10.7 cm radio flux at the actual distance of the Earth C from the Sun rather than the radio flux at 1 AU. The following C site provides both classes of values: C ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SOLAR_RADIO/FLUX/ C C F107, F107A, and AP effects are neither large nor well C established below 80 km and these parameters should be set to C 150., 150., and 4. respectively. C C OUTPUT VARIABLES: C D(1) - HE NUMBER DENSITY(CM-3) C D(2) - O NUMBER DENSITY(CM-3) C D(3) - N2 NUMBER DENSITY(CM-3) C D(4) - O2 NUMBER DENSITY(CM-3) C D(5) - AR NUMBER DENSITY(CM-3) C D(6) - TOTAL MASS DENSITY(GM/CM3) [Anomalous O NOT included] C D(7) - H NUMBER DENSITY(CM-3) C D(8) - N NUMBER DENSITY(CM-3) C D(9) - Anomalous oxygen NUMBER DENSITY(CM-3) C T(1) - EXOSPHERIC TEMPERATURE C T(2) - TEMPERATURE AT ALT C DIMENSION ZN1(5),ALPHA(9) COMMON/GTS3C/TLB,S,DB04,DB16,DB28,DB32,DB40,DB48,DB01,ZA,T0,Z0 & ,G0,RL,DD,DB14,TR12 COMMON/MESO7/TN1(5),TN2(4),TN3(5),TGN1(2),TGN2(2),TGN3(2) DIMENSION D(9),T(2),MT(11),AP(1),ALTL(8) COMMON/LOWER7/PTM(10),PDM(10,8) COMMON/PARM7/PT(150),PD(150,9),PS(150),PDL(25,2),PTL(100,4), $ PMA(100,10),SAM(100) COMMON/CSW/SW(25),ISW,SWC(25) COMMON/TTEST/TINFG,GB,ROUT,TT(15) COMMON/DMIX/DM04,DM16,DM28,DM32,DM40,DM01,DM14 COMMON/METSEL/IMR SAVE DATA MT/48,0,4,16,28,32,40,1,49,14,17/ DATA ALTL/200.,300.,160.,250.,240.,450.,320.,450./ DATA MN1/5/,ZN1/120.,110.,100.,90.,72.5/ DATA DGTR/1.74533E-2/,DR/1.72142E-2/,ALAST/-999./ DATA ALPHA/-0.38,0.,0.,0.,0.17,0.,-0.38,0.,0./ C Test for changed input V2=VTST7(IYD,SEC,GLAT,GLONG,STL,F107A,F107,AP,2) C YRD=IYD ZA=PDL(16,2) ZN1(1)=ZA DO 2 J=1,9 D(J)=0. 2 CONTINUE C TINF VARIATIONS NOT IMPORTANT BELOW ZA OR ZN1(1) IF(ALT.GT.ZN1(1)) THEN IF(V2.EQ.1..OR.ALAST.LE.ZN1(1)) TINF=PTM(1)*PT(1) $ *(1.+SW(16)*GLOBE7(YRD,SEC,GLAT,GLONG,STL,F107A,F107,AP,PT)) ELSE TINF=PTM(1)*PT(1) ENDIF T(1)=TINF C GRADIENT VARIATIONS NOT IMPORTANT BELOW ZN1(5) IF(ALT.GT.ZN1(5)) THEN IF(V2.EQ.1.OR.ALAST.LE.ZN1(5)) G0=PTM(4)*PS(1) $ *(1.+SW(19)*GLOBE7(YRD,SEC,GLAT,GLONG,STL,F107A,F107,AP,PS)) ELSE G0=PTM(4)*PS(1) ENDIF C Calculate these temperatures only if input changed IF(V2.EQ.1. .OR. ALT.LT.300.) $ TLB=PTM(2)*(1.+SW(17)*GLOBE7(YRD,SEC,GLAT,GLONG,STL, $ F107A,F107,AP,PD(1,4)))*PD(1,4) S=G0/(TINF-TLB) C Lower thermosphere temp variations not significant for C density above 300 km IF(ALT.LT.300.) THEN IF(V2.EQ.1..OR.ALAST.GE.300.) THEN TN1(2)=PTM(7)*PTL(1,1)/(1.-SW(18)*GLOB7S(PTL(1,1))) TN1(3)=PTM(3)*PTL(1,2)/(1.-SW(18)*GLOB7S(PTL(1,2))) TN1(4)=PTM(8)*PTL(1,3)/(1.-SW(18)*GLOB7S(PTL(1,3))) TN1(5)=PTM(5)*PTL(1,4)/(1.-SW(18)*SW(20)*GLOB7S(PTL(1,4))) TGN1(2)=PTM(9)*PMA(1,9)*(1.+SW(18)*SW(20)*GLOB7S(PMA(1,9))) $ *TN1(5)*TN1(5)/(PTM(5)*PTL(1,4))**2 ENDIF ELSE TN1(2)=PTM(7)*PTL(1,1) TN1(3)=PTM(3)*PTL(1,2) TN1(4)=PTM(8)*PTL(1,3) TN1(5)=PTM(5)*PTL(1,4) TGN1(2)=PTM(9)*PMA(1,9) $ *TN1(5)*TN1(5)/(PTM(5)*PTL(1,4))**2 ENDIF C Z0=ZN1(4) T0=TN1(4) TR12=1. C IF(MASS.EQ.0) GO TO 50 C N2 variation factor at Zlb G28=SW(21)*GLOBE7(YRD,SEC,GLAT,GLONG,STL,F107A,F107, & AP,PD(1,3)) DAY=AMOD(YRD,1000.) C VARIATION OF TURBOPAUSE HEIGHT ZHF=PDL(25,2) $ *(1.+SW(5)*PDL(25,1)*SIN(DGTR*GLAT)*COS(DR*(DAY-PT(14)))) YRD=IYD T(1)=TINF XMM=PDM(5,3) Z=ALT C DO 10 J = 1,11 IF(MASS.EQ.MT(J)) GO TO 15 10 CONTINUE WRITE(6,100) MASS GO TO 90 15 IF(Z.GT.ALTL(6).AND.MASS.NE.28.AND.MASS.NE.48) GO TO 17 C C **** N2 DENSITY **** C C Diffusive density at Zlb DB28 = PDM(1,3)*EXP(G28)*PD(1,3) C Diffusive density at Alt D(3)=DENSU(Z,DB28,TINF,TLB, 28.,ALPHA(3),T(2),PTM(6),S,MN1,ZN1, & TN1,TGN1) DD=D(3) C Turbopause ZH28=PDM(3,3)*ZHF ZHM28=PDM(4,3)*PDL(6,2) XMD=28.-XMM C Mixed density at Zlb B28=DENSU(ZH28,DB28,TINF,TLB,XMD,ALPHA(3)-1.,TZ,PTM(6),S,MN1, & ZN1,TN1,TGN1) IF(Z.GT.ALTL(3).OR.SW(15).EQ.0.) GO TO 17 C Mixed density at Alt DM28=DENSU(Z,B28,TINF,TLB,XMM,ALPHA(3),TZ,PTM(6),S,MN1, & ZN1,TN1,TGN1) C Net density at Alt D(3)=DNET(D(3),DM28,ZHM28,XMM,28.) 17 CONTINUE GO TO (20,50,20,25,90,35,40,45,25,48,46), J 20 CONTINUE C C **** HE DENSITY **** C C Density variation factor at Zlb G4 = SW(21)*GLOBE7(YRD,SEC,GLAT,GLONG,STL,F107A,F107,AP,PD(1,1)) C Diffusive density at Zlb DB04 = PDM(1,1)*EXP(G4)*PD(1,1) C Diffusive density at Alt D(1)=DENSU(Z,DB04,TINF,TLB, 4.,ALPHA(1),T(2),PTM(6),S,MN1,ZN1, & TN1,TGN1) DD=D(1) IF(Z.GT.ALTL(1).OR.SW(15).EQ.0.) GO TO 24 C Turbopause ZH04=PDM(3,1) C Mixed density at Zlb B04=DENSU(ZH04,DB04,TINF,TLB,4.-XMM,ALPHA(1)-1., $ T(2),PTM(6),S,MN1,ZN1,TN1,TGN1) C Mixed density at Alt DM04=DENSU(Z,B04,TINF,TLB,XMM,0.,T(2),PTM(6),S,MN1,ZN1,TN1,TGN1) ZHM04=ZHM28 C Net density at Alt D(1)=DNET(D(1),DM04,ZHM04,XMM,4.) C Correction to specified mixing ratio at ground RL=ALOG(B28*PDM(2,1)/B04) ZC04=PDM(5,1)*PDL(1,2) HC04=PDM(6,1)*PDL(2,2) C Net density corrected at Alt D(1)=D(1)*CCOR(Z,RL,HC04,ZC04) 24 CONTINUE IF(MASS.NE.48) GO TO 90 25 CONTINUE C C **** O DENSITY **** C C Density variation factor at Zlb G16= SW(21)*GLOBE7(YRD,SEC,GLAT,GLONG,STL,F107A,F107,AP,PD(1,2)) C Diffusive density at Zlb DB16 = PDM(1,2)*EXP(G16)*PD(1,2) C Diffusive density at Alt D(2)=DENSU(Z,DB16,TINF,TLB, 16.,ALPHA(2),T(2),PTM(6),S,MN1, $ ZN1,TN1,TGN1) DD=D(2) IF(Z.GT.ALTL(2).OR.SW(15).EQ.0.) GO TO 34 C Corrected from PDM(3,1) to PDM(3,2) 12/2/85 C Turbopause ZH16=PDM(3,2) C Mixed density at Zlb B16=DENSU(ZH16,DB16,TINF,TLB,16-XMM,ALPHA(2)-1., $ T(2),PTM(6),S,MN1,ZN1,TN1,TGN1) C Mixed density at Alt DM16=DENSU(Z,B16,TINF,TLB,XMM,0.,T(2),PTM(6),S,MN1,ZN1,TN1,TGN1) ZHM16=ZHM28 C Net density at Alt D(2)=DNET(D(2),DM16,ZHM16,XMM,16.) C 3/16/99 Change form to match O2 departure from diff equil near 150 C km and add dependence on F10.7 C RL=ALOG(B28*PDM(2,2)*ABS(PDL(17,2))/B16) RL=PDM(2,2)*PDL(17,2)*(1.+SW(1)*PDL(24,1)*(F107A-150.)) HC16=PDM(6,2)*PDL(4,2) ZC16=PDM(5,2)*PDL(3,2) D(2)=D(2)*CCOR(Z,RL,HC16,ZC16) C Chemistry correction HCC16=PDM(8,2)*PDL(14,2) ZCC16=PDM(7,2)*PDL(13,2) RC16=PDM(4,2)*PDL(15,2) C Net density corrected at Alt D(2)=D(2)*CCOR(Z,RC16,HCC16,ZCC16) 34 CONTINUE IF(MASS.NE.48.AND.MASS.NE.49) GO TO 90 35 CONTINUE C C **** O2 DENSITY **** C C Density variation factor at Zlb G32= SW(21)*GLOBE7(YRD,SEC,GLAT,GLONG,STL,F107A,F107,AP,PD(1,5)) C Diffusive density at Zlb DB32 = PDM(1,4)*EXP(G32)*PD(1,5) C Diffusive density at Alt D(4)=DENSU(Z,DB32,TINF,TLB, 32.,ALPHA(4),T(2),PTM(6),S,MN1, $ ZN1,TN1,TGN1) IF(MASS.EQ.49) THEN DD=DD+2.*D(4) ELSE DD=D(4) ENDIF IF(SW(15).EQ.0.) GO TO 39 IF(Z.GT.ALTL(4)) GO TO 38 C Turbopause ZH32=PDM(3,4) C Mixed density at Zlb B32=DENSU(ZH32,DB32,TINF,TLB,32.-XMM,ALPHA(4)-1., $ T(2),PTM(6),S,MN1,ZN1,TN1,TGN1) C Mixed density at Alt DM32=DENSU(Z,B32,TINF,TLB,XMM,0.,T(2),PTM(6),S,MN1,ZN1,TN1,TGN1) ZHM32=ZHM28 C Net density at Alt D(4)=DNET(D(4),DM32,ZHM32,XMM,32.) C Correction to specified mixing ratio at ground RL=ALOG(B28*PDM(2,4)/B32) HC32=PDM(6,4)*PDL(8,2) ZC32=PDM(5,4)*PDL(7,2) D(4)=D(4)*CCOR(Z,RL,HC32,ZC32) 38 CONTINUE C Correction for general departure from diffusive equilibrium above Zlb HCC32=PDM(8,4)*PDL(23,2) ZCC32=PDM(7,4)*PDL(22,2) RC32=PDM(4,4)*PDL(24,2)*(1.+SW(1)*PDL(24,1)*(F107A-150.)) C Net density corrected at Alt D(4)=D(4)*CCOR(Z,RC32,HCC32,ZCC32) 39 CONTINUE IF(MASS.NE.48) GO TO 90 40 CONTINUE C C **** AR DENSITY **** C C Density variation factor at Zlb G40= SW(21)*GLOBE7(YRD,SEC,GLAT,GLONG,STL,F107A,F107,AP,PD(1,6)) C Diffusive density at Zlb DB40 = PDM(1,5)*EXP(G40)*PD(1,6) C Diffusive density at Alt D(5)=DENSU(Z,DB40,TINF,TLB, 40.,ALPHA(5),T(2),PTM(6),S,MN1, $ ZN1,TN1,TGN1) DD=D(5) IF(Z.GT.ALTL(5).OR.SW(15).EQ.0.) GO TO 44 C Turbopause ZH40=PDM(3,5) C Mixed density at Zlb B40=DENSU(ZH40,DB40,TINF,TLB,40.-XMM,ALPHA(5)-1., $ T(2),PTM(6),S,MN1,ZN1,TN1,TGN1) C Mixed density at Alt DM40=DENSU(Z,B40,TINF,TLB,XMM,0.,T(2),PTM(6),S,MN1,ZN1,TN1,TGN1) ZHM40=ZHM28 C Net density at Alt D(5)=DNET(D(5),DM40,ZHM40,XMM,40.) C Correction to specified mixing ratio at ground RL=ALOG(B28*PDM(2,5)/B40) HC40=PDM(6,5)*PDL(10,2) ZC40=PDM(5,5)*PDL(9,2) C Net density corrected at Alt D(5)=D(5)*CCOR(Z,RL,HC40,ZC40) 44 CONTINUE IF(MASS.NE.48) GO TO 90 45 CONTINUE C C **** HYDROGEN DENSITY **** C C Density variation factor at Zlb G1 = SW(21)*GLOBE7(YRD,SEC,GLAT,GLONG,STL,F107A,F107,AP,PD(1,7)) C Diffusive density at Zlb DB01 = PDM(1,6)*EXP(G1)*PD(1,7) C Diffusive density at Alt D(7)=DENSU(Z,DB01,TINF,TLB,1.,ALPHA(7),T(2),PTM(6),S,MN1, $ ZN1,TN1,TGN1) DD=D(7) IF(Z.GT.ALTL(7).OR.SW(15).EQ.0.) GO TO 47 C Turbopause ZH01=PDM(3,6) C Mixed density at Zlb B01=DENSU(ZH01,DB01,TINF,TLB,1.-XMM,ALPHA(7)-1., $ T(2),PTM(6),S,MN1,ZN1,TN1,TGN1) C Mixed density at Alt DM01=DENSU(Z,B01,TINF,TLB,XMM,0.,T(2),PTM(6),S,MN1,ZN1,TN1,TGN1) ZHM01=ZHM28 C Net density at Alt D(7)=DNET(D(7),DM01,ZHM01,XMM,1.) C Correction to specified mixing ratio at ground RL=ALOG(B28*PDM(2,6)*ABS(PDL(18,2))/B01) HC01=PDM(6,6)*PDL(12,2) ZC01=PDM(5,6)*PDL(11,2) D(7)=D(7)*CCOR(Z,RL,HC01,ZC01) C Chemistry correction HCC01=PDM(8,6)*PDL(20,2) ZCC01=PDM(7,6)*PDL(19,2) RC01=PDM(4,6)*PDL(21,2) C Net density corrected at Alt D(7)=D(7)*CCOR(Z,RC01,HCC01,ZCC01) 47 CONTINUE IF(MASS.NE.48) GO TO 90 48 CONTINUE C C **** ATOMIC NITROGEN DENSITY **** C C Density variation factor at Zlb G14 = SW(21)*GLOBE7(YRD,SEC,GLAT,GLONG,STL,F107A,F107,AP,PD(1,8)) C Diffusive density at Zlb DB14 = PDM(1,7)*EXP(G14)*PD(1,8) C Diffusive density at Alt D(8)=DENSU(Z,DB14,TINF,TLB,14.,ALPHA(8),T(2),PTM(6),S,MN1, $ ZN1,TN1,TGN1) DD=D(8) IF(Z.GT.ALTL(8).OR.SW(15).EQ.0.) GO TO 49 C Turbopause ZH14=PDM(3,7) C Mixed density at Zlb B14=DENSU(ZH14,DB14,TINF,TLB,14.-XMM,ALPHA(8)-1., $ T(2),PTM(6),S,MN1,ZN1,TN1,TGN1) C Mixed density at Alt DM14=DENSU(Z,B14,TINF,TLB,XMM,0.,T(2),PTM(6),S,MN1,ZN1,TN1,TGN1) ZHM14=ZHM28 C Net density at Alt D(8)=DNET(D(8),DM14,ZHM14,XMM,14.) C Correction to specified mixing ratio at ground RL=ALOG(B28*PDM(2,7)*ABS(PDL(3,1))/B14) HC14=PDM(6,7)*PDL(2,1) ZC14=PDM(5,7)*PDL(1,1) D(8)=D(8)*CCOR(Z,RL,HC14,ZC14) C Chemistry correction HCC14=PDM(8,7)*PDL(5,1) ZCC14=PDM(7,7)*PDL(4,1) RC14=PDM(4,7)*PDL(6,1) C Net density corrected at Alt D(8)=D(8)*CCOR(Z,RC14,HCC14,ZCC14) 49 CONTINUE IF(MASS.NE.48) GO TO 90 46 CONTINUE C C **** Anomalous OXYGEN DENSITY **** C G16H = SW(21)*GLOBE7(YRD,SEC,GLAT,GLONG,STL,F107A,F107,AP,PD(1,9)) DB16H = PDM(1,8)*EXP(G16H)*PD(1,9) THO=PDM(10,8)*PDL(7,1) DD=DENSU(Z,DB16H,THO,THO,16.,ALPHA(9),T2,PTM(6),S,MN1, $ ZN1,TN1,TGN1) ZSHT=PDM(6,8) ZMHO=PDM(5,8) ZSHO=SCALH(ZMHO,16.,THO) D(9)=DD*EXP(-ZSHT/ZSHO*(EXP(-(Z-ZMHO)/ZSHT)-1.)) IF(MASS.NE.48) GO TO 90 C C TOTAL MASS DENSITY C D(6) = 1.66E-24*(4.*D(1)+16.*D(2)+28.*D(3)+32.*D(4)+40.*D(5)+ & D(7)+14.*D(8)) DB48=1.66E-24*(4.*DB04+16.*DB16+28.*DB28+32.*DB32+40.*DB40+DB01+ & 14.*DB14) GO TO 90 C TEMPERATURE AT ALTITUDE 50 CONTINUE Z=ABS(ALT) DDUM = DENSU(Z,1., TINF,TLB,0.,0.,T(2),PTM(6),S,MN1,ZN1,TN1,TGN1) 90 CONTINUE C ADJUST DENSITIES FROM CGS TO KGM IF(IMR.EQ.1) THEN DO 95 I=1,9 D(I)=D(I)*1.E6 95 CONTINUE D(6)=D(6)/1000. ENDIF ALAST=ALT RETURN 100 FORMAT(1X,'MASS', I5, ' NOT VALID') END C----------------------------------------------------------------------- SUBROUTINE METERS(METER) C Convert outputs to Kg & Meters if METER true LOGICAL METER COMMON/METSEL/IMR SAVE IMR=0 IF(METER) IMR=1 END C----------------------------------------------------------------------- FUNCTION SCALH(ALT,XM,TEMP) C Calculate scale height (km) COMMON/PARMB/GSURF,RE SAVE DATA RGAS/831.4/ G=GSURF/(1.+ALT/RE)**2 SCALH=RGAS*TEMP/(G*XM) RETURN END C----------------------------------------------------------------------- FUNCTION GLOBE7(YRD,SEC,LAT,LONG,TLOC,F107A,F107,AP,P) C CALCULATE G(L) FUNCTION C Upper Thermosphere Parameters REAL LAT, LONG DIMENSION P(1),SV(25),AP(1) COMMON/TTEST/TINF,GB,ROUT,T(15) COMMON/CSW/SW(25),ISW,SWC(25) COMMON/LPOLY/PLG(9,4),CTLOC,STLOC,C2TLOC,S2TLOC,C3TLOC,S3TLOC, $ IYR,DAY,DF,DFA,APD,APDF,APT(4),XLONG SAVE DATA DGTR/1.74533E-2/,DR/1.72142E-2/, XL/1000./,TLL/1000./ DATA SW9/1./,DAYL/-1./,P14/-1000./,P18/-1000./,P32/-1000./ DATA HR/.2618/,SR/7.2722E-5/,SV/25*1./,NSW/14/,P39/-1000./ C 3hr Magnetic activity functions C Eq. A24d G0(A)=(A-4.+(P(26)-1.)*(A-4.+(EXP(-ABS(P(25))*(A-4.))-1.)/ABS(P(25 *)))) C Eq. A24c SUMEX(EX)=1.+(1.-EX**19)/(1.-EX)*EX**(.5) C Eq. A24a SG0(EX)=(G0(AP(2))+(G0(AP(3))*EX+G0(AP(4))*EX*EX+G0(AP(5))*EX**3 $ +(G0(AP(6))*EX**4+G0(AP(7))*EX**12)*(1.-EX**8)/(1.-EX)) $ )/SUMEX(EX) C IF(ISW.NE.64999) CALL TSELEC(SV) DO 10 J=1,14 T(J)=0 10 CONTINUE IF(SW(9).GT.0) SW9=1. IF(SW(9).LT.0) SW9=-1. IYR = YRD/1000. DAY = YRD - IYR*1000. XLONG=LONG C Eq. A22 (remainder of code) IF(XL.EQ.LAT) GO TO 15 C CALCULATE LEGENDRE POLYNOMIALS C = SIN(LAT*DGTR) S = COS(LAT*DGTR) C2 = C*C C4 = C2*C2 S2 = S*S PLG(2,1) = C PLG(3,1) = 0.5*(3.*C2 -1.) PLG(4,1) = 0.5*(5.*C*C2-3.*C) PLG(5,1) = (35.*C4 - 30.*C2 + 3.)/8. PLG(6,1) = (63.*C2*C2*C - 70.*C2*C + 15.*C)/8. PLG(7,1) = (11.*C*PLG(6,1) - 5.*PLG(5,1))/6. C PLG(8,1) = (13.*C*PLG(7,1) - 6.*PLG(6,1))/7. PLG(2,2) = S PLG(3,2) = 3.*C*S PLG(4,2) = 1.5*(5.*C2-1.)*S PLG(5,2) = 2.5*(7.*C2*C-3.*C)*S PLG(6,2) = 1.875*(21.*C4 - 14.*C2 +1.)*S PLG(7,2) = (11.*C*PLG(6,2)-6.*PLG(5,2))/5. C PLG(8,2) = (13.*C*PLG(7,2)-7.*PLG(6,2))/6. C PLG(9,2) = (15.*C*PLG(8,2)-8.*PLG(7,2))/7. PLG(3,3) = 3.*S2 PLG(4,3) = 15.*S2*C PLG(5,3) = 7.5*(7.*C2 -1.)*S2 PLG(6,3) = 3.*C*PLG(5,3)-2.*PLG(4,3) PLG(7,3)=(11.*C*PLG(6,3)-7.*PLG(5,3))/4. PLG(8,3)=(13.*C*PLG(7,3)-8.*PLG(6,3))/5. PLG(4,4) = 15.*S2*S PLG(5,4) = 105.*S2*S*C PLG(6,4)=(9.*C*PLG(5,4)-7.*PLG(4,4))/2. PLG(7,4)=(11.*C*PLG(6,4)-8.*PLG(5,4))/3. XL=LAT 15 CONTINUE IF(TLL.EQ.TLOC) GO TO 16 IF(SW(7).EQ.0.AND.SW(8).EQ.0.AND.SW(14).EQ.0) GOTO 16 STLOC = SIN(HR*TLOC) CTLOC = COS(HR*TLOC) S2TLOC = SIN(2.*HR*TLOC) C2TLOC = COS(2.*HR*TLOC) S3TLOC = SIN(3.*HR*TLOC) C3TLOC = COS(3.*HR*TLOC) TLL = TLOC 16 CONTINUE IF(DAY.NE.DAYL.OR.P(14).NE.P14) CD14=COS(DR*(DAY-P(14))) IF(DAY.NE.DAYL.OR.P(18).NE.P18) CD18=COS(2.*DR*(DAY-P(18))) IF(DAY.NE.DAYL.OR.P(32).NE.P32) CD32=COS(DR*(DAY-P(32))) IF(DAY.NE.DAYL.OR.P(39).NE.P39) CD39=COS(2.*DR*(DAY-P(39))) DAYL = DAY P14 = P(14) P18 = P(18) P32 = P(32) P39 = P(39) C F10.7 EFFECT DF = F107 - F107A DFA=F107A-150. T(1) = P(20)*DF*(1.+P(60)*DFA) + P(21)*DF*DF + P(22)*DFA & + P(30)*DFA**2 F1 = 1. + (P(48)*DFA +P(20)*DF+P(21)*DF*DF)*SWC(1) F2 = 1. + (P(50)*DFA+P(20)*DF+P(21)*DF*DF)*SWC(1) C TIME INDEPENDENT T(2) = 1 (P(2)*PLG(3,1) + P(3)*PLG(5,1)+P(23)*PLG(7,1)) & +(P(15)*PLG(3,1))*DFA*SWC(1) 2 +P(27)*PLG(2,1) C SYMMETRICAL ANNUAL T(3) = 1 (P(19) )*CD32 C SYMMETRICAL SEMIANNUAL T(4) = 1 (P(16)+P(17)*PLG(3,1))*CD18 C ASYMMETRICAL ANNUAL T(5) = F1* 1 (P(10)*PLG(2,1)+P(11)*PLG(4,1))*CD14 C ASYMMETRICAL SEMIANNUAL T(6) = P(38)*PLG(2,1)*CD39 C DIURNAL IF(SW(7).EQ.0) GOTO 200 T71 = (P(12)*PLG(3,2))*CD14*SWC(5) T72 = (P(13)*PLG(3,2))*CD14*SWC(5) T(7) = F2* 1 ((P(4)*PLG(2,2) + P(5)*PLG(4,2) + P(28)*PLG(6,2) 2 + T71)*CTLOC 4 + (P(7)*PLG(2,2) + P(8)*PLG(4,2) +P(29)*PLG(6,2) 5 + T72)*STLOC) 200 CONTINUE C SEMIDIURNAL IF(SW(8).EQ.0) GOTO 210 T81 = (P(24)*PLG(4,3)+P(36)*PLG(6,3))*CD14*SWC(5) T82 = (P(34)*PLG(4,3)+P(37)*PLG(6,3))*CD14*SWC(5) T(8) = F2* 1 ((P(6)*PLG(3,3) + P(42)*PLG(5,3) + T81)*C2TLOC 3 +(P(9)*PLG(3,3) + P(43)*PLG(5,3) + T82)*S2TLOC) 210 CONTINUE C TERDIURNAL IF(SW(14).EQ.0) GOTO 220 T(14) = F2* 1 ((P(40)*PLG(4,4)+(P(94)*PLG(5,4)+P(47)*PLG(7,4))*CD14*SWC(5))* $ S3TLOC 2 +(P(41)*PLG(4,4)+(P(95)*PLG(5,4)+P(49)*PLG(7,4))*CD14*SWC(5))* $ C3TLOC) 220 CONTINUE C MAGNETIC ACTIVITY BASED ON DAILY AP IF(SW9.EQ.-1.) GO TO 30 APD=(AP(1)-4.) P44=P(44) P45=P(45) IF(P44.LT.0) P44=1.E-5 APDF = APD+(P45-1.)*(APD+(EXP(-P44 *APD)-1.)/P44) IF(SW(9).EQ.0) GOTO 40 T(9)=APDF*(P(33)+P(46)*PLG(3,1)+P(35)*PLG(5,1)+ $ (P(101)*PLG(2,1)+P(102)*PLG(4,1)+P(103)*PLG(6,1))*CD14*SWC(5)+ $ (P(122)*PLG(2,2)+P(123)*PLG(4,2)+P(124)*PLG(6,2))*SWC(7)* $ COS(HR*(TLOC-P(125)))) GO TO 40 30 CONTINUE IF(P(52).EQ.0) GO TO 40 EXP1 = EXP(-10800.*ABS(P(52))/(1.+P(139)*(45.-ABS(LAT)))) IF(EXP1.GT..99999) EXP1=.99999 IF(P(25).LT.1.E-4) P(25)=1.E-4 APT(1)=SG0(EXP1) C APT(2)=SG2(EXP1) c APT(3)=SG0(EXP2) C APT(4)=SG2(EXP2) IF(SW(9).EQ.0) GOTO 40 T(9) = APT(1)*(P(51)+P(97)*PLG(3,1)+P(55)*PLG(5,1)+ $ (P(126)*PLG(2,1)+P(127)*PLG(4,1)+P(128)*PLG(6,1))*CD14*SWC(5)+ $ (P(129)*PLG(2,2)+P(130)*PLG(4,2)+P(131)*PLG(6,2))*SWC(7)* $ COS(HR*(TLOC-P(132)))) 40 CONTINUE IF(SW(10).EQ.0.OR.LONG.LE.-1000.) GO TO 49 C LONGITUDINAL IF(SW(11).EQ.0) GOTO 230 T(11)= (1.+P(81)*DFA*SWC(1))* $((P(65)*PLG(3,2)+P(66)*PLG(5,2)+P(67)*PLG(7,2) $ +P(104)*PLG(2,2)+P(105)*PLG(4,2)+P(106)*PLG(6,2) $ +SWC(5)*(P(110)*PLG(2,2)+P(111)*PLG(4,2)+P(112)*PLG(6,2))*CD14)* $ COS(DGTR*LONG) $ +(P(91)*PLG(3,2)+P(92)*PLG(5,2)+P(93)*PLG(7,2) $ +P(107)*PLG(2,2)+P(108)*PLG(4,2)+P(109)*PLG(6,2) $ +SWC(5)*(P(113)*PLG(2,2)+P(114)*PLG(4,2)+P(115)*PLG(6,2))*CD14)* $ SIN(DGTR*LONG)) 230 CONTINUE C UT AND MIXED UT,LONGITUDE IF(SW(12).EQ.0) GOTO 240 T(12)=(1.+P(96)*PLG(2,1))*(1.+P(82)*DFA*SWC(1))* $(1.+P(120)*PLG(2,1)*SWC(5)*CD14)* $((P(69)*PLG(2,1)+P(70)*PLG(4,1)+P(71)*PLG(6,1))* $ COS(SR*(SEC-P(72)))) T(12)=T(12)+SWC(11)* $ (P(77)*PLG(4,3)+P(78)*PLG(6,3)+P(79)*PLG(8,3))* $ COS(SR*(SEC-P(80))+2.*DGTR*LONG)*(1.+P(138)*DFA*SWC(1)) 240 CONTINUE C UT,LONGITUDE MAGNETIC ACTIVITY IF(SW(13).EQ.0) GOTO 48 IF(SW9.EQ.-1.) GO TO 45 T(13)= APDF*SWC(11)*(1.+P(121)*PLG(2,1))* $((P( 61)*PLG(3,2)+P( 62)*PLG(5,2)+P( 63)*PLG(7,2))* $ COS(DGTR*(LONG-P( 64)))) $ +APDF*SWC(11)*SWC(5)* $ (P(116)*PLG(2,2)+P(117)*PLG(4,2)+P(118)*PLG(6,2))* $ CD14*COS(DGTR*(LONG-P(119))) $ + APDF*SWC(12)* $ (P( 84)*PLG(2,1)+P( 85)*PLG(4,1)+P( 86)*PLG(6,1))* $ COS(SR*(SEC-P( 76))) GOTO 48 45 CONTINUE IF(P(52).EQ.0) GOTO 48 T(13)=APT(1)*SWC(11)*(1.+P(133)*PLG(2,1))* $((P(53)*PLG(3,2)+P(99)*PLG(5,2)+P(68)*PLG(7,2))* $ COS(DGTR*(LONG-P(98)))) $ +APT(1)*SWC(11)*SWC(5)* $ (P(134)*PLG(2,2)+P(135)*PLG(4,2)+P(136)*PLG(6,2))* $ CD14*COS(DGTR*(LONG-P(137))) $ +APT(1)*SWC(12)* $ (P(56)*PLG(2,1)+P(57)*PLG(4,1)+P(58)*PLG(6,1))* $ COS(SR*(SEC-P(59))) 48 CONTINUE C PARMS NOT USED: 83, 90,100,140-150 49 CONTINUE TINF=P(31) DO 50 I = 1,NSW 50 TINF = TINF + ABS(SW(I))*T(I) GLOBE7 = TINF RETURN END C----------------------------------------------------------------------- SUBROUTINE TSELEC(SV) C SET SWITCHES C Output in COMMON/CSW/SW(25),ISW,SWC(25) C SW FOR MAIN TERMS, SWC FOR CROSS TERMS C C TO TURN ON AND OFF PARTICULAR VARIATIONS CALL TSELEC(SV), C WHERE SV IS A 25 ELEMENT ARRAY CONTAINING 0. FOR OFF, 1. C FOR ON, OR 2. FOR MAIN EFFECTS OFF BUT CROSS TERMS ON C C To get current values of SW: CALL TRETRV(SW) C DIMENSION SV(1),SAV(25),SVV(1) COMMON/CSW/SW(25),ISW,SWC(25) SAVE DO 100 I = 1,25 SAV(I)=SV(I) SW(I)=AMOD(SV(I),2.) IF(ABS(SV(I)).EQ.1.OR.ABS(SV(I)).EQ.2.) THEN SWC(I)=1. ELSE SWC(I)=0. ENDIF 100 CONTINUE ISW=64999 RETURN ENTRY TRETRV(SVV) DO 200 I=1,25 SVV(I)=SAV(I) 200 CONTINUE END C----------------------------------------------------------------------- FUNCTION GLOB7S(P) C VERSION OF GLOBE FOR LOWER ATMOSPHERE 10/26/99 REAL LONG COMMON/LPOLY/PLG(9,4),CTLOC,STLOC,C2TLOC,S2TLOC,C3TLOC,S3TLOC, $ IYR,DAY,DF,DFA,APD,APDF,APT(4),LONG COMMON/CSW/SW(25),ISW,SWC(25) DIMENSION P(1),T(14) SAVE DATA DR/1.72142E-2/,DGTR/1.74533E-2/,PSET/2./ DATA DAYL/-1./,P32,P18,P14,P39/4*-1000./ C CONFIRM PARAMETER SET IF(P(100).EQ.0) P(100)=PSET IF(P(100).NE.PSET) THEN WRITE(6,900) PSET,P(100) 900 FORMAT(1X,'WRONG PARAMETER SET FOR GLOB7S',3F10.1) STOP ENDIF DO 10 J=1,14 T(J)=0. 10 CONTINUE IF(DAY.NE.DAYL.OR.P32.NE.P(32)) CD32=COS(DR*(DAY-P(32))) IF(DAY.NE.DAYL.OR.P18.NE.P(18)) CD18=COS(2.*DR*(DAY-P(18))) IF(DAY.NE.DAYL.OR.P14.NE.P(14)) CD14=COS(DR*(DAY-P(14))) IF(DAY.NE.DAYL.OR.P39.NE.P(39)) CD39=COS(2.*DR*(DAY-P(39))) DAYL=DAY P32=P(32) P18=P(18) P14=P(14) P39=P(39) C C F10.7 T(1)=P(22)*DFA C TIME INDEPENDENT T(2)=P(2)*PLG(3,1)+P(3)*PLG(5,1)+P(23)*PLG(7,1) $ +P(27)*PLG(2,1)+P(15)*PLG(4,1)+P(60)*PLG(6,1) C SYMMETRICAL ANNUAL T(3)=(P(19)+P(48)*PLG(3,1)+P(30)*PLG(5,1))*CD32 C SYMMETRICAL SEMIANNUAL T(4)=(P(16)+P(17)*PLG(3,1)+P(31)*PLG(5,1))*CD18 C ASYMMETRICAL ANNUAL T(5)=(P(10)*PLG(2,1)+P(11)*PLG(4,1)+P(21)*PLG(6,1))*CD14 C ASYMMETRICAL SEMIANNUAL T(6)=(P(38)*PLG(2,1))*CD39 C DIURNAL IF(SW(7).EQ.0) GOTO 200 T71 = P(12)*PLG(3,2)*CD14*SWC(5) T72 = P(13)*PLG(3,2)*CD14*SWC(5) T(7) = 1 ((P(4)*PLG(2,2) + P(5)*PLG(4,2) 2 + T71)*CTLOC 4 + (P(7)*PLG(2,2) + P(8)*PLG(4,2) 5 + T72)*STLOC) 200 CONTINUE C SEMIDIURNAL IF(SW(8).EQ.0) GOTO 210 T81 = (P(24)*PLG(4,3)+P(36)*PLG(6,3))*CD14*SWC(5) T82 = (P(34)*PLG(4,3)+P(37)*PLG(6,3))*CD14*SWC(5) T(8) = 1 ((P(6)*PLG(3,3) + P(42)*PLG(5,3) + T81)*C2TLOC 3 +(P(9)*PLG(3,3) + P(43)*PLG(5,3) + T82)*S2TLOC) 210 CONTINUE C TERDIURNAL IF(SW(14).EQ.0) GOTO 220 T(14) = P(40)*PLG(4,4)*S3TLOC $ +P(41)*PLG(4,4)*C3TLOC 220 CONTINUE C MAGNETIC ACTIVITY IF(SW(9).EQ.0) GOTO 40 IF(SW(9).EQ.1) $ T(9)=APDF*(P(33)+P(46)*PLG(3,1)*SWC(2)) IF(SW(9).EQ.-1) $ T(9)=(P(51)*APT(1)+P(97)*PLG(3,1)*APT(1)*SWC(2)) 40 CONTINUE IF(SW(10).EQ.0.OR.SW(11).EQ.0.OR.LONG.LE.-1000.) GO TO 49 C LONGITUDINAL T(11)= (1.+PLG(2,1)*(P(81)*SWC(5)*COS(DR*(DAY-P(82))) $ +P(86)*SWC(6)*COS(2.*DR*(DAY-P(87)))) $ +P(84)*SWC(3)*COS(DR*(DAY-P(85))) $ +P(88)*SWC(4)*COS(2.*DR*(DAY-P(89)))) $ *((P(65)*PLG(3,2)+P(66)*PLG(5,2)+P(67)*PLG(7,2) $ +P(75)*PLG(2,2)+P(76)*PLG(4,2)+P(77)*PLG(6,2) $ )*COS(DGTR*LONG) $ +(P(91)*PLG(3,2)+P(92)*PLG(5,2)+P(93)*PLG(7,2) $ +P(78)*PLG(2,2)+P(79)*PLG(4,2)+P(80)*PLG(6,2) $ )*SIN(DGTR*LONG)) 49 CONTINUE TT=0. DO 50 I=1,14 50 TT=TT+ABS(SW(I))*T(I) GLOB7S=TT RETURN END C-------------------------------------------------------------------- FUNCTION DENSU(ALT,DLB,TINF,TLB,XM,ALPHA,TZ,ZLB,S2, $ MN1,ZN1,TN1,TGN1) C Calculate Temperature and Density Profiles for MSIS models C New lower thermo polynomial 10/30/89 DIMENSION ZN1(MN1),TN1(MN1),TGN1(2),XS(5),YS(5),Y2OUT(5) COMMON/PARMB/GSURF,RE COMMON/LSQV/MP,II,JG,LT,QPB(50),IERR,IFUN,N,J,DV(60) SAVE DATA RGAS/831.4/ ZETA(ZZ,ZL)=(ZZ-ZL)*(RE+ZL)/(RE+ZZ) CCCCCCWRITE(6,*) 'DB',ALT,DLB,TINF,TLB,XM,ALPHA,ZLB,S2,MN1,ZN1,TN1 DENSU=1. C Joining altitude of Bates and spline ZA=ZN1(1) Z=AMAX1(ALT,ZA) C Geopotential altitude difference from ZLB ZG2=ZETA(Z,ZLB) C Bates temperature TT=TINF-(TINF-TLB)*EXP(-S2*ZG2) TA=TT TZ=TT DENSU=TZ IF(ALT.GE.ZA) GO TO 10 C C CALCULATE TEMPERATURE BELOW ZA C Temperature gradient at ZA from Bates profile DTA=(TINF-TA)*S2*((RE+ZLB)/(RE+ZA))**2 TGN1(1)=DTA TN1(1)=TA Z=AMAX1(ALT,ZN1(MN1)) MN=MN1 Z1=ZN1(1) Z2=ZN1(MN) T1=TN1(1) T2=TN1(MN) C Geopotental difference from Z1 ZG=ZETA(Z,Z1) ZGDIF=ZETA(Z2,Z1) C Set up spline nodes DO 20 K=1,MN XS(K)=ZETA(ZN1(K),Z1)/ZGDIF YS(K)=1./TN1(K) 20 CONTINUE C End node derivatives YD1=-TGN1(1)/(T1*T1)*ZGDIF YD2=-TGN1(2)/(T2*T2)*ZGDIF*((RE+Z2)/(RE+Z1))**2 C Calculate spline coefficients CALL SPLINE(XS,YS,MN,YD1,YD2,Y2OUT) X=ZG/ZGDIF CALL SPLINT(XS,YS,Y2OUT,MN,X,Y) C temperature at altitude TZ=1./Y DENSU=TZ 10 IF(XM.EQ.0.) GO TO 50 C C CALCULATE DENSITY ABOVE ZA GLB=GSURF/(1.+ZLB/RE)**2 GAMMA=XM*GLB/(S2*RGAS*TINF) EXPL=EXP(-S2*GAMMA*ZG2) IF(EXPL.GT.50.OR.TT.LE.0.) THEN EXPL=50. ENDIF C Density at altitude DENSA=DLB*(TLB/TT)**(1.+ALPHA+GAMMA)*EXPL DENSU=DENSA IF(ALT.GE.ZA) GO TO 50 C C CALCULATE DENSITY BELOW ZA GLB=GSURF/(1.+Z1/RE)**2 GAMM=XM*GLB*ZGDIF/RGAS C integrate spline temperatures CALL SPLINI(XS,YS,Y2OUT,MN,X,YI) EXPL=GAMM*YI IF(EXPL.GT.50..OR.TZ.LE.0.) THEN EXPL=50. ENDIF C Density at altitude DENSU=DENSU*(T1/TZ)**(1.+ALPHA)*EXP(-EXPL) 50 CONTINUE RETURN END C-------------------------------------------------------------------- FUNCTION DENSM(ALT,D0,XM,TZ,MN3,ZN3,TN3,TGN3,MN2,ZN2,TN2,TGN2) C Calculate Temperature and Density Profiles for lower atmos. DIMENSION ZN3(MN3),TN3(MN3),TGN3(2),XS(10),YS(10),Y2OUT(10) DIMENSION ZN2(MN2),TN2(MN2),TGN2(2) COMMON/PARMB/GSURF,RE COMMON/FIT/TAF COMMON/LSQV/MP,II,JG,LT,QPB(50),IERR,IFUN,N,J,DV(60) SAVE DATA RGAS/831.4/ ZETA(ZZ,ZL)=(ZZ-ZL)*(RE+ZL)/(RE+ZZ) DENSM=D0 IF(ALT.GT.ZN2(1)) GOTO 50 C STRATOSPHERE/MESOSPHERE TEMPERATURE Z=AMAX1(ALT,ZN2(MN2)) MN=MN2 Z1=ZN2(1) Z2=ZN2(MN) T1=TN2(1) T2=TN2(MN) ZG=ZETA(Z,Z1) ZGDIF=ZETA(Z2,Z1) C Set up spline nodes DO 210 K=1,MN XS(K)=ZETA(ZN2(K),Z1)/ZGDIF YS(K)=1./TN2(K) 210 CONTINUE YD1=-TGN2(1)/(T1*T1)*ZGDIF YD2=-TGN2(2)/(T2*T2)*ZGDIF*((RE+Z2)/(RE+Z1))**2 C Calculate spline coefficients CALL SPLINE(XS,YS,MN,YD1,YD2,Y2OUT) X=ZG/ZGDIF CALL SPLINT(XS,YS,Y2OUT,MN,X,Y) C Temperature at altitude TZ=1./Y IF(XM.EQ.0.) GO TO 20 C C CALCULATE STRATOSPHERE/MESOSPHERE DENSITY GLB=GSURF/(1.+Z1/RE)**2 GAMM=XM*GLB*ZGDIF/RGAS C Integrate temperature profile CALL SPLINI(XS,YS,Y2OUT,MN,X,YI) EXPL=GAMM*YI IF(EXPL.GT.50.) EXPL=50. C Density at altitude DENSM=DENSM*(T1/TZ)*EXP(-EXPL) 20 CONTINUE IF(ALT.GT.ZN3(1)) GOTO 50 C C TROPOSPHERE/STRATOSPHERE TEMPERATURE Z=ALT MN=MN3 Z1=ZN3(1) Z2=ZN3(MN) T1=TN3(1) T2=TN3(MN) ZG=ZETA(Z,Z1) ZGDIF=ZETA(Z2,Z1) C Set up spline nodes DO 220 K=1,MN XS(K)=ZETA(ZN3(K),Z1)/ZGDIF YS(K)=1./TN3(K) 220 CONTINUE YD1=-TGN3(1)/(T1*T1)*ZGDIF YD2=-TGN3(2)/(T2*T2)*ZGDIF*((RE+Z2)/(RE+Z1))**2 C Calculate spline coefficients CALL SPLINE(XS,YS,MN,YD1,YD2,Y2OUT) X=ZG/ZGDIF CALL SPLINT(XS,YS,Y2OUT,MN,X,Y) C temperature at altitude TZ=1./Y IF(XM.EQ.0.) GO TO 30 C C CALCULATE TROPOSPHERIC/STRATOSPHERE DENSITY C GLB=GSURF/(1.+Z1/RE)**2 GAMM=XM*GLB*ZGDIF/RGAS C Integrate temperature profile CALL SPLINI(XS,YS,Y2OUT,MN,X,YI) EXPL=GAMM*YI IF(EXPL.GT.50.) EXPL=50. C Density at altitude DENSM=DENSM*(T1/TZ)*EXP(-EXPL) 30 CONTINUE 50 CONTINUE IF(XM.EQ.0) DENSM=TZ RETURN END C----------------------------------------------------------------------- SUBROUTINE SPLINE(X,Y,N,YP1,YPN,Y2) C CALCULATE 2ND DERIVATIVES OF CUBIC SPLINE INTERP FUNCTION C ADAPTED FROM NUMERICAL RECIPES BY PRESS ET AL C X,Y: ARRAYS OF TABULATED FUNCTION IN ASCENDING ORDER BY X C N: SIZE OF ARRAYS X,Y C YP1,YPN: SPECIFIED DERIVATIVES AT X(1) AND X(N); VALUES C >= 1E30 SIGNAL SIGNAL SECOND DERIVATIVE ZERO C Y2: OUTPUT ARRAY OF SECOND DERIVATIVES PARAMETER (NMAX=100) DIMENSION X(N),Y(N),Y2(N),U(NMAX) SAVE IF(YP1.GT..99E30) THEN Y2(1)=0 U(1)=0 ELSE Y2(1)=-.5 U(1)=(3./(X(2)-X(1)))*((Y(2)-Y(1))/(X(2)-X(1))-YP1) ENDIF DO 11 I=2,N-1 SIG=(X(I)-X(I-1))/(X(I+1)-X(I-1)) P=SIG*Y2(I-1)+2. Y2(I)=(SIG-1.)/P U(I)=(6.*((Y(I+1)-Y(I))/(X(I+1)-X(I))-(Y(I)-Y(I-1)) $ /(X(I)-X(I-1)))/(X(I+1)-X(I-1))-SIG*U(I-1))/P 11 CONTINUE IF(YPN.GT..99E30) THEN QN=0 UN=0 ELSE QN=.5 UN=(3./(X(N)-X(N-1)))*(YPN-(Y(N)-Y(N-1))/(X(N)-X(N-1))) ENDIF Y2(N)=(UN-QN*U(N-1))/(QN*Y2(N-1)+1.) DO 12 K=N-1,1,-1 Y2(K)=Y2(K)*Y2(K+1)+U(K) 12 CONTINUE RETURN END C----------------------------------------------------------------------- SUBROUTINE SPLINT(XA,YA,Y2A,N,X,Y) C CALCULATE CUBIC SPLINE INTERP VALUE C ADAPTED FROM NUMERICAL RECIPES BY PRESS ET AL. C XA,YA: ARRAYS OF TABULATED FUNCTION IN ASCENDING ORDER BY X C Y2A: ARRAY OF SECOND DERIVATIVES C N: SIZE OF ARRAYS XA,YA,Y2A C X: ABSCISSA FOR INTERPOLATION C Y: OUTPUT VALUE DIMENSION XA(N),YA(N),Y2A(N) SAVE KLO=1 KHI=N 1 CONTINUE IF(KHI-KLO.GT.1) THEN K=(KHI+KLO)/2 IF(XA(K).GT.X) THEN KHI=K ELSE KLO=K ENDIF GOTO 1 ENDIF H=XA(KHI)-XA(KLO) IF(H.EQ.0) WRITE(6,*) 'BAD XA INPUT TO SPLINT' A=(XA(KHI)-X)/H B=(X-XA(KLO))/H Y=A*YA(KLO)+B*YA(KHI)+ $ ((A*A*A-A)*Y2A(KLO)+(B*B*B-B)*Y2A(KHI))*H*H/6. RETURN END C----------------------------------------------------------------------- SUBROUTINE SPLINI(XA,YA,Y2A,N,X,YI) C INTEGRATE CUBIC SPLINE FUNCTION FROM XA(1) TO X C XA,YA: ARRAYS OF TABULATED FUNCTION IN ASCENDING ORDER BY X C Y2A: ARRAY OF SECOND DERIVATIVES C N: SIZE OF ARRAYS XA,YA,Y2A C X: ABSCISSA ENDPOINT FOR INTEGRATION C Y: OUTPUT VALUE DIMENSION XA(N),YA(N),Y2A(N) SAVE YI=0 KLO=1 KHI=2 1 CONTINUE IF(X.GT.XA(KLO).AND.KHI.LE.N) THEN XX=X IF(KHI.LT.N) XX=AMIN1(X,XA(KHI)) H=XA(KHI)-XA(KLO) A=(XA(KHI)-XX)/H B=(XX-XA(KLO))/H A2=A*A B2=B*B YI=YI+((1.-A2)*YA(KLO)/2.+B2*YA(KHI)/2.+ $ ((-(1.+A2*A2)/4.+A2/2.)*Y2A(KLO)+ $ (B2*B2/4.-B2/2.)*Y2A(KHI))*H*H/6.)*H KLO=KLO+1 KHI=KHI+1 GOTO 1 ENDIF RETURN END C----------------------------------------------------------------------- FUNCTION DNET(DD,DM,ZHM,XMM,XM) C TURBOPAUSE CORRECTION FOR MSIS MODELS C Root mean density C 8/20/80 C DD - diffusive density C DM - full mixed density C ZHM - transition scale length C XMM - full mixed molecular weight C XM - species molecular weight C DNET - combined density SAVE A=ZHM/(XMM-XM) IF(DM.GT.0.AND.DD.GT.0) GOTO 5 WRITE(6,*) 'DNET LOG ERROR',DM,DD,XM IF(DD.EQ.0.AND.DM.EQ.0) DD=1. IF(DM.EQ.0) GOTO 10 IF(DD.EQ.0) GOTO 20 5 CONTINUE YLOG=A*ALOG(DM/DD) IF(YLOG.LT.-10.) GO TO 10 IF(YLOG.GT.10.) GO TO 20 DNET=DD*(1.+EXP(YLOG))**(1/A) GO TO 50 10 CONTINUE DNET=DD GO TO 50 20 CONTINUE DNET=DM GO TO 50 50 CONTINUE RETURN END C----------------------------------------------------------------------- FUNCTION CCOR(ALT, R,H1,ZH) C CHEMISTRY/DISSOCIATION CORRECTION FOR MSIS MODELS C ALT - altitude C R - target ratio C H1 - transition scale length C ZH - altitude of 1/2 R SAVE E=(ALT-ZH)/H1 IF(E.GT.70.) GO TO 20 IF(E.LT.-70.) GO TO 10 EX=EXP(E) CCOR=R/(1.+EX) GO TO 50 10 CCOR=R GO TO 50 20 CCOR=0. GO TO 50 50 CONTINUE CCOR=EXP(CCOR) RETURN END C----------------------------------------------------------------------- BLOCK DATA GTD7BK C NRLMSISE-00 13-APR-00 COMMON/PARM7/PT1(50),PT2(50),PT3(50),PA1(50),PA2(50),PA3(50), $ PB1(50),PB2(50),PB3(50),PC1(50),PC2(50),PC3(50), $ PD1(50),PD2(50),PD3(50),PE1(50),PE2(50),PE3(50), $ PF1(50),PF2(50),PF3(50),PG1(50),PG2(50),PG3(50), $ PH1(50),PH2(50),PH3(50),PI1(50),PI2(50),PI3(50), $ PJ1(50),PJ2(50),PJ3(50),PK1(50),PL1(50),PL2(50), $ PM1(50),PM2(50),PN1(50),PN2(50),PO1(50),PO2(50), $ PP1(50),PP2(50),PQ1(50),PQ2(50),PR1(50),PR2(50), $ PS1(50),PS2(50),PU1(50),PU2(50),PV1(50),PV2(50), $ PW1(50),PW2(50),PX1(50),PX2(50),PY1(50),PY2(50), $ PZ1(50),PZ2(50),PAA1(50),PAA2(50) COMMON/LOWER7/PTM(10),PDM(10,8) COMMON/MAVG7/PAVGM(10) COMMON/DATIM7/ISDATE(3),ISTIME(2),NAME(2) COMMON/METSEL/IMR DATA IMR/0/ DATA ISDATE/'13-A','PR-0','0 '/,ISTIME/'17:4','6:08'/ DATA NAME/'MSIS','E-00'/ C TEMPERATURE DATA PT1/ * 9.86573E-01, 1.62228E-02, 1.55270E-02,-1.04323E-01,-3.75801E-03, * -1.18538E-03,-1.24043E-01, 4.56820E-03, 8.76018E-03,-1.36235E-01, * -3.52427E-02, 8.84181E-03,-5.92127E-03,-8.61650E+00, 0.00000E+00, * 1.28492E-02, 0.00000E+00, 1.30096E+02, 1.04567E-02, 1.65686E-03, * -5.53887E-06, 2.97810E-03, 0.00000E+00, 5.13122E-03, 8.66784E-02, * 1.58727E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00,-7.27026E-06, * 0.00000E+00, 6.74494E+00, 4.93933E-03, 2.21656E-03, 2.50802E-03, * 0.00000E+00, 0.00000E+00,-2.08841E-02,-1.79873E+00, 1.45103E-03, * 2.81769E-04,-1.44703E-03,-5.16394E-05, 8.47001E-02, 1.70147E-01, * 5.72562E-03, 5.07493E-05, 4.36148E-03, 1.17863E-04, 4.74364E-03/ DATA PT2/ * 6.61278E-03, 4.34292E-05, 1.44373E-03, 2.41470E-05, 2.84426E-03, * 8.56560E-04, 2.04028E-03, 0.00000E+00,-3.15994E+03,-2.46423E-03, * 1.13843E-03, 4.20512E-04, 0.00000E+00,-9.77214E+01, 6.77794E-03, * 5.27499E-03, 1.14936E-03, 0.00000E+00,-6.61311E-03,-1.84255E-02, * -1.96259E-02, 2.98618E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 6.44574E+02, 8.84668E-04, 5.05066E-04, 0.00000E+00, 4.02881E+03, * -1.89503E-03, 0.00000E+00, 0.00000E+00, 8.21407E-04, 2.06780E-03, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * -1.20410E-02,-3.63963E-03, 9.92070E-05,-1.15284E-04,-6.33059E-05, * -6.05545E-01, 8.34218E-03,-9.13036E+01, 3.71042E-04, 0.00000E+00/ DATA PT3/ * 4.19000E-04, 2.70928E-03, 3.31507E-03,-4.44508E-03,-4.96334E-03, * -1.60449E-03, 3.95119E-03, 2.48924E-03, 5.09815E-04, 4.05302E-03, * 2.24076E-03, 0.00000E+00, 6.84256E-03, 4.66354E-04, 0.00000E+00, * -3.68328E-04, 0.00000E+00, 0.00000E+00,-1.46870E+02, 0.00000E+00, * 0.00000E+00, 1.09501E-03, 4.65156E-04, 5.62583E-04, 3.21596E+00, * 6.43168E-04, 3.14860E-03, 3.40738E-03, 1.78481E-03, 9.62532E-04, * 5.58171E-04, 3.43731E+00,-2.33195E-01, 5.10289E-04, 0.00000E+00, * 0.00000E+00,-9.25347E+04, 0.00000E+00,-1.99639E-03, 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/ C HE DENSITY DATA PA1/ * 1.09979E+00,-4.88060E-02,-1.97501E-01,-9.10280E-02,-6.96558E-03, * 2.42136E-02, 3.91333E-01,-7.20068E-03,-3.22718E-02, 1.41508E+00, * 1.68194E-01, 1.85282E-02, 1.09384E-01,-7.24282E+00, 0.00000E+00, * 2.96377E-01,-4.97210E-02, 1.04114E+02,-8.61108E-02,-7.29177E-04, * 1.48998E-06, 1.08629E-03, 0.00000E+00, 0.00000E+00, 8.31090E-02, * 1.12818E-01,-5.75005E-02,-1.29919E-02,-1.78849E-02,-2.86343E-06, * 0.00000E+00,-1.51187E+02,-6.65902E-03, 0.00000E+00,-2.02069E-03, * 0.00000E+00, 0.00000E+00, 4.32264E-02,-2.80444E+01,-3.26789E-03, * 2.47461E-03, 0.00000E+00, 0.00000E+00, 9.82100E-02, 1.22714E-01, * -3.96450E-02, 0.00000E+00,-2.76489E-03, 0.00000E+00, 1.87723E-03/ DATA PA2/ * -8.09813E-03, 4.34428E-05,-7.70932E-03, 0.00000E+00,-2.28894E-03, * -5.69070E-03,-5.22193E-03, 6.00692E-03,-7.80434E+03,-3.48336E-03, * -6.38362E-03,-1.82190E-03, 0.00000E+00,-7.58976E+01,-2.17875E-02, * -1.72524E-02,-9.06287E-03, 0.00000E+00, 2.44725E-02, 8.66040E-02, * 1.05712E-01, 3.02543E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, * -6.01364E+03,-5.64668E-03,-2.54157E-03, 0.00000E+00, 3.15611E+02, * -5.69158E-03, 0.00000E+00, 0.00000E+00,-4.47216E-03,-4.49523E-03, * 4.64428E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 4.51236E-02, 2.46520E-02, 6.17794E-03, 0.00000E+00, 0.00000E+00, * -3.62944E-01,-4.80022E-02,-7.57230E+01,-1.99656E-03, 0.00000E+00/ DATA PA3/ * -5.18780E-03,-1.73990E-02,-9.03485E-03, 7.48465E-03, 1.53267E-02, * 1.06296E-02, 1.18655E-02, 2.55569E-03, 1.69020E-03, 3.51936E-02, * -1.81242E-02, 0.00000E+00,-1.00529E-01,-5.10574E-03, 0.00000E+00, * 2.10228E-03, 0.00000E+00, 0.00000E+00,-1.73255E+02, 5.07833E-01, * -2.41408E-01, 8.75414E-03, 2.77527E-03,-8.90353E-05,-5.25148E+00, * -5.83899E-03,-2.09122E-02,-9.63530E-03, 9.77164E-03, 4.07051E-03, * 2.53555E-04,-5.52875E+00,-3.55993E-01,-2.49231E-03, 0.00000E+00, * 0.00000E+00, 2.86026E+01, 0.00000E+00, 3.42722E-04, 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/ C O DENSITY DATA PB1/ * 1.02315E+00,-1.59710E-01,-1.06630E-01,-1.77074E-02,-4.42726E-03, * 3.44803E-02, 4.45613E-02,-3.33751E-02,-5.73598E-02, 3.50360E-01, * 6.33053E-02, 2.16221E-02, 5.42577E-02,-5.74193E+00, 0.00000E+00, * 1.90891E-01,-1.39194E-02, 1.01102E+02, 8.16363E-02, 1.33717E-04, * 6.54403E-06, 3.10295E-03, 0.00000E+00, 0.00000E+00, 5.38205E-02, * 1.23910E-01,-1.39831E-02, 0.00000E+00, 0.00000E+00,-3.95915E-06, * 0.00000E+00,-7.14651E-01,-5.01027E-03, 0.00000E+00,-3.24756E-03, * 0.00000E+00, 0.00000E+00, 4.42173E-02,-1.31598E+01,-3.15626E-03, * 1.24574E-03,-1.47626E-03,-1.55461E-03, 6.40682E-02, 1.34898E-01, * -2.42415E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00, 6.13666E-04/ DATA PB2/ * -5.40373E-03, 2.61635E-05,-3.33012E-03, 0.00000E+00,-3.08101E-03, * -2.42679E-03,-3.36086E-03, 0.00000E+00,-1.18979E+03,-5.04738E-02, * -2.61547E-03,-1.03132E-03, 1.91583E-04,-8.38132E+01,-1.40517E-02, * -1.14167E-02,-4.08012E-03, 1.73522E-04,-1.39644E-02,-6.64128E-02, * -6.85152E-02,-1.34414E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 6.07916E+02,-4.12220E-03,-2.20996E-03, 0.00000E+00, 1.70277E+03, * -4.63015E-03, 0.00000E+00, 0.00000E+00,-2.25360E-03,-2.96204E-03, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 3.92786E-02, 1.31186E-02,-1.78086E-03, 0.00000E+00, 0.00000E+00, * -3.90083E-01,-2.84741E-02,-7.78400E+01,-1.02601E-03, 0.00000E+00/ DATA PB3/ * -7.26485E-04,-5.42181E-03,-5.59305E-03, 1.22825E-02, 1.23868E-02, * 6.68835E-03,-1.03303E-02,-9.51903E-03, 2.70021E-04,-2.57084E-02, * -1.32430E-02, 0.00000E+00,-3.81000E-02,-3.16810E-03, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00,-9.05762E-04,-2.14590E-03,-1.17824E-03, 3.66732E+00, * -3.79729E-04,-6.13966E-03,-5.09082E-03,-1.96332E-03,-3.08280E-03, * -9.75222E-04, 4.03315E+00,-2.52710E-01, 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/ C N2 DENSITY DATA PC1/ * 1.16112E+00, 0.00000E+00, 0.00000E+00, 3.33725E-02, 0.00000E+00, * 3.48637E-02,-5.44368E-03, 0.00000E+00,-6.73940E-02, 1.74754E-01, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.74712E+02, 0.00000E+00, * 1.26733E-01, 0.00000E+00, 1.03154E+02, 5.52075E-02, 0.00000E+00, * 0.00000E+00, 8.13525E-04, 0.00000E+00, 0.00000E+00, 8.66784E-02, * 1.58727E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00,-2.50482E+01, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-2.48894E-03, * 6.16053E-04,-5.79716E-04, 2.95482E-03, 8.47001E-02, 1.70147E-01, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ DATA PC2/ * 0.00000E+00, 2.47425E-05, 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.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.00000E+00, 0.00000E+00/ DATA PC3/ * 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.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.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ C TLB DATA PD1/ * 9.44846E-01, 0.00000E+00, 0.00000E+00,-3.08617E-02, 0.00000E+00, * -2.44019E-02, 6.48607E-03, 0.00000E+00, 3.08181E-02, 4.59392E-02, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.74712E+02, 0.00000E+00, * 2.13260E-02, 0.00000E+00,-3.56958E+02, 0.00000E+00, 1.82278E-04, * 0.00000E+00, 3.07472E-04, 0.00000E+00, 0.00000E+00, 8.66784E-02, * 1.58727E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 3.83054E-03, 0.00000E+00, 0.00000E+00, * -1.93065E-03,-1.45090E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00,-1.23493E-03, 1.36736E-03, 8.47001E-02, 1.70147E-01, * 3.71469E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ DATA PD2/ * 5.10250E-03, 2.47425E-05, 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.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, 3.68756E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00/ DATA PD3/ * 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.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.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ C O2 DENSITY DATA PE1/ * 1.38720E+00, 1.44816E-01, 0.00000E+00, 6.07767E-02, 0.00000E+00, * 2.94777E-02, 7.46900E-02, 0.00000E+00,-9.23822E-02, 8.57342E-02, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.38636E+01, 0.00000E+00, * 7.71653E-02, 0.00000E+00, 8.18751E+01, 1.87736E-02, 0.00000E+00, * 0.00000E+00, 1.49667E-02, 0.00000E+00, 0.00000E+00, 8.66784E-02, * 1.58727E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00,-3.67874E+02, 5.48158E-03, 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, 8.47001E-02, 1.70147E-01, * 1.22631E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ DATA PE2/ * 8.17187E-03, 3.71617E-05, 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,-2.10826E-03, * -3.13640E-03, 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.00000E+00, * -7.35742E-02,-5.00266E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 1.94965E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00/ DATA PE3/ * 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.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.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ C AR DENSITY DATA PF1/ * 1.04761E+00, 2.00165E-01, 2.37697E-01, 3.68552E-02, 0.00000E+00, * 3.57202E-02,-2.14075E-01, 0.00000E+00,-1.08018E-01,-3.73981E-01, * 0.00000E+00, 3.10022E-02,-1.16305E-03,-2.07596E+01, 0.00000E+00, * 8.64502E-02, 0.00000E+00, 9.74908E+01, 5.16707E-02, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 8.66784E-02, * 1.58727E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 3.46193E+02, 1.34297E-02, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-3.48509E-03, * -1.54689E-04, 0.00000E+00, 0.00000E+00, 8.47001E-02, 1.70147E-01, * 1.47753E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ DATA PF2/ * 1.89320E-02, 3.68181E-05, 1.32570E-02, 0.00000E+00, 0.00000E+00, * 3.59719E-03, 7.44328E-03,-1.00023E-03,-6.50528E+03, 0.00000E+00, * 1.03485E-02,-1.00983E-03,-4.06916E-03,-6.60864E+01,-1.71533E-02, * 1.10605E-02, 1.20300E-02,-5.20034E-03, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * -2.62769E+03, 7.13755E-03, 4.17999E-03, 0.00000E+00, 1.25910E+04, * 0.00000E+00, 0.00000E+00, 0.00000E+00,-2.23595E-03, 4.60217E-03, * 5.71794E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * -3.18353E-02,-2.35526E-02,-1.36189E-02, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 2.03522E-02,-6.67837E+01,-1.09724E-03, 0.00000E+00/ DATA PF3/ * -1.38821E-02, 1.60468E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.51574E-02, * -5.44470E-04, 0.00000E+00, 7.28224E-02, 6.59413E-02, 0.00000E+00, * -5.15692E-03, 0.00000E+00, 0.00000E+00,-3.70367E+03, 0.00000E+00, * 0.00000E+00, 1.36131E-02, 5.38153E-03, 0.00000E+00, 4.76285E+00, * -1.75677E-02, 2.26301E-02, 0.00000E+00, 1.76631E-02, 4.77162E-03, * 0.00000E+00, 5.39354E+00, 0.00000E+00,-7.51710E-03, 0.00000E+00, * 0.00000E+00,-8.82736E+01, 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/ C H DENSITY DATA PG1/ * 1.26376E+00,-2.14304E-01,-1.49984E-01, 2.30404E-01, 2.98237E-02, * 2.68673E-02, 2.96228E-01, 2.21900E-02,-2.07655E-02, 4.52506E-01, * 1.20105E-01, 3.24420E-02, 4.24816E-02,-9.14313E+00, 0.00000E+00, * 2.47178E-02,-2.88229E-02, 8.12805E+01, 5.10380E-02,-5.80611E-03, * 2.51236E-05,-1.24083E-02, 0.00000E+00, 0.00000E+00, 8.66784E-02, * 1.58727E-01,-3.48190E-02, 0.00000E+00, 0.00000E+00, 2.89885E-05, * 0.00000E+00, 1.53595E+02,-1.68604E-02, 0.00000E+00, 1.01015E-02, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.84552E-04, * -1.22181E-03, 0.00000E+00, 0.00000E+00, 8.47001E-02, 1.70147E-01, * -1.04927E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00,-5.91313E-03/ DATA PG2/ * -2.30501E-02, 3.14758E-05, 0.00000E+00, 0.00000E+00, 1.26956E-02, * 8.35489E-03, 3.10513E-04, 0.00000E+00, 3.42119E+03,-2.45017E-03, * -4.27154E-04, 5.45152E-04, 1.89896E-03, 2.89121E+01,-6.49973E-03, * -1.93855E-02,-1.48492E-02, 0.00000E+00,-5.10576E-02, 7.87306E-02, * 9.51981E-02,-1.49422E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 2.65503E+02, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 6.37110E-03, 3.24789E-04, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 6.14274E-02, 1.00376E-02,-8.41083E-04, 0.00000E+00, 0.00000E+00, * 0.00000E+00,-1.27099E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00/ DATA PG3/ * -3.94077E-03,-1.28601E-02,-7.97616E-03, 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,-6.71465E-03,-1.69799E-03, 1.93772E-03, 3.81140E+00, * -7.79290E-03,-1.82589E-02,-1.25860E-02,-1.04311E-02,-3.02465E-03, * 2.43063E-03, 3.63237E+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/ C N DENSITY DATA PH1/ * 7.09557E+01,-3.26740E-01, 0.00000E+00,-5.16829E-01,-1.71664E-03, * 9.09310E-02,-6.71500E-01,-1.47771E-01,-9.27471E-02,-2.30862E-01, * -1.56410E-01, 1.34455E-02,-1.19717E-01, 2.52151E+00, 0.00000E+00, * -2.41582E-01, 5.92939E-02, 4.39756E+00, 9.15280E-02, 4.41292E-03, * 0.00000E+00, 8.66807E-03, 0.00000E+00, 0.00000E+00, 8.66784E-02, * 1.58727E-01, 9.74701E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 6.70217E+01,-1.31660E-03, 0.00000E+00,-1.65317E-02, * 0.00000E+00, 0.00000E+00, 8.50247E-02, 2.77428E+01, 4.98658E-03, * 6.15115E-03, 9.50156E-03,-2.12723E-02, 8.47001E-02, 1.70147E-01, * -2.38645E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.37380E-03/ DATA PH2/ * -8.41918E-03, 2.80145E-05, 7.12383E-03, 0.00000E+00,-1.66209E-02, * 1.03533E-04,-1.68898E-02, 0.00000E+00, 3.64526E+03, 0.00000E+00, * 6.54077E-03, 3.69130E-04, 9.94419E-04, 8.42803E+01,-1.16124E-02, * -7.74414E-03,-1.68844E-03, 1.42809E-03,-1.92955E-03, 1.17225E-01, * -2.41512E-02, 1.50521E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 1.60261E+03, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00,-3.54403E-04,-1.87270E-02, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 2.76439E-02, 6.43207E-03,-3.54300E-02, 0.00000E+00, 0.00000E+00, * 0.00000E+00,-2.80221E-02, 8.11228E+01,-6.75255E-04, 0.00000E+00/ DATA PH3/ * -1.05162E-02,-3.48292E-03,-6.97321E-03, 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,-1.45546E-03,-1.31970E-02,-3.57751E-03,-1.09021E+00, * -1.50181E-02,-7.12841E-03,-6.64590E-03,-3.52610E-03,-1.87773E-02, * -2.22432E-03,-3.93895E-01, 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/ C HOT O DENSITY DATA PI1/ * 6.04050E-02, 1.57034E+00, 2.99387E-02, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-1.51018E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00,-8.61650E+00, 1.26454E-02, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 5.50878E-03, 0.00000E+00, 0.00000E+00, 8.66784E-02, * 1.58727E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 6.23881E-02, 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, 8.47001E-02, 1.70147E-01, * -9.45934E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ DATA PI2/ * 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.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.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ DATA PI3/ * 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.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.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ C S PARAM DATA PJ1/ * 9.56827E-01, 6.20637E-02, 3.18433E-02, 0.00000E+00, 0.00000E+00, * 3.94900E-02, 0.00000E+00, 0.00000E+00,-9.24882E-03,-7.94023E-03, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.74712E+02, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 2.74677E-03, 0.00000E+00, 1.54951E-02, 8.66784E-02, * 1.58727E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00,-6.99007E-04, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 1.24362E-02,-5.28756E-03, 8.47001E-02, 1.70147E-01, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ DATA PJ2/ * 0.00000E+00, 2.47425E-05, 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.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.00000E+00, 0.00000E+00/ DATA PJ3/ * 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.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.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ C TURBO DATA PK1/ * 1.09930E+00, 3.90631E+00, 3.07165E+00, 9.86161E-01, 1.63536E+01, * 4.63830E+00, 1.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, 2.93318E-02, 1.18339E-01, * 1.22732E+00, 1.02669E-01, 1.17681E+00, 2.12185E+00, 1.00000E+00, * 1.00000E+00, 1.08607E+00, 1.34836E+00, 1.10016E+00, 7.34129E-01, * 1.15241E+00, 2.22784E+00, 7.95907E-01, 4.03601E+00, 4.39732E+00, * 1.23435E+02,-4.52411E-02, 1.68986E-06, 7.44294E-01, 1.03604E+00, * 1.72783E+02, 1.17681E+00, 2.12185E+00,-7.83697E-01, 9.49154E-01/ C LOWER BOUNDARY DATA PTM/ L 1.04130E+03, 3.86000E+02, 1.95000E+02, 1.66728E+01, 2.13000E+02, L 1.20000E+02, 2.40000E+02, 1.87000E+02,-2.00000E+00, 0.00000E+00/ DATA PDM/ L 2.45600E+07, 6.71072E-06, 1.00000E+02, 0.00000E+00, 1.10000E+02, L 1.00000E+01, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, C L 8.59400E+10, 1.00000E+00, 1.05000E+02,-8.00000E+00, 1.10000E+02, L 1.00000E+01, 9.00000E+01, 2.00000E+00, 0.00000E+00, 0.00000E+00, C L 2.81000E+11, 0.00000E+00, 1.05000E+02, 2.80000E+01, 2.89500E+01, L 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, C L 3.30000E+10, 2.68270E-01, 1.05000E+02, 1.00000E+00, 1.10000E+02, L 1.00000E+01, 1.10000E+02,-1.00000E+01, 0.00000E+00, 0.00000E+00, C L 1.33000E+09, 1.19615E-02, 1.05000E+02, 0.00000E+00, 1.10000E+02, L 1.00000E+01, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, C L 1.76100E+05, 1.00000E+00, 9.50000E+01,-8.00000E+00, 1.10000E+02, L 1.00000E+01, 9.00000E+01, 2.00000E+00, 0.00000E+00, 0.00000E+00, C L 1.00000E+07, 1.00000E+00, 1.05000E+02,-8.00000E+00, 1.10000E+02, L 1.00000E+01, 9.00000E+01, 2.00000E+00, 0.00000E+00, 0.00000E+00, C L 1.00000E+06, 1.00000E+00, 1.05000E+02,-8.00000E+00, 5.50000E+02, L 7.60000E+01, 9.00000E+01, 2.00000E+00, 0.00000E+00, 4.00000E+03/ C TN1(2) DATA PL1/ * 1.00858E+00, 4.56011E-02,-2.22972E-02,-5.44388E-02, 5.23136E-04, * -1.88849E-02, 5.23707E-02,-9.43646E-03, 6.31707E-03,-7.80460E-02, * -4.88430E-02, 0.00000E+00, 0.00000E+00,-7.60250E+00, 0.00000E+00, * -1.44635E-02,-1.76843E-02,-1.21517E+02, 2.85647E-02, 0.00000E+00, * 0.00000E+00, 6.31792E-04, 0.00000E+00, 5.77197E-03, 8.66784E-02, * 1.58727E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00,-8.90272E+03, 3.30611E-03, 3.02172E-03, 0.00000E+00, * -2.13673E-03,-3.20910E-04, 0.00000E+00, 0.00000E+00, 2.76034E-03, * 2.82487E-03,-2.97592E-04,-4.21534E-03, 8.47001E-02, 1.70147E-01, * 8.96456E-03, 0.00000E+00,-1.08596E-02, 0.00000E+00, 0.00000E+00/ DATA PL2/ * 5.57917E-03, 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.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, 9.65405E-03, 0.00000E+00, 0.00000E+00, 2.00000E+00/ C TN1(3) DATA PM1/ * 9.39664E-01, 8.56514E-02,-6.79989E-03, 2.65929E-02,-4.74283E-03, * 1.21855E-02,-2.14905E-02, 6.49651E-03,-2.05477E-02,-4.24952E-02, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.19148E+01, 0.00000E+00, * 1.18777E-02,-7.28230E-02,-8.15965E+01, 1.73887E-02, 0.00000E+00, * 0.00000E+00, 0.00000E+00,-1.44691E-02, 2.80259E-04, 8.66784E-02, * 1.58727E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 2.16584E+02, 3.18713E-03, 7.37479E-03, 0.00000E+00, * -2.55018E-03,-3.92806E-03, 0.00000E+00, 0.00000E+00,-2.89757E-03, * -1.33549E-03, 1.02661E-03, 3.53775E-04, 8.47001E-02, 1.70147E-01, * -9.17497E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ DATA PM2/ * 3.56082E-03, 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.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,-1.00902E-02, 0.00000E+00, 0.00000E+00, 2.00000E+00/ C TN1(4) DATA PN1/ * 9.85982E-01,-4.55435E-02, 1.21106E-02, 2.04127E-02,-2.40836E-03, * 1.11383E-02,-4.51926E-02, 1.35074E-02,-6.54139E-03, 1.15275E-01, * 1.28247E-01, 0.00000E+00, 0.00000E+00,-5.30705E+00, 0.00000E+00, * -3.79332E-02,-6.24741E-02, 7.71062E-01, 2.96315E-02, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 6.81051E-03,-4.34767E-03, 8.66784E-02, * 1.58727E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 1.07003E+01,-2.76907E-03, 4.32474E-04, 0.00000E+00, * 1.31497E-03,-6.47517E-04, 0.00000E+00,-2.20621E+01,-1.10804E-03, * -8.09338E-04, 4.18184E-04, 4.29650E-03, 8.47001E-02, 1.70147E-01, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ DATA PN2/ * -4.04337E-03, 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,-9.52550E-04, * 8.56253E-04, 4.33114E-04, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.21223E-03, * 2.38694E-04, 9.15245E-04, 1.28385E-03, 8.67668E-04,-5.61425E-06, * 1.04445E+00, 3.41112E+01, 0.00000E+00,-8.40704E-01,-2.39639E+02, * 7.06668E-01,-2.05873E+01,-3.63696E-01, 2.39245E+01, 0.00000E+00, * -1.06657E-03,-7.67292E-04, 1.54534E-04, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.00000E+00/ C TN1(5) TN2(1) DATA PO1/ * 1.00320E+00, 3.83501E-02,-2.38983E-03, 2.83950E-03, 4.20956E-03, * 5.86619E-04, 2.19054E-02,-1.00946E-02,-3.50259E-03, 4.17392E-02, * -8.44404E-03, 0.00000E+00, 0.00000E+00, 4.96949E+00, 0.00000E+00, * -7.06478E-03,-1.46494E-02, 3.13258E+01,-1.86493E-03, 0.00000E+00, * -1.67499E-02, 0.00000E+00, 0.00000E+00, 5.12686E-04, 8.66784E-02, * 1.58727E-01,-4.64167E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 4.37353E-03,-1.99069E+02, 0.00000E+00,-5.34884E-03, 0.00000E+00, * 1.62458E-03, 2.93016E-03, 2.67926E-03, 5.90449E+02, 0.00000E+00, * 0.00000E+00,-1.17266E-03,-3.58890E-04, 8.47001E-02, 1.70147E-01, * 0.00000E+00, 0.00000E+00, 1.38673E-02, 0.00000E+00, 0.00000E+00/ DATA PO2/ * 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, 1.60571E-03, * 6.28078E-04, 5.05469E-05, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-1.57829E-03, * -4.00855E-04, 5.04077E-05,-1.39001E-03,-2.33406E-03,-4.81197E-04, * 1.46758E+00, 6.20332E+00, 0.00000E+00, 3.66476E-01,-6.19760E+01, * 3.09198E-01,-1.98999E+01, 0.00000E+00,-3.29933E+02, 0.00000E+00, * -1.10080E-03,-9.39310E-05, 1.39638E-04, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.00000E+00/ C TN2(2) DATA PP1/ * 9.81637E-01,-1.41317E-03, 3.87323E-02, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-3.58707E-02, * -8.63658E-03, 0.00000E+00, 0.00000E+00,-2.02226E+00, 0.00000E+00, * -8.69424E-03,-1.91397E-02, 8.76779E+01, 4.52188E-03, 0.00000E+00, * 2.23760E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00,-7.07572E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, * -4.11210E-03, 3.50060E+01, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00,-8.36657E-03, 1.61347E+01, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00,-1.45130E-02, 0.00000E+00, 0.00000E+00/ DATA PP2/ * 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, 1.24152E-03, * 6.43365E-04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.33255E-03, * 2.42657E-03, 1.60666E-03,-1.85728E-03,-1.46874E-03,-4.79163E-06, * 1.22464E+00, 3.53510E+01, 0.00000E+00, 4.49223E-01,-4.77466E+01, * 4.70681E-01, 8.41861E+00,-2.88198E-01, 1.67854E+02, 0.00000E+00, * 7.11493E-04, 6.05601E-04, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.00000E+00/ C TN2(3) DATA PQ1/ * 1.00422E+00,-7.11212E-03, 5.24480E-03, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-5.28914E-02, * -2.41301E-02, 0.00000E+00, 0.00000E+00,-2.12219E+01,-1.03830E-02, * -3.28077E-03, 1.65727E-02, 1.68564E+00,-6.68154E-03, 0.00000E+00, * 1.45155E-02, 0.00000E+00, 8.42365E-03, 0.00000E+00, 0.00000E+00, * 0.00000E+00,-4.34645E-03, 0.00000E+00, 0.00000E+00, 2.16780E-02, * 0.00000E+00,-1.38459E+02, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 7.04573E-03,-4.73204E+01, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 1.08767E-02, 0.00000E+00, 0.00000E+00/ DATA PQ2/ * 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,-8.08279E-03, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 5.21769E-04, * -2.27387E-04, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 3.26769E-03, * 3.16901E-03, 4.60316E-04,-1.01431E-04, 1.02131E-03, 9.96601E-04, * 1.25707E+00, 2.50114E+01, 0.00000E+00, 4.24472E-01,-2.77655E+01, * 3.44625E-01, 2.75412E+01, 0.00000E+00, 7.94251E+02, 0.00000E+00, * 2.45835E-03, 1.38871E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.00000E+00/ C TN2(4) TN3(1) DATA PR1/ * 1.01890E+00,-2.46603E-02, 1.00078E-02, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-6.70977E-02, * -4.02286E-02, 0.00000E+00, 0.00000E+00,-2.29466E+01,-7.47019E-03, * 2.26580E-03, 2.63931E-02, 3.72625E+01,-6.39041E-03, 0.00000E+00, * 9.58383E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00,-1.85291E-03, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 1.39717E+02, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 9.19771E-03,-3.69121E+02, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00,-1.57067E-02, 0.00000E+00, 0.00000E+00/ DATA PR2/ * 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,-7.07265E-03, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-2.92953E-03, * -2.77739E-03,-4.40092E-04, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.47280E-03, * 2.95035E-04,-1.81246E-03, 2.81945E-03, 4.27296E-03, 9.78863E-04, * 1.40545E+00,-6.19173E+00, 0.00000E+00, 0.00000E+00,-7.93632E+01, * 4.44643E-01,-4.03085E+02, 0.00000E+00, 1.15603E+01, 0.00000E+00, * 2.25068E-03, 8.48557E-04,-2.98493E-04, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.00000E+00/ C TN3(2) DATA PS1/ * 9.75801E-01, 3.80680E-02,-3.05198E-02, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 3.85575E-02, * 5.04057E-02, 0.00000E+00, 0.00000E+00,-1.76046E+02, 1.44594E-02, * -1.48297E-03,-3.68560E-03, 3.02185E+01,-3.23338E-03, 0.00000E+00, * 1.53569E-02, 0.00000E+00,-1.15558E-02, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 4.89620E-03, 0.00000E+00, 0.00000E+00,-1.00616E-02, * -8.21324E-03,-1.57757E+02, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 6.63564E-03, 4.58410E+01, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00,-2.51280E-02, 0.00000E+00, 0.00000E+00/ DATA PS2/ * 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, 9.91215E-03, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-8.73148E-04, * -1.29648E-03,-7.32026E-05, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-4.68110E-03, * -4.66003E-03,-1.31567E-03,-7.39390E-04, 6.32499E-04,-4.65588E-04, * -1.29785E+00,-1.57139E+02, 0.00000E+00, 2.58350E-01,-3.69453E+01, * 4.10672E-01, 9.78196E+00,-1.52064E-01,-3.85084E+03, 0.00000E+00, * -8.52706E-04,-1.40945E-03,-7.26786E-04, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.00000E+00/ C TN3(3) DATA PU1/ * 9.60722E-01, 7.03757E-02,-3.00266E-02, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.22671E-02, * 4.10423E-02, 0.00000E+00, 0.00000E+00,-1.63070E+02, 1.06073E-02, * 5.40747E-04, 7.79481E-03, 1.44908E+02, 1.51484E-04, 0.00000E+00, * 1.97547E-02, 0.00000E+00,-1.41844E-02, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 5.77884E-03, 0.00000E+00, 0.00000E+00, 9.74319E-03, * 0.00000E+00,-2.88015E+03, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00,-4.44902E-03,-2.92760E+01, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 2.34419E-02, 0.00000E+00, 0.00000E+00/ DATA PU2/ * 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, 5.36685E-03, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-4.65325E-04, * -5.50628E-04, 3.31465E-04, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-2.06179E-03, * -3.08575E-03,-7.93589E-04,-1.08629E-04, 5.95511E-04,-9.05050E-04, * 1.18997E+00, 4.15924E+01, 0.00000E+00,-4.72064E-01,-9.47150E+02, * 3.98723E-01, 1.98304E+01, 0.00000E+00, 3.73219E+03, 0.00000E+00, * -1.50040E-03,-1.14933E-03,-1.56769E-04, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.00000E+00/ C TN3(4) DATA PV1/ * 1.03123E+00,-7.05124E-02, 8.71615E-03, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-3.82621E-02, * -9.80975E-03, 0.00000E+00, 0.00000E+00, 2.89286E+01, 9.57341E-03, * 0.00000E+00, 0.00000E+00, 8.66153E+01, 7.91938E-04, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 4.68917E-03, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 7.86638E-03, 0.00000E+00, 0.00000E+00, 9.90827E-03, * 0.00000E+00, 6.55573E+01, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00,-4.00200E+01, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 7.07457E-03, 0.00000E+00, 0.00000E+00/ DATA PV2/ * 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, 5.72268E-03, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-2.04970E-04, * 1.21560E-03,-8.05579E-06, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-2.49941E-03, * -4.57256E-04,-1.59311E-04, 2.96481E-04,-1.77318E-03,-6.37918E-04, * 1.02395E+00, 1.28172E+01, 0.00000E+00, 1.49903E-01,-2.63818E+01, * 0.00000E+00, 4.70628E+01,-2.22139E-01, 4.82292E-02, 0.00000E+00, * -8.67075E-04,-5.86479E-04, 5.32462E-04, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.00000E+00/ C TN3(5) SURFACE TEMP TSL DATA PW1/ * 1.00828E+00,-9.10404E-02,-2.26549E-02, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-2.32420E-02, * -9.08925E-03, 0.00000E+00, 0.00000E+00, 3.36105E+01, 0.00000E+00, * 0.00000E+00, 0.00000E+00,-1.24957E+01,-5.87939E-03, 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, 2.79765E+01, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.01237E+03, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00,-1.75553E-02, 0.00000E+00, 0.00000E+00/ DATA PW2/ * 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, 3.29699E-03, * 1.26659E-03, 2.68402E-04, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.17894E-03, * 1.48746E-03, 1.06478E-04, 1.34743E-04,-2.20939E-03,-6.23523E-04, * 6.36539E-01, 1.13621E+01, 0.00000E+00,-3.93777E-01, 2.38687E+03, * 0.00000E+00, 6.61865E+02,-1.21434E-01, 9.27608E+00, 0.00000E+00, * 1.68478E-04, 1.24892E-03, 1.71345E-03, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.00000E+00/ C TGN3(2) SURFACE GRAD TSLG DATA PX1/ * 1.57293E+00,-6.78400E-01, 6.47500E-01, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-7.62974E-02, * -3.60423E-01, 0.00000E+00, 0.00000E+00, 1.28358E+02, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 4.68038E+01, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00,-1.67898E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 2.90994E+04, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 3.15706E+01, 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/ DATA PX2/ * 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.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.00000E+00, 0.00000E+00, 0.00000E+00, 2.00000E+00/ C TGN2(1) TGN1(2) DATA PY1/ * 8.60028E-01, 3.77052E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,-1.17570E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 7.77757E-03, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 1.01024E+02, 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, 6.54251E+02, 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/ DATA PY2/ * 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.00000E+00,-1.56959E-02, * 1.91001E-02, 3.15971E-02, 1.00982E-02,-6.71565E-03, 2.57693E-03, * 1.38692E+00, 2.82132E-01, 0.00000E+00, 0.00000E+00, 3.81511E+02, * 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, 2.00000E+00/ C TGN3(1) TGN2(2) DATA PZ1/ * 1.06029E+00,-5.25231E-02, 3.73034E-01, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 3.31072E-02, * -3.88409E-01, 0.00000E+00, 0.00000E+00,-1.65295E+02,-2.13801E-01, * -4.38916E-02,-3.22716E-01,-8.82393E+01, 1.18458E-01, 0.00000E+00, * -4.35863E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00,-1.19782E-01, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 2.62229E+01, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00,-5.37443E+01, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00,-4.55788E-01, 0.00000E+00, 0.00000E+00/ DATA PZ2/ * 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, 3.84009E-02, * 3.96733E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 5.05494E-02, * 7.39617E-02, 1.92200E-02,-8.46151E-03,-1.34244E-02, 1.96338E-02, * 1.50421E+00, 1.88368E+01, 0.00000E+00, 0.00000E+00,-5.13114E+01, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 5.11923E-02, 3.61225E-02, 0.00000E+00, 0.00000E+00, 0.00000E+00, * 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2.00000E+00/ C SEMIANNUAL MULT SAM DATA PAA1/ * 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, * 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, * 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, * 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, * 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, * 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, * 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, * 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, * 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, * 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00/ DATA PAA2/ * 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.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.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ C MIDDLE ATMOSPHERE AVERAGES DATA PAVGM/ M 2.61000E+02, 2.64000E+02, 2.29000E+02, 2.17000E+02, 2.17000E+02, M 2.23000E+02, 2.86760E+02,-2.93940E+00, 2.50000E+00, 0.00000E+00/ END $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ C TEST DRIVER FOR GTD7 (ATMOSPHERIC MODEL) DIMENSION D(9,16),T(2,16),SW(25),APH(7) DIMENSION IDAY(15),UT(15),ALT(15),XLAT(15),XLONG(15),XLST(15), & F107A(15),F107(15),AP(15) COMMON/GTS3C/DL(16) COMMON/DATIME/ISDATE(3),ISTIME(2),NAME(2) DATA IDAY/172,81,13*172/ DATA UT/29000.,29000.,75000.,12*29000./ DATA ALT/400.,400.,1000.,100.,6*400.,0,10.,30.,50.,70./ DATA XLAT/4*60.,0.,10*60./ DATA XLONG/5*-70.,0.,9*-70./ DATA XLST/6*16.,4.,8*16./ DATA F107A/7*150.,70.,7*150./ DATA F107/8*150.,180.,6*150./ DATA AP/9*4.,40.,5*4./ DATA APH/7*100./,SW/8*1.,-1.,16*1./ DO I=1,15 CALL GTD7(IDAY(I),UT(I),ALT(I),XLAT(I),XLONG(I),XLST(I), & F107A(I),F107(I),AP(I),48,D(1,I),T(1,I)) WRITE(6,100) (D(J,I),J=1,9),T(1,I),T(2,I),DL ENDDO CALL TSELEC(SW) I=16 CALL GTD7(IDAY(1),UT(1),ALT(1),XLAT(1),XLONG(1),XLST(1), & F107A(1),F107(1),APH,48,D(1,I),T(1,I)) WRITE(6,100) (D(J,I),J=1,9),T(1,I),T(2,I),DL CALL GTD7(IDAY(1),UT(1),ALT(4),XLAT(1),XLONG(1),XLST(1), & F107A(1),F107(1),APH,48,D(1,I),T(1,I)) WRITE(6,100) (D(J,I),J=1,9),T(1,I),T(2,I),DL WRITE(6,300) NAME,ISDATE,ISTIME WRITE(6,200) (IDAY(I),I=1,5) WRITE(6,201) (UT(I),I=1,5) WRITE(6,202) (ALT(I),I=1,5) WRITE(6,203) (XLAT(I),I=1,5) WRITE(6,204) (XLONG(I),I=1,5) WRITE(6,205) (XLST(I),I=1,5) WRITE(6,206) (F107A(I),I=1,5) WRITE(6,207) (F107(I),I=1,5) WRITE(6,210) (T(1,I),I=1,5) WRITE(6,211) (T(2,I),I=1,5) WRITE(6,212) (D(1,I),I=1,5) WRITE(6,213) (D(2,I),I=1,5) WRITE(6,214) (D(3,I),I=1,5) WRITE(6,215) (D(4,I),I=1,5) WRITE(6,216) (D(5,I),I=1,5) WRITE(6,217) (D(7,I),I=1,5) WRITE(6,219) (D(8,I),I=1,5) WRITE(6,220) (D(9,I),I=1,5) WRITE(6,218) (D(6,I),I=1,5) WRITE(6,200) (IDAY(I),I=6,10) WRITE(6,201) (UT(I),I=6,10) WRITE(6,202) (ALT(I),I=6,10) WRITE(6,203) (XLAT(I),I=6,10) WRITE(6,204) (XLONG(I),I=6,10) WRITE(6,205) (XLST(I),I=6,10) WRITE(6,206) (F107A(I),I=6,10) WRITE(6,207) (F107(I),I=6,10) WRITE(6,210) (T(1,I),I=6,10) WRITE(6,211) (T(2,I),I=6,10) WRITE(6,212) (D(1,I),I=6,10) WRITE(6,213) (D(2,I),I=6,10) WRITE(6,214) (D(3,I),I=6,10) WRITE(6,215) (D(4,I),I=6,10) WRITE(6,216) (D(5,I),I=6,10) WRITE(6,217) (D(7,I),I=6,10) WRITE(6,219) (D(8,I),I=6,10) WRITE(6,220) (D(9,I),I=6,10) WRITE(6,218) (D(6,I),I=6,10) WRITE(6,200) (IDAY(I),I=11,15) WRITE(6,201) (UT(I),I=11,15) WRITE(6,202) (ALT(I),I=11,15) WRITE(6,203) (XLAT(I),I=11,15) WRITE(6,204) (XLONG(I),I=11,15) WRITE(6,205) (XLST(I),I=11,15) WRITE(6,206) (F107A(I),I=11,15) WRITE(6,207) (F107(I),I=11,15) WRITE(6,210) (T(1,I),I=11,15) WRITE(6,211) (T(2,I),I=11,15) WRITE(6,212) (D(1,I),I=11,15) WRITE(6,213) (D(2,I),I=11,15) WRITE(6,214) (D(3,I),I=11,15) WRITE(6,215) (D(4,I),I=11,15) WRITE(6,216) (D(5,I),I=11,15) WRITE(6,217) (D(7,I),I=11,15) WRITE(6,219) (D(8,I),I=11,15) WRITE(6,220) (D(9,I),I=11,15) WRITE(6,218) (D(6,I),I=11,15) 100 FORMAT(1X,1P8E9.2/4X,1PE9.2,2E10.3/4X,8E9.2/4X,8E9.2/) 200 FORMAT(//' DAY ',5I12) 201 FORMAT(' UT ',5F12.0) 202 FORMAT(' ALT ',5F12.0) 203 FORMAT(' LAT ',5F12.0) 204 FORMAT(' LONG ',5F12.0) 205 FORMAT(' LST ',5F12.0) 206 FORMAT(' F107A',5F12.0) 207 FORMAT(' F107 ',5F12.0) 210 FORMAT(/' TINF ',5F12.2) 211 FORMAT(' TG ',5F12.2) 212 FORMAT(' HE ',1P5E12.3) 213 FORMAT(' O ',1P5E12.3) 214 FORMAT(' N2 ',1P5E12.3) 215 FORMAT(' O2 ',1P5E12.3) 216 FORMAT(' AR ',1P5E12.3) 217 FORMAT(' H ',1P5E12.3) 219 FORMAT(' N ',1P5E12.3) 220 FORMAT(' ANM O',1P5E12.3) 218 FORMAT(' RHO ',1P5E12.3) 300 FORMAT(1X,2A4,2X,3A4,2X,2A4) STOP END $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 6.67E+05 1.14E+08 2.00E+07 3.96E+05 3.56E+03 4.07E-15 3.48E+04 4.10E+06 2.67E+04 1.251E+03 1.241E+03 4.02E+02 1.69E-02 4.54E+06 4.63E+10 3.16E+11 4.65E+10 2.65E+09 1.86E-11 9.49E+04 1.23E+02 1.84E+02 9.00E+01 1.43E+01 7.78E+00 1.76E+05 9.02E+08 3.41E+06 1.59E+08 1.39E+07 3.21E+05 1.56E+03 5.00E-15 4.85E+04 4.38E+06 6.96E+03 1.167E+03 1.162E+03 3.76E+02 1.88E-02 2.47E+07 8.24E+10 3.37E+11 6.14E+10 2.14E+09 2.13E-11 1.35E+05 1.23E+02 1.84E+02 9.00E+01 1.49E+01 7.63E+00 4.58E+04 1.20E+09 1.12E+05 6.93E+04 4.25E+01 1.30E-01 2.62E-05 2.76E-18 2.02E+04 5.74E+03 2.37E+04 1.240E+03 1.240E+03 4.02E+02 1.71E-02 4.95E+06 4.92E+10 3.16E+11 4.65E+10 2.48E+09 1.86E-11 8.78E+04 1.23E+02 1.84E+02 9.00E+01 1.43E+01 7.63E+00 1.75E+04 8.70E+08 7.08E+07 1.95E+11 6.12E+12 1.18E+12 6.02E+10 3.56E-10 1.06E+07 2.62E+05 0.00E+00 1.027E+03 2.069E+02 4.02E+02 2.29E-02 4.54E+06 4.63E+10 3.16E+11 4.65E+10 2.65E+09 1.86E-11 9.49E+04 1.23E+02 1.47E+02 9.00E+01 1.43E+01 7.68E+00 3.73E+07 9.02E+08 1.85E+06 1.48E+08 1.58E+07 2.60E+05 1.59E+03 4.81E-15 5.82E+04 5.48E+06 1.26E+03 1.212E+03 1.208E+03 4.12E+02 1.95E-02 1.21E+07 5.46E+10 2.20E+11 2.65E+10 9.98E+08 1.32E-11 1.53E+05 1.23E+02 1.84E+02 9.00E+01 1.56E+01 7.27E+00 8.26E+03 1.10E+09 8.67E+05 1.28E+08 1.82E+07 2.88E+05 2.40E+03 4.36E-15 3.69E+04 3.90E+06 2.67E+04 1.220E+03 1.213E+03 4.02E+02 1.75E-02 5.93E+06 5.43E+10 3.16E+11 3.76E+10 2.05E+09 1.83E-11 9.96E+04 1.23E+02 1.84E+02 9.00E+01 1.43E+01 7.79E+00 1.76E+05 8.89E+08 5.78E+05 6.98E+07 1.24E+07 2.46E+05 1.41E+03 2.47E-15 5.29E+04 1.07E+06 2.67E+04 1.116E+03 1.113E+03 3.98E+02 2.00E-02 4.03E+06 3.60E+10 3.19E+11 5.11E+10 2.16E+09 1.86E-11 1.39E+05 1.23E+02 1.84E+02 9.00E+01 1.43E+01 8.94E+00 1.76E+05 2.86E+08 3.74E+05 4.78E+07 5.24E+06 1.71E+05 5.50E+02 1.57E-15 8.90E+04 1.98E+06 9.12E+03 1.031E+03 1.025E+03 3.98E+02 1.71E-02 2.80E+06 3.72E+10 2.96E+11 1.40E+10 2.65E+09 1.57E-11 2.28E+05 1.23E+02 1.84E+02 9.00E+01 1.08E+01 7.90E+00 6.01E+04 7.52E+08 6.75E+05 1.25E+08 2.37E+07 4.84E+05 4.58E+03 4.56E-15 3.24E+04 5.37E+06 2.67E+04 1.306E+03 1.293E+03 4.04E+02 1.59E-02 4.57E+06 4.67E+10 3.16E+11 4.65E+10 2.65E+09 1.86E-11 8.99E+04 1.23E+02 1.84E+02 9.00E+01 1.43E+01 7.57E+00 1.76E+05 1.11E+09 5.53E+05 1.20E+08 3.50E+07 9.20E+05 1.10E+04 4.97E-15 2.69E+04 4.89E+06 2.81E+04 1.362E+03 1.347E+03 4.36E+02 1.55E-02 3.48E+06 3.54E+10 3.16E+11 5.70E+10 3.67E+09 1.89E-11 7.19E+04 1.23E+02 1.84E+02 9.00E+01 1.43E+01 7.87E+00 1.85E+05 8.14E+08 1.38E+14 0.00E+00 2.05E+19 5.50E+18 2.45E+17 1.26E-03 0.00E+00 0.00E+00 0.00E+00 1.027E+03 2.815E+02 4.02E+02 2.55E-02 3.48E+06 3.54E+10 3.16E+11 5.70E+10 3.67E+09 1.89E-11 7.19E+04 1.23E+02 1.47E+02 9.00E+01 1.60E+01 7.87E+00 1.19E+15 8.14E+08 4.43E+13 0.00E+00 6.60E+18 1.77E+18 7.89E+16 4.06E-04 0.00E+00 0.00E+00 0.00E+00 1.027E+03 2.274E+02 4.02E+02 2.55E-02 3.48E+06 3.54E+10 3.16E+11 5.70E+10 3.67E+09 1.89E-11 7.19E+04 1.23E+02 1.47E+02 9.00E+01 1.60E+01 7.87E+00 1.19E+15 8.14E+08 2.13E+12 0.00E+00 3.17E+17 8.51E+16 3.79E+15 1.95E-05 0.00E+00 0.00E+00 0.00E+00 1.027E+03 2.374E+02 4.02E+02 2.55E-02 3.48E+06 3.54E+10 3.16E+11 5.70E+10 3.67E+09 1.89E-11 7.19E+04 1.23E+02 1.47E+02 9.00E+01 1.60E+01 7.87E+00 1.19E+15 8.14E+08 1.41E+11 0.00E+00 2.10E+16 5.65E+15 2.52E+14 1.29E-06 0.00E+00 0.00E+00 0.00E+00 1.027E+03 2.796E+02 4.02E+02 2.55E-02 3.48E+06 3.54E+10 3.16E+11 5.70E+10 3.67E+09 1.89E-11 7.19E+04 1.23E+02 1.47E+02 9.00E+01 1.60E+01 7.87E+00 1.19E+15 8.14E+08 1.26E+10 0.00E+00 1.87E+15 4.81E+14 2.24E+13 1.14E-07 0.00E+00 0.00E+00 0.00E+00 1.027E+03 2.191E+02 4.02E+02 2.55E-02 4.54E+06 4.63E+10 3.16E+11 4.65E+10 2.65E+09 1.86E-11 9.49E+04 1.23E+02 1.47E+02 9.00E+01 1.60E+01 7.68E+00 7.53E+08 9.02E+08 5.20E+05 1.27E+08 4.85E+07 1.70E+06 2.35E+04 5.88E-15 2.50E+04 6.28E+06 2.67E+04 1.426E+03 1.409E+03 4.62E+02 1.49E-02 3.10E+06 3.10E+10 3.16E+11 7.23E+10 4.96E+09 1.97E-11 6.56E+04 1.23E+02 1.84E+02 9.00E+01 1.43E+01 7.78E+00 1.76E+05 8.82E+08 5.67E+07 1.26E+11 4.93E+12 1.02E+12 4.99E+10 2.90E-10 8.83E+06 2.25E+05 0.00E+00 1.027E+03 1.934E+02 4.62E+02 2.54E-02 3.10E+06 3.10E+10 3.16E+11 7.23E+10 4.96E+09 1.97E-11 6.56E+04 1.23E+02 1.37E+02 9.00E+01 1.43E+01 7.49E+00 3.19E+07 8.82E+08 MSISE-00 13-APR-00 17:46:08 DAY 172 81 172 172 172 UT 29000. 29000. 75000. 29000. 29000. ALT 400. 400. 1000. 100. 400. LAT 60. 60. 60. 60. 0. LONG -70. -70. -70. -70. -70. LST 16. 16. 16. 16. 16. F107A 150. 150. 150. 150. 150. F107 150. 150. 150. 150. 150. TINF 1250.54 1166.75 1239.89 1027.32 1212.40 TG 1241.42 1161.71 1239.89 206.89 1208.14 HE 6.665E+05 3.407E+06 1.124E+05 7.078E+07 1.851E+06 O 1.139E+08 1.586E+08 6.934E+04 1.946E+11 1.477E+08 N2 1.998E+07 1.391E+07 4.247E+01 6.116E+12 1.579E+07 O2 3.964E+05 3.215E+05 1.304E-01 1.184E+12 2.595E+05 AR 3.557E+03 1.560E+03 2.619E-05 6.023E+10 1.589E+03 H 3.475E+04 4.854E+04 2.017E+04 1.060E+07 5.816E+04 N 4.096E+06 4.381E+06 5.741E+03 2.616E+05 5.479E+06 ANM O 2.667E+04 6.957E+03 2.374E+04 0.000E+00 1.264E+03 RHO 4.074E-15 5.002E-15 2.757E-18 3.563E-10 4.809E-15 DAY 172 172 172 172 172 UT 29000. 29000. 29000. 29000. 29000. ALT 400. 400. 400. 400. 400. LAT 60. 60. 60. 60. 60. LONG 0. -70. -70. -70. -70. LST 16. 4. 16. 16. 16. F107A 150. 150. 70. 150. 150. F107 150. 150. 150. 180. 150. TINF 1220.15 1116.39 1031.25 1306.05 1361.87 TG 1212.71 1113.00 1024.85 1293.37 1347.39 HE 8.673E+05 5.776E+05 3.740E+05 6.748E+05 5.529E+05 O 1.279E+08 6.979E+07 4.783E+07 1.245E+08 1.198E+08 N2 1.823E+07 1.237E+07 5.240E+06 2.369E+07 3.496E+07 O2 2.880E+05 2.456E+05 1.712E+05 4.840E+05 9.203E+05 AR 2.403E+03 1.406E+03 5.502E+02 4.579E+03 1.096E+04 H 3.686E+04 5.292E+04 8.897E+04 3.245E+04 2.686E+04 N 3.897E+06 1.070E+06 1.980E+06 5.371E+06 4.890E+06 ANM O 2.667E+04 2.667E+04 9.122E+03 2.667E+04 2.805E+04 RHO 4.356E-15 2.470E-15 1.572E-15 4.564E-15 4.974E-15 DAY 172 172 172 172 172 UT 29000. 29000. 29000. 29000. 29000. ALT 0. 10. 30. 50. 70. LAT 60. 60. 60. 60. 60. LONG -70. -70. -70. -70. -70. LST 16. 16. 16. 16. 16. F107A 150. 150. 150. 150. 150. F107 150. 150. 150. 150. 150. TINF 1027.32 1027.32 1027.32 1027.32 1027.32 TG 281.46 227.42 237.44 279.56 219.07 HE 1.375E+14 4.427E+13 2.128E+12 1.412E+11 1.261E+10 O 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 N2 2.050E+19 6.598E+18 3.171E+17 2.104E+16 1.875E+15 O2 5.499E+18 1.770E+18 8.506E+16 5.645E+15 4.809E+14 AR 2.452E+17 7.892E+16 3.793E+15 2.517E+14 2.240E+13 H 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 N 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 ANM O 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 RHO 1.261E-03 4.059E-04 1.951E-05 1.295E-06 1.142E-07 $$$$$$$$$$$$$ end of NRLMSISE-00 distribution package $$$$$$$$$$$$$$