Changeset 4131 for trunk/LMDZ.PLUTO/libf/muphypluto/mp2m_clouds.F90

Timestamp:

Mar 16, 2026, 3:12:16 PM (6 days ago)

Author:

debatzbr

Message:

Pluto PCM: Fixes for debug mode.
BBT

File:

• trunk/LMDZ.PLUTO/libf/muphypluto/mp2m_clouds.F90 (modified) (5 diffs)

trunk/LMDZ.PLUTO/libf/muphypluto/mp2m_clouds.F90

@@ -581 +581 @@
         ! Equilibrium saturation near the drop:
         !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-        S_eq = exp(min((2._mm_wp*sig*xESP%masmol) / (mm_rgas*mm_temp*xESP%rho_s*drad),30._mm_wp))
+        where(drad .gt. mm_drad_min)
+          S_eq = exp(min((2._mm_wp*sig*xESP%masmol) / (mm_rgas*mm_temp*xESP%rho_s*drad),30._mm_wp))
+        elsewhere
+          S_eq = exp(min((2._mm_wp*sig*xESP%masmol) / (mm_rgas*mm_temp*xESP%rho_s*mm_drad_min),30._mm_wp))
+        endwhere
 
         ! Gets growth rate:
@@ -765 +769 @@
 
         ! The knudsen number of the drop:
-        !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-        knu(:) = L(:) / drad(:)
+        !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~        
+        where(drad .gt. mm_drad_min)
+          knu = L / drad
+        elsewhere
+          knu = L / mm_drad_min
+        endwhere
 
         ! Slip flow correction:
@@ -1095 +1103 @@
 
         nz = SIZE(rad)
+        alpha = 0._mm_wp
 
         ! Linear extrapolation of velocity:
@@ -1132 +1141 @@
           wi = wsettle(mm_play(idx),mm_temp(idx),zi,rhoair_mid,rho(idx),rad(idx))
           ! ws & wi must be different !
-          IF(dabs(wi-ws)/wi <= 1.e-3_mm_wp)  wi=ws/1.001_mm_wp
+          IF(wi/=0._mm_wp.AND.dabs(wi-ws)/wi <= 1.e-3_mm_wp)  wi=ws/1.001_mm_wp
           IF (wi/=0._mm_wp.AND.(ws/wi)>0._mm_wp.AND.(zs-zi)/=0._mm_wp) alpha = dlog(ws/wi)/(zs-zi) ! Alpha < 0 if wi > ws
           !   -es < argexp < es
@@ -1160 +1169 @@
             ENDIF
           ENDIF
-          zni = mm_zlev(idx)-dlog(arg1)/alpha
-          zns = mm_zlev(idx)-mm_dzlev(idx)-dlog(arg2)/alpha
+          IF (alpha /= 0._mm_wp) THEN
+            zni = mm_zlev(idx)-dlog(arg1)/alpha
+            zns = mm_zlev(idx)-mm_dzlev(idx)-dlog(arg2)/alpha
+          ELSE
+            zni = mm_zlev(idx)
+            zns = mm_zlev(idx)-mm_dzlev(idx)
+          ENDIF
           RETURN
         ENDIF ! end of ichx