Changeset 1554 for trunk/NEMO/OPA_SRC/SBC

Timestamp:

2009-07-29T12:30:10+02:00 (15 years ago)

Author:

ctlod

Message:

add the possibility to bound the sea surface salinity damping term erp(:,:), see ticket: #499

File:

• trunk/NEMO/OPA_SRC/SBC/sbcssr.F90 (modified) (5 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcssr.F90

@@ -35 +35 @@
    INTEGER ::   nn_sstr, nn_sssr   ! SST/SSS restoring indicator
    REAL(wp) ::  dqdt   , deds      ! restoring factor on SST and SSS
+   LOGICAL ::   ln_sssr_bnd        ! flag to bound erp term 
+   REAL(wp) ::  rn_sssr_bnd        ! ABS(Max./Min.) value of erp term [mm/day]
 
    REAL(wp) , ALLOCATABLE, DIMENSION(:) ::   buffer   ! Temporary buffer for exchange
@@ -70 +72 @@
       REAL(wp) ::   zqrp     ! local scalar for heat flux damping
       REAL(wp) ::   zsrp     ! local scalar for unit conversion of deds factor
+      REAL(wp) ::   zerp_bnd ! local scalar for unit conversion of rn_epr_max factor
       INTEGER  ::   ierror   ! return error code
       !!
       CHARACTER(len=100) ::  cn_dir          ! Root directory for location of ssr files
       TYPE(FLD_N) ::   sn_sst, sn_sss        ! informations about the fields to be read
-      NAMELIST/namsbc_ssr/ cn_dir, nn_sstr, nn_sssr, dqdt, deds, sn_sst, sn_sss
+      NAMELIST/namsbc_ssr/ cn_dir, nn_sstr, nn_sssr, dqdt, deds, sn_sst, sn_sss, ln_sssr_bnd, rn_sssr_bnd
       !!----------------------------------------------------------------------
       !
@@ -103 +106 @@
             WRITE(numout,*) '          dQ/dT (restoring magnitude on SST)     dqdt    = ', dqdt, ' W/m2/K'
             WRITE(numout,*) '          dE/dS (restoring magnitude on SST)     deds    = ', deds, ' mm/day'
+            WRITE(numout,*) '          flag to bound erp term             ln_sssr_bnd = ', ln_sssr_bnd
+            WRITE(numout,*) '          ABS(Max./Min.) erp threshold       rn_sssr_bnd = ', rn_sssr_bnd, ' mm/day'
          ENDIF
 
@@ -172 +177 @@
             ELSEIF( nn_sssr == 2 ) THEN               !* Salinity damping term (volume flux, emp and emps)
                zsrp = deds / rday                                     ! from [mm/day] to [kg/m2/s]
+               zerp_bnd = rn_sssr_bnd / rday                          !       -              -    
 !CDIR COLLAPSE
                DO jj = 1, jpj
@@ -178 +184 @@
                         &        * ( sss_m(ji,jj) - sf_sss(1)%fnow(ji,jj) )   &
                         &        / ( sss_m(ji,jj) + 1.e-20   )
+                     IF( ln_sssr_bnd )   zerp = SIGN( 1., zerp ) * MIN( zerp_bnd, ABS(zerp) )
                      emp (ji,jj) = emp (ji,jj) + zerp
                      emps(ji,jj) = emps(ji,jj) + zerp