Changeset 5239 for LMDZ6/trunk/libf

Timestamp:

Oct 17, 2024, 11:20:55 AM (6 weeks ago)

Author:

lebasn

Message:

StratAer?: bugfix adding sanity check to avoid a floating point exception in debug mode.

File:

• LMDZ6/trunk/libf/phylmd/StratAer/traccoag_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/StratAer/traccoag_mod.F90

@@ -406 +406 @@
            ENDDO
 !     reff_sulfate: effective radius of WET sulfate aerosols (cm)
-           reff_sulfate(i,j) = (samoment3 / samoment2) &
-                & *1.e2                                              ! conversion from m to cm
+           IF(samoment2 > 1.e-30) THEN
+              reff_sulfate(i,j) = (samoment3 / samoment2) &
+                   & *1.e2                                              ! conversion from m to cm
+              
+              ! Sanity check
+              IF(reff_sulfate(i,j) > 5.e-4) reff_sulfate(i,j) = 5.e-4   ! reff_sulfate max = 5 micron
+              IF(reff_sulfate(i,j) < 1.e-6) reff_sulfate(i,j) = 1.e-6   ! reff_sulfate min = 10 nm
+           ELSE
+              reff_sulfate(i,j) = 1.e-5                                 ! reff_sulfate N part = nul (ref = 100 nm)
+           ENDIF
         ENDDO
      ENDDO