term@box:~/Documents/temp> diff -u hm_orig.f hm.f
--- hm_orig.f 2012-03-04 00:22:55.000000000 +0100
+++ hm.f 2012-03-04 00:18:56.000000000 +0100
@@ -134,6 +134,7 @@
C
DIMENSION D(9),T(2),AP(7),DS(9),TS(2)
DIMENSION ZN3(5),ZN2(4),SV(25)
+ CHARACTER(len=4) :: ISDATE(3),ISTIME(2),NAME(2)
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)
@@ -141,7 +142,7 @@
COMMON/PARM7/PT(150),PD(150,9),PS(150),PDL(25,2),PTL(100,4),
$ PMA(100,10),SAM(100)
COMMON/DATIM7/ISD(3),IST(2),NAM(2)
- COMMON/DATIME/ISDATE(3),ISTIME(2),NAME(2)
+ COMMON/DATIME/ISDATE,ISTIME,NAME
COMMON/CSW/SW(25),ISW,SWC(25)
COMMON/MAVG7/PAVGM(10)
COMMON/DMIX/DM04,DM16,DM28,DM32,DM40,DM01,DM14
@@ -156,11 +157,14 @@
IF(ISW.NE.64999) CALL TSELEC(SV)
C Put identification data into common/datime/
DO 1 I=1,3
- ISDATE(I)=ISD(I)
+C ISDATE(I)=ISD(I)
+ write(ISDATE(I),'(A)') ISD(I)
1 CONTINUE
DO 2 I=1,2
- ISTIME(I)=IST(I)
- NAME(I)=NAM(I)
+C ISTIME(I)=IST(I)
+ write(ISTIME(I),'(A)') IST(I)
+C NAME(I)=NAM(I)
+ write(NAME(I),'(A)') NAM(I)
2 CONTINUE
C
C Test for changed input
@@ -574,7 +578,7 @@
COMMON/GTS3C/TLB,S,DB04,DB16,DB28,DB32,DB40,DB48,DB01,ZA,T0,Z0
& ,G0,RL,DD,DB14,TR12
COMMON/MESO7/TN1(5),TN2(4),TN3(5),TGN1(2),TGN2(2),TGN3(2)
- DIMENSION D(9),T(2),MT(11),AP(1),ALTL(8)
+ DIMENSION D(9),T(2),MT(11),AP(200),ALTL(8)
COMMON/LOWER7/PTM(10),PDM(10,8)
COMMON/PARM7/PT(150),PD(150,9),PS(150),PDL(25,2),PTL(100,4),
$ PMA(100,10),SAM(100)
@@ -963,7 +967,7 @@
C CALCULATE G(L) FUNCTION
C Upper Thermosphere Parameters
REAL LAT, LONG
- DIMENSION P(1),SV(25),AP(1)
+ DIMENSION P(200),SV(25),AP(200)
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,
@@ -1217,7 +1221,7 @@
COMMON/LPOLY/PLG(9,4),CTLOC,STLOC,C2TLOC,S2TLOC,C3TLOC,S3TLOC,
$ IYR,DAY,DF,DFA,APD,APDF,APT(4),LONG
COMMON/CSW/SW(25),ISW,SWC(25)
- DIMENSION P(1),T(14)
+ DIMENSION P(200),T(14)
SAVE
DATA DR/1.72142E-2/,DGTR/1.74533E-2/,PSET/2./
DATA DAYL/-1./,P32,P18,P14,P39/4*-1000./
@@ -1651,6 +1655,7 @@
C-----------------------------------------------------------------------
BLOCK DATA GTD7BK
C MSISE-00 01-FEB-02
+ CHARACTER(len=4) :: ISDATE,ISTIME,NAME
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),
@@ -2434,7 +2439,7 @@
M 2.23000E+02, 2.86760E+02,-2.93940E+00, 2.50000E+00, 0.00000E+00/
END
-$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
+
C TEST DRIVER FOR GTD7 (ATMOSPHERIC MODEL)
DIMENSION D(9,16),T(2,16),SW(25),APH(7)
|
