C Neutral temperature and density model used in IRI C C Includes the following subroutines and functions: C GTD7, GTD7D, GHP7, GLATF, VTST7, GTS7, METERS, SCALH, GLOBE7, TSELEC, C GLOBE7S, DENSU, DENSM, SPLINEM, SPLINTM, SPLINI,DNET, CCOR, CCOR2, C BLOCKDATA GTD7BK C C 2011.00 09/27/10 GLOBE7: AP(1),P(1)->AP(7),P(150) C 2011.00 09/27/10 GLOB7S: P(1)->P(150) C 2011.00 09/27/10 GTS7: AP(1)->AP(7), PD(150,9)->PDA1(150),..,PDA9(150) C 2012.00 10/05/11 IRI-2012: bottomside B0 B1 model (SHAMDB0D, SHAB1D), C 2012.00 10/05/11 bottomside Ni model (iriflip.for), auroral foE C 2012.00 10/05/11 storm model (storme_ap), Te with PF10.7 (elteik), C 2012.00 10/05/11 oval kp model (auroral_boundary), IGRF-11(igrf.for), C 2012.00 10/05/11 NRLMSIS00 (cira.for), CGM coordinates, F10.7 daily C 2012.00 10/05/11 81-day 365-day indices (apf107.dat), ap->kp (ckp), C 2012.00 10/05/11 array size change jf(50) outf(20,1000), oarr(100). C 2014.02 07/24/14 COMMON/iounit/konsol, mess turn off messages C 2014.05 12/22/14 TSELEC: SV(25),SVV(25) C 2015.01 02/18/15 switches set (TSELEC) only in main not in GTD7,GLOBE7 C 2015.01 07/20/15 GLOB7S: P(150)->P(100) C SUBROUTINE GTD7(IYD,SEC,ALT,GLAT,GLONG,STL,F107A,F107,AP,MASS,D,T) C----------------------------------------------------------------------- C C NRLMSISE-00 C ----------- C Neutral Atmosphere Empirical Model from the surface to lower C exosphere C J.M. Picone, A.E. Hedin, D.P. Drob, and A.C. Aikin, NRLMSISE-00 C empirical model of the atmosphere: Statistical comparisons C and scientific issues, J. Geophys. Res., 107(A12), 1468, C doi:10.1029/2002JA009430, 2002. 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 NOTES: The default should be sw(1:25)=1.0; Uses the full Ap history C if sw(9) = -1.0; sw(9)=-1 uses the full ap history and is supposed C to give the best results for geomagentic storms; sw(9) = 0 turns off C all ap variations. C C To get current values of SW: CALL TRETRV(SW) C C Simplified version for F77 compilers and IRI (Sep 27,2010, dbilitza): C GTS7, GLOBE7: AP(1), P(1) -> AP(7), P(150) C PD(150,9) -> PDA1(150),..,PDA9(150) C----------------------------------------------------------------------- CHARACTER*4 ISDATE,ISTIME,NAME,ISD,IST,NAM 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),PDA1(150),PDA2(150),PDA3(150), $ PDA4(150),PDA5(150),PDA6(150),PDA7(150),PDA8(150),PDA9(150), $ 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/ c DATA SV/25*1./ c 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 DS(1)=D(1) !.. switches on PGR mod 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.0,TZ,MN3,ZN3,TN3,TGN3,MN2,ZN2,TN2,TGN2) T(2)=TZ 90 CONTINUE ALAST=ALT RETURN END C C SUBROUTINE GTD7D(IYD,SEC,ALT,GLAT,GLONG,STL,F107A,F107,AP,MASS, $ D,T) C----------------------------------------------------------------------- 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 C SUBROUTINE GHP7(IYD,SEC,ALT,GLAT,GLONG,STL,F107A,F107,AP, $ D,T,PRESS) C----------------------------------------------------------------------- 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----------------------------------------------------------------------- LOGICAL mess COMMON/PARMB/GSURF,RE COMMON/METSEL/IMR COMMON/iounit/konsol,mess 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(mess.and.(L.EQ.LTEST)) WRITE(konsol,100) PRESS,DIFF 100 FORMAT(1X,29HGHP7 NOT CONVERGING FOR PRESS, 1PE12.2,E12.2) ALT=Z RETURN END C C SUBROUTINE GLATF(LAT,GV,REFF) C----------------------------------------------------------------------- C CALCULATE LATITUDE VARIABLE GRAVITY (GV) AND EFFECTIVE C RADIUS (REFF) C----------------------------------------------------------------------- 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 C FUNCTION VTST7(IYD,SEC,GLAT,GLONG,STL,F107A,F107,AP,IC) C----------------------------------------------------------------------- C Test if geophysical variables or switches changed and save C Return 0 if unchanged and 1 if changed C----------------------------------------------------------------------- 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 C SUBROUTINE GTS7(IYD,SEC,ALT,GLAT,GLONG,STL,F107A,F107,AP,MASS,D,T) C----------------------------------------------------------------------- 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),APLOW(7) LOGICAL mess COMMON/iounit/konsol,mess 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) C DIMENSION D(9),T(2),MT(11),AP(1),ALTL(8) DIMENSION D(9),T(2),MT(11),AP(7),ALTL(8) COMMON/LOWER7/PTM(10),PDM(10,8) C COMMON/PARM7/PT(150),PD(150,9),PS(150),PDL(25,2),PTL(100,4), C $ PMA(100,10),SAM(100) COMMON/PARM7/PT(150),PDA1(150),PDA2(150),PDA3(150),PDA4(150), $ PDA5(150),PDA6(150),PDA7(150),PDA8(150),PDA9(150), $ 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./ TNMOD=0 !.. for switching on mod MSIS IF(D(1).LT.0) TNMOD=-D(1) !.. PGR 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 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 IF(TNMOD.GT.600) TINF=TNMOD !.. PGR 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,PDA4))*PDA4(1) 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,PDA3) 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 if(mess) WRITE(konsol,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)*PDA3(1) 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,PDA1) C Diffusive density at Zlb DB04 = PDM(1,1)*EXP(G4)*PDA1(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 C G16= SW(21)*GLOBE7(YRD,SEC,GLAT,GLONG,STL,F107A,F107,AP,PDA2) G16= SW(21)*GLOBE7(YRD,SEC,GLAT,GLONG,STL,F107A,F107,AP,PDA2) C Diffusive density at Zlb C DB16 = PDM(1,2)*EXP(G16)*PDA2(1) DB16 = PDM(1,2)*EXP(G16)*PDA2(1) 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) HC216=PDM(6,2)*PDL(5,2) D(2)=D(2)*CCOR2(Z,RL,HC16,ZC16,HC216) 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,PDA5) C Diffusive density at Zlb DB32 = PDM(1,4)*EXP(G32)*PDA5(1) 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) HCC232=PDM(8,4)*PDL(23,1) 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)*CCOR2(Z,RC32,HCC32,ZCC32,HCC232) 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,PDA6) C Diffusive density at Zlb DB40 = PDM(1,5)*EXP(G40)*PDA6(1) 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,PDA7) C Diffusive density at Zlb DB01 = PDM(1,6)*EXP(G1)*PDA7(1) 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,PDA8) C Diffusive density at Zlb DB14 = PDM(1,7)*EXP(G14)*PDA8(1) 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,PDA9) DB16H = PDM(1,8)*EXP(G16H)*PDA9(1) 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 C SUBROUTINE METERS(METER) C----------------------------------------------------------------------- C Convert outputs to Kg & Meters if METER true C----------------------------------------------------------------------- LOGICAL METER COMMON/METSEL/IMR SAVE IMR=0 IF(METER) IMR=1 END C C FUNCTION SCALH(ALT,XM,TEMP) C----------------------------------------------------------------------- C Calculate scale height (km) C----------------------------------------------------------------------- COMMON/PARMB/GSURF,RE SAVE DATA RGAS/831.4/ G=GSURF/(1.+ALT/RE)**2 SCALH=RGAS*TEMP/(G*XM) RETURN END C C FUNCTION GLOBE7(YRD,SEC,LAT,LONG,TLOC,F107A,F107,AP,P) C----------------------------------------------------------------------- C CALCULATE G(L) FUNCTION C Upper Thermosphere Parameters C----------------------------------------------------------------------- REAL LAT, LONG C DIMENSION P(1),SV(25),AP(1) DIMENSION P(150),SV(25),AP(7) 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 C SUBROUTINE TSELEC(SV) C----------------------------------------------------------------------- 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----------------------------------------------------------------------- c DIMENSION SV(1),SAV(25),SVV(1) DIMENSION SV(25),SAV(25),SVV(25) 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 C FUNCTION GLOB7S(P) C----------------------------------------------------------------------- C VERSION OF GLOBE FOR LOWER ATMOSPHERE 10/26/99 C----------------------------------------------------------------------- REAL LONG LOGICAL mess COMMON/iounit/konsol,mess 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) C DIMENSION P(1),T(14) C DIMENSION P(150),T(14) DIMENSION P(100),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 if(mess) WRITE(konsol,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 C FUNCTION DENSU(ALT,DLB,TINF,TLB,XM,ALPHA,TZ,ZLB,S2, $ MN1,ZN1,TN1,TGN1) C-------------------------------------------------------------------- C Calculate Temperature and Density Profiles for MSIS models C New lower thermo polynomial 10/30/89 C-------------------------------------------------------------------- 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) 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 SPLINEM(XS,YS,MN,YD1,YD2,Y2OUT) X=ZG/ZGDIF CALL SPLINTM(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 C FUNCTION DENSM(ALT,D0,XM,TZ,MN3,ZN3,TN3,TGN3,MN2,ZN2,TN2,TGN2) C-------------------------------------------------------------------- C Calculate Temperature and Density Profiles for lower atmos. C-------------------------------------------------------------------- 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 SPLINEM(XS,YS,MN,YD1,YD2,Y2OUT) X=ZG/ZGDIF CALL SPLINTM(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 SPLINEM(XS,YS,MN,YD1,YD2,Y2OUT) X=ZG/ZGDIF CALL SPLINTM(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 C SUBROUTINE SPLINEM(X,Y,N,YP1,YPN,Y2) C----------------------------------------------------------------------- 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 C----------------------------------------------------------------------- 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 C SUBROUTINE SPLINTM(XA,YA,Y2A,N,X,Y) C----------------------------------------------------------------------- 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 C----------------------------------------------------------------------- DIMENSION XA(N),YA(N),Y2A(N) LOGICAL mess COMMON/iounit/konsol,mess 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(mess.and.(H.EQ.0)) WRITE(konsol,*) '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 C SUBROUTINE SPLINI(XA,YA,Y2A,N,X,YI) C----------------------------------------------------------------------- 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 C----------------------------------------------------------------------- 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 C FUNCTION DNET(DD,DM,ZHM,XMM,XM) C----------------------------------------------------------------------- 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 C----------------------------------------------------------------------- LOGICAL mess COMMON/iounit/konsol,mess SAVE A=ZHM/(XMM-XM) IF(DM.GT.0.AND.DD.GT.0) GOTO 5 if(mess) WRITE(konsol,*) '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 C FUNCTION CCOR(ALT, R,H1,ZH) C----------------------------------------------------------------------- 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 C----------------------------------------------------------------------- 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 C FUNCTION CCOR2(ALT, R,H1,ZH,H2) C----------------------------------------------------------------------- C O&O2 CHEMISTRY/DISSOCIATION CORRECTION FOR MSIS MODELS C----------------------------------------------------------------------- E1=(ALT-ZH)/H1 E2=(ALT-ZH)/H2 IF(E1.GT.70. .OR. E2.GT.70.) GO TO 20 IF(E1.LT.-70. .AND. E2.LT.-70) GO TO 10 EX1=EXP(E1) EX2=EXP(E2) CCOR2=R/(1.+.5*(EX1+EX2)) GO TO 50 10 CCOR2=R GO TO 50 20 CCOR2=0. GO TO 50 50 CONTINUE CCOR2=EXP(CCOR2) RETURN END C C BLOCK DATA GTD7BK C----------------------------------------------------------------------- C MSISE-00 01-FEB-02 C----------------------------------------------------------------------- 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) CHARACTER*4 ISDATE,ISTIME,NAME 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/'01-F','EB-0','2 '/,ISTIME/'15:4','9:27'/ 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.35580E+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, 1.28840E+00, 3.10302E-02, 1.18339E-01, * 1.00000E+00, 7.00000E-01, 1.15020E+00, 3.44689E+00, 1.28840E+00, * 1.00000E+00, 1.08738E+00, 1.22947E+00, 1.10016E+00, 7.34129E-01, * 1.15241E+00, 2.22784E+00, 7.95046E-01, 4.01612E+00, 4.47749E+00, * 1.23435E+02,-7.60535E-02, 1.68986E-06, 7.44294E-01, 1.03604E+00, * 1.72783E+02, 1.15020E+00, 3.44689E+00,-7.46230E-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