Changeset 1516 for LMDZ5/trunk/libf/phylmd/physiq.F

Timestamp:

May 6, 2011, 4:45:19 PM (14 years ago)

Author:

idelkadi

Message:

• Modification dans physiq.F pour que les pluies stratiformes n'alimentent les pochers que si le spoches occupent plus de 10% de la surface.
• Dans le cas iflag_wake=2, on ne prend en compte les tendances de fisrtilp pour forcer les poches que si les poches existent deja (en controlant si elles couvrent deja 10% de la maille)
• On autorise differentes formulations pour la vitesse a la base de la convection WB (plus precisement au niveau de convection libre LFC, en haut de la zone d'inhibition) en fonction de flag_wb, lu dans conv_param.data

flag_wb=0 : wb=wbmax, wbmax etant lu dans conv_param.data
flag_wb=1 : wb=wbmax/[1+ 500hPa / (ps-p_LFC) ]
Concerne cv3_routines.F cv3param.h cv3p1_closure.F

File:

• LMDZ5/trunk/libf/phylmd/physiq.F (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -615 +615 @@
       real, save :: alp_max=2.
 
+      real, save :: wake_s_min_lsp=0.1
+
 c$OMP THREADPRIVATE(alp_bl_prescr,ale_bl_prescr)
 c$OMP THREADPRIVATE(ale_max,alp_max)
+c$OMP THREADPRIVATE(wake_s_min_lsp)
 
       real ale_wake(klon)
       real alp_wake(klon)
+
+      real ok_wk_lsp(klon)
+
 cRC
 c Variables liées à la poche froide (jyg et rr)
@@ -2363 +2369 @@
 
       if (iflag_wake==2) then
-         dt_dwn(:,:)= dt_dwn(:,:)+(d_t_eva(:,:)+d_t_lsc(:,:))/dtime
-         dq_dwn(:,:)= dq_dwn(:,:)+(d_q_eva(:,:)+d_q_lsc(:,:))/dtime
+        ok_wk_lsp(:)=max(sign(1.,wake_s(:)-wake_s_min_lsp),0.)
+        DO k = 1,klev
+         dt_dwn(:,k)= dt_dwn(:,k)+
+     :            ok_wk_lsp(:)*(d_t_eva(:,k)+d_t_lsc(:,k))/dtime
+         dq_dwn(:,k)= dq_dwn(:,k)+
+     :            ok_wk_lsp(:)*(d_q_eva(:,k)+d_q_lsc(:,k))/dtime
+        ENDDO
       endif
 c