trunk/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
  ! ======================================================================
  ! 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
  ! ======================================================================
  ! Arguments:
  ! nbtr------input-I- nombre de vrais traceurs (sans l'eau)
  ! ftsol-------input-R- Temperature du sol (Kelvin)
  ! pctsrf-----input-R-  Nature de sol (pourcentage de sol)
  ! masktr---output-R- Masque reservoir de sol traceur (1 = reservoir)
  ! fshtr----output-R- Flux surfacique de production dans le sol
  ! hsoltr---output-R- Epaisseur du reservoir de sol
  ! tautr----output-R- Constante de decroissance du traceur
  ! vdeptr---output-R- Vitesse de depot sec dans la couche Brownienne
  ! scavtr---output-R- Coefficient de lessivage
  ! ======================================================================
  ! ======================================================================

  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
  ! ! le cas echeant
  REAL masktr(klon, nbtr) ! Masque de l'echange avec la surface
  ! (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

  PRINT *, 'Call sequence information: initrrnpb'
  PRINT *, 'nbtr= ', nbtr
  PRINT *, 'nbsrf= ', nbsrf
  PRINT *, 'klon= ', klon

  ! Puis les initialisation specifiques a chaque traceur (pour le moment,
  ! Rn222)


  ! Radon it = 1

  IF (nbtr<=0) STOP 'initrrnpb pas glop pas glop'
  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
  ! AA
  ! tautr(it) = 4.765E55  ! Decroissance du radon,infinie
  ! AA
  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))==1) masktr(i, it) = 1.
    fshtr(i, it) = s*masktr(i, it)

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

  END DO

  ! 210Pb it = 2

  IF (nbtr<=1) STOP 'initrrnpb pas glop pas glop'
  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
  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)

  WRITE (*, *) 'initialisation rnpb ok'

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