Changeset 467 for LMDZ.3.3/branches/rel-LF/libf/phylmd/clcdrag.F90

Timestamp:

Aug 6, 2003, 4:50:49 PM (21 years ago)

Author:

lmdzadmin

Message:

Modifs sur les seuils (cdrag etc...), inclusion des diagnostics ISCCP par Ionela
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/clcdrag.F90 (modified) (5 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/clcdrag.F90

@@ -40 +40 @@
 !
 ! Quelques constantes et options:
-      REAL, PARAMETER :: ckap=0.35, cb=5.0, cc=5.0, cd=5.0, cepdu2=(0.1)**2
+!!$PB      REAL, PARAMETER :: ckap=0.35, cb=5.0, cc=5.0, cd=5.0, cepdu2=(0.1)**2
+      REAL, PARAMETER :: ckap=0.40, cb=5.0, cc=5.0, cd=5.0, cepdu2=(0.1)**2
 !
 ! Variables locales :
@@ -68 +69 @@
         zcdn(i) = (ckap/log(1.+zgeop(i)/(RG*rugos(i))))**2
 !
-        IF (zri(i) .ge. 0.) THEN      ! situation stable
+!!$        IF (zri(i) .ge. 0.) THEN      ! situation stable
+        IF (zri(i) .gt. 0.) THEN      ! situation stable
           zri(i) = min(20.,zri(i))
           IF (.NOT.zxli) THEN
@@ -75 +77 @@
             zcfm1(i) = zcdn(i) * FRIV
             FRIH = AMAX1(1./ (1.+3.*CB*zri(i)*ZSCF), 0.1 )
-            zcfh1(i) = zcdn(i) * FRIH
+!!$  PB          zcfh1(i) = zcdn(i) * FRIH
+            zcfh1(i) = 0.8 * zcdn(i) * FRIH
             pcfm(i) = zcfm1(i)
             pcfh(i) = zcfh1(i)
@@ -87 +90 @@
                  *(1.0+zgeop(i)/(RG*rugos(i)))))
             zcfm2(i) = zcdn(i)*amax1((1.-2.0*cb*zri(i)*zucf),0.1)
-            zcfh2(i) = zcdn(i)*amax1((1.-3.0*cb*zri(i)*zucf),0.1)
+!!$PB            zcfh2(i) = zcdn(i)*amax1((1.-3.0*cb*zri(i)*zucf),0.1)
+            zcfh2(i) = 0.8 * zcdn(i)*amax1((1.-3.0*cb*zri(i)*zucf),0.1)
             pcfm(i) = zcfm2(i)
             pcfh(i) = zcfh2(i)
@@ -94 +98 @@
             pcfh(i) = zcdn(i)* fins(zri(i))
           ENDIF
-          zcr = (0.0016/(zcdn(i)*SQRT(zdu2)))*ABS(ztvd-ztsolv)**(1./3.)
-          IF(nsrf.EQ.is_oce) pcfh(i) = zcdn(i)*(1.0+zcr**1.25)**(1./1.25)
+                  zcr = (0.0016/(zcdn(i)*SQRT(zdu2)))*ABS(ztvd-ztsolv)**(1./3.)
+          IF(nsrf.EQ.is_oce) pcfh(i) =0.8* zcdn(i)*(1.0+zcr**1.25)**(1./1.25)
         ENDIF
       END DO