Changeset 5652 for LMDZ6

Timestamp:

May 14, 2025, 2:59:06 PM (3 months ago)

Author:

lebasn

Message:

StratAer?: Bug fix on class sulfur number to avoid problem when using Reprobus or other outside code.

File:

• LMDZ6/trunk/libf/phylmd/StratAer/strataer_local_var_mod.f90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/StratAer/strataer_local_var_mod.f90

@@ -259 +259 @@
     
     !--initialising dry diameters to geometrically spaced mass/volume (see Jacobson 1994)
-    mdw(1)=mdwmin
-    IF (V_rat.LT.1.62) THEN ! compensate for dip in second bin for lower volume ratio
-       mdw(2)=mdw(1)*2.**(1./3.)
-       DO it=3, nbtr_bin
-          mdw(it)=mdw(it-1)*V_rat**(1./3.)
-       ENDDO
+    IF(nbtr_bin < 3) THEN
+       WRITE(lunout,*) 'WARNING: There are less than 3 sulfur aerosol class, it could be a problem for StratAer usage !'
+       WRITE(lunout,*) 'NBTR_BIN=',nbtr_bin
     ELSE
-       DO it=2, nbtr_bin
-          mdw(it)=mdw(it-1)*V_rat**(1./3.)
-       ENDDO
-    ENDIF
-    IF (is_master) WRITE(lunout,*) 'init mdw=', mdw
+       mdw(1)=mdwmin
+       IF (V_rat.LT.1.62) THEN ! compensate for dip in second bin for lower volume ratio
+          mdw(2)=mdw(1)*2.**(1./3.)
+          DO it=3, nbtr_bin
+             mdw(it)=mdw(it-1)*V_rat**(1./3.)
+          ENDDO
+       ELSE
+          DO it=2, nbtr_bin
+             mdw(it)=mdw(it-1)*V_rat**(1./3.)
+          ENDDO
+       ENDIF
+       IF (is_master) WRITE(lunout,*) 'init mdw=', mdw
+    ENDIF
     
     !   compute particle radius RRSI [cm] and volume Vbin [m3] from diameter mdw [m]