Changeset 1409 for LMDZ4/trunk/libf/phylmd/radio_decay.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/radio_decay.F90 (modified) (3 diffs)

LMDZ4/trunk/libf/phylmd/radio_decay.F90

@@ -8 +8 @@
   USE dimphy
   USE infotrac, ONLY : nbtr
+  USE traclmdz_mod, ONLY : id_rn, id_pb
   IMPLICIT NONE
 !-----------------------------------------------------------------------
 ! Auteur(s): AA + CG (LGGE/CNRS) Date 24-06-94
 ! Objet: Calcul de la tendance radioactive des traceurs type radioelements
-!CG240694 : Pour un traceur, le radon
-!CG161294 : Plus un 2eme traceur, le 210Pb. Le radon decroit en plomb.
+!        Cas particulier pour le couple radon-plomb : Le radon decroit en plomb
 !-----------------------------------------------------------------------
 !
@@ -34 +34 @@
 
   DO it = 1,nbtr
+     d_tr(:,:,it) = 0.
      IF ( radio(it) ) THEN
         IF (tautr(it) .GT. 0.) THEN
@@ -41 +42 @@
               END DO
            END DO
-        ELSE 
-           d_tr(:,:,it) = 0.
         END IF
-     ELSE
-        d_tr(:,:,it) = 0.
      END IF
   END DO
+
 !-------------------------------------------------------
-!CG161294 : Cas particulier radon [it=1] => plomb [it=2]
+! Cas particulier radon (id_rn) => plomb (id_pb)
 !-------------------------------------------------------
   IF ( rnpb ) THEN
      DO k = 1,klev
         DO i = 1,klon
-           d_tr(i,k,2) = d_tr(i,k,2) - d_tr(i,k,1)
+           d_tr(i,k,id_pb) = d_tr(i,k,id_pb) - d_tr(i,k,id_rn)
         ENDDO
      ENDDO