Changeset 4911 for LMDZ6/trunk/libf/phylmd/ecrad/lmdz/setup_config_from_lmdz.F90

Timestamp:

Apr 19, 2024, 11:12:57 AM (4 weeks ago)

Author:

idelkadi

Message:

Introduction of 2 routines for LMDZ-ECRAD :

• One for LMDZ (calcul_cloud_overlap_decorr_len.F90) to calculate the decorrelation length (Ld) in the case of "Exp-Ran" cloud overlap.
  1. Ld=constant
  2. Ld=linear function of the absolute value of latitude (Shonk and Hogan 2010)
  3. Ld=function of vertical level
• One for Ecrad (set_overlap_param_var2D in radiation_cloud.F90) : Compute and store the overlap parameter from the provided overlap decorrelation length, which depends on vertical level.

Translated with DeepL.com (free version)

File:

• LMDZ6/trunk/libf/phylmd/ecrad/lmdz/setup_config_from_lmdz.F90 (modified) (4 diffs)

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

@@ -19 +19 @@
      real(jprb) :: frac_std = 0.75_jprb
      real(jprb) :: overlap_decorr_length = 2000.0_jprb
-
+     ! KDECOLAT : 0 > cste, 1 > Shonk et al. (2010) 2 > 1 lisse a l'equateur
+     !            3 : 1+variation sur la verticale
+     !            4 : 2+variation sur la verticale
+     integer    :: kdecolat = 0
+     real(jprb) :: low_decorrelation_length = 2000.0_jprb 
+     real(jprb) :: mid_decorrelation_length = 2000.0_jprb
+     real(jprb) :: high_decorrelation_length = 2000.0_jprb
 
  contains
@@ -47 +53 @@
     real(jprb) :: frac_std
     real(jprb) :: overlap_decorr_length
+    integer    :: kdecolat 
+    real(jprb) :: low_decorrelation_length
+    real(jprb) :: mid_decorrelation_length 
+    real(jprb) :: high_decorrelation_length 
 
     namelist /radiation_driver/ ok_effective_size, ok_separation, &
-         &  frac_std, overlap_decorr_length, &
+         &  frac_std, overlap_decorr_length, kdecolat, &
+         &  low_decorrelation_length, mid_decorrelation_length, high_decorrelation_length, &
          &  high_inv_effective_size, middle_inv_effective_size, low_inv_effective_size, &
          &  cloud_inhom_separation_factor, cloud_separation_scale_surface, &
@@ -65 +76 @@
     frac_std = 0.75_jprb
     overlap_decorr_length = 2000.0_jprb
+    kdecolat = 0
+    low_decorrelation_length = 2000.0_jprb
+    mid_decorrelation_length = 2000.0_jprb
+    high_decorrelation_length = 2000.0_jprb
 
     ! Open the namelist file and read the radiation_driver namelist
@@ -92 +107 @@
     this%frac_std = frac_std
     this%overlap_decorr_length = overlap_decorr_length
+    this%kdecolat = kdecolat
+    this%low_decorrelation_length = low_decorrelation_length
+    this%mid_decorrelation_length = mid_decorrelation_length
+    this%high_decorrelation_length = high_decorrelation_length
     this%cloud_inhom_separation_factor = cloud_inhom_separation_factor
     this%cloud_separation_scale_surface = cloud_separation_scale_surface