source: LMDZ5/trunk/libf/phylmd/initrrnpb.F90 @ 1705

Last change on this file since 1705 was 1409, checked in by jghattas, 14 years ago
  • 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.
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 3.6 KB
Line 
1!
2! $Id: initrrnpb.F90 1409 2010-07-08 11:46:29Z emillour $
3!
4SUBROUTINE  initrrnpb(ftsol,pctsrf,masktr,fshtr,hsoltr,tautr,vdeptr,scavtr)
5  USE dimphy
6  USE infotrac, ONLY : nbtr
7  USE traclmdz_mod, ONLY : id_rn, id_pb
8  IMPLICIT NONE
9!======================================================================
10! Auteur(s): AA + CG (LGGE/CNRS) Date 24-06-94
11! Objet: initialisation des constantes des traceurs
12! id_rn : identificateur du traceur radon
13! id_pb : identificateur du traceur plomb
14!======================================================================
15! Arguments:
16! nbtr.............. nombre de vrais traceurs (sans l'eau)
17! ftsol....input-R-  Temperature du sol (Kelvin)
18! pctsrf...input-R-  Nature de sol (pourcentage de sol)
19! masktr...output-R- Masque reservoir de sol traceur (1 = reservoir)
20! fshtr....output-R- Flux surfacique de production dans le reservoir de sol
21! hsoltr...output-R- Epaisseur equivalente du reservoir de sol
22! tautr....output-R- Constante de decroissance radioactive du traceur
23! vdeptr...output-R- Vitesse de depot sec dans la couche Brownienne
24! scavtr...output-R- Coefficient de lessivage
25!======================================================================
26  INCLUDE "indicesol.h"
27!======================================================================
28
29  REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: pctsrf
30  REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: ftsol
31  REAL,DIMENSION(klon,nbtr),INTENT(OUT) :: masktr
32  REAL,DIMENSION(klon,nbtr),INTENT(OUT) :: fshtr
33  REAL,DIMENSION(nbtr),INTENT(OUT)      :: hsoltr
34  REAL,DIMENSION(nbtr),INTENT(OUT)      :: tautr
35  REAL,DIMENSION(nbtr),INTENT(OUT)      :: vdeptr
36  REAL,DIMENSION(nbtr),INTENT(OUT)      :: scavtr
37  INTEGER                               :: i, it
38  REAL                                  :: s
39
40  CHARACTER (LEN=20) :: modname='initrrnpb'
41  CHARACTER (LEN=80) :: abort_message
42
43!
44! Radon it = id_rn
45!----------------
46  IF (id_rn /= 0) THEN
47     it = id_rn
48     s = 1.E4             ! Source: atome par m2
49     hsoltr(it) = 0.1     ! Hauteur equivalente du reservoir :
50                          ! 1 m * porosite 0.1
51     tautr(it) = 4.765E5  ! Decroissance du radon, secondes
52     vdeptr(it) = 0.      ! Pas de depot sec pour le radon
53     scavtr(it) = 0.      ! Pas de lessivage pour le radon
54     
55     WRITE(*,*)'-------------- SOURCE DU RADON ------------------------ '
56     WRITE(*,*)'it = ',it
57     WRITE(*,*)'Source : ', s
58     WRITE(*,*)'Hauteur equivalente du reservoir de sol: ',hsoltr(it)
59     WRITE(*,*)'Decroissance (s): ', tautr(it)
60     WRITE(*,*)'Vitesse de depot sec: ',vdeptr(it)
61     WRITE(*,*)'Facteur de lessivage: ',scavtr(it)
62
63     DO i = 1,klon
64        masktr(i,it) = 0.
65        IF ( NINT(pctsrf(i,1)) .EQ. 1 ) masktr(i,it) = 1.
66        fshtr(i,it) = s * masktr(i,it)
67     END DO
68
69  END IF ! id_rn /= 0
70
71!
72! 210Pb it = id_pb
73!----------------
74  IF (id_pb /= 0) THEN
75     it = id_pb
76     s = 0.                ! Pas de source
77     hsoltr(it) = 10.      ! Hauteur equivalente du reservoir
78                           ! a partir duquel le depot Brownien a lieu
79     tautr(it) = 1.028E9   ! Decroissance du Pb210, secondes
80     vdeptr(it) = 1.E-3    ! 1 mm/s pour le 210Pb
81     scavtr(it) =  .5      ! Lessivage du Pb210
82     DO i = 1,klon
83        masktr(i,it) = 1.  ! Le depot sec peut avoir lieu partout
84        fshtr(i,it) = s * masktr(i,it)
85     END DO
86     WRITE(*,*)'-------------- SOURCE DU PLOMB ------------------------ '
87     WRITE(*,*)'it = ',it
88     WRITE(*,*)'Source : ', s
89     WRITE(*,*)'Hauteur equivalente du reservoir : ',hsoltr(it)
90     WRITE(*,*)'Decroissance (s): ', tautr(it)
91     WRITE(*,*)'Vitesse de depot sec: ',vdeptr(it)
92     WRITE(*,*)'Facteur de lessivage: ',scavtr(it)
93     
94  END IF
95     
96END SUBROUTINE initrrnpb
Note: See TracBrowser for help on using the repository browser.