Changeset 1191 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/initrrnpb.F90

Timestamp:

Jun 24, 2009, 11:56:13 AM (15 years ago)

Author:

jghattas

Message:

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.

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/initrrnpb.F90 (moved) (moved from LMDZ4/branches/LMDZ4-dev/libf/phylmd/initrrnpb.F) (1 diff)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/initrrnpb.F90

@@ -1 +1 @@
 !
-! $Header$
+! $Id $
 !
-      SUBROUTINE  initrrnpb(ftsol,pctsrf,masktr,fshtr,hsoltr,tautr
-     .                   ,vdeptr,scavtr)
-      USE dimphy
-      USE infotrac, ONLY : nbtr
-      IMPLICIT none
-c======================================================================
-c Auteur(s): AA + CG (LGGE/CNRS) Date 24-06-94
-c Objet: initialisation des constantes des traceurs
-CAA Revison pour le controle avec la temperature du sol
-cAA
-CAA   it = 1 radon ss controle de ts
-cAA   it = 2 plomb ss controle de ts  
-c======================================================================
-c Arguments:
-c nbtr------input-I- nombre de vrais traceurs (sans l'eau)
-c ftsol-------input-R- Temperature du sol (Kelvin)
-c pctsrf-----input-R-  Nature de sol (pourcentage de sol)
-c masktr---output-R- Masque reservoir de sol traceur (1 = reservoir)
-c fshtr----output-R- Flux surfacique de production dans le sol
-c hsoltr---output-R- Epaisseur du reservoir de sol
-c tautr----output-R- Constante de decroissance du traceur
-c vdeptr---output-R- Vitesse de depot sec dans la couche Brownienne
-c scavtr---output-R- Coefficient de lessivage
-c======================================================================
-cym#include "dimensions.h"
-cym#include "dimphy.h"
-#include "indicesol.h"
-c======================================================================
-C
-      INTEGER i, it
-      REAL pctsrf(klon,nbsrf) !Pourcentage de sol (f(nature du sol))
-      REAL ftsol(klon,nbsrf)  ! Temperature du sol pour le controle Rn
-c                             ! le cas echeant
-      REAL masktr(klon,nbtr)  ! Masque de l'echange avec la surface
-c                                 (possible => 1 )
-      REAL fshtr(klon,nbtr)  ! Flux surfacique dans le reservoir de sol
-      REAL hsoltr(nbtr)      ! Epaisseur equivalente du reservoir de sol
-      REAL tautr(nbtr)       ! Constante de decroissance radioactive
-      REAL vdeptr(nbtr)      ! Vitesse de depot sec dans la couche Brownienne
-      REAL scavtr(nbtr)      ! Coefficient de lessivage
-      REAL s
-C
-      WRITE(*,'(''PASSAGE initrrnpb ...'',$)')
-      print*,'nbtr= ',nbtr
-      print*,'nbsrf= ',nbsrf
-      print*,'klon= ',klon
-C
-C Puis les initialisation specifiques a chaque traceur (pour le moment, Rn222)
-C
-C
-C Radon it = 1
-c
-      IF ( nbtr .LE. 0 ) STOP 'initrrnpb pas glop pas glop' 
-      it = 1
-      s = 1.E4  !  Source: atome par m2
-      hsoltr(it) = 0.1      ! Hauteur equivalente du reservoir : 
-c                              1 m * porosite 0.1
-      tautr(it) = 4.765E5  ! Decroissance du radon, secondes
-cAA
-c      tautr(it) = 4.765E55  ! Decroissance du radon,infinie
-cAA
-      vdeptr(it) = 0. ! Pas de depot sec pour le radon
-      scavtr(it) = 0. ! Pas de lessivage pour le radon
+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"
+!======================================================================
 
-      print*, '-------------- SOURCE DU RADON ------------------------ '
-      print*,'it = ',it
-      print*,'Source : ', s
-      print*,'Hauteur equivalente du reservoir de sol: ',hsoltr(it) 
-      print*,'Decroissance (s): ', tautr(it)
-      print*,'Vitesse de depot sec: ',vdeptr(it) 
-      print*,'Facteur de lessivage: ',scavtr(it)
+  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
 
-      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)
+  WRITE(*,*)'PASSAGE initrrnpb ...'
+!
+! Radon it = 1
+!----------------
+  IF ( nbtr .LE. 0 ) STOP '**PHYTRAC:initrrnpb:** nbtr < 0; verifier RN dans traceur.def' 
+  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)
 
-cAA
-cAA quelques tests
-cAA POur l'instant le pctsrf(i,3) = 1.0 
-cAA lorsqu'il ya de la terre mias ne prend aucune autre valeur
-cAA il n'est donc pas necessaire de multiplier fshtr par pctsrf
-cAA 
-c       print*, '------------------------------------------ '
-c        print*, 'masktr(',i,it,')= ',masktr(i,it)
-c        print*, 'fshtr(',i,it,')= ',fshtr(i,it)
-c        print*, 'pctsrf(',i,',1)= ',pctsrf(i,1)
-c        print*, 'pctsrf(',i,',2)= ',pctsrf(i,2)
-c        print*, 'pctsrf(',i,',3)= ',pctsrf(i,3)
-c        print*, 'pctsrf(',i,',4)= ',pctsrf(i,4)
-c        print*, 's = ',s
-c        print*, '------------------------------------------ '
+  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 ) STOP '**PHYTRAC**:initrrnpb:** nbtr <= 1; verifier PB dans traceur.def' 
+  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)
 
-      END DO
-C
-C 210Pb it = 2
-C
-      IF ( nbtr .LE. 1 ) STOP 'initrrnpb pas glop pas glop' 
-      it = 2
-      s = 0. !  Pas de source !!!
-      hsoltr(it) = 10.     ! Hauteur equivalente du reservoir 
-c                              a partir duquel le
-c                              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
-      print*, '-------------- SOURCE DU PLOMB ------------------------ '
-      print*,'it = ',it
-      print*,'Source : ', s
-      print*,'Hauteur equivalente du reservoir : ',hsoltr(it) 
-      print*,'Decroissance (s): ', tautr(it)
-      print*,'Vitesse de depot sec: ',vdeptr(it) 
-      print*,'Facteur de lessivage: ',scavtr(it)
-c
-      WRITE(*,*) 'initialisation rnpb ok'
-c
-      RETURN
-      END
+  WRITE(*,*) 'Initialisation RN et PB ok'
+
+END SUBROUTINE initrrnpb