source: LMDZ.3.3/trunk/libf/phylmd/initrrnpb.F @ 5392

Last change on this file since 5392 was 2, checked in by lmdz, 25 years ago

Initial revision

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