Changeset 4046 for LMDZ6/trunk/libf/phylmd/Dust

Timestamp:

Dec 15, 2021, 11:18:49 PM (4 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".

Location:

LMDZ6/trunk/libf/phylmd/Dust

Files:

• phytracr_spl_mod.F90 (modified) (1 diff)
• splaeropt_5wv_rrtm.F90 (modified) (1 diff)
• splaeropt_6bands_rrtm.F90 (modified) (1 diff)
• splaeropt_lw_rrtm.F90 (modified) (1 diff)
• splaerosol_optic_rrtm.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90

@@ -1416 +1416 @@
        !print *,nbtr
        do it=1,nbtr
-        print *, it, tname(it+nqo)
-        if (tname(it+nqo) == 'PREC' ) then
-            id_prec=it
-        endif
-        if (tname(it+nqo) == 'FINE' ) then
-            id_fine=it
-        endif
-        if (tname(it+nqo) == 'COSS' ) then
-            id_coss=it
-        endif
-        if (tname(it+nqo) == 'CODU' ) then
-            id_codu=it
-        endif
-        if (tname(it+nqo) == 'SCDU' ) then
-            id_scdu=it
-        endif
+        print *, it, tracers(it+nqo)%name
+        SELECT CASE(tracers(it+nqo)%name)
+          CASE('PREC'); id_prec=it
+          CASE('FINE'); id_fine=it
+          CASE('COSS'); id_coss=it
+          CASE('CODU'); id_codu=it
+          CASE('SCDU'); id_scdu=it
+        END SELECT
        enddo
        ! check consistency with dust emission scheme:

LMDZ6/trunk/libf/phylmd/Dust/splaeropt_5wv_rrtm.F90

@@ -113 +113 @@
 
   DO itr=1,nbtr    !--loop over tracers  
-
-    IF (tname(itr+nqo)=='PREC') THEN       !--fine mode accumulation mode
-      CYCLE
-    ELSE IF (tname(itr+nqo)=='FINE') THEN  !--fine mode accumulation mode
-      soluble=.TRUE.
-      spsol=1
-      aerindex=1
-    ELSE IF (tname(itr+nqo)=='COSS') THEN  !--coarse mode sea salt
-      soluble=.TRUE.
-      spsol=2
-      aerindex=2
-    ELSE IF (tname(itr+nqo)=='CODU') THEN  !--coarse mode dust
-      soluble=.FALSE.
-      spinsol=1
-      aerindex=3
-    ELSE IF (tname(itr+nqo)=='SCDU') THEN  !--super coarse mode dust
-      soluble=.FALSE.
-      spinsol=2
-      aerindex=4
-    ELSE
-       CALL abort_physic(modname,'I cannot do aerosol optics for '//tname(itr+nqo),1)
-    ENDIF
+    SELECT CASE(tracers(itr+nqo)%name)
+      CASE('PREC'); CYCLE                                  !--precursor
+      CASE('FINE'); soluble=.TRUE.;  spsol=1; aerindex=1   !--fine mode accumulation mode
+      CASE('COSS'); soluble=.TRUE.;  spsol=2; aerindex=2   !--coarse mode sea salt
+      CASE('CODU'); soluble=.FALSE.; spsol=1; aerindex=3   !--coarse mode dust
+      CASE('SCDU'); soluble=.FALSE.; spsol=2; aerindex=4   !--super coarse mode dust
+      CASE DEFAULT; CALL abort_physic(modname,'I cannot do aerosol optics for '//tracers(itr+nqo)%name,1)
+    END SELECT
 
     DO la=1,las

LMDZ6/trunk/libf/phylmd/Dust/splaeropt_6bands_rrtm.F90

@@ -165 +165 @@
   cg_ae(:,:,:,:)=0.
     
-  DO itr=1, nbtr
-
-    IF (tname(itr+nqo)=='PREC') THEN       !--precursor
-      CYCLE
-    ELSE IF (tname(itr+nqo)=='FINE') THEN  !--fine mode accumulation mode
-      soluble=.TRUE.
-      spsol=1
-      aerindex=1
-    ELSE IF (tname(itr+nqo)=='COSS') THEN  !--coarse mode sea salt
-      soluble=.TRUE.
-      spsol=2
-      aerindex=2
-    ELSE IF (tname(itr+nqo)=='CODU') THEN  !--coarse mode dust
-      soluble=.FALSE.
-      spinsol=1
-      aerindex=3
-    ELSE IF (tname(itr+nqo)=='SCDU') THEN  !--super coarse mode dust
-      soluble=.FALSE.
-      spinsol=2
-      aerindex=4
-    ELSE
-       CALL abort_physic(modname,'I cannot do aerosol optics for '//tname(itr+nqo),1)
-    ENDIF
+  DO itr=1,nbtr    !--loop over tracers  
+    SELECT CASE(tracers(itr+nqo)%name)
+      CASE('PREC'); CYCLE                                  !--precursor
+      CASE('FINE'); soluble=.TRUE.;  spsol=1; aerindex=1   !--fine mode accumulation mode
+      CASE('COSS'); soluble=.TRUE.;  spsol=2; aerindex=2   !--coarse mode sea salt
+      CASE('CODU'); soluble=.FALSE.; spsol=1; aerindex=3   !--coarse mode dust
+      CASE('SCDU'); soluble=.FALSE.; spsol=2; aerindex=4   !--super coarse mode dust
+      CASE DEFAULT; CALL abort_physic(modname,'I cannot do aerosol optics for '//tracers(itr+nqo)%name,1)
+    END SELECT
 
     IF (soluble) THEN ! For aerosol soluble components

LMDZ6/trunk/libf/phylmd/Dust/splaeropt_lw_rrtm.F90

@@ -55 +55 @@
     !
     DO itr=1,nbtr
-      !
-      IF (tname(itr+nqo)=='PREC') THEN       !--precursor
-        CYCLE
-      ELSE IF (tname(itr+nqo)=='FINE') THEN  !--fine mode accumulation mode
-        CYCLE
-      ELSE IF (tname(itr+nqo)=='COSS') THEN  !--coarse mode sea salt
-        CYCLE
-      ELSE IF (tname(itr+nqo)=='CODU') THEN  !--coarse mode dust
-        spinsol=1
-      ELSE IF (tname(itr+nqo)=='SCDU') THEN  !--super coarse mode dust
-        spinsol=2
-      ELSE
-         CALL abort_physic(modname,'I cannot do aerosol optics for '//tname(itr+nqo),1)
-      ENDIF
+      SELECT CASE(tracers(itr+nqo)%name)
+        CASE('PREC','FINE''COSS'); CYCLE                   !--precursor or fine/coarde accumulation mode
+        CASE('CODU'); spinsol=1                            !--coarse mode dust
+        CASE('SCDU'); spinsol=2                            !--super coarse mode dust
+        CASE DEFAULT; CALL abort_physic(modname,'I cannot do aerosol optics for '//tracers(itr+nqo)%name,1)
+      END SELECT
       !
       DO inu=1,NLW

LMDZ6/trunk/libf/phylmd/Dust/splaerosol_optic_rrtm.F90

@@ -51 +51 @@
   !
   DO itr=1,nbtr
-    IF (tname(itr+nqo)=='FINE') THEN
+    IF (tracers(itr+nqo)%name=='FINE') THEN
       mass_solu_aero(:,:)    = tr_seri(:,:,itr)
       mass_solu_aero_pi(:,:) = tr_seri(:,:,itr)