libf/phylmd/initrrnpb.f

!
! $Header: /home/cvsroot/LMDZ4/libf/phylmd/initrrnpb.F,v 1.2 2004/06/22 11:45:33 lmdzadmin Exp $
!
      SUBROUTINE  initrrnpb(ftsol,pctsrf,masktr,fshtr,hsoltr,tautr
     .                   ,vdeptr,scavtr)
      use dimens_m
      use indicesol
      use dimphy
      IMPLICIT none
c======================================================================
c Auteur(s): AA + CG (LGGE/CNRS) Date 24-06-94
c Objet: initialisation des constantes des traceurs
CAA Revison pour le controle avec la temperature du sol
cAA
CAA   it = 1 radon ss controle de ts
cAA   it = 2 plomb ss controle de ts  
c======================================================================
c Arguments:
c nbtr------input-I- nombre de vrais traceurs (sans l'eau)
c ftsol-------input-R- Temperature du sol (Kelvin)
c pctsrf-----input-R-  Nature de sol (pourcentage de sol)
c masktr---output-R- Masque reservoir de sol traceur (1 = reservoir)
c fshtr----output-R- Flux surfacique de production dans le sol
c hsoltr---output-R- Epaisseur du reservoir de sol
c tautr----output-R- Constante de decroissance du traceur
c vdeptr---output-R- Vitesse de depot sec dans la couche Brownienne
c scavtr---output-R- Coefficient de lessivage
c======================================================================
c======================================================================
C
      INTEGER i, it
      REAL pctsrf(klon,nbsrf) !Pourcentage de sol (f(nature du sol))
      REAL ftsol(klon,nbsrf)  ! Temperature du sol pour le controle Rn
c                             ! le cas echeant
      REAL masktr(klon,nbtr)  ! Masque de l'echange avec la surface
c                                 (possible => 1 )
      REAL fshtr(klon,nbtr)  ! Flux surfacique dans le reservoir de sol
      REAL hsoltr(nbtr)      ! Epaisseur equivalente du reservoir de sol
      REAL tautr(nbtr)       ! Constante de decroissance radioactive
      REAL vdeptr(nbtr)      ! Vitesse de depot sec dans la couche Brownienne
      REAL scavtr(nbtr)      ! Coefficient de lessivage
      REAL s
C
      print *, "Call sequence information: initrrnpb"
      print*,'nbtr= ',nbtr
      print*,'nbsrf= ',nbsrf
      print*,'klon= ',klon
C
C Puis les initialisation specifiques a chaque traceur (pour le moment, Rn222)
C
C
C Radon it = 1
c
      IF ( nbtr .LE. 0 ) STOP 'initrrnpb pas glop pas glop' 
      it = 1
      s = 1.E4  !  Source: atome par m2
      hsoltr(it) = 0.1      ! Hauteur equivalente du reservoir : 
c                              1 m * porosite 0.1
      tautr(it) = 4.765E5  ! Decroissance du radon, secondes
cAA
c      tautr(it) = 4.765E55  ! Decroissance du radon,infinie
cAA
      vdeptr(it) = 0. ! Pas de depot sec pour le radon
      scavtr(it) = 0. ! Pas de lessivage pour le radon

      print*, '-------------- SOURCE DU RADON ------------------------ '
      print*,'it = ',it
      print*,'Source : ', s
      print*,'Hauteur equivalente du reservoir de sol: ',hsoltr(it) 
      print*,'Decroissance (s): ', tautr(it)
      print*,'Vitesse de depot sec: ',vdeptr(it) 
      print*,'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)

cAA
cAA POur l'instant le pctsrf(i,3) = 1.0 
cAA lorsqu'il ya de la terre mias ne prend aucune autre valeur
cAA il n'est donc pas necessaire de multiplier fshtr par pctsrf
cAA 

      END DO
C
C 210Pb it = 2
C
      IF ( nbtr .LE. 1 ) STOP 'initrrnpb pas glop pas glop' 
      it = 2
      s = 0. !  Pas de source !!!
      hsoltr(it) = 10.     ! Hauteur equivalente du reservoir 
c                              a partir duquel le
c                              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
      print*, '-------------- SOURCE DU PLOMB ------------------------ '
      print*,'it = ',it
      print*,'Source : ', s
      print*,'Hauteur equivalente du reservoir : ',hsoltr(it) 
      print*,'Decroissance (s): ', tautr(it)
      print*,'Vitesse de depot sec: ',vdeptr(it) 
      print*,'Facteur de lessivage: ',scavtr(it)
c
      WRITE(*,*) 'initialisation rnpb ok'
c
      RETURN
      END
1	!
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