source: LMDZ4/trunk/libf/phylmd/initrrnpb.F90 @ 1404

Last change on this file since 1404 was 1403, checked in by Laurent Fairhead, 14 years ago

Merged LMDZ4V5.0-dev branch changes r1292:r1399 to trunk.

Validation:
Validation consisted in compiling the HEAD revision of the trunk,
LMDZ4V5.0-dev branch and the merged sources and running different
configurations on local and SX8 machines comparing results.

Local machine: bench configuration, 32x24x11, gfortran

  • IPSLCM5A configuration (comparison between trunk and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent
  • MH07 configuration, new physics package (comparison between LMDZ4V5.0-dev branch and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent

SX8 machine (brodie), 96x95x39 on 4 processors:

  • IPSLCM5A configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent
  • MH07 configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent

Changes to the makegcm and create_make_gcm scripts to take into account
main programs in F90 files


Fusion de la branche LMDZ4V5.0-dev (r1292:r1399) au tronc principal

Validation:
La validation a consisté à compiler la HEAD de le trunk et de la banche
LMDZ4V5.0-dev et les sources fusionnées et de faire tourner le modéle selon
différentes configurations en local et sur SX8 et de comparer les résultats

En local: 32x24x11, config bench/gfortran

  • pour une config IPSLCM5A (comparaison tronc/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux (à part sur RN et Pb)
    • fichiers histoire égaux
  • pour une config nlle physique (MH07) (comparaison LMDZ4v5.0-dev/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux
    • fichiers histoire équivalents

Sur brodie, 96x95x39 sur 4 proc:

  • pour une config IPSLCM5A:
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc
  • pour une config MH07
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc

Changement sur makegcm et create_make-gcm pour pouvoir prendre en compte des
programmes principaux en *F90

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 3.8 KB
Line 
1!
2! $Id: initrrnpb.F90 1403 2010-07-01 09:02:53Z fairhead $
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  CHARACTER (LEN=20) :: modname='initrrnpb'
42  CHARACTER (LEN=80) :: abort_message
43
44
45  WRITE(*,*)'PASSAGE initrrnpb ...'
46!
47! Radon it = 1
48!----------------
49  IF ( nbtr .LE. 0 ) then
50    abort_message = '**PHYTRAC:initrrnpb:** nbtr < 0; verifier RN dans traceur.def'
51    CALL abort_gcm (modname,abort_message,1)
52  ENDIF
53  it = 1
54  s = 1.E4             ! Source: atome par m2
55  hsoltr(it) = 0.1     ! Hauteur equivalente du reservoir :
56                       ! 1 m * porosite 0.1
57  tautr(it) = 4.765E5  ! Decroissance du radon, secondes
58  vdeptr(it) = 0.      ! Pas de depot sec pour le radon
59  scavtr(it) = 0.      ! Pas de lessivage pour le radon
60 
61  WRITE(*,*)'-------------- SOURCE DU RADON ------------------------ '
62  WRITE(*,*)'it = ',it
63  WRITE(*,*)'Source : ', s
64  WRITE(*,*)'Hauteur equivalente du reservoir de sol: ',hsoltr(it)
65  WRITE(*,*)'Decroissance (s): ', tautr(it)
66  WRITE(*,*)'Vitesse de depot sec: ',vdeptr(it)
67  WRITE(*,*)'Facteur de lessivage: ',scavtr(it)
68
69  DO i = 1,klon
70     masktr(i,it) = 0.
71     IF ( NINT(pctsrf(i,1)) .EQ. 1 ) masktr(i,it) = 1.
72     fshtr(i,it) = s * masktr(i,it)
73  END DO
74!
75! 210Pb it = 2
76!----------------
77  IF ( nbtr .LE. 1 ) THEN
78    abort_message='**PHYTRAC**:initrrnpb:** nbtr <= 1; verifier PB dans traceur.def'
79    CALL abort_gcm (modname,abort_message,1)
80  ENDIF
81  it = 2
82  s = 0.                ! Pas de source
83  hsoltr(it) = 10.      ! Hauteur equivalente du reservoir
84                        ! a partir duquel le depot Brownien a lieu
85  tautr(it) = 1.028E9   ! Decroissance du Pb210, secondes
86  vdeptr(it) = 1.E-3    ! 1 mm/s pour le 210Pb
87  scavtr(it) =  .5      ! Lessivage du Pb210
88  DO i = 1,klon
89     masktr(i,it) = 1.  ! Le depot sec peut avoir lieu partout
90     fshtr(i,it) = s * masktr(i,it)
91  END DO
92  WRITE(*,*)'-------------- SOURCE DU PLOMB ------------------------ '
93  WRITE(*,*)'it = ',it
94  WRITE(*,*)'Source : ', s
95  WRITE(*,*)'Hauteur equivalente du reservoir : ',hsoltr(it)
96  WRITE(*,*)'Decroissance (s): ', tautr(it)
97  WRITE(*,*)'Vitesse de depot sec: ',vdeptr(it)
98  WRITE(*,*)'Facteur de lessivage: ',scavtr(it)
99
100  WRITE(*,*) 'Initialisation RN et PB ok'
101
102END SUBROUTINE initrrnpb
Note: See TracBrowser for help on using the repository browser.