Changeset 1409 for LMDZ4/trunk/libf/phylmd/cltracrn.F90

Timestamp:

Jul 8, 2010, 1:46:29 PM (14 years ago)

Author:

jghattas

Message:

• Added tracer "Age of stratospheric air", activated by putting Aga in tracer.def.
• The logical rnpb is now controled during execution. rnpb is true only if both tracers RN and PB existe in tracer.def. RN and PB can now be removed from tracer.def
• In tracer.def, the 2 water traceurs (H2Ov and H2Ol) must still be the first 2 tracers. The following tracers have no specific order(RN and PB can now change places). Still a minimum of 3 tracers are required.

File:

• LMDZ4/trunk/libf/phylmd/cltracrn.F90 (modified) (6 diffs)

LMDZ4/trunk/libf/phylmd/cltracrn.F90

@@ -8 +8 @@
   
   USE dimphy
+  USE traclmdz_mod, ONLY : id_rn, id_pb
   IMPLICIT NONE
 !======================================================================
@@ -20 +21 @@
 !---------------------------------------------------------------------
 ! Arguments:
-! itr......input-R-  le type de traceur 1- Rn 2 - Pb
+! itr......input-R-  le type de traceur : id_rn(radon), id_pb(plomb)
 ! dtime....input-R-  intervalle du temps (en secondes) ~ pdtphys
 ! u1lay....input-R-  vent u de la premiere couche (m/s)
@@ -195 +196 @@
 !--------------------------------------------------------
 
-     IF ( (itr.eq.1.AND.NINT(masktr(i)).EQ.1.AND.lat(i).GE.60..AND.lat(i).LE.70.).OR.      &
-          (itr.eq.2.AND.NINT(masktr(i)).EQ.1.AND.lat(i).GE.60..AND.lat(i).LE.70.) ) THEN
+     IF ( (itr.eq.id_rn.AND.NINT(masktr(i)).EQ.1.AND.lat(i).GE.60..AND.lat(i).LE.70.).OR.      &
+          (itr.eq.id_pb.AND.NINT(masktr(i)).EQ.1.AND.lat(i).GE.60..AND.lat(i).LE.70.) ) THEN
         zx_trs(i) = local_trs(i)
         zx_a = zx_trs(i)                                           &
@@ -221 +222 @@
 !---------------------------------------------------------------
 
-     IF ( (itr.EQ.1.AND.NINT(masktr(i)).EQ.0).OR.       &
-          (itr.EQ.1.AND.NINT(masktr(i)).EQ.1.AND.lat(i).LT.-60.)) THEN
+     IF ( (itr.EQ.id_rn.AND.NINT(masktr(i)).EQ.0).OR.       &
+          (itr.EQ.id_rn.AND.NINT(masktr(i)).EQ.1.AND.lat(i).LT.-60.)) THEN
         zx_trs(i) = 0.
         local_trs(i) = 0.
@@ -230 +231 @@
 ! des oceans et des continents
 !--------------------------------------------------------------
-     IF ( (itr.EQ.1.AND.NINT(masktr(i)).EQ.0).OR.    &
-          (itr.EQ.1.AND.NINT(masktr(i)).EQ.1.AND.lat(i).GT.70.)) THEN
+     IF ( (itr.EQ.id_rn.AND.NINT(masktr(i)).EQ.0).OR.    &
+          (itr.EQ.id_rn.AND.NINT(masktr(i)).EQ.1.AND.lat(i).GT.70.)) THEN
         zx_trs(i) = 0.
         local_trs(i) = 0.
@@ -239 +240 @@
 !--------------------------------------------
 
-     IF (itr.eq.1.AND.NINT(masktr(i)).EQ.0) THEN
+     IF (itr.eq.id_rn.AND.NINT(masktr(i)).EQ.0) THEN
         zx_trs(i) = 0.
         local_trs(i) = 0.