C********************************************************************* C********************************************************************* C SUBROUTINE SETGRD O ( ZGRID, FDZ, FDZM, DX, D NZ, NZB ) C C********************************************************************* IMPLICIT REAL*8 ( A-H, O-Z ) C********************************************************************* C DIMENSION ZLV(0:100) , ZGN(-10:100) C DIMENSION ZGRID(-2:NZB), FDZ(-2:NZB), FDZM(-2:NZB) C C--------------------------------------------------------------------- C DATA ( ZLV(ILV), ILV=0,39 ) 1 / 0.00, 0.10, 0.22, 0.36, 0.52, 0.70, 2 0.90, 1.14, 1.42, 1.74, 2.16, 2.58, 3 3.06, 3.60, 4.20, 4.80, 5.40, 2 6.0, 6.6, 7.2, 7.8, 8.4, 3 9.0, 9.6, 10.2, 10.8, 11.4, 4 12.0, 12.7, 13.4, 14.2, 15.1, 5 16.0, 17.1, 18.3, 19.6, 21.0, 6 22.6, 24.2, 26.2 / C C--------------------------------------------------------------------- C DX = 2000.0D0 C C DO 11 IZ = 0, NZ ZGN ( IZ ) = 1000.0D0 * ZLV(IZ) 11 CONTINUE C ZGN(NZ+1) = 2.0D0 * ZGN(NZ) - ZGN(NZ-1) ZGN(NZ+2) = 2.0D0 * ZGN(NZ) - ZGN(NZ-2) ZGN(NZ+3) = 2.0D0 * ZGN(NZ) - ZGN(NZ-3) ZGN( -1) = - ZGN( 1) ZGN( -2) = - ZGN( 2) C DO 12 IZ = 0, NZ+2 ZGRID(IZ) = ( ZGN(IZ) + ZGN(IZ+1) ) / 2.0D0 12 CONTINUE C ZGRID(-1) = - ZGRID(0) ZGRID(-2) = - ZGRID(1) C DO 21 IZ = 0, NZ-1 FDZ(IZ) = ( ZGN(IZ+1) - ZGN(IZ) ) / DX 21 CONTINUE C FDZ(-1) = FDZ(0) FDZ(-2) = FDZ(1) FDZ(NZ) = FDZ(NZ-1) FDZ(NZ+1) = FDZ(NZ-2) FDZ(NZ+2) = FDZ(NZ-3) C C DO 22 IZ = 0, NZ FDZM(IZ) = ( FDZ(IZ-1) + FDZ(IZ) ) / 2.0D0 22 CONTINUE C FDZM(-2) = FDZM(2) FDZM(-1) = FDZM(1) FDZM(NZ+1) = FDZM(NZ-1) FDZM(NZ+2) = FDZM(NZ-2) C C===================================================================== C RETURN E N D