Changeset 1376 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phytrac.F90

Timestamp:

May 5, 2010, 3:23:18 PM (14 years ago)

Author:

musat

Message:

Output all tracers defined in .def in hist files with dynamic
declaration of LMDZ atmospheric tracers' output levels
Add 6 pseudo-water tracers with and without transport by boundary layer
IM

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phytrac.F90 (modified) (5 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phytrac.F90

@@ +1 @@
+
 !$Id $
 
@@ -213 +214 @@
      SELECT CASE(type_trac)
      CASE('lmdz')
-        CALL traclmdz_init(pctsrf, ftsol, tr_seri, aerosol, lessivage)
+!IM ajout t_seri, pplay, sh    CALL traclmdz_init(pctsrf, ftsol, tr_seri, aerosol, lessivage)
+        CALL traclmdz_init(pctsrf, ftsol, tr_seri, t_seri, pplay, sh, aerosol, lessivage)
      CASE('inca')
         source(:,:)=0.
@@ -239 +241 @@
           paprs,    pplay,        cdragh,  coefh,   &
           yu1,      yv1,          ftsol,   pctsrf,  &
-          xlat,     couchelimite,                   &
+          xlat,     couchelimite, sh,               &
           tr_seri,  source,       solsym,  d_tr_cl)
      
@@ -277 +279 @@
         END IF
 
+!IM ajout traceurs RR 
+!      print*,'phytrac it,nseuil=',it,nseuil
+       IF (it.lt.nseuil) THEN
         DO k = 1, klev
            DO i = 1, klon        
@@ -282 +287 @@
            END DO
         END DO
+       END IF !(it.lt.nseuil) then
 
         CALL minmaxqfi(tr_seri(:,:,it),0.,1.e33,'convection it = '//solsym(it))