Changeset 1061 for LMDZ4/trunk/libf


Ignore:
Timestamp:
Dec 16, 2008, 7:22:11 PM (16 years ago)
Author:
Laurent Fairhead
Message:

Ajout limite inferieure vent calcul diagnostic cdrag
FH/IM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/coefcdrag.F90

    r793 r1061  
    1 !
    2 ! $Header$
    3 !
    41!
    52!
     
    5451      include "indicesol.h"
    5552! Quelques constantes :
    56       REAL, parameter :: RKAR=0.40, CB=5.0, CC=5.0, CD=5.0
     53      REAL, parameter :: RKAR=0.40, CB=5.0, CC=5.0, CD=5.0, cepdu2=(0.1)**2
    5754!
    5855! Variables locales :
     
    7269!
    7370       zdphi(i) = zgeop(i)
    74        zdu2(i) = speed(i)**2
     71       zdu2(i) = max(cepdu2,speed(i)**2)
    7572       pref(i) = exp(log(psol(i)) - zdphi(i)/(RD*t(i)* &
    7673                 (1.+ RETV * max(q(i),0.0))))
Note: See TracChangeset for help on using the changeset viewer.