Changeset 250

Timestamp:

2005-03-30T12:31:27+02:00 (19 years ago)

Author:

opalod

Message:

CT : UPDATE173 : Allow the possibility to have constant thickness levels for the vertical grid

File:

• trunk/NEMO/OPA_SRC/DOM/domzgr.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -173 +173 @@
          WRITE(numout,*) '    zgr_z   : Reference vertical z-coordinates'
          WRITE(numout,*) '    ~~~~~~~'
-         IF ( ppa1 == 0. .AND. ppa0 == 0. .AND. ppsur == 0. ) THEN
-            WRITE(numout,*) '         zsur, za0, za1 computed from '
-            WRITE(numout,*) '                 zdzmin = ', zdzmin
-            WRITE(numout,*) '                 zhmax  = ', zhmax
-         ENDIF
-         WRITE(numout,*) '              Namelist namzgr : value of coefficients for vertical mesh:'
-         WRITE(numout,*) '                 zsur = ', zsur
-         WRITE(numout,*) '                 za0  = ', za0
-         WRITE(numout,*) '                 za1  = ', za1
-         WRITE(numout,*) '                 zkth = ', zkth
-         WRITE(numout,*) '                 zacr = ', zacr
+         IF (  ppkth == 0. ) THEN              
+              WRITE(numout,*) '            Uniform grid with ',jpk-1,' layers'
+              WRITE(numout,*) '            Total depth    :', zhmax
+              WRITE(numout,*) '            Layer thickness:', zhmax/(jpk-1)
+         ELSE
+            IF ( ppa1 == 0. .AND. ppa0 == 0. .AND. ppsur == 0. ) THEN
+               WRITE(numout,*) '         zsur, za0, za1 computed from '
+               WRITE(numout,*) '                 zdzmin = ', zdzmin
+               WRITE(numout,*) '                 zhmax  = ', zhmax
+            ENDIF
+            WRITE(numout,*) '           Value of coefficients for vertical mesh:'
+            WRITE(numout,*) '                 zsur = ', zsur
+            WRITE(numout,*) '                 za0  = ', za0
+            WRITE(numout,*) '                 za1  = ', za1
+            WRITE(numout,*) '                 zkth = ', zkth
+            WRITE(numout,*) '                 zacr = ', zacr
+         ENDIF
       ENDIF
 
@@ -189 +195 @@
       ! Reference z-coordinate (depth - scale factor at T- and W-points)
       ! ======================
-
-      DO jk = 1, jpk
-         zw = FLOAT( jk )
-         zt = FLOAT( jk ) + 0.5
-         gdepw(jk) = ( zsur + za0 * zw + za1 * zacr * LOG( COSH( (zw-zkth)/zacr ) )  )
-         gdept(jk) = ( zsur + za0 * zt + za1 * zacr * LOG( COSH( (zt-zkth)/zacr ) )  )
-         e3w  (jk) =          za0      + za1        * TANH(      (zw-zkth)/zacr   )
-         e3t  (jk) =          za0      + za1        * TANH(      (zt-zkth)/zacr   )
-      END DO
-      gdepw(1) = 0.e0   ! force first w-level to be exactly at zero
-
+      IF (  ppkth == 0. ) THEN            !  uniform vertical grid       
+
+         za1 = zhmax/FLOAT(jpk-1) 
+         DO jk = 1, jpk
+            zw = FLOAT( jk )
+            zt = FLOAT( jk ) + 0.5
+            gdepw(jk) = ( zw - 1 ) * za1
+            gdept(jk) = ( zt - 1 ) * za1
+            e3w  (jk) =  za1
+            e3t  (jk) =  za1
+         END DO
+
+      ELSE
+
+         DO jk = 1, jpk
+            zw = FLOAT( jk )
+            zt = FLOAT( jk ) + 0.5
+            gdepw(jk) = ( zsur + za0 * zw + za1 * zacr * LOG( COSH( (zw-zkth)/zacr ) )  )
+            gdept(jk) = ( zsur + za0 * zt + za1 * zacr * LOG( COSH( (zt-zkth)/zacr ) )  )
+            e3w  (jk) =          za0      + za1        * TANH(      (zw-zkth)/zacr   )
+            e3t  (jk) =          za0      + za1        * TANH(      (zt-zkth)/zacr   )
+         END DO
+         gdepw(1) = 0.e0   ! force first w-level to be exactly at zero
+
+      ENDIF
 
       ! Control and  print