Changeset 2452

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

Location:

LMDZ5/trunk/libf/phylmd

Files:

• pbl_surface_mod.F90 (modified) (1 diff)
• screenc.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -1915 +1915 @@
 
           DO j = 1, knon
-            yt1_new=(1./RCPD)*(AcoefH(j)+BcoefH(j)*yfluxsens(j)*dtime)
-            ytsurf_new(j)=yt1_new-yfluxsens(j)/(Kech_h(j)*RCPD)
+            yt1_new=(1./RCPD)*(AcoefH(j)+BcoefH(j)*y_flux_t1(j)*dtime)
+            ytsurf_new(j)=yt1_new-y_flux_t1(j)/(Kech_h(j)*RCPD)
           ENDDO
 

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)