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
