Changeset 104 for trunk/libf/phytitan/phytrac.F

Timestamp:

Mar 25, 2011, 4:58:38 PM (14 years ago)

Author:

slebonnois

Message:

SLebonnois: modification de makelmdz et create_make_gcm pour pouvoir
compiler la chimie titan. Pas de raison que ca gene les autres.
Dans cette version, les compilations de Venus et Titan fonctionnent.

Phytitan: modifications pour pouvoir compiler correctement.
Il ne manque plus que physiq.F a faire.

File:

• trunk/libf/phytitan/phytrac.F (modified) (4 diffs)

trunk/libf/phytitan/phytrac.F

@@ -2 +2 @@
      .                   nqmax,nmicro,ptimestep,appkim,dtkim,
      .                   pplev,pplay,delp,ptemp,pmu0,pfract,pdecli,
-     .                   lonsol,tr_seri,d_tr_mph,d_tr_kim)
+     .                   lonsol,tr_seri,qaer,d_tr_mph,d_tr_kim)
 
 c======================================================================
@@ -45 +45 @@
       REAL pmu0(klon), pfract(klon), pdecli, lonsol
       REAL tr_seri(klon,klev,nqmax)
+      REAL qaer(klon,klev,nqmax)
       REAL d_tr_mph(klon,klev,nqmax),d_tr_kim(klon,klev,nqmax)
 
 c======================================================================
 c Local variables
-
-* common relatifs aux aerosols 
-      REAL qaer(klon,klev,nqtot)
-      common/traceurs/qaer
 
 c grandeurs en moyennes zonales
@@ -62 +59 @@
       REAL pdqmfi(jjm+1,klev,nqmax)
       REAL ychim(jjm+1,klev,nqmax-nmicro) 
-      REAL qysat(klev,nqtot)        ! dim nqtot, mais en fait nqmax-nmicro (save...)
-      REAL pdyfi(jjm+1,klev,nqtot)  ! dim nqtot, mais en fait nqmax-nmicro (save...)
-      character*10 nomqy(nqmax-nmicro+1)
-      integer      i,j,l,iq,ig0
-
 c La saturation n est calculee qu une seule fois: sauvegarde qysat
 c La chimie n est pas calculee tous les pas, il faut donc
 c                      sauvegarder les sorties de la chimie
-
-      SAVE pdyfi,qysat  
-      
-c======================================================================
-c======================================================================
+      REAL,save,allocatable :: qysat(:,:),pdyfi(:,:,:) 
+      
+      character*10 nomqy(nqmax-nmicro+1)
+      integer      i,j,l,iq,ig0
+      
+c======================================================================
+c======================================================================
+
+      if (firstcall) then
+       allocate(qysat(klev,nqmax-nmicro),pdyfi(jjm+1,klev,nqmax-nmicro))
+      endif
 
 c-----------------------------------------------------------------------
@@ -349 +347 @@
          DO l=1,llm
           DO i = 1, klon
-c incrementation de la tendance sur qaer (pour sorties dans physiq.F)
-            qaer(i,l,iq) = qaer(i,l,iq) + d_tr_mph(i,l,iq)*ptimestep
 c  Traceurs microphysiques: passage en intensif: n/m^2 --> n/kg
             d_tr_mph(i,l,iq) = d_tr_mph(i,l,iq)*RG/delp(i,l)