Changeset 4046 for LMDZ6/trunk/libf/phylmd/traclmdz_mod.F90

Timestamp:

Dec 15, 2021, 11:18:49 PM (2 years ago)

Author:

dcugnet

Message:

First commit for new tracers.

• parser routines readTracFiles, strings_mod and tracer_types added in misc using revision 4 of ​https://svn.lmd.jussieu.fr/tracers-parser
• tested in sequential and parallel mode using ioipsl.
• for now, only two fields of "tracers(:)" derived type vector are used: "name" and "longName".

File:

• LMDZ6/trunk/libf/phylmd/traclmdz_mod.F90 (modified) (3 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
@@ -175 +176 @@
 !!       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
-! RomP >>> profil initial de PB210 
-     open (ilesfil2,file='prof.pb210',status='old',iostat=irr2)
-     IF (irr2 == 0) THEN 
-      read(ilesfil2,*) kradio2
-      print*,'number of levels for pb210 profile ',kradio2
-      do k=kradio2,1,-1
-       read (ilesfil2,*) plomb(:,k)
-      enddo
-      close(ilesfil2)
-      do k=1,klev
-       do i=1,klon
-         tr_seri(i,k,id_pb)=plomb(i,k)
-!!        print*, 'tr_seri',i,k,tr_seri(i,k,id_pb)
-        enddo
-      enddo
-     ELSE
-       print *, 'Prof.pb210 does not exist: use restart values'
-     ENDIF
-! RomP <<<
-       ELSE IF ( tname(iiq) == "Aga" .OR. tname(iiq) == "AGA" ) THEN
-          ! Age of stratospheric air
-          id_aga=it
-          radio(id_aga) = .FALSE.
-          aerosol(id_aga) = .FALSE.
-          pbl_flg(id_aga) = 0 
-          
-          ! Find the first model layer above 1.5km from the surface
-          IF (klev>=30) THEN
-             lev_1p5km=6   ! NB! This value is for klev=39
-          ELSE IF (klev>=10) THEN
-             lev_1p5km=5   ! NB! This value is for klev=19
-          ELSE
-             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
-          ALLOCATE( srcbe(klon,klev) )
-          radio(id_be) = .TRUE.
-          aerosol(id_be) = .TRUE. ! le Be est un aerosol
-!jyg le 13/03/2013 ; ajout de pplay en argument de init_be
-!!!          CALL init_be(pctsrf,masktr(:,id_be),tautr(id_be),vdeptr(id_be),scavtr(id_be),srcbe)
-          CALL init_be(pctsrf,pplay,masktr(:,id_be),tautr(id_be),vdeptr(id_be),scavtr(id_be),srcbe)
-          WRITE(lunout,*) 'Initialisation srcBe: OK'
-! RomP >>> profil initial de Be7
-      open (ilesfil,file='prof.be7',status='old',iostat=irr)
-      IF (irr == 0) THEN 
-       read(ilesfil,*) kradio
-       print*,'number of levels for Be7 profile ',kradio
-       do k=kradio,1,-1
-        read (ilesfil,*) beryllium(:,k)
-       enddo
-       close(ilesfil)
-       do k=1,klev
-         do i=1,klon
-          tr_seri(i,k,id_be)=beryllium(i,k)
-!!        print*, 'tr_seri',i,k,tr_seri(i,k,id_be)
-         enddo
-       enddo
-     ELSE
-       print *, 'Prof.Be7 does not exist: use restart values'
-     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
+print*,'###'//TRIM(tracers(iiq)%name)//'###'
+print*,'###'//TRIM(strLower(tracers(iiq)%name))//'###'
+       SELECT CASE(strLower(tracers(iiq)%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
+         CASE DEFAULT
+           WRITE(lunout,*) 'This is an unknown tracer in LMDZ : ', trim(tracers(iiq)%name)
+       END SELECT
+
+       SELECT CASE(strLower(tracers(iiq)%name))
+         CASE("pb")                        !--- RomP >>> profil initial de PB210 
+           OPEN(ilesfil2,file='prof.pb210',status='old',iostat=irr2)
+           IF(irr2 == 0) THEN 
+             READ(ilesfil2,*) kradio2
+             WRITE(lunout,*)'number of levels for pb210 profile ',kradio2
+             DO k=kradio2,1,-1; READ (ilesfil2,*) plomb(:,k); END DO
+             CLOSE(ilesfil2)
+             tr_seri(:,:,id_pb) = plomb(:,:)
+           ELSE
+             WRITE(lunout,*)'Prof. Pb210 does not exist: use restart values'
+           END IF
+         CASE("aga")
+           radio(id_aga) = .FALSE.
+           aerosol(id_aga) = .FALSE.
+           pbl_flg(id_aga) = 0 
+           ! Find the first model layer above 1.5km from the surface
+           IF (klev>=30) THEN
+              lev_1p5km=6                  !--- NB: This value is for klev=39
+           ELSE IF (klev>=10) THEN
+              lev_1p5km=5                  !--- NB: This value is for klev=19
+           ELSE
+              lev_1p5km=klev/2
+           END IF
+         CASE("be","be7")
+           ALLOCATE( srcbe(klon,klev) )
+           radio(id_be) = .TRUE.
+           aerosol(id_be) = .TRUE.         !--- Le Be est un aerosol
+           CALL init_be(pctsrf,pplay,masktr(:,id_be),tautr(id_be),vdeptr(id_be),scavtr(id_be),srcbe)
+           WRITE(lunout,*) 'Initialisation srcBe: OK'
+                                           !--- RomP >>> profil initial de Be7
+           OPEN(ilesfil,file='prof.be7',status='old',iostat=irr)
+           IF(irr == 0) THEN
+             READ(ilesfil,*) kradio
+             WRITE(lunout,*)'number of levels for Be7 profile ',kradio
+             DO k=kradio,1,-1; READ(ilesfil,*) beryllium(:,k); END DO
+             CLOSE(ilesfil)
+             tr_seri(:,:,id_be)=beryllium(:,:)
+           ELSE
+             WRITE(lunout,*)'Prof. Be7 does not exist: use restart values'
+           END IF
+         CASE("o3")                         !--- Parametrisation par la chimie de Cariolle
+           CALL alloc_coefoz                !--- Allocate ozone coefficients
+           CALL press_coefoz                !--- Read input pressure levels
+         CASE("pcs0","pcos0","pcq0")
+           conv_flg(it)=0                   !--- No transport by convection for this tracer
+       END SELECT
     END DO
 
@@ -309 +284 @@
        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