source: LMDZ4/branches/LMDZ4-dev-20091210/libf/phylmd/initrrnpb.F90 @ 2370

Last change on this file since 2370 was 1191, checked in by jghattas, 15 years ago

Reecriture de phytrac et les routines concernes (Anthony Jamelot)

  • les suffix change de F -> F90 (nflxtr.F90,cltracrn.F90,initrrnpb.F90,cvltr.F90,minmaxqfi.F90,cltrac.F90,phytrac.F90)

Traitement d'un nouveau traceur berelium (optionel, toujours pour des
tests)(Anthony Jamelot)

  • radiornpb.F change du nom pour radio_decay.F90 car il traite maintenant tout les traceurs radioactives
  • ajoute init_be.F90

Nouveau interface dans phytrac pour serparer les calculs et appels
specifique a INCA avec les traitements des traceurs specifiques au LMDZ
(JG)

  • ajoute tracinca_mod.F90 pour les appeles a INCA
  • ajoute traclmdz_mod.F90 pour les calculs des traceurs specifiques a LMDZ
  • enleve fichier restartrac et ajoute la variable trs dans restartphy.nc

La convergence numerique a etait rompue uniquement pour les traceurs
LMDZ RN et PB.

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 3.6 KB
Line 
1!
2! $Id $
3!
4SUBROUTINE  initrrnpb(ftsol,pctsrf,masktr,fshtr,hsoltr,tautr,vdeptr,scavtr)
5  USE dimphy
6  USE infotrac, ONLY : nbtr
7  IMPLICIT NONE
8!======================================================================
9! Auteur(s): AA + CG (LGGE/CNRS) Date 24-06-94
10! Objet: initialisation des constantes des traceurs
11!AA Revison pour le controle avec la temperature du sol
12!AA
13!AA   it = 1 radon ss controle de ts
14!AA   it = 2 plomb ss controle de ts 
15!======================================================================
16! Arguments:
17! nbtr.............. nombre de vrais traceurs (sans l'eau)
18! ftsol....input-R-  Temperature du sol (Kelvin)
19! pctsrf...input-R-  Nature de sol (pourcentage de sol)
20! masktr...output-R- Masque reservoir de sol traceur (1 = reservoir)
21! fshtr....output-R- Flux surfacique de production dans le reservoir de sol
22! hsoltr...output-R- Epaisseur equivalente du reservoir de sol
23! tautr....output-R- Constante de decroissance radioactive du traceur
24! vdeptr...output-R- Vitesse de depot sec dans la couche Brownienne
25! scavtr...output-R- Coefficient de lessivage
26!======================================================================
27  INCLUDE "indicesol.h"
28!======================================================================
29
30  REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: pctsrf
31  REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: ftsol
32  REAL,DIMENSION(klon,nbtr),INTENT(OUT) :: masktr
33  REAL,DIMENSION(klon,nbtr),INTENT(OUT) :: fshtr
34  REAL,DIMENSION(nbtr),INTENT(OUT)      :: hsoltr
35  REAL,DIMENSION(nbtr),INTENT(OUT)      :: tautr
36  REAL,DIMENSION(nbtr),INTENT(OUT)      :: vdeptr
37  REAL,DIMENSION(nbtr),INTENT(OUT)      :: scavtr
38  INTEGER                               :: i, it
39  REAL                                  :: s
40
41  WRITE(*,*)'PASSAGE initrrnpb ...'
42!
43! Radon it = 1
44!----------------
45  IF ( nbtr .LE. 0 ) STOP '**PHYTRAC:initrrnpb:** nbtr < 0; verifier RN dans traceur.def'
46  it = 1
47  s = 1.E4             ! Source: atome par m2
48  hsoltr(it) = 0.1     ! Hauteur equivalente du reservoir :
49                       ! 1 m * porosite 0.1
50  tautr(it) = 4.765E5  ! Decroissance du radon, secondes
51  vdeptr(it) = 0.      ! Pas de depot sec pour le radon
52  scavtr(it) = 0.      ! Pas de lessivage pour le radon
53 
54  WRITE(*,*)'-------------- SOURCE DU RADON ------------------------ '
55  WRITE(*,*)'it = ',it
56  WRITE(*,*)'Source : ', s
57  WRITE(*,*)'Hauteur equivalente du reservoir de sol: ',hsoltr(it)
58  WRITE(*,*)'Decroissance (s): ', tautr(it)
59  WRITE(*,*)'Vitesse de depot sec: ',vdeptr(it)
60  WRITE(*,*)'Facteur de lessivage: ',scavtr(it)
61
62  DO i = 1,klon
63     masktr(i,it) = 0.
64     IF ( NINT(pctsrf(i,1)) .EQ. 1 ) masktr(i,it) = 1.
65     fshtr(i,it) = s * masktr(i,it)
66  END DO
67!
68! 210Pb it = 2
69!----------------
70  IF ( nbtr .LE. 1 ) STOP '**PHYTRAC**:initrrnpb:** nbtr <= 1; verifier PB dans traceur.def'
71  it = 2
72  s = 0.                ! Pas de source
73  hsoltr(it) = 10.      ! Hauteur equivalente du reservoir
74                        ! a partir duquel le depot Brownien a lieu
75  tautr(it) = 1.028E9   ! Decroissance du Pb210, secondes
76  vdeptr(it) = 1.E-3    ! 1 mm/s pour le 210Pb
77  scavtr(it) =  .5      ! Lessivage du Pb210
78  DO i = 1,klon
79     masktr(i,it) = 1.  ! Le depot sec peut avoir lieu partout
80     fshtr(i,it) = s * masktr(i,it)
81  END DO
82  WRITE(*,*)'-------------- SOURCE DU PLOMB ------------------------ '
83  WRITE(*,*)'it = ',it
84  WRITE(*,*)'Source : ', s
85  WRITE(*,*)'Hauteur equivalente du reservoir : ',hsoltr(it)
86  WRITE(*,*)'Decroissance (s): ', tautr(it)
87  WRITE(*,*)'Vitesse de depot sec: ',vdeptr(it)
88  WRITE(*,*)'Facteur de lessivage: ',scavtr(it)
89
90  WRITE(*,*) 'Initialisation RN et PB ok'
91
92END SUBROUTINE initrrnpb
Note: See TracBrowser for help on using the repository browser.