Changeset 4355 for branches/2013/dev_r4050_NOC_WaD/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

Timestamp:

2014-01-17T18:55:02+01:00 (10 years ago)

Author:

hliu

Message:

update Wetting/drying branch. Will keep commit in the next a few days.

File:

• branches/2013/dev_r4050_NOC_WaD/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90 (modified) (2 diffs)

branches/2013/dev_r4050_NOC_WaD/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

@@ -81 +81 @@
       !
                              CALL dom_nam      ! read namelist ( namrun, namdom, namcla )
+                             CALL dom_wad      ! read namelist ( namwad ), wetting/drying
                              CALL dom_clo      ! Closed seas and lake
                              CALL dom_hgr      ! Horizontal mesh
@@ -269 +270 @@
       !
    END SUBROUTINE dom_nam
+
+   SUBROUTINE dom_wad
+      !!----------------------------------------------------------------------
+      !!                     ***  ROUTINE dom_wad  ***
+      !!                    
+      !! ** Purpose :   read wetting/drying namelists and print the variables.
+      !!
+      !! ** input   : - namwad namelist
+      !!----------------------------------------------------------------------
+      USE ioipsl
+      NAMELIST/namwad/ ln_wad , rn_wadmin, rn_wadmine, rn_landele
+      !!----------------------------------------------------------------------
+
+      REWIND( numnam )              ! Namelist namrun : parameters of the run
+      READ  ( numnam, namwad )
+      !
+      IF(lwp) THEN                  ! control print
+         WRITE(numout,*)
+         WRITE(numout,*) 'dom_wad  : set up wetting/drying '
+         WRITE(numout,*) '~~~~~~~ '
+         WRITE(numout,*) '   Namelist namwad'
+         WRITE(numout,*) '      turn on(off) wetting/drying     ln_wad        = ', ln_wad
+         WRITE(numout,*) '      minimal dry cell depth          rn_ wadmin    = ', rn_wadmin 
+         WRITE(numout,*) '      tollerance of minimual depth    rn_wadmine    = ', rn_wadmine
+         WRITE(numout,*) ' threshold elevation for land masking rn_landele    = ', rn_landele
+      ENDIF
+
+   END SUBROUTINE dom_wad