Changeset 4368 for LMDZ6/branches/Ocean_skin/libf/phylmd/Dust/splaeropt_lw_rrtm.F90

Timestamp:

Dec 6, 2022, 12:01:16 AM (2 years ago)

Author:

lguez

Message:

Sync latest trunk changes to Ocean_skin

Location:

LMDZ6/branches/Ocean_skin

Files:

• . (modified) (1 prop)
• libf/phylmd/Dust/splaeropt_lw_rrtm.F90 (modified) (3 diffs)

LMDZ6/branches/Ocean_skin/libf/phylmd/Dust/splaeropt_lw_rrtm.F90

@@ -10 +10 @@
   USE dimphy
   USE aero_mod
-  USE infotrac_phy
+  USE infotrac_phy, ONLY: nqtot, nbtr, tracers
   USE phys_state_var_mod, ONLY : tau_aero_lw_rrtm
   USE YOERAD, ONLY : NLW
@@ -30 +30 @@
   INTEGER, PARAMETER :: naero=naero_soluble+naero_insoluble
   !
-  INTEGER inu, itr, spinsol
+  INTEGER inu, itr, iq, spinsol
   CHARACTER*20 modname
   !
@@ -54 +54 @@
     tau_aero_lw_rrtm = 0.0
     !
-    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
+    
+    itr = 0
+    DO iq = 1, nqtot
+      IF(.NOT.tracers(iq)%isInPhysics) CYCLE
+      itr = itr+1
+      SELECT CASE(tracers(iq)%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(iq)%name,1)
+      END SELECT
       !
       DO inu=1,NLW