Changeset 3805 for trunk

Timestamp:

Jun 13, 2025, 5:33:46 PM (4 days ago)

Author:

debatzbr

Message:

Pluto PCM: Fix reffrad and nueffrad calculation
BBT

File:

• trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90 (modified) (2 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90

@@ -459 +459 @@
             ! Spherical aerosols
             sig = 0.2
-            where (mp2m_rc_sph(:,:) > 1e-9)
+            where (mp2m_rc_sph(:,:) >= 5e-9)
                reffrad(:,:,1) = mp2m_rc_sph(:,:) * exp(5.*sig**2 / 2.)
             elsewhere
                reffrad(:,:,1) = 0d0
             endwhere
-            if (exp(sig**2) - 1 > 0.1) then
-               nueffrad(:,:,1) = exp(sig**2) - 1
-            else
-               nueffrad(:,:,1) = 0.1
-            endif
+            nueffrad(:,:,1) = exp(sig**2) - 1
+            !if (exp(sig**2) - 1 > 0.1) then
+            !   nueffrad(:,:,1) = exp(sig**2) - 1
+            !else
+            !   nueffrad(:,:,1) = 0.1
+            !endif
             ! Fractal aerosols
             sig = 0.35
-            where (mp2m_rc_fra(:,:) > 1e-8)
+            where (mp2m_rc_fra(:,:) >= 1e-8)
                reffrad(:,:,2) = mp2m_rc_fra(:,:) * exp(5.*sig**2 / 2.)
-               elsewhere
+            elsewhere
                reffrad(:,:,2) = 0d0
-               endwhere
-            if (exp(sig**2) - 1 > 0.1) then
-               nueffrad(:,:,2) = exp(sig**2) - 1
-            else
-               nueffrad(:,:,2) = 0.1
-            endif
+            endwhere
+            nueffrad(:,:,2) = exp(sig**2) - 1
+            !if (exp(sig**2) - 1 > 0.1) then
+            !   nueffrad(:,:,2) = exp(sig**2) - 1
+            !else
+            !   nueffrad(:,:,2) = 0.1
+            !endif
 
          else
@@ -511 +513 @@
             QREFvis3d,QREFir3d)
 
-         ! Get aerosol optical depths.
+         ! Get aerosol optical depths dtau_aer and diagnostic tau_col.
          call aeropacity(ngrid,nlayer,nq,pplay,pplev,zzlev,pt,pq,dtau_aer,      &
             reffrad,nueffrad,QREFvis3d,QREFir3d,                             &