Changeset 1742 for LMDZ5/trunk/libf/phylmd/concvl.F

Timestamp:

Apr 5, 2013, 1:49:35 PM (11 years ago)

Author:

idelkadi

Message:

1- Inclusion des developpements de la these de Romain Pilon sur le
lessivage des aerosols :

a/ par les pluies convectives (modifs cv30_routines et cv3_routines pour

sortir les champs nécessaires au calcul off-line ; modif cvltr)

b/ par les pluies stratiformes (modifs phytrac et introduction

lsc_scav).

2- Choix entre plusieurs schemas pour les pluies stratiformes, commande
par iflag_lscav.

3- Quelques corrections dans la convection "Nouvelle Physique" pour
assurer la conservation des traceurs (cv3p1_mixing et cva_driver) (travail
de Robin Locatelli).

File:

• LMDZ5/trunk/libf/phylmd/concvl.F (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/concvl.F

@@ -10 +10 @@
      .             pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr,
      .             qcondc,wd,pmflxr,pmflxs,
-     .             da,phi,mp,dd_t,dd_q,lalim_conv,wght_th)
+! RomP >>>
+!!     .             da,phi,mp,dd_t,dd_q,lalim_conv,wght_th)
+     .             da,phi,mp,phi2,d1a,dam,sij,clw,elij,        ! RomP
+     .             dd_t,dd_q,lalim_conv,wght_th,               ! RomP
+     .             evap, ep, epmlmMm,eplaMm,                   ! RomP
+     .             wdtrainA,wdtrainM)                          ! RomP
+! RomP <<<
 ***************************************************************
 *                                                             *
@@ -91 +97 @@
 
        real da(klon,klev),phi(klon,klev,klev),mp(klon,klev)
+! RomP >>>
+       real phi2(klon,klev,klev)                                    
+       real d1a(klon,klev),dam(klon,klev)
+       real sij(klon,klev,klev),clw(klon,klev),elij(klon,klev,klev)
+       REAL wdtrainA(klon,klev),wdtrainM(klon,klev)
+       REAL evap(klon,klev),ep(klon,klev)
+       REAL epmlmMm(klon,klev,klev),eplaMm(klon,klev)
+! RomP <<<
        REAL cape(klon),cin(klon),tvp(klon,klev)
        REAL Tconv(klon,klev)
@@ -248 +262 @@
          DO i = 1, klon
           cbmf(i) = 0.
-          plcl(i) = 0.
+!!          plcl(i) = 0.
           sigd(i) = 0.
          ENDDO
@@ -256 +270 @@
       plfc(:)  = 0.
       wbeff(:) = 100.
+      plcl(:) = 0.
 
       DO k = 1, klev+1
@@ -339 +354 @@
       if (iflag_con.eq.30) then
 
-      CALL cv_driver(klon,klev,klev+1,ntra,iflag_con,
+      print *, '-> cv_driver'      !jyg
+      CALL cv_driver(klon,klev,klevp1,ntra,iflag_con,
      :              t,q,qs,u,v,tra,
      $              em_p,em_ph,iflag,
      $              d_t,d_q,d_u,d_v,d_tra,rain,
-!!     $              pmflxr,cbmf,work1,work2,           !jyg
-     $              Vprecip,cbmf,work1,work2,            !jyg
+     $              Vprecip,cbmf,work1,work2,                  !jyg
      $              kbas,ktop,
      $              dtime,Ma,upwd,dnwd,dnwdbis,qcondc,wd,cape,
-     $              da,phi,mp)
+     $              da,phi,mp,phi2,d1a,dam,sij,clw,elij,       !RomP
+     $              evap,ep,epmlmMm,eplaMm,                    !RomP
+     $              wdtrainA,wdtrainM)                         !RomP
+      print *, 'cv_driver ->'      !jyg
 c
       DO i = 1,klon
@@ -369 +387 @@
      $              dd_t,dd_q,Plim1,Plim2,asupmax,supmax0,
      $              asupmaxmin,lalim_conv,
-!AC!
-     $              da,phi)
-!AC!
+!AC!+!RomP
+     $              da,phi,mp,phi2,d1a,dam,sij,clw,  ! RomP
+     $              elij,evap,ep,wdtrainA,wdtrainM)  ! RomP
+!AC!+!RomP
       endif  
 C------------------------------------------------------------------