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

Timestamp:

Mar 19, 2025, 3:53:17 PM (3 months ago)

Author:

aborella

Message:

Added resuspension of ice crystals in poprecip + contrails do not precipitate anymore + added different options to control the behavior of contrails

File:

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

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

@@ -245 +245 @@
   REAL, SAVE, PROTECTED :: initial_width_contrails=200.      ! [m] initial width of the linear contrails formed
   !$OMP THREADPRIVATE(initial_width_contrails)
+
+  REAL, SAVE, PROTECTED :: initial_height_contrails=200.     ! [m] initial height of the linear contrails formed
+  !$OMP THREADPRIVATE(initial_height_contrails)
+
+  REAL, SAVE, PROTECTED :: aviation_coef=1.                  ! [-] scaling factor for aviation emissions and flown distance
+  !$OMP THREADPRIVATE(aviation_coef)
   !--End of the parameters for aviation
 
@@ -337 +343 @@
   REAL, SAVE, PROTECTED :: snow_fallspeed_cld               ! Snow fall velocity in cloudy sky [m/s]
   !$OMP THREADPRIVATE(snow_fallspeed_cld)
+
+  INTEGER, SAVE, PROTECTED :: ok_precip_resuspension=.FALSE.! Flag to activate the resuspension of ice crystals
+  !$OMP THREADPRIVATE(ok_precip_resuspension)
+
+  REAL, SAVE, PROTECTED :: snow_thresh_resuspension=1.e-5   ! [kg/m2/s] Threshold flux below which part of the snow flux will be resuspended
+  !$OMP THREADPRIVATE(snow_thresh_resuspension)
   !--End of the parameters for poprecip
 
@@ -455 +467 @@
     CALL getin_p('snow_fallspeed_clr',snow_fallspeed_clr)
     CALL getin_p('snow_fallspeed_cld',snow_fallspeed_cld)
+    CALL getin_p('ok_precip_resuspension',ok_precip_resuspension)
+    CALL getin_p('snow_thresh_resuspension',snow_thresh_resuspension)
     ! for condensation and ice supersaturation
     CALL getin_p('ok_unadjusted_clouds',ok_unadjusted_clouds)
@@ -488 +502 @@
     CALL getin_p('linear_contrails_lifetime',linear_contrails_lifetime)
     CALL getin_p('initial_width_contrails',initial_width_contrails)
+    CALL getin_p('initial_height_contrails',initial_height_contrails)
+    CALL getin_p('aviation_coef',aviation_coef)
 
 
@@ -550 +566 @@
     WRITE(lunout,*) 'lscp_ini, snow_fallspeed_clr:', snow_fallspeed_clr
     WRITE(lunout,*) 'lscp_ini, snow_fallspeed_cld:', snow_fallspeed_cld
+    WRITE(lunout,*) 'lscp_ini, ok_precip_resuspension:', ok_precip_resuspension
+    WRITE(lunout,*) 'lscp_ini, snow_thresh_resuspension:', snow_thresh_resuspension
     ! for condensation and ice supersaturation
     WRITE(lunout,*) 'lscp_ini, ok_ice_supersat:', ok_ice_supersat
@@ -581 +599 @@
     WRITE(lunout,*) 'lscp_ini, linear_contrails_lifetime:', linear_contrails_lifetime
     WRITE(lunout,*) 'lscp_ini, initial_width_contrails:', initial_width_contrails
+    WRITE(lunout,*) 'lscp_ini, initial_height_contrails:', initial_height_contrails
+    WRITE(lunout,*) 'lscp_ini, aviation_coef:', aviation_coef