Changeset 4170 for LMDZ6/trunk/libf/dyn3d_common

Timestamp:

Jun 16, 2022, 8:16:59 PM (3 years ago)

Author:

dcugnet

Message:

The variable "types_trac" is the equivalent of "type_trac" in case multiple sections must be read
and used in "tracer.def" file.
Tests on the "type_trac" were replaced with tests on the vector "types_trac".
Most of the time, there are two components: 'lmdz' and a second one. The later has priority on 'lmdz'
and must be used for the tests. For more components, care must be taken to execute specific parts
of the code on the right tracers ; the tracers(:)%component has been created in that respect.

File:

• LMDZ6/trunk/libf/dyn3d_common/infotrac.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/dyn3d_common/infotrac.F90

@@ -14 +14 @@
    !=== FOR TRACERS:
    PUBLIC :: infotrac_init                                 !--- Initialization of the tracers
-   PUBLIC :: tracers, type_trac                            !--- Full tracers database, tracers type keyword
+   PUBLIC :: tracers, type_trac, types_trac                !--- Full tracers database, tracers type keyword
    PUBLIC :: nqtot,   nbtr,   nqo,   nqCO2,   nqtottr      !--- Main dimensions
    PUBLIC :: conv_flg, pbl_flg                             !--- Convection & boundary layer activation keys
@@ -108 +108 @@
                                     nqtottr, &                  !--- Number of tracers passed to phytrac (TO BE DELETED ?)
                                     nqCO2                       !--- Number of tracers of CO2  (ThL)
-   CHARACTER(LEN=maxlen),   SAVE :: type_trac                   !--- Keyword for tracers type
+   CHARACTER(LEN=maxlen),   SAVE :: type_trac                   !--- Keyword for tracers type(s)
+   CHARACTER(LEN=maxlen),   SAVE, ALLOCATABLE :: types_trac(:)  !--- Keyword for tracers type(s), parsed version
 
    !=== DERIVED TYPES EMBEDDING MOST INFORMATIONS ABOUT TRACERS AND ISOTOPES FAMILIES
@@ -175 +176 @@
    CHARACTER(LEN=3)      :: descrq(30)                               !--- Advection scheme description tags
    CHARACTER(LEN=maxlen) :: msg1                                     !--- String for messages
-   CHARACTER(LEN=maxlen), ALLOCATABLE :: str(:)                      !--- Temporary storage
    INTEGER :: fType                                                  !--- Tracers description file type ; 0: none
                                                                      !--- 1/2/3: "traceur.def"/"tracer.def"/"tracer_*.def"
@@ -201 +201 @@
     
    CALL msg('type_trac = "'//TRIM(type_trac)//'"', modname)
-   IF(strParse(type_trac, '|', str, n=nt)) CALL abort_gcm(modname,'can''t parse "type_trac = '//TRIM(type_trac)//'"',1)
+   IF(strParse(type_trac, '|', types_trac, n=nt)) CALL abort_gcm(modname,'can''t parse "type_trac = '//TRIM(type_trac)//'"',1)
 
    !---------------------------------------------------------------------------------------------------------------------------
@@ -207 +207 @@
    !---------------------------------------------------------------------------------------------------------------------------
       !--- MESSAGE ABOUT THE CHOSEN CONFIGURATION
-      msg1 = 'For type_trac = "'//TRIM(str(it))//'":'
-      SELECT CASE(str(it))
+      msg1 = 'For type_trac = "'//TRIM(types_trac(it))//'":'
+      SELECT CASE(types_trac(it))
          CASE('inca'); CALL msg(TRIM(msg1)//' coupling with INCA chemistry model, config_inca='//config_inca, modname)
          CASE('inco'); CALL msg(TRIM(msg1)//' coupling jointly with INCA and CO2 cycle',  modname)
@@ -215 +215 @@
          CASE('coag'); CALL msg(TRIM(msg1)//' tracers are treated for COAGULATION tests', modname)
          CASE('lmdz'); CALL msg(TRIM(msg1)//' tracers are treated in LMDZ only',          modname)
-         CASE DEFAULT; CALL abort_gcm(modname,'type_trac='//TRIM(str(it))//' not possible yet.',1)
+         CASE DEFAULT; CALL abort_gcm(modname,'type_trac='//TRIM(types_trac(it))//' not possible yet.',1)
       END SELECT
 
       !--- COHERENCE TEST BETWEEN "type_trac" AND "config_inca"
-      IF(ANY(['inca', 'inco'] == str(it)) .AND. ALL(['aero', 'aeNP', 'chem'] /= config_inca)) &
+      IF(ANY(['inca', 'inco'] == types_trac(it)) .AND. ALL(['aero', 'aeNP', 'chem'] /= config_inca)) &
          CALL abort_gcm(modname, 'Incoherence between type_trac and config_inca. Please modify "run.def"',1)
 
       !--- COHERENCE TEST BETWEEN "type_trac" AND PREPROCESSING KEYS
-      SELECT CASE(str(it))
+      SELECT CASE(types_trac(it))
          CASE('inca', 'inco')
 #ifndef INCA
@@ -244 +244 @@
    !--- DISABLE "config_inca" OPTION FOR A RUN WITHOUT "INCA" IF IT DIFFERS FROM "none"
    IF(fmsg('Setting config_inca="none" as you do not couple with INCA model', &
-           modname, ALL(str /= 'inca') .AND. ALL(str /= 'inco') .AND. config_inca /= 'none')) config_inca = 'none'
-
-   nqCO2 = 0; IF(ANY(str == 'inco')) nqCO2 = 1
+         modname, ALL(types_trac /= 'inca') .AND. ALL(types_trac /= 'inco') .AND. config_inca /= 'none')) config_inca = 'none'
+
+   nqCO2 = 0; IF(ANY(types_trac == 'inco')) nqCO2 = 1
 
 !==============================================================================================================================