Changeset 1044

Timestamp:

Nov 6, 2008, 5:29:35 PM (16 years ago)

Author:

lmdzadmin

Message:

Correction oscillations verticalles de la pluie lorsque celle-ci s'evapore en totalite (correction a la
programmation originale de Emanuel)
JYG/IM

Location:

LMDZ4/trunk/libf/phylmd

Files:

• cv3_routines.F (modified) (4 diffs)
• cva_driver.F (modified) (1 diff)

LMDZ4/trunk/libf/phylmd/cv3_routines.F

@@ -55 +55 @@
       pbcrit = 150.0
       ptcrit = 500.0
-      epmax  = 0.993
+cIM lu dans physiq.def via conf_phys.F90     epmax  = 0.993
 
       omtrain = 45.0 ! used also for snow (no disctinction rain/snow)
@@ -2093 +2093 @@
        water(il,i)=revap*revap
       else
-       evap(il,i)=-evap(il,i+1)
-     :            +(wdtrain(il)+sigd(il)*wt(il,i)*water(il,i+1))
-     :                 /(sigd(il)*(ph(il,i)-ph(il,i+1))*50.)
+cJYG/IM : ci-dessous formulation originale de KE
+c      evap(il,i)=-evap(il,i+1)
+c    :            +(wdtrain(il)+sigd(il)*wt(il,i)*water(il,i+1))
+c    :                 /(sigd(il)*(ph(il,i)-ph(il,i+1))*50.)
+c
+cJYG/IM : ci-dessous modification formulation originale de KE
+c        pour eliminer oscillations verticales de pluie se produisant
+c        lorsqu'il y a evaporation totale de la pluie
+c
+       evap(il,i)= +(wdtrain(il)+sigd(il)*wt(il,i)*water(il,i+1))
+     :                 /(sigd(il)*(ph(il,i)-ph(il,i+1))*100.)
       end if
 ccc
@@ -2249 +2257 @@
      :                    ,icb,inb,delt
      :                    ,t,rr,t_wake,rr_wake,u,v,tra
-     :                    ,gz,p,ph,h,hp,lv,cpn,th
+     :                    ,gz,p,ph,h,hp,lv,cpn,th,th_wake
      :                    ,ep,clw,m,tp,mp,rp,up,vp,trap
      :                    ,wt,water,evap,b,sigd
@@ -3016 +3024 @@
       enddo
 c
-      do i=1,nl
-       do il=1,ncum
-          th_wake(il,i)=t_wake(il,i)*(1000.0/p(il,i))**rdcp
-       enddo
-      enddo
+c     do i=1,nl
+c      do il=1,ncum
+c         th_wake(il,i)=t_wake(il,i)*(1000.0/p(il,i))**rdcp
+c      enddo
+c     enddo
 c
       do i=1,nl

LMDZ4/trunk/libf/phylmd/cva_driver.F

@@ -857 +857 @@
      :                     ,icb,inb,delt
      :                     ,t,q,t_wake,q_wake,u,v,tra
-     :                     ,gz,p,ph,h,hp,lv,cpn,th
+     :                     ,gz,p,ph,h,hp,lv,cpn,th,th_wake
      :                     ,ep,clw,m,tp,mp,qp,up,vp,trap
      :                     ,wt,water,evap,b,sigd