Changeset 4124 for LMDZ6/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90

Timestamp:

Apr 8, 2022, 4:47:04 PM (2 years ago)

Author:

dcugnet

Message:

Remove solsym, ok_isotopes (=niso>0), ok_isotrac (=nzone>0)

File:

• LMDZ6/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90

@@ -20 +20 @@
        load_tmp8,load_tmp9,load_tmp10
 
-  USE infotrac_phy
+  USE infotrac_phy, ONLY: tracers, nqtot, nbtr
   USE YOMCST
 
@@ -91 +91 @@
   integer :: id_ASBCM, id_ASPOMM, id_ASSO4M, id_ASMSAM, id_CSSO4M, id_CSMSAM, id_SSSSM
   integer :: id_CSSSM, id_ASSSM, id_CIDUSTM, id_AIBCM, id_AIPOMM, id_ASNO3M, id_CSNO3M, id_CINO3M
-  INTEGER :: k, i
+  INTEGER :: k, i, iq, itr
 
   !--air density
@@ -107 +107 @@
      !--convert to ug m-3 unit for consistency with offline fields
      !
-     DO i=1,nbtr
-        SELECT CASE(trim(solsym(i)))
-           CASE ("ASBCM")
-              id_ASBCM = i
-           CASE ("ASPOMM") 
-              id_ASPOMM = i 
-           CASE ("ASSO4M")
-              id_ASSO4M = i 
-           CASE ("ASMSAM")
-              id_ASMSAM = i 
-           CASE ("CSSO4M")
-              id_CSSO4M = i 
-           CASE ("CSMSAM")
-              id_CSMSAM = i 
-           CASE ("SSSSM")
-              id_SSSSM = i 
-           CASE ("CSSSM")
-              id_CSSSM = i 
-           CASE ("ASSSM")
-              id_ASSSM = i 
-           CASE ("CIDUSTM")
-              id_CIDUSTM = i 
-           CASE ("AIBCM")
-              id_AIBCM = i 
-           CASE ("AIPOMM")
-              id_AIPOMM = i 
-           CASE ("ASNO3M")
-              id_ASNO3M = i 
-           CASE ("CSNO3M")
-              id_CSNO3M = i 
-           CASE ("CINO3M")
-              id_CINO3M = i 
-           END SELECT
-     ENDDO
+     itr = 0
+     DO iq = 1,nqtot
+        IF(.NOT. tracers(iq)%isInPhysics) CYCLE
+        itr = itr+1
+        SELECT CASE(trim(tracers(iq)%name))
+           CASE ("ASBCM");  id_ASBCM  = itr
+           CASE ("ASPOMM"); id_ASPOMM = itr
+           CASE ("ASSO4M"); id_ASSO4M = itr
+           CASE ("ASMSAM"); id_ASMSAM = itr
+           CASE ("CSSO4M"); id_CSSO4M = itr
+           CASE ("CSMSAM"); id_CSMSAM = itr
+           CASE ("SSSSM");  id_SSSSM  = itr
+           CASE ("CSSSM");  id_CSSSM  = itr
+           CASE ("ASSSM");  id_ASSSM  = itr
+           CASE ("CIDUSTM");id_CIDUSTM= itr
+           CASE ("AIBCM");  id_AIBCM  = itr
+           CASE ("AIPOMM"); id_AIPOMM = itr
+           CASE ("ASNO3M"); id_ASNO3M = itr
+           CASE ("CSNO3M"); id_CSNO3M = itr
+           CASE ("CINO3M"); id_CINO3M = itr
+        END SELECT
+     END DO
 
      bcsol(:,:)        =   tr_seri(:,:,id_ASBCM)                         *zrho(:,:)*1.e9  ! ASBCM