Changeset 10499 for NEMO/trunk/src/OCE/DYN

Timestamp:

2019-01-10T16:12:24+01:00 (5 years ago)

Author:

deazer

Message:

Fix ticket #2154

File:

• NEMO/trunk/src/OCE/DYN/wet_dry.F90 (modified) (4 diffs)

NEMO/trunk/src/OCE/DYN/wet_dry.F90

@@ -45 +45 @@
    REAL(wp), PUBLIC  ::   r_rn_wdmin1 !: 1/minimum water depth on dried cells 
    REAL(wp), PUBLIC  ::   rn_wdmin2   !: tolerance of minimum water depth on dried cells
+   REAL(wp), PUBLIC  ::   rn_wd_sbcdep   !: Depth at which to taper sbc fluxes
+   REAL(wp), PUBLIC  ::   rn_wd_sbcfra   !: Fraction of SBC at taper depth
    REAL(wp), PUBLIC  ::   rn_wdld     !: land elevation below which wetting/drying will be considered
    INTEGER , PUBLIC  ::   nn_wdit     !: maximum number of iteration for W/D limiter
@@ -74 +76 @@
       !!
       NAMELIST/namwad/ ln_wd_il, ln_wd_dl   , rn_wdmin0, rn_wdmin1, rn_wdmin2, rn_wdld,   &
-         &             nn_wdit , ln_wd_dl_bc, ln_wd_dl_rmp
+         &             nn_wdit , ln_wd_dl_bc, ln_wd_dl_rmp, rn_wd_sbcdep,rn_wd_sbcfra
       !!----------------------------------------------------------------------
       !
@@ -85 +87 @@
       IF(lwm) WRITE ( numond, namwad )
       !
+      IF( rn_wd_sbcfra>=1 )   CALL ctl_stop( 'STOP', 'rn_wd_sbcfra >=1 : must be < 1' )
       IF(lwp) THEN                  ! control print
          WRITE(numout,*)
@@ -99 +102 @@
          WRITE(numout,*) '      T => baroclinic u,v=0 at dry pts: ln_wd_dl_bc = ', ln_wd_dl_bc     
          WRITE(numout,*) '      use a ramp for rwd limiter:  ln_wd_dl_rwd_rmp = ', ln_wd_dl_rmp
+         WRITE(numout,*) '      cut off depth sbc for wd   rn_wd_sbcdep       = ', rn_wd_sbcdep
+         WRITE(numout,*) '      fraction to start sbc wgt rn_wd_sbcfra        = ', rn_wd_sbcfra
       ENDIF
       IF( .NOT. ln_read_cfg ) THEN