Changeset 10456 for NEMO/branches/UKMO/dev_10448_WAD_SBC_BUGFIX/src/OCE/DYN/wet_dry.F90

Timestamp:

2019-01-06T17:17:13+01:00 (5 years ago)

Author:

deazer

Message:

Added option to taper sbc fluxes near very shallow water when using WAD
Corrected some IO bugs in dia25h, diatmb for WAD case.
User has control of the tapering. At what depth to start it, and at what fraction to start
the tanh tapering. At the WAD limit SBC is turned off completely.
Dry cells do not have any communication with the atmosphere
To DO: Documentation update.
Although not all sette tests are passed (AGRIF etc.)
it does no worse than the trunk at the revision the branch is made

File:

• NEMO/branches/UKMO/dev_10448_WAD_SBC_BUGFIX/src/OCE/DYN/wet_dry.F90 (modified) (4 diffs)

NEMO/branches/UKMO/dev_10448_WAD_SBC_BUGFIX/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