!
! $Id: initrrnpb.F90 1299 2010-01-20 14:27:21Z jyg $
!
SUBROUTINE  initrrnpb(ftsol,pctsrf,masktr,fshtr,hsoltr,tautr,vdeptr,scavtr)
  USE dimphy
  USE infotrac, ONLY : nbtr
  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.............. 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 reservoir de sol
! hsoltr...output-R- Epaisseur equivalente du reservoir de sol
! tautr....output-R- Constante de decroissance radioactive du traceur
! vdeptr...output-R- Vitesse de depot sec dans la couche Brownienne
! scavtr...output-R- Coefficient de lessivage
!======================================================================
  INCLUDE "indicesol.h"
!======================================================================

  REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: pctsrf
  REAL,DIMENSION(klon,nbsrf),INTENT(IN) :: ftsol
  REAL,DIMENSION(klon,nbtr),INTENT(OUT) :: masktr
  REAL,DIMENSION(klon,nbtr),INTENT(OUT) :: fshtr
  REAL,DIMENSION(nbtr),INTENT(OUT)      :: hsoltr
  REAL,DIMENSION(nbtr),INTENT(OUT)      :: tautr
  REAL,DIMENSION(nbtr),INTENT(OUT)      :: vdeptr
  REAL,DIMENSION(nbtr),INTENT(OUT)      :: scavtr
  INTEGER                               :: i, it
  REAL                                  :: s

  CHARACTER (LEN=20) :: modname='initrrnpb'
  CHARACTER (LEN=80) :: abort_message


  WRITE(*,*)'PASSAGE initrrnpb ...'
!
! Radon it = 1
!----------------
  IF ( nbtr .LE. 0 ) then
    abort_message = '**PHYTRAC:initrrnpb:** nbtr < 0; verifier RN dans traceur.def' 
    CALL abort_gcm (modname,abort_message,1)
  ENDIF
  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
  vdeptr(it) = 0.      ! Pas de depot sec pour le radon
  scavtr(it) = 0.      ! Pas de lessivage pour le radon
  
  WRITE(*,*)'-------------- SOURCE DU RADON ------------------------ '
  WRITE(*,*)'it = ',it
  WRITE(*,*)'Source : ', s
  WRITE(*,*)'Hauteur equivalente du reservoir de sol: ',hsoltr(it) 
  WRITE(*,*)'Decroissance (s): ', tautr(it)
  WRITE(*,*)'Vitesse de depot sec: ',vdeptr(it) 
  WRITE(*,*)'Facteur de lessivage: ',scavtr(it)

  DO i = 1,klon
     masktr(i,it) = 0.
     IF ( NINT(pctsrf(i,1)) .EQ. 1 ) masktr(i,it) = 1.
     fshtr(i,it) = s * masktr(i,it)
  END DO
!
! 210Pb it = 2
!----------------
  IF ( nbtr .LE. 1 ) THEN
    abort_message='**PHYTRAC**:initrrnpb:** nbtr <= 1; verifier PB dans traceur.def' 
    CALL abort_gcm (modname,abort_message,1)
  ENDIF
  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
  WRITE(*,*)'-------------- SOURCE DU PLOMB ------------------------ '
  WRITE(*,*)'it = ',it
  WRITE(*,*)'Source : ', s
  WRITE(*,*)'Hauteur equivalente du reservoir : ',hsoltr(it) 
  WRITE(*,*)'Decroissance (s): ', tautr(it)
  WRITE(*,*)'Vitesse de depot sec: ',vdeptr(it) 
  WRITE(*,*)'Facteur de lessivage: ',scavtr(it)

  WRITE(*,*) 'Initialisation RN et PB ok'

END SUBROUTINE initrrnpb