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

Timestamp:

May 6, 2011, 4:45:19 PM (13 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/cv3p1_closure.F (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/cv3p1_closure.F

@@ -70 +70 @@
       real cbmflim(nloc),cbmf1(nloc),cbmfmax(nloc),cbmf(nloc)
       real cbmflast(nloc)
+      real wbeff(nloc)
+
       real coef(nloc)
       real xp(nloc),xq(nloc),xr(nloc),discr(nloc),b3(nloc),b4(nloc)
@@ -515 +517 @@
       ENDDO
 c
+      if (flag_wb==0) then
+         wbeff(:)=wbmax
+      else if (flag_wb==1) then
+         wbeff(1:ncum)=wbmax/(1.+500./(ph(1:ncum,1)-plfc(1:ncum)))
+      else
+         stop'cas non prevu iflag_wb'
+! ligne a suprimer un jour ...
+      endif
+!     if (ncum==1) print*,'WB ',plfc,wbeff,alp2,cin
+
       DO il = 1,ncum
-cjyg    Modification du coef de wb*wb pour conformite avec papier Wake
-cc       cbmf1(il) = alp2(il)/(0.5*wb*wb-Cin(il))
-       cbmf1(il) = alp2(il)/(2.*wb*wb-Cin(il))
+
+!jyg    Modification du coef de wb*wb pour conformite avec papier Wake
+!       cbmf1(il) = alp2(il)/(0.5*wb*wb-Cin(il))
+! A partir de mai 2011 on revient au 2 wb^2, conformement au papier.
+! les versions d'avant (NPv2 et al) tournaient avec 0.5 wb^2 et wb=1m/s
+! equivalent  a wb=1 m/s pour la nouvelle formulation
+! Puis passage de wb constant a wbeff variable suivant flag_wb
+
+       cbmf1(il) = alp2(il)/(2.*wbeff(il)*wbeff(il)-Cin(il))
+
        if(cbmf1(il).EQ.0.AND.alp2(il).NE.0.) THEN
         write(lunout,*)