source: LMDZ4/trunk/libf/phylmd/initrrnpb.F @ 1167

Last change on this file since 1167 was 1146, checked in by Laurent Fairhead, 16 years ago

Réintegration dans le tronc des modifications issues de la branche LMDZ-dev
comprises entre la révision 1074 et 1145
Validation: une simulation de 1 jour en séquentiel sur PC donne les mêmes
résultats entre la trunk et la dev
LF

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