//LEVITUS JOB NOTIFY=6461 /*SETUP SLOT=P5271 VOLUME=NC5271 NOWRITE /*SERVICE NONPRIME /*JOBPARM S=349,L=999,R=1024,PRINT=ALL //STEP1 EXEC FORTVCLG,REGION=1024K //FORT.SYSIN DD * C----------------------------------------------------------------------- C*********************************************************************** C * C PROGRAM NAME : LEVITUS * C FUNCTION : GENERATION OF A PARTICULAR MAP OF DYNAMIC * C TOPOGRAPHY OUT OF THE 32 AVAILABLE * C REFERENCE LEVELS. * C COMPILER : VS FORTRAN VERSION 2 * C AUTHOR : NICK CHRISTOU * C HISTORY : OCTOBER 3, 1985 - VERSION 1.0 * C*********************************************************************** C----------------------------------------------------------------------- C * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * C* * C* SPECIFY BELOW WHICH LEVEL YOU WANT * C* * C* TO BE PLOTTED OUT OF THE 32 * C* * C* I.E.: LEVEL = 26 * C * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * IMPLICIT REAL*8 (A-H,O-Z) REAL*4 DTOP DIMENSION D(64800) DATA LEVEL/26/ C--------------------------------------------------------------- C FIND THE LAST AND FIRST RECORD NUMBER OF THE LEVEL SPECIFIED C ABOVE IN THE FILE CONTAINING THE DYNAMIC TOPOGRAPHY FOR ALL C 32 REFERENCE LEVELS. C--------------------------------------------------------------- LL = LEVEL * 64802 N =(( LEVEL - 1 ) * 64802) + 2 C--------------------------------------------------------------- C PUT THE DYNAMIC TOPOGRAPHY 1X1 DEGREE VALUES OF THE C SPECIFIED LEVEL IN THE DOUBLE PRECISION ARRAY D(K) C READ FROM UNIT 12 C--------------------------------------------------------------- DO 10 L=1,LL READ(21,1000) DTOP IF(L.LE.N) GO TO 10 K = L - N IF(DTOP.LT.-1000.0) DTOP=-9.9 D(K)=DBLE(DTOP) 10 CONTINUE C--------------------------------------------------------------- C GENERATE THE GRID COORDINATES AND ASSOCIATE THE ELEMENTS C OF THE D(K) ARRAY TO THE CORRESPONDING GRID NODE C--------------------------------------------------------------- DO 30 J=1,180 DLAT=-89.5D0+DFLOAT((J-1)) DO 20 I=1,360 DLON=0.5D0+DFLOAT((I-1)) M=I+((J-1)*360) C-------------------------------------------------------------- C IF YOU WISH TO OBTAIN THE DYNAMIC TOPOGRAPHY OF A SMALLER C AREA REMOVE COMMENT CARDS AND SPECIFY LATITUDE AND C LONGITUDE BOUNDARIES IN THE TWO STATEMENTS BELOW C-------------------------------------------------------------- C--------- C IF(DLAT.LT.35.0D0.OR.DLAT.GT.70.0D0) GO TO 20 C IF(DLON.LT.260.0D0.OR.DLON.GT.350.0D0) GO TO 20 C--------- WRITE(6,9000) DLAT,DLON,D(M) 20 CONTINUE 30 CONTINUE C-------------------------------------------------------------- 1000 FORMAT(E12.4) 9000 FORMAT(3F12.4) STOP END //GO.FT21F001 DD UNIT=3480,DSN=LEV32.GRID1X1.DEG,VOL=SER=NC5271, // LABEL=(5,SL),DISP=OLD, // DCB=(RECFM=FB,LRECL=12,BLKSIZE=32640) //GO.SYSIN DD * //