Changeset 5821 for LMDZ6/trunk/libf/phylmd/ecrad/lmdz

Timestamp:

Sep 18, 2025, 10:29:26 AM (4 months ago)

Author:

idelkadi

Message:

Added a formulation to prescribe effective cloud size as a hyperbolic tangent function of pressure for calculating radiative fluxes related to 3D cloud effects.
Activation is controlled in namelist_ecrad by the logical key ok_separation_tanh

Location:

LMDZ6/trunk/libf/phylmd/ecrad/lmdz

Files:

• radiation_scheme_mod.f90 (modified) (2 diffs)
• setup_config_from_lmdz.f90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/ecrad/lmdz/radiation_scheme_mod.f90

@@ -428 +428 @@
                &  driver_config%high_inv_effective_size, 0.8_jprb, 0.45_jprb, &
                &  KIDIA, KFDIA)
-   else if (driver_config%ok_separation) then
+   else if (driver_config%ok_separation_eta) then
      call cloud%param_cloud_effective_separation_eta(klon, klev, &
+               &  thermodynamics%pressure_hl, &
+               &  driver_config%cloud_separation_scale_surface, &
+               &  driver_config%cloud_separation_scale_toa, &
+               &  driver_config%cloud_separation_scale_power, &
+               &  driver_config%cloud_inhom_separation_factor, &
+               &  KIDIA, KFDIA)
+   else if (driver_config%ok_separation_tanh) then
+     call cloud%param_cloud_effective_separation_tanh(klon, klev, &
                &  thermodynamics%pressure_hl, &
                &  driver_config%cloud_separation_scale_surface, &
@@ -1134 +1142 @@
                &  driver_config%high_inv_effective_size, 0.8_jprb, 0.45_jprb, &
                &  KIDIA, KFDIA)
-   else if (driver_config%ok_separation) then
+   else if (driver_config%ok_separation_eta) then
      call cloud%param_cloud_effective_separation_eta(klon, klev, &
+               &  thermodynamics%pressure_hl, &
+               &  driver_config%cloud_separation_scale_surface, &
+               &  driver_config%cloud_separation_scale_toa, &
+               &  driver_config%cloud_separation_scale_power, &
+               &  driver_config%cloud_inhom_separation_factor, &
+               &  KIDIA, KFDIA)
+   else if (driver_config%ok_separation_tanh) then
+     call cloud%param_cloud_effective_separation_tanh(klon, klev, &
                &  thermodynamics%pressure_hl, &
                &  driver_config%cloud_separation_scale_surface, &

LMDZ6/trunk/libf/phylmd/ecrad/lmdz/setup_config_from_lmdz.f90

@@ -9 +9 @@
   type driver_config_type
      logical    :: ok_effective_size  = .true.
-     logical    :: ok_separation = .false.   
+     logical    :: ok_separation_eta = .false.   
+     logical    :: ok_separation_tanh = .false.
      real(jprb) :: high_inv_effective_size   = -1.0_jprb ! m-1
      real(jprb) :: middle_inv_effective_size = -1.0_jprb ! m-1
@@ -44 +45 @@
     character(*), intent(in)          :: file_name
     logical, intent(out), optional    :: is_success
-    logical    :: ok_effective_size, ok_separation
+    logical    :: ok_effective_size, ok_separation_eta, ok_separation_tanh
     integer    :: iosopen ! Status after calling open
     real(jprb) :: high_inv_effective_size 
@@ -61 +62 @@
     logical :: do_save_inputs
 
-    namelist /radiation_driver/ ok_effective_size, ok_separation, &
+    namelist /radiation_driver/ ok_effective_size, ok_separation_eta, ok_separation_tanh, &
          &  frac_std, overlap_decorr_length, kdecolat, &
          &  low_decorrelation_length, mid_decorrelation_length, high_decorrelation_length, &
@@ -70 +71 @@
 
     ok_effective_size = .false.
-    ok_separation = .false.
+    ok_separation_eta = .false.
+    ok_separation_tanh = .false.
     high_inv_effective_size   = -1.0_jprb 
     middle_inv_effective_size = -1.0_jprb 
@@ -109 +111 @@
     ! Copy namelist data into configuration object
     this%ok_effective_size = ok_effective_size 
-    this%ok_separation = ok_separation
+    this%ok_separation_eta = ok_separation_eta
+    this%ok_separation_tanh = ok_separation_tanh
     this%frac_std = frac_std
     this%overlap_decorr_length = overlap_decorr_length