1 |
guez |
3 |
! |
2 |
|
|
! $Header: /home/cvsroot/LMDZ4/libf/phylmd/initrrnpb.F,v 1.2 2004/06/22 11:45:33 lmdzadmin Exp $ |
3 |
|
|
! |
4 |
|
|
SUBROUTINE initrrnpb(ftsol,pctsrf,masktr,fshtr,hsoltr,tautr |
5 |
|
|
. ,vdeptr,scavtr) |
6 |
|
|
use dimens_m |
7 |
|
|
use indicesol |
8 |
|
|
use dimphy |
9 |
|
|
IMPLICIT none |
10 |
|
|
c====================================================================== |
11 |
|
|
c Auteur(s): AA + CG (LGGE/CNRS) Date 24-06-94 |
12 |
|
|
c Objet: initialisation des constantes des traceurs |
13 |
|
|
CAA Revison pour le controle avec la temperature du sol |
14 |
|
|
cAA |
15 |
|
|
CAA it = 1 radon ss controle de ts |
16 |
|
|
cAA it = 2 plomb ss controle de ts |
17 |
|
|
c====================================================================== |
18 |
|
|
c Arguments: |
19 |
|
|
c nbtr------input-I- nombre de vrais traceurs (sans l'eau) |
20 |
|
|
c ftsol-------input-R- Temperature du sol (Kelvin) |
21 |
|
|
c pctsrf-----input-R- Nature de sol (pourcentage de sol) |
22 |
|
|
c masktr---output-R- Masque reservoir de sol traceur (1 = reservoir) |
23 |
|
|
c fshtr----output-R- Flux surfacique de production dans le sol |
24 |
|
|
c hsoltr---output-R- Epaisseur du reservoir de sol |
25 |
|
|
c tautr----output-R- Constante de decroissance du traceur |
26 |
|
|
c vdeptr---output-R- Vitesse de depot sec dans la couche Brownienne |
27 |
|
|
c scavtr---output-R- Coefficient de lessivage |
28 |
|
|
c====================================================================== |
29 |
|
|
c====================================================================== |
30 |
|
|
C |
31 |
|
|
INTEGER i, it |
32 |
|
|
REAL pctsrf(klon,nbsrf) !Pourcentage de sol (f(nature du sol)) |
33 |
|
|
REAL ftsol(klon,nbsrf) ! Temperature du sol pour le controle Rn |
34 |
|
|
c ! le cas echeant |
35 |
|
|
REAL masktr(klon,nbtr) ! Masque de l'echange avec la surface |
36 |
|
|
c (possible => 1 ) |
37 |
|
|
REAL fshtr(klon,nbtr) ! Flux surfacique dans le reservoir de sol |
38 |
|
|
REAL hsoltr(nbtr) ! Epaisseur equivalente du reservoir de sol |
39 |
|
|
REAL tautr(nbtr) ! Constante de decroissance radioactive |
40 |
|
|
REAL vdeptr(nbtr) ! Vitesse de depot sec dans la couche Brownienne |
41 |
|
|
REAL scavtr(nbtr) ! Coefficient de lessivage |
42 |
|
|
REAL s |
43 |
|
|
C |
44 |
|
|
print *, "Call sequence information: initrrnpb" |
45 |
|
|
print*,'nbtr= ',nbtr |
46 |
|
|
print*,'nbsrf= ',nbsrf |
47 |
|
|
print*,'klon= ',klon |
48 |
|
|
C |
49 |
|
|
C Puis les initialisation specifiques a chaque traceur (pour le moment, Rn222) |
50 |
|
|
C |
51 |
|
|
C |
52 |
|
|
C Radon it = 1 |
53 |
|
|
c |
54 |
|
|
IF ( nbtr .LE. 0 ) STOP 'initrrnpb pas glop pas glop' |
55 |
|
|
it = 1 |
56 |
|
|
s = 1.E4 ! Source: atome par m2 |
57 |
|
|
hsoltr(it) = 0.1 ! Hauteur equivalente du reservoir : |
58 |
|
|
c 1 m * porosite 0.1 |
59 |
|
|
tautr(it) = 4.765E5 ! Decroissance du radon, secondes |
60 |
|
|
cAA |
61 |
|
|
c tautr(it) = 4.765E55 ! Decroissance du radon,infinie |
62 |
|
|
cAA |
63 |
|
|
vdeptr(it) = 0. ! Pas de depot sec pour le radon |
64 |
|
|
scavtr(it) = 0. ! Pas de lessivage pour le radon |
65 |
|
|
|
66 |
|
|
print*, '-------------- SOURCE DU RADON ------------------------ ' |
67 |
|
|
print*,'it = ',it |
68 |
|
|
print*,'Source : ', s |
69 |
|
|
print*,'Hauteur equivalente du reservoir de sol: ',hsoltr(it) |
70 |
|
|
print*,'Decroissance (s): ', tautr(it) |
71 |
|
|
print*,'Vitesse de depot sec: ',vdeptr(it) |
72 |
|
|
print*,'Facteur de lessivage: ',scavtr(it) |
73 |
|
|
|
74 |
|
|
DO i = 1,klon |
75 |
|
|
masktr(i,it) = 0. |
76 |
|
|
IF ( NINT(pctsrf(i,1)) .EQ. 1 ) masktr(i,it) = 1. |
77 |
|
|
fshtr(i,it) = s * masktr(i,it) |
78 |
|
|
|
79 |
|
|
cAA |
80 |
|
|
cAA POur l'instant le pctsrf(i,3) = 1.0 |
81 |
|
|
cAA lorsqu'il ya de la terre mias ne prend aucune autre valeur |
82 |
|
|
cAA il n'est donc pas necessaire de multiplier fshtr par pctsrf |
83 |
|
|
cAA |
84 |
|
|
|
85 |
|
|
END DO |
86 |
|
|
C |
87 |
|
|
C 210Pb it = 2 |
88 |
|
|
C |
89 |
|
|
IF ( nbtr .LE. 1 ) STOP 'initrrnpb pas glop pas glop' |
90 |
|
|
it = 2 |
91 |
|
|
s = 0. ! Pas de source !!! |
92 |
|
|
hsoltr(it) = 10. ! Hauteur equivalente du reservoir |
93 |
|
|
c a partir duquel le |
94 |
|
|
c depot Brownien a lieu |
95 |
|
|
tautr(it) = 1.028E9 ! Decroissance du Pb210, secondes |
96 |
|
|
vdeptr(it) = 1.E-3 ! 1 mm/s pour le 210Pb |
97 |
|
|
scavtr(it) = .5 ! Lessivage du Pb210 |
98 |
|
|
DO i = 1,klon |
99 |
|
|
masktr(i,it) = 1. ! Le depot sec peut avoir lieu partout |
100 |
|
|
fshtr(i,it) = s * masktr(i,it) |
101 |
|
|
END DO |
102 |
|
|
print*, '-------------- SOURCE DU PLOMB ------------------------ ' |
103 |
|
|
print*,'it = ',it |
104 |
|
|
print*,'Source : ', s |
105 |
|
|
print*,'Hauteur equivalente du reservoir : ',hsoltr(it) |
106 |
|
|
print*,'Decroissance (s): ', tautr(it) |
107 |
|
|
print*,'Vitesse de depot sec: ',vdeptr(it) |
108 |
|
|
print*,'Facteur de lessivage: ',scavtr(it) |
109 |
|
|
c |
110 |
|
|
WRITE(*,*) 'initialisation rnpb ok' |
111 |
|
|
c |
112 |
|
|
RETURN |
113 |
|
|
END |