Changeset 1742 for LMDZ5/trunk/libf/phylmd/physiq.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/physiq.F (modified) (9 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -227 +227 @@
       REAL d_qx(klon,klev,nqtot)
       REAL d_ps(klon)
+! Variables pour le transport convectif
       real da(klon,klev),phi(klon,klev,klev),mp(klon,klev)
+! Variables pour le lessivage convectif
+! RomP >>> 
+      real phi2(klon,klev,klev)
+      real d1a(klon,klev),dam(klon,klev)
+      real ev(klon,klev),ep(klon,klev)
+      real clw(klon,klev),elij(klon,klev,klev)
+      real epmlmMm(klon,klev,klev),eplaMm(klon,klev)
+      real wdtrainA(klon,klev),wdtrainM(klon,klev)
+! RomP <<<
 !IM definition dynamique o_trac dans phys_output_open
 !      type(ctrl_out) :: o_trac(nqtot)
@@ -697 +707 @@
       REAL frac_impa(klon,klev) ! fractions d'aerosols lessivees (impaction)
       REAL frac_nucl(klon,klev) ! idem (nucleation)
+! RomP >>> 
+      REAL beta_prec_fisrt(klon,klev) ! taux de conv de l'eau cond (fisrt)
+      REAL beta_prec(klon,klev)       ! taux de conv de l'eau cond (utilise)
+! RomP <<<
       INTEGER       :: iii
       REAL          :: calday
@@ -1747 +1761 @@
       mp(:,:)=0.
       phi(:,:,:)=0.
+! RomP >>>
+      phi2(:,:,:)=0.
+      beta_prec_fisrt(:,:)=0.
+      beta_prec(:,:)=0.
+      epmlmMm(:,:,:)=0.
+      eplaMm(:,:)=0.
+      d1a(:,:)=0.
+      dam(:,:)=0.
+! RomP <<<
+
 c
 c Ne pas affecter les valeurs entrees de u, v, h, et q
@@ -1812 +1836 @@
          ENDDO
          ENDDO
+!!! RomP >>>   td dyn traceur
+       IF (nqtot.GE.3) THEN
+          DO iq = 3, nqtot
+          DO k = 1, klev
+          DO i = 1, klon
+            d_tr_dyn(i,k,iq-2)=
+     $       (tr_seri(i,k,iq-2)-tr_ancien(i,k,iq-2))/dtime
+!         iiq=niadv(iq)
+!         print*,i,k," d_tr_dyn",d_tr_dyn(i,k,iq-2),"tra:",iq,tname(iiq)
+          ENDDO
+          ENDDO
+          ENDDO
+       ENDIF
+!!! RomP <<<
       ELSE
          DO k = 1, klev
@@ -1821 +1859 @@
          ENDDO
          ENDDO
+!!! RomP >>>   td dyn traceur
+        IF (nqtot.GE.3) THEN
+          DO iq = 3, nqtot
+          DO k = 1, klev
+          DO i = 1, klon
+            d_tr_dyn(i,k,iq-2)= 0.0
+          ENDDO
+          ENDDO
+          ENDDO
+       ENDIF
+!!! RomP <<<
          ancien_ok = .TRUE.
       ENDIF
@@ -2300 +2349 @@
      .        Ma,mip,Vprecip,cape,cin,tvp,Tconv,iflagctrl,
      .        pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr,qcondc,wd,
-     .        pmflxr,pmflxs,da,phi,mp,
-     .        ftd,fqd,lalim_conv,wght_th)
+! RomP >>>
+!!     .        pmflxr,pmflxs,da,phi,mp,
+!!     .        ftd,fqd,lalim_conv,wght_th)
+     .        pmflxr,pmflxs,da,phi,mp,phi2,d1a,dam,sij,clw,elij,
+     .        ftd,fqd,lalim_conv,wght_th,
+     .        ev, ep,epmlmMm,eplaMm,
+     .        wdtrainA,wdtrainM)
+! RomP <<<
 
 cIM begin
@@ -2783 +2838 @@
      .           rain_lsc, snow_lsc,
      .           pfrac_impa, pfrac_nucl, pfrac_1nucl,
-     .           frac_impa, frac_nucl,
+     .           frac_impa, frac_nucl, beta_prec_fisrt,
      .           prfl, psfl, rhcl, 
      .           zqasc, fraca,ztv,zpspsk,ztla,zthl,iflag_cldcon )
@@ -3638 +3693 @@
      I     itap,     days_elapsed+1,    jH_cur,   debut,
      I     lafin,    dtime,     u, v,     t,
-     I     paprs,    pplay,     pmfu,     pmfd, 
+     I     paprs,    pplay,     pmfu,     pmfd,
      I     pen_u,    pde_u,     pen_d,    pde_d,
      I     cdragh,   coefh,     fm_therm, entr_therm,
      I     u1,       v1,        ftsol,    pctsrf,
      I     ustar,     u10m,      v10m,
-     I     rlat,     frac_impa, frac_nucl,rlon,
+     I     rlat,     rlon,
+     I     frac_impa,frac_nucl, beta_prec_fisrt,beta_prec,
      I     presnivs, pphis,     pphi,     albsol1,
-     I     sh_in,    rhcl,      cldfra,   rneb, 
+     I     sh_in,    rhcl,      cldfra,   rneb,
      I     diafra,   cldliq,    itop_con, ibas_con,
      I     pmflxr,   pmflxs,    prfl,     psfl,
-     I     da,       phi,       mp,       upwd,     
+     I     da,       phi,       mp,       upwd,
+     I     phi2,     d1a,       dam,      sij,         !<<RomP
+     I     wdtrainA, wdtrainM,  sigd,     clw,elij,    !<<RomP
+     I     ev,       ep,        epmlmMm,  eplaMm,      !<<RomP
      I     dnwd,     aerosol_couple,      flxmass_w,
      I     tau_aero, piz_aero,  cg_aero,  ccm,
      I     rfname,
+     I     d_tr_dyn,                                   !<<RomP
      O     tr_seri)
 
@@ -3823 +3883 @@
       ENDDO
       ENDDO
-c
+
+!!! RomP >>>
+      IF (nqtot.GE.3) THEN
+        DO iq = 3, nqtot
+        DO k = 1, klev
+        DO i = 1, klon
+           tr_ancien(i,k,iq-2) = tr_seri(i,k,iq-2)
+        ENDDO
+        ENDDO
+        ENDDO
+      ENDIF
+!!! RomP <<<
 !==========================================================================
 ! Sorties des tendances pour un point particulier