Changeset 1269

Timestamp:

2009-01-14T17:16:05+01:00 (15 years ago)

Author:

ctlod

Message:

limit the heat transfer coefficient zchcm value, see ticket: #285

File:

• trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90

@@ -256 +256 @@
       REAL(wp) ::   zdtetar, ztvmoyr, zlxins, zchcm, zclcm      !    -         -
       REAL(wp) ::   zmt1, zmt2, zmt3, ztatm3, ztamr, ztaevbk    !    -         -
-      REAL(wp) ::   zsst, ztatm, zcco1, zpatm                   !    -         -
+      REAL(wp) ::   zsst, ztatm, zcco1, zpatm, zcmax, zrmax     !    -         -
       REAL(wp) ::   zrhoa, zev, zes, zeso, zqatm, zevsqr        !    -         -
       !!
@@ -364 +364 @@
             zcln           = 0.0346 * zcmn
             zcmcmn         = 1. / ( 1. - zcmn * zpsim / vkarmn )
-            zchcm          = zcmcmn / ( 1. - zchn * zpsih / ( vkarmn * zcmn ) )
+            ! sometimes the ratio zchn * zpsih / ( vkarmn * zcmn ) is too close to 1 and zchcm becomes very very big
+            zcmax = 0.1               ! choice for maximum value of the heat transfer coefficient, guided by my intuition
+            zrmax = 1 - 3.e-4 / zcmax ! maximum value of the ratio
+            zchcm = zcmcmn / ( 1. - MIN ( zchn * zpsih / ( vkarmn * zcmn ) , zrmax ) )
             zclcm          = zchcm
             !                                                          ! transfert coef. (Large and Pond 1981,1982)