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

LMDZ4/trunk/libf/phylmd/initrrnpb.F90

@@ -5 +5 @@
   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: initialisation des constantes des traceurs
-!AA Revison pour le controle avec la temperature du sol
-!AA
-!AA   it = 1 radon ss controle de ts
-!AA   it = 2 plomb ss controle de ts  
+! id_rn : identificateur du traceur radon
+! id_pb : identificateur du traceur plomb
 !======================================================================
 ! Arguments:
@@ -42 +41 @@
   CHARACTER (LEN=80) :: abort_message
 
+!
+! Radon it = id_rn
+!----------------
+  IF (id_rn /= 0) THEN
+     it = id_rn
+     s = 1.E4             ! Source: atome par m2
+     hsoltr(it) = 0.1     ! Hauteur equivalente du reservoir : 
+                          ! 1 m * porosite 0.1
+     tautr(it) = 4.765E5  ! Decroissance du radon, secondes
+     vdeptr(it) = 0.      ! Pas de depot sec pour le radon
+     scavtr(it) = 0.      ! Pas de lessivage pour le radon
+     
+     WRITE(*,*)'-------------- SOURCE DU RADON ------------------------ '
+     WRITE(*,*)'it = ',it
+     WRITE(*,*)'Source : ', s
+     WRITE(*,*)'Hauteur equivalente du reservoir de sol: ',hsoltr(it) 
+     WRITE(*,*)'Decroissance (s): ', tautr(it)
+     WRITE(*,*)'Vitesse de depot sec: ',vdeptr(it) 
+     WRITE(*,*)'Facteur de lessivage: ',scavtr(it)
 
-  WRITE(*,*)'PASSAGE initrrnpb ...'
+     DO i = 1,klon
+        masktr(i,it) = 0.
+        IF ( NINT(pctsrf(i,1)) .EQ. 1 ) masktr(i,it) = 1.
+        fshtr(i,it) = s * masktr(i,it)
+     END DO
+
+  END IF ! id_rn /= 0
+
 !
-! Radon it = 1
+! 210Pb it = id_pb
 !----------------
-  IF ( nbtr .LE. 0 ) then
-    abort_message = '**PHYTRAC:initrrnpb:** nbtr < 0; verifier RN dans traceur.def' 
-    CALL abort_gcm (modname,abort_message,1)
-  ENDIF
-  it = 1
-  s = 1.E4             ! Source: atome par m2
-  hsoltr(it) = 0.1     ! Hauteur equivalente du reservoir : 
-                       ! 1 m * porosite 0.1
-  tautr(it) = 4.765E5  ! Decroissance du radon, secondes
-  vdeptr(it) = 0.      ! Pas de depot sec pour le radon
-  scavtr(it) = 0.      ! Pas de lessivage pour le radon
-  
-  WRITE(*,*)'-------------- SOURCE DU RADON ------------------------ '
-  WRITE(*,*)'it = ',it
-  WRITE(*,*)'Source : ', s
-  WRITE(*,*)'Hauteur equivalente du reservoir de sol: ',hsoltr(it) 
-  WRITE(*,*)'Decroissance (s): ', tautr(it)
-  WRITE(*,*)'Vitesse de depot sec: ',vdeptr(it) 
-  WRITE(*,*)'Facteur de lessivage: ',scavtr(it)
-
-  DO i = 1,klon
-     masktr(i,it) = 0.
-     IF ( NINT(pctsrf(i,1)) .EQ. 1 ) masktr(i,it) = 1.
-     fshtr(i,it) = s * masktr(i,it)
-  END DO
-!
-! 210Pb it = 2
-!----------------
-  IF ( nbtr .LE. 1 ) THEN
-    abort_message='**PHYTRAC**:initrrnpb:** nbtr <= 1; verifier PB dans traceur.def' 
-    CALL abort_gcm (modname,abort_message,1)
-  ENDIF
-  it = 2
-  s = 0.                ! Pas de source 
-  hsoltr(it) = 10.      ! Hauteur equivalente du reservoir 
-                        ! a partir duquel le depot Brownien a lieu
-  tautr(it) = 1.028E9   ! Decroissance du Pb210, secondes
-  vdeptr(it) = 1.E-3    ! 1 mm/s pour le 210Pb
-  scavtr(it) =  .5      ! Lessivage du Pb210
-  DO i = 1,klon
-     masktr(i,it) = 1.  ! Le depot sec peut avoir lieu partout
-     fshtr(i,it) = s * masktr(i,it)
-  END DO
-  WRITE(*,*)'-------------- SOURCE DU PLOMB ------------------------ '
-  WRITE(*,*)'it = ',it
-  WRITE(*,*)'Source : ', s
-  WRITE(*,*)'Hauteur equivalente du reservoir : ',hsoltr(it) 
-  WRITE(*,*)'Decroissance (s): ', tautr(it)
-  WRITE(*,*)'Vitesse de depot sec: ',vdeptr(it) 
-  WRITE(*,*)'Facteur de lessivage: ',scavtr(it)
-
-  WRITE(*,*) 'Initialisation RN et PB ok'
-
+  IF (id_pb /= 0) THEN
+     it = id_pb
+     s = 0.                ! Pas de source 
+     hsoltr(it) = 10.      ! Hauteur equivalente du reservoir 
+                           ! a partir duquel le depot Brownien a lieu
+     tautr(it) = 1.028E9   ! Decroissance du Pb210, secondes
+     vdeptr(it) = 1.E-3    ! 1 mm/s pour le 210Pb
+     scavtr(it) =  .5      ! Lessivage du Pb210
+     DO i = 1,klon
+        masktr(i,it) = 1.  ! Le depot sec peut avoir lieu partout
+        fshtr(i,it) = s * masktr(i,it)
+     END DO
+     WRITE(*,*)'-------------- SOURCE DU PLOMB ------------------------ '
+     WRITE(*,*)'it = ',it
+     WRITE(*,*)'Source : ', s
+     WRITE(*,*)'Hauteur equivalente du reservoir : ',hsoltr(it) 
+     WRITE(*,*)'Decroissance (s): ', tautr(it)
+     WRITE(*,*)'Vitesse de depot sec: ',vdeptr(it) 
+     WRITE(*,*)'Facteur de lessivage: ',scavtr(it)
+     
+  END IF
+     
 END SUBROUTINE initrrnpb