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

Timestamp:

Jun 6, 2025, 11:55:21 AM (3 days ago)

Author:

aborella

Message:

Adaptations made to contrails radiative transfer and to ice sedimentation, following discussion with Bernd Karcher

File:

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

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

@@ -228 +228 @@
   !$OMP THREADPRIVATE(ok_plane_contrail)
 
-  LOGICAL, SAVE, PROTECTED :: ok_precip_lincontrails=.TRUE.  ! if True, linear contrails can precipitate
-  !$OMP THREADPRIVATE(ok_precip_lincontrails)
+  LOGICAL, SAVE, PROTECTED :: ok_precip_contrails=.TRUE.     ! if True, contrails can be autoconverted to snow
+  !$OMP THREADPRIVATE(ok_precip_contrails)
 
   REAL, SAVE, PROTECTED :: aspect_ratio_lincontrails=.1      ! [-] aspect ratio of linear contrails
@@ -260 +260 @@
   REAL, SAVE, PROTECTED :: initial_height_contrails=200.     ! [m] initial height of the linear contrails formed
   !$OMP THREADPRIVATE(initial_height_contrails)
+
+  REAL, SAVE, PROTECTED :: fallice_linear_contrails=1.       ! [m/s] Ice fallspeed velocity in linear contrails
+  !$OMP THREADPRIVATE(fallice_linear_contrails)
+
+  REAL, SAVE, PROTECTED :: fallice_cirrus_contrails=1.       ! [m/s] Ice fallspeed velocity in cirrus contrails
+  !$OMP THREADPRIVATE(fallice_cirrus_contrails)
 
   REAL, SAVE, PROTECTED :: aviation_coef=1.                  ! [-] scaling factor for aviation emissions and flown distance
@@ -512 +518 @@
     CALL getin_p('coef_shear_lscp',coef_shear_lscp)
     ! for aviation
-    CALL getin_p('ok_precip_lincontrails',ok_precip_lincontrails)
+    CALL getin_p('ok_precip_contrails',ok_precip_contrails)
     CALL getin_p('aspect_ratio_lincontrails',aspect_ratio_lincontrails)
     coef_mixing_lincontrails=coef_mixing_lscp
@@ -525 +531 @@
     CALL getin_p('initial_width_contrails',initial_width_contrails)
     CALL getin_p('initial_height_contrails',initial_height_contrails)
+    CALL getin_p('fallice_linear_contrails',fallice_linear_contrails)
+    CALL getin_p('fallice_cirrus_contrails',fallice_cirrus_contrails)
     CALL getin_p('aviation_coef',aviation_coef)
 
@@ -615 +623 @@
     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, ok_precip_contrails:', ok_precip_contrails
     WRITE(lunout,*) 'lscp_ini, aspect_ratio_lincontrails:', aspect_ratio_lincontrails
     WRITE(lunout,*) 'lscp_ini, coef_mixing_lincontrails:', coef_mixing_lincontrails
@@ -626 +634 @@
     WRITE(lunout,*) 'lscp_ini, initial_width_contrails:', initial_width_contrails
     WRITE(lunout,*) 'lscp_ini, initial_height_contrails:', initial_height_contrails
+    WRITE(lunout,*) 'lscp_ini, fallice_linear_contrails:', fallice_linear_contrails
+    WRITE(lunout,*) 'lscp_ini, fallice_cirrus_contrails:', fallice_cirrus_contrails
     WRITE(lunout,*) 'lscp_ini, aviation_coef:', aviation_coef