Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/infotrac.F90

    r1180 r1191  
    2121  INTEGER, ALLOCATABLE, DIMENSION(:), SAVE    :: niadv ! equivalent dyn / physique
    2222
    23 ! Variables for INCA
     23! conv_flg(it)=0 : convection desactivated for tracer number it
    2424  INTEGER, ALLOCATABLE, DIMENSION(:), SAVE  :: conv_flg
     25! pbl_flg(it)=0  : boundary layer diffusion desactivaded for tracer number it
    2526  INTEGER, ALLOCATABLE, DIMENSION(:), SAVE  :: pbl_flg
    2627
     28  CHARACTER(len=4),SAVE :: type_trac
     29 
    2730CONTAINS
    2831
     
    8083    descrq(20)='SLP'
    8184    descrq(30)='PRA'
     85   
     86
     87    IF (config_inca=='none') THEN
     88       type_trac='lmdz'
     89    ELSE
     90       type_trac='inca'
     91    END IF
    8292
    8393!-----------------------------------------------------------------------
     
    8797!
    8898!-----------------------------------------------------------------------
    89     IF (config_inca == 'none') THEN
     99    IF (type_trac == 'lmdz') THEN
    90100       OPEN(90,file='traceur.def',form='formatted',status='old', iostat=ierr)
    91101       IF(ierr.EQ.0) THEN
     
    109119    END IF
    110120!
    111 ! Allocate variables depending on nqtrue
     121! Allocate variables depending on nqtrue and nbtr
    112122!
    113123    ALLOCATE(tnom_0(nqtrue), hadv(nqtrue), vadv(nqtrue))
    114 
    115     IF (config_inca /= 'none') THEN
    116        ! Varaibles only needed in case of INCA
    117        ALLOCATE(conv_flg(nbtr), pbl_flg(nbtr), tracnam(nbtr))
    118     END IF
    119        
     124    ALLOCATE(conv_flg(nbtr), pbl_flg(nbtr), tracnam(nbtr))
     125    conv_flg(:) = 1 ! convection activated for all tracers
     126    pbl_flg(:)  = 1 ! boundary layer activated for all tracers
     127
    120128!-----------------------------------------------------------------------
    121129! 2)     Choix  des schemas d'advection pour l'eau et les traceurs
     
    144152!    Get choice of advection schema from file tracer.def or from INCA
    145153!---------------------------------------------------------------------
    146     IF (config_inca == 'none') THEN
     154    IF (type_trac == 'lmdz') THEN
    147155       IF(ierr.EQ.0) THEN
    148156          ! Continue to read tracer.def
     
    172180       END DO
    173181
    174     ELSE  ! config_inca='aero' ou 'chem'
     182    ELSE  ! type_trac=inca : config_inca='aero' ou 'chem'
    175183! le module de chimie fournit les noms des traceurs
    176184! et les schemas d'advection associes.
     
    191199       END DO
    192200
    193     END IF ! config_inca
     201    END IF ! type_trac
    194202
    195203!-----------------------------------------------------------------------
     
    319327!
    320328    DEALLOCATE(tnom_0, hadv, vadv)
    321     IF (config_inca /= 'none') DEALLOCATE(tracnam)
     329    DEALLOCATE(tracnam)
    322330
    323331999 FORMAT (i2,1x,i2,1x,a8)
Note: See TracChangeset for help on using the changeset viewer.