Changeset 5643 for LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_ini.f90

Timestamp:

May 5, 2025, 3:02:32 PM (5 weeks ago)

Author:

aborella

Message:

Various minor corrections to ISSR and contrails param

File:

• LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_ini.f90 (modified) (6 diffs)

LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_ini.f90

@@ -162 +162 @@
   !$OMP THREADPRIVATE(ok_weibull_warm_clouds)
 
+  LOGICAL, SAVE, PROTECTED :: ok_nodeep_lscp=.FALSE.         ! if True, the deep convection clouds are removed from the lscp calculations
+  !$OMP THREADPRIVATE(ok_nodeep_lscp)
+
+  LOGICAL, SAVE, PROTECTED :: ok_nodeep_lscp_rad=.FALSE.     ! if True, the deep convection clouds are not accounted two times in radiative transfer
+  !$OMP THREADPRIVATE(ok_nodeep_lscp_rad)
+
   REAL, SAVE, PROTECTED :: ffallv_issr                       ! tuning coefficient crystal fall velocity, cirrus clouds (with ISSR)
   !$OMP THREADPRIVATE(ffallv_issr)
@@ -221 +227 @@
   LOGICAL, SAVE, PROTECTED :: ok_plane_contrail              ! activates the contrails parameterisation
   !$OMP THREADPRIVATE(ok_plane_contrail)
+
+  LOGICAL, SAVE, PROTECTED :: ok_precip_lincontrails=.TRUE.  ! if True, linear contrails can precipitate
+  !$OMP THREADPRIVATE(ok_precip_lincontrails)
 
   REAL, SAVE, PROTECTED :: aspect_ratio_lincontrails=.1      ! [-] aspect ratio of linear contrails
@@ -482 +491 @@
     CALL getin_p('ok_unadjusted_clouds',ok_unadjusted_clouds)
     CALL getin_p('ok_weibull_warm_clouds',ok_weibull_warm_clouds)
+    CALL getin_p('ok_nodeep_lscp',ok_nodeep_lscp)
+    CALL getin_p('ok_nodeep_lscp_rad',ok_nodeep_lscp_rad)
     ffallv_issr=ffallv_lsc
     CALL getin_p('ffallv_issr',ffallv_issr)
@@ -501 +512 @@
     CALL getin_p('coef_shear_lscp',coef_shear_lscp)
     ! for aviation
+    CALL getin_p('ok_precip_lincontrails',ok_precip_lincontrails)
     CALL getin_p('aspect_ratio_lincontrails',aspect_ratio_lincontrails)
     coef_mixing_lincontrails=coef_mixing_lscp
@@ -583 +595 @@
     WRITE(lunout,*) 'lscp_ini, ok_unadjusted_clouds:', ok_unadjusted_clouds
     WRITE(lunout,*) 'lscp_ini, ok_weibull_warm_clouds:', ok_weibull_warm_clouds
+    WRITE(lunout,*) 'lscp_ini, ok_nodeep_lscp:', ok_nodeep_lscp
+    WRITE(lunout,*) 'lscp_ini, ok_nodeep_lscp_rad:', ok_nodeep_lscp_rad
     WRITE(lunout,*) 'lscp_ini, ffallv_issr', ffallv_issr
     WRITE(lunout,*) 'lscp_ini, cooling_rate_ice_thresh', cooling_rate_ice_thresh
@@ -601 +615 @@
     WRITE(lunout,*) 'lscp_ini, coef_shear_lscp:', coef_shear_lscp
     ! for aviation
+    WRITE(lunout,*) 'lscp_ini, ok_precip_lincontrails:', ok_precip_lincontrails
     WRITE(lunout,*) 'lscp_ini, aspect_ratio_lincontrails:', aspect_ratio_lincontrails
     WRITE(lunout,*) 'lscp_ini, coef_mixing_lincontrails:', coef_mixing_lincontrails