Changeset 497 for trunk/NEMO/OFF_SRC/DOM/domain.F90

Timestamp:

2006-09-12T13:03:53+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_064:CE:re-organization of coordinate definition and scale factors

File:

• trunk/NEMO/OFF_SRC/DOM/domain.F90 (modified) (5 diffs)

trunk/NEMO/OFF_SRC/DOM/domain.F90

@@ -29 +29 @@
    PUBLIC dom_init       ! called by opa.F90
 
+   !! * Module variables
+      REAL(wp) ::          & !!: Namelist nam_zgr_sco
+      sbot_min =  300.  ,  &  !: minimum depth of s-bottom surface (>0) (m)
+      sbot_max = 5250.  ,  &  !: maximum depth of s-bottom surface (= ocean depth) (>0) (m)
+      theta    =    6.0 ,  &  !: surface control parameter (0<=theta<=20)
+      thetb    =    0.75,  &  !: bottom control parameter  (0<=thetb<= 1)
+      r_max    =    0.15      !: maximum cut-off r-value allowed (0<r_max<1)
+
+
    !! * Substitutions
 #  include "domzgr_substitute.h90"
@@ -98 +107 @@
       !! * Modules used
       USE ioipsl
-      NAMELIST/namrun/ no    , cexper   , ln_rstart , nrstdt , nit000,          &
+      INTEGER ::   ioptio = 0      ! temporary integer
+
+      NAMELIST/nam_run/ no    , cexper   , ln_rstart , nrstdt , nit000,          &
          &             nitend, ndate0   , nleapy   , ninist , nstock,           &
          &             nprint, nwrite   , nrunoff  , ln_ctl , nictls, nictle,   &
          &             njctls, njctle   , nbench   , isplt  , jsplt
 
-      NAMELIST/namdom/ e3zps_min, e3zps_rat, nmsh  ,   &
+      NAMELIST/nam_zgr/ ln_zco, ln_zps, ln_sco
+
+      NAMELIST/nam_dom/ e3zps_min, e3zps_rat, nmsh  ,   &
          &             nacc  , atfp     , rdt      , rdtmin , rdtmax,   &
          &             rdth 
 
-      NAMELIST/namcla/ n_cla
+      NAMELIST/nam_cla/ n_cla
       !!----------------------------------------------------------------------
 
@@ -118 +131 @@
       ! Namelist namrun : parameters of the run
       REWIND( numnam )
-      READ  ( numnam, namrun )
+      READ  ( numnam, nam_run )
 
       IF(lwp) THEN
@@ -256 +269 @@
       ENDIF
 
+      ! Read Namelist nam_zgr : vertical coordinate'
+      ! ---------------------
+      REWIND ( numnam )
+      READ   ( numnam, nam_zgr )
+
+      ! Parameter control and print
+      ! ---------------------------
+      ! Control print
+      IF(lwp) THEN
+         WRITE(numout,*)
+         WRITE(numout,*) 'Namelist namzgr : vertical coordinate'
+         WRITE(numout,*) '~~~~~~~'
+         WRITE(numout,*) '          Namelist nam_zgr : set vertical coordinate'
+         WRITE(numout,*) '             z-coordinate - full steps      ln_zco = ', ln_zco
+         WRITE(numout,*) '             z-coordinate - partial steps   ln_zps = ', ln_zps
+         WRITE(numout,*) '             s- or hybrid z-s-coordinate    ln_sco = ', ln_sco
+      ENDIF
+
+      ! Check Vertical coordinate options
+      ioptio = 0
+      IF( ln_zco ) ioptio = ioptio + 1
+      IF( ln_zps ) ioptio = ioptio + 1
+      IF( ln_sco ) ioptio = ioptio + 1
+      IF ( ioptio /= 1 ) CALL ctl_stop( ' none or several vertical coordinate options used' )
+
+      IF( ln_zco ) THEN
+          IF(lwp) WRITE(numout,*) '          z-coordinate with reduced incore memory requirement'
+          IF( ln_zps .OR. ln_sco ) CALL ctl_stop( ' reduced memory with zps or sco option is impossible' )
+      ENDIF
+
+
       ! Namelist namdom : space/time domain (bathymetry, mesh, timestep)
       REWIND( numnam )
-      READ  ( numnam, namdom )
+      READ  ( numnam, nam_dom )
 
       IF(lwp) THEN
@@ -285 +329 @@
       ! Namelist cross land advection
       REWIND( numnam )
-      READ  ( numnam, namcla )
+      READ  ( numnam, nam_cla )
       IF(lwp) THEN
          WRITE(numout,*)