Changeset 1191 for LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar

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/dyn3dpar/infotrac.F90 (modified) (8 diffs)

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/infotrac.F90

@@ -21 +21 @@
   INTEGER, ALLOCATABLE, DIMENSION(:), SAVE    :: niadv ! equivalent dyn / physique
 
-! Variables for INCA
+! conv_flg(it)=0 : convection desactivated for tracer number it 
   INTEGER, ALLOCATABLE, DIMENSION(:), SAVE  :: conv_flg
+! pbl_flg(it)=0  : boundary layer diffusion desactivaded for tracer number it 
   INTEGER, ALLOCATABLE, DIMENSION(:), SAVE  :: pbl_flg
 
+  CHARACTER(len=4),SAVE :: type_trac
+ 
 CONTAINS
 
@@ -80 +83 @@
     descrq(20)='SLP'
     descrq(30)='PRA'
+    
+
+    IF (config_inca=='none') THEN
+       type_trac='lmdz'
+    ELSE
+       type_trac='inca'
+    END IF
 
 !-----------------------------------------------------------------------
@@ -87 +97 @@
 !
 !-----------------------------------------------------------------------
-    IF (config_inca == 'none') THEN
+    IF (type_trac == 'lmdz') THEN
        OPEN(90,file='traceur.def',form='formatted',status='old', iostat=ierr)
        IF(ierr.EQ.0) THEN
@@ -109 +119 @@
     END IF
 !
-! Allocate variables depending on nqtrue
+! Allocate variables depending on nqtrue and nbtr
 !
     ALLOCATE(tnom_0(nqtrue), hadv(nqtrue), vadv(nqtrue))
-
-    IF (config_inca /= 'none') THEN
-       ! Varaibles only needed in case of INCA
-       ALLOCATE(conv_flg(nbtr), pbl_flg(nbtr), tracnam(nbtr))
-    END IF
-       
+    ALLOCATE(conv_flg(nbtr), pbl_flg(nbtr), tracnam(nbtr))
+    conv_flg(:) = 1 ! convection activated for all tracers
+    pbl_flg(:)  = 1 ! boundary layer activated for all tracers
+
 !-----------------------------------------------------------------------
 ! 2)     Choix  des schemas d'advection pour l'eau et les traceurs
@@ -144 +152 @@
 !    Get choice of advection schema from file tracer.def or from INCA
 !---------------------------------------------------------------------
-    IF (config_inca == 'none') THEN
+    IF (type_trac == 'lmdz') THEN
        IF(ierr.EQ.0) THEN
           ! Continue to read tracer.def
@@ -172 +180 @@
        END DO
 
-    ELSE  ! config_inca='aero' ou 'chem'
+    ELSE  ! type_trac=inca : config_inca='aero' ou 'chem'
 ! le module de chimie fournit les noms des traceurs
 ! et les schemas d'advection associes.
@@ -191 +199 @@
        END DO
 
-    END IF ! config_inca
+    END IF ! type_trac
 
 !-----------------------------------------------------------------------
@@ -319 +327 @@
 !
     DEALLOCATE(tnom_0, hadv, vadv)
-    IF (config_inca /= 'none') DEALLOCATE(tracnam)
+    DEALLOCATE(tracnam)
 
 999 FORMAT (i2,1x,i2,1x,a8)