Changeset 3852 for LMDZ6/branches/LMDZ-tracers/libf/phylmd/traclmdz_mod.F90

Timestamp:

Feb 22, 2021, 5:28:31 PM (3 years ago)

Author:

dcugnet

Message:

Extension of the tracers management.

The tracers files can be:

1) "traceur.def": old format, with:

• the number of tracers on the first line
• one line for each tracer: <tracer name> <hadv> <vadv> [<parent name>]

2) "tracer.def": new format with one section each model component.
3) "tracer_<name>.def": new format with a single section.

The formats 2 and 3 reading is driven by the "type_trac" key, which can be a

coma-separated list of components.

• Format 2: read the sections from the "tracer.def" file.
• format 3: read one section each "tracer_<section name>.def" file.
• the first line of a section is "&<section name>
• the other lines start with a tracer name followed by <key>=<val> pairs.
• the "default" tracer name is reserved ; the other tracers of the section inherit its <key>=<val>, except for the keys that are redefined locally.

This format helps keeping the tracers files compact, thanks to the "default"
special tracer and the three levels of factorization:

• on the tracers names: a tracer name can be a coma-separated list of tracers => all the tracers of the list have the same <key>=<val> properties
• on the parents names: the value of the "parent" property can be a coma-separated list of tracers => only possible for geographic tagging tracers
• on the phases: the property "phases" is [g](l][s] (gas/liquid/solid)

Read information is stored in the vector "tracers(:)", of derived type "tra".

"isotopes_params.def" is a similar file, with one section each isotopes family.
It contains a database of isotopes properties ; if there are second generation
tracers (isotopes), the corresponding sections are read.

Read information is stored in the vector "isotopes(:)", of derived type "iso".

The "getKey" function helps to get the values of the parameters stored in
"tracers" or "isotopes".

File:

• LMDZ6/branches/LMDZ-tracers/libf/phylmd/traclmdz_mod.F90 (modified) (7 diffs)

LMDZ6/branches/LMDZ-tracers/libf/phylmd/traclmdz_mod.F90

@@ -67 +67 @@
     
     USE dimphy
-    USE infotrac_phy
+    USE infotrac_phy, ONLY: nbtr
     
     ! Input argument
@@ -89 +89 @@
     ! Initialization of the tracers should be done here only for those not found in the restart file.
     USE dimphy
-    USE infotrac_phy
+    USE infotrac_phy, ONLY: tracers, nqo, nbtr, niadv, pbl_flg, conv_flg
     USE regr_pr_comb_coefoz_m, ONLY: alloc_coefoz
     USE press_coefoz_m, ONLY: press_coefoz
@@ -175 +175 @@
 !!       iiq=niadv(it+2)                                                            ! jyg
        iiq=niadv(it+nqo)                                                            ! jyg
-       IF ( tname(iiq) == "RN" ) THEN                                               
-          id_rn=it ! radon
-       ELSE IF ( tname(iiq) == "PB") THEN
-          id_pb=it ! plomb
+       !-----------------------------------------------------------------------
+       SELECT CASE(tracers(iiq)%name)
+       !-----------------------------------------------------------------------
+         CASE("RN"); id_rn=it ! radon
+       !-----------------------------------------------------------------------
+         CASE("PB"); id_pb=it ! plomb
 ! RomP >>> profil initial de PB210 
      open (ilesfil2,file='prof.pb210',status='old',iostat=irr2)
@@ -198 +200 @@
      ENDIF
 ! RomP <<<
-       ELSE IF ( tname(iiq) == "Aga" .OR. tname(iiq) == "AGA" ) THEN
-          ! Age of stratospheric air
-          id_aga=it
+       !-----------------------------------------------------------------------
+         CASE("Aga","AGA");           id_aga   = it ! Age of stratospheric air
           radio(id_aga) = .FALSE.
           aerosol(id_aga) = .FALSE.
@@ -213 +214 @@
              lev_1p5km=klev/2
           END IF
-       ELSE IF ( tname(iiq) == "BE" .OR. tname(iiq) == "Be" .OR.  &
-            tname(iiq) == "BE7" .OR. tname(iiq) == "Be7" ) THEN  
-          ! Recherche du Beryllium 7
-          id_be=it
+       !-----------------------------------------------------------------------
+         CASE("BE","Be","BE7","Be7"); id_be    = it ! Recherche du Beryllium 7
           ALLOCATE( srcbe(klon,klev) )
           radio(id_be) = .TRUE.
@@ -243 +242 @@
      ENDIF
 ! RomP <<<
-       ELSE IF (tname(iiq)=="O3" .OR. tname(iiq)=="o3") THEN
-          ! Recherche de l'ozone : parametrization de la chimie par Cariolle
-          id_o3=it
-          CALL alloc_coefoz   ! allocate ozone coefficients
-          CALL press_coefoz   ! read input pressure levels
-       ELSE IF ( tname(iiq) == "pcsat" .OR. tname(iiq) == "Pcsat" ) THEN
-          id_pcsat=it
-       ELSE IF ( tname(iiq) == "pcocsat" .OR. tname(iiq) == "Pcocsat" ) THEN
-          id_pcocsat=it
-       ELSE IF ( tname(iiq) == "pcq" .OR. tname(iiq) == "Pcq" ) THEN
-          id_pcq=it
-       ELSE IF ( tname(iiq) == "pcs0" .OR. tname(iiq) == "Pcs0" ) THEN
-          id_pcs0=it
-          conv_flg(it)=0 ! No transport by convection for this tracer
-       ELSE IF ( tname(iiq) == "pcos0" .OR. tname(iiq) == "Pcos0" ) THEN
-          id_pcos0=it
-          conv_flg(it)=0 ! No transport by convection for this tracer
-       ELSE IF ( tname(iiq) == "pcq0" .OR. tname(iiq) == "Pcq0" ) THEN
-          id_pcq0=it
-          conv_flg(it)=0 ! No transport by convection for this tracer
-       ELSE
-          WRITE(lunout,*) 'This is an unknown tracer in LMDZ : ', trim(tname(iiq))
-       END IF
+       !-----------------------------------------------------------------------
+         CASE("O3","o3");           id_o3      = it
+           ! Recherche de l'ozone : parametrization de la chimie par Cariolle
+           CALL alloc_coefoz   ! allocate ozone coefficients
+           CALL press_coefoz   ! read input pressure levels
+       !-----------------------------------------------------------------------
+         CASE("pcsat"  ,"Pcsat");   id_pcsat   = it
+       !-----------------------------------------------------------------------
+         CASE("pcocsat","Pcocsat"); id_pcocsat = it
+       !-----------------------------------------------------------------------
+         CASE("pcq"    ,"Pcq");     id_pcq     = it
+       !-----------------------------------------------------------------------
+         CASE("pcs0"   ,"Pcs0");    id_pcs0    = it
+           conv_flg(it)=0 ! No transport by convection for this tracer
+       !-----------------------------------------------------------------------
+         CASE("pcos0"  ,"Pcos0");   id_pcos0   = it
+           conv_flg(it)=0 ! No transport by convection for this tracer
+       !-----------------------------------------------------------------------
+         CASE("pcq0"   ,"Pcq0");    id_pcq0    = it
+           conv_flg(it)=0 ! No transport by convection for this tracer
+       !-----------------------------------------------------------------------
+         CASE DEFAULT
+           WRITE(lunout,*) 'This is an unknown tracer in LMDZ : ', trim(tracers(iiq)%name)
+       !-----------------------------------------------------------------------
+       END SELECT
+       !-----------------------------------------------------------------------
     END DO
 
@@ -309 +311 @@
        IF (zero) THEN
           ! The tracer was not found in restart file or it was equal zero everywhere.
-          WRITE(lunout,*) "The tracer ",trim(tname(iiq))," will be initialized"
+          WRITE(lunout,*) "The tracer ",trim(tracers(iiq)%name)," will be initialized"
           IF (it==id_pcsat .OR. it==id_pcq .OR. &
                it==id_pcs0 .OR. it==id_pcq0) THEN