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

Timestamp:

Apr 5, 2013, 1:49:35 PM (12 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/cv3a_uncompress.F (modified) (7 diffs)

LMDZ5/trunk/libf/phylmd/cv3a_uncompress.F

@@ -9 +9 @@
      :         ,Plim1,Plim2,asupmax,supmax0
      :         ,asupmaxmin
-!AC!
-     :         ,da,phi
-!AC!
+!
+     :         ,da,phi                               !AC!
+     :         ,mp,phi2,d1a,dam,sigij                  !RomP
+     :         ,wdtrainA,wdtrainM,elij,clw           !RomP
+     :         ,evap,ep                              !RomP
+!
      o         ,iflag1,kbas1,ktop1
      :         ,precip1,cbmf1,plcl1,plfc1,wbeff1,sig1,w01,ptop21
@@ -21 +24 @@
      :         ,Plim11,Plim21,asupmax1,supmax01
      :         ,asupmaxmin1     
-!AC!
-     :         ,da1,phi1  )
-!AC!
+!
+     o         ,da1,phi1                             !AC!
+     o         ,mp1,phi21,d1a1,dam1,sigij1             !RomP
+     o         ,wdtrainA1,wdtrainM1,elij1,clw1       !RomP
+     o         ,evap1,ep1)                           !RomP
+!
 ***************************************************************
 *                                                             *
@@ -56 +62 @@
       real asupmax(nloc,nd),supmax0(nloc)
       real asupmaxmin(nloc)
-!AC!
-      real da(nloc,nd),phi(nloc,nd,nd)
-!AC!
+!
+      real da(nloc,nd),phi(nloc,nd,nd)                    !AC!
+      real mp(nloc,nd)                                    !RomP
+      real phi2(nloc,nd,nd)                               !RomP
+      real d1a(nloc,nd),dam(nloc,nd)                      !RomP
+      real wdtrainA(nloc,nd), wdtrainM(nloc,nd)           !RomP
+      real sigij(nloc,nd,nd)                                !RomP
+      real elij(nloc,nd,nd),clw(nloc,nd)                  !RomP
+      real evap(nloc,nd),ep(nloc,nd)                      !RomP
+!
 c outputs:
       integer iflag1(len),kbas1(len),ktop1(len)
@@ -76 +89 @@
       real asupmax1(len,nd),supmax01(len)
       real asupmaxmin1(len)
-!AC!
-      real da1(nloc,nd),phi1(nloc,nd,nd)
-!AC!
+!
+      real da1(nloc,nd),phi1(nloc,nd,nd)                  !AC!
+      real mp1(nloc,nd)                                   !RomP
+      real phi21(nloc,nd,nd)                              !RomP
+      real d1a1(nloc,nd),dam1(nloc,nd)                    !RomP
+      real wdtrainA1(len,nd), wdtrainM1(len,nd)           !RomP
+      real sigij1(len,nd,nd)                                !RomP
+      real elij1(len,nd,nd),clw1(len,nd)                  !RomP
+      real evap1(len,nd),ep1(len,nd)                      !RomP
+!
 c
 c local variables:
@@ -122 +142 @@
             fqd1(idcum(i),k)=fqd(i,k)
             asupmax1(idcum(i),k)=asupmax(i,k)
-!AC!
-            da1(idcum(i),k)=da(i,k)
-!AC!
+!
+            da1(idcum(i),k)=da(i,k)                       !AC!
+            mp1(idcum(i),k)      = mp(i,k)                !RomP
+            d1a1(idcum(i),k)     = d1a(i,k)               !RomP
+            dam1(idcum(i),k)     = dam(i,k)               !RomP
+            wdtrainA1(idcum(i),k)= wdtrainA(i,k)          !RomP
+            wdtrainM1(idcum(i),k)= wdtrainM(i,k)          !RomP
+            clw1(idcum(i),k)     = clw(i,k)               !RomP
+            evap1(idcum(i),k)    = evap(i,k)              !RomP
+            ep1(idcum(i),k)      = ep(i,k)                !RomP
+!
  2010    continue
  2020   continue
@@ -146 +174 @@
          do k=1,nd 
           do i=1,ncum
-            phi1(idcum(i),k,j)=phi(i,k,j)
+            phi1(idcum(i),k,j)=phi(i,k,j)                 !AC!
+            phi21(idcum(i),k1,k2)= phi2(idcum(i),k1,k2)   !RomP
+            sigij1(idcum(i),k1,k2) = sigij(idcum(i),k1,k2)    !RomP
+            elij1(idcum(i),k1,k2)= elij(idcum(i),k1,k2)   !RomP
           end do
          end do
@@ -157 +188 @@
 c          do 2200 i=1,ncum
 c            ment1(idcum(i),k1,k2) = ment(i,k1,k2)
-c            sij1(idcum(i),k1,k2) = sij(i,k1,k2)
+c            sigij1(idcum(i),k1,k2) = sigij(i,k1,k2)
 c2200      enddo
 c2210     enddo