Changeset 1061

Timestamp:

Dec 16, 2008, 7:22:11 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Ajout limite inferieure vent calcul diagnostic cdrag
FH/IM

File:

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

LMDZ4/trunk/libf/phylmd/coefcdrag.F90

@@ -1 @@
-!
-! $Header$
-!
 !
 !
@@ -54 +51 @@
       include "indicesol.h"
 ! Quelques constantes :
-      REAL, parameter :: RKAR=0.40, CB=5.0, CC=5.0, CD=5.0
+      REAL, parameter :: RKAR=0.40, CB=5.0, CC=5.0, CD=5.0, cepdu2=(0.1)**2
 !
 ! Variables locales :
@@ -72 +69 @@
 !
        zdphi(i) = zgeop(i)
-       zdu2(i) = speed(i)**2
+       zdu2(i) = max(cepdu2,speed(i)**2)
        pref(i) = exp(log(psol(i)) - zdphi(i)/(RD*t(i)* &
                  (1.+ RETV * max(q(i),0.0))))