Changeset 3687 for trunk/LMDZ.PLUTO

Timestamp:

Mar 17, 2025, 10:00:19 AM (5 months ago)

Author:

debatzbr

Message:

Initialize dtlw_hcn_c2h2 if cooling = false
BBT

File:

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

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

@@ -450 +450 @@
 !--------------------------------------------------
 !     Effective radius and variance of the aerosols
+!     Madeleine's PhD (eq. 2.3-2.4):
+!     --> r_eff = <r^3> / <r^2>
+!     --> nu_eff = <r^4>*<r^2> / <r^3>^2 - 1
 !--------------------------------------------------
       ! Radiative Hazes
@@ -456 +459 @@
             ! Spherical aerosols
             sig = 0.2
-            where (mp2m_rc_sph(:,:) > 1e-10)
+            where (mp2m_rc_sph(:,:) > 1e-9)
                reffrad(:,:,1) = mp2m_rc_sph(:,:) * exp(5.*sig**2 / 2.)
             elsewhere
@@ -468 +471 @@
             ! Fractal aerosols
             sig = 0.35
-            where (mp2m_rc_fra(:,:) > 1e-10)
+            where (mp2m_rc_fra(:,:) > 1e-8)
                reffrad(:,:,2) = mp2m_rc_fra(:,:) * exp(5.*sig**2 / 2.)
                elsewhere
@@ -543 +546 @@
 !     Net atmospheric radiative cooling rate from C2H2 (K.s-1):
 !     ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-    !   dtlw_hcn_c2h2=0.
+      dtlw_hcn_c2h2=0.
       if (cooling) then
          CALL cooling_hcn_c2h2(ngrid,nlayer,pplay,&
@@ -1103 +1106 @@
         end do
         ! adding c2h2 if cooling active
-        dtlw(ig,l)=dtlw(ig,l)+dtlw_hcn_c2h2(ig,l)
+        if (cooling) then
+         dtlw(ig,l)=dtlw(ig,l)+dtlw_hcn_c2h2(ig,l)
+        endif
 
         !SW