Changeset 1516

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

Location:

LMDZ5/trunk/libf/phylmd

Files:

• cv3_routines.F (modified) (4 diffs)
• cv3p1_closure.F (modified) (2 diffs)
• cv3param.h (modified) (3 diffs)
• physiq.F (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/cv3_routines.F

@@ -36 +36 @@
       real delt ! timestep (seconds)
 
+
       CHARACTER (LEN=20) :: modname='cv3_param'
       CHARACTER (LEN=80) :: abort_message
@@ -74 +75 @@
 ccc      dttrig = 5.   ! (loose) condition for triggering
       dttrig = 10.   ! (loose) condition for triggering
-      wb = 1.   ! (m/s) adiab updraught speed at LFC (used in cv3p1_closure)
+      flag_wb=0
+      wbmax = 1.   ! (m/s) adiab updraught speed at LFC (used in cv3p1_closure)
 
 c -- rate of approach to quasi-equilibrium:
@@ -97 +99 @@
       READ(99,*,end=9998) spfac
       READ(99,*,end=9998) tau
-      READ(99,*,end=9998) wb
+      READ(99,*,end=9998) flag_wb
+      READ(99,*,end=9998) wbmax
 9998  Continue
       CLOSE(99)
@@ -107 +110 @@
         WRITE(*,*)'spfac=',spfac
         WRITE(*,*)'tau=',tau
-        WRITE(*,*)'wb =',wb
+        WRITE(*,*)'flag_wb =',flag_wb
+        WRITE(*,*)'wbmax =',wbmax
 
 cIM Lecture du fichier ep_param.data

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,*)

LMDZ5/trunk/libf/phylmd/cv3param.h

@@ -9 +9 @@
       integer noff, minorig, nl, nlp, nlm
       real sigdz, spfac
-      integer flag_epKEorig
+      integer flag_epKEorig,flag_wb
       real pbcrit, ptcrit
       real elcrit, tlcrit
@@ -15 +15 @@
       real dtovsh, dpbase, dttrig
       real dtcrit, tau, beta, alpha, alpha1
-      real wb
+      real wbmax
       real delta
       real betad
@@ -27 +27 @@
      :                ,dtovsh, dpbase, dttrig
      :                ,dtcrit, tau, beta, alpha, alpha1
-     :                ,wb
+     :                ,flag_wb,wbmax
      :                ,delta, betad
 !$OMP THREADPRIVATE(/cv3param/)

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