Changeset 1577 for trunk/NEMO/OPA_SRC/DOM/domzgr.F90

Timestamp:

2009-08-04T16:56:59+02:00 (15 years ago)

Author:

smasson

Message:

update diahth and zdfmxl, see ticket:468

File:

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

trunk/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -160 +160 @@
       REAL(wp) ::   zsur, za0, za1, zkth   ! Values set from parameters in
       REAL(wp) ::   zacr, zdzmin, zhmax    ! par_CONFIG_Rxx.h90
+      REAL(wp) ::   zrefdep                ! depth of the reference level (~10m)
       !!----------------------------------------------------------------------
 
@@ -230 +231 @@
          gdepw_0(1) = 0.e0                     ! force first w-level to be exactly at zero
       ENDIF
+
+      ! deepest/shallowest W level Above/Bellow ~10m
+      zrefdep = 10. - ( 0.1*MINVAL(e3w_0) )                          ! ref. depth with tolerance (10% of minimum layer thickness)
+      nlb10 = MINLOC( gdepw_0, mask = gdepw_0 > zrefdep, dim = 1 )   ! shallowest W level Bellow ~10m
+      nla10 = nlb10 - 1                                              ! deepest    W level Above  ~10m
 
       IF(lwp) THEN                        ! control print