Changeset 4393 for LMDZ6/trunk/libf/phylmd

Timestamp:

Jan 23, 2023, 5:39:48 PM (23 months ago)

Author:

dcugnet

Message:

No longer use "strLower" for the SELECT CASE statements in traclmdz.

File:

• LMDZ6/trunk/libf/phylmd/traclmdz_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/traclmdz_mod.F90

@@ -96 +96 @@
     USE indice_sol_mod
     USE print_control_mod, ONLY: lunout
-    USE strings_mod, ONLY: strLower
 
 ! Input variables
@@ -177 +176 @@
        IF(.NOT.(tracers(iq)%isInPhysics)) CYCLE
        it = it+1
-       SELECT CASE(strLower(tracers(iq)%name))
-         CASE("rn");      id_rn     = it ! radon
-         CASE("pb");      id_pb     = it ! plomb
-         CASE("aga");     id_aga    = it ! Age of stratospheric air
-         CASE("be","be7");id_be     = it ! Recherche du Beryllium 7
-         CASE("o3");      id_o3     = it ! Recherche de l'ozone
-         CASE("pcsat");   id_pcsat  = it
-         CASE("pcocsat"); id_pcocsat= it
-         CASE("pcq");     id_pcq    = it
-         CASE("pcs0");    id_pcs0   = it
-         CASE("pcos0");   id_pcos0  = it
-         CASE("pcq0");    id_pcq0   = it
+       SELECT CASE(tracers(iq)%name)
+         CASE("RN");                  id_rn     = it ! radon
+         CASE("PB");                  id_pb     = it ! plomb
+         CASE("Aga","AGA");           id_aga    = it ! Age of stratospheric air
+         CASE("Be","BE","Be7","BE7"); id_be     = it ! Recherche du Beryllium 7
+         CASE("o3","O3");             id_o3     = it ! Recherche de l'ozone
+         CASE("pcsat",  "Pcsat");     id_pcsat  = it
+         CASE("pcocsat","Pcocsat");   id_pcocsat= it
+         CASE("pcq",    "Pcq");       id_pcq    = it
+         CASE("pcs0",   "Pcs0");      id_pcs0   = it
+         CASE("pcos0",  "Pcos0");     id_pcos0  = it
+         CASE("pcq0",   "Pcq0");      id_pcq0   = it
          CASE DEFAULT
            WRITE(lunout,*) 'This is an unknown tracer in LMDZ : ', trim(tracers(iq)%name)
        END SELECT
 
-       SELECT CASE(strLower(tracers(iq)%name))
-         CASE("pb")                        !--- RomP >>> profil initial de PB210 
+       SELECT CASE(tracers(iq)%name)
+         CASE("PB")                        !--- RomP >>> profil initial de PB210 
            OPEN(ilesfil2,file='prof.pb210',status='old',iostat=irr2)
            IF(irr2 == 0) THEN 
@@ -205 +204 @@
              WRITE(lunout,*)'prof.pb210 does not exist: use restart values'
            END IF
-         CASE("aga")
+         CASE("Aga","AGA")
            radio(id_aga) = .FALSE.
            aerosol(id_aga) = .FALSE.
@@ -217 +216 @@
               lev_1p5km=klev/2
            END IF
-         CASE("be","be7")
+         CASE("Be","BE","Be7","BE7")
            ALLOCATE( srcbe(klon,klev) )
            radio(id_be) = .TRUE.
@@ -234 +233 @@
              WRITE(lunout,*)'Prof. Be7 does not exist: use restart values'
            END IF
-         CASE("o3")                         !--- Parametrisation par la chimie de Cariolle
+         CASE("o3","O3")                    !--- Parametrisation par la chimie de Cariolle
            CALL alloc_coefoz                !--- Allocate ozone coefficients
            CALL press_coefoz                !--- Read input pressure levels
-         CASE("pcs0","pcos0","pcq0")
+         CASE("pcs0","Pcs0", "pcos0","Pcos0", "pcq0","Pcq0")
            conv_flg(it)=0                   !--- No transport by convection for this tracer
        END SELECT