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