Changeset 2452 for LMDZ5/trunk/libf/phylmd/screenc.F90

Timestamp:

Feb 24, 2016, 6:21:32 PM (9 years ago)

Author:

musat

Message:

Correction bug cas Dice:

• remplace yfluxsens par y_flux_t1 (pbl_surface_mod)
• ajout IF(ok_prescr_ust) pour calcul cdram (screenc)

MPL/JYG/FCheruy

File:

• LMDZ5/trunk/libf/phylmd/screenc.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/screenc.F90

@@ -54 +54 @@
 !-----------------------------------------------------------------------
       include "YOMCST.h"
+      include "flux_arp.h"
 !
 ! Variables locales  
       INTEGER :: i 
-      REAL, dimension(klon) :: cdram, cdrah, cdran, zri1, gref
+      REAL, dimension(klon) :: cdram, cdrah, cdran, zri1, gref,ycdragm
 !
 !------------------------------------------------------------------------- 
@@ -78 +79 @@
                     cdram, cdrah, zri1, pref)
       DO i = 1, knon
+        IF(ok_prescr_ust) THEN
+! La aussi il faut forcer avec ust (FC + MPL 20160210)
+        ycdragm(i) = ust*ust/(1.+speed(i))/speed(i)
+        cdram=ycdragm
+        delu(i) = ust/sqrt(cdram(i))
+        ELSE
         delu(i) = ustar(i)/sqrt(cdram(i))
+        ENDIF
         delte(i)= (testar(i)* sqrt(cdram(i)))/ &
                    cdrah(i)