Changeset 2126 for LMDZ5/trunk/libf/phylmd/coefcdrag.F90

Timestamp:

Oct 15, 2014, 2:03:57 AM (10 years ago)

Author:

fhourdin

Message:

Introduction du cas Dice couplé atmosphere/surface
+ nouveau paramètre de contrôle f_ri_cd_min, seuil minimum
sur la fonction f(Ri) en facteur du coefficient de traîné neutre.
Par défaut : f_ri_cd_min=0.1 (comme avant)

Introduction of the coupled atmosphere/surface dice case.
+ new parameter f_ri_cd_min, minimum threshold on the f(Ri) fonction
that multiply the neutral drag coefficient.

File:

• LMDZ5/trunk/libf/phylmd/coefcdrag.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/coefcdrag.F90

@@ -63 +63 @@
       REAL, dimension(klon) :: zcfm2, zcfh2
       REAL, dimension(klon) :: trm0, trm1
+
 !-------------------------------------------------------------------------
       REAL :: fsta, fins, x
@@ -110 +111 @@
          IF (.NOT.zxli) THEN
            zscf(i) = SQRT(1.+CD*ABS(zri1(i)))
-           friv(i) = max(1. / (1.+2.*CB*zri1(i)/ zscf(i)), 0.1)
+           friv(i) = max(1. / (1.+2.*CB*zri1(i)/ zscf(i)), f_ri_cd_min)
            zcfm1(i) = cdran(i) * friv(i)
-           frih(i) = max(1./ (1.+3.*CB*zri1(i)*zscf(i)), 0.1 )
+           frih(i) = max(1./ (1.+3.*CB*zri1(i)*zscf(i)), f_ri_cd_min )
 !           zcfh1(i) = cdran(i) * frih(i)
            zcfh1(i) = f_cdrag_ter*cdran(i) * frih(i)
@@ -130 +131 @@
            zucf(i) = 1./(1.+3.0*CB*CC*cdran(i)*SQRT(ABS(zri1(i)) &
                  *(1.0+zdphi(i)/(RG*rugos(i)))))
-           zcfm2(i) = cdran(i)*max((1.-2.0*CB*zri1(i)*zucf(i)),0.1)
-!           zcfh2(i) = cdran(i)*max((1.-3.0*CB*zri1(i)*zucf(i)),0.1)
-           zcfh2(i) = f_cdrag_ter*cdran(i)*max((1.-3.0*CB*zri1(i)*zucf(i)),0.1)
+           zcfm2(i) = cdran(i)*max((1.-2.0*CB*zri1(i)*zucf(i)),f_ri_cd_min)
+!           zcfh2(i) = cdran(i)*max((1.-3.0*CB*zri1(i)*zucf(i)),f_ri_cd_min)
+           zcfh2(i) = f_cdrag_ter*cdran(i)*max((1.-3.0*CB*zri1(i)*zucf(i)),f_ri_cd_min)
            cdram(i) = zcfm2(i)
            cdrah(i) = zcfh2(i)