Changeset 5790 for LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_tools.f90

Timestamp:

Jul 28, 2025, 6:44:28 PM (8 days ago)

Author:

aborella

Message:

Major modifs to treatment of contrails (from 2 classes to 2 moments) + diagnostics. Increased numerical efficiency

File:

• LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_tools.f90 (modified) (1 diff)

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

@@ -108 +108 @@
 
 END SUBROUTINE FALLICE_VELOCITY
+!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+
+!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+FUNCTION ICECRYSTAL_VELO(mice, temp, pres)
+
+    ! Ref:
+    ! Spichtinger and Gierens (2009)
+
+    USE lmdz_lscp_ini, ONLY: RPI, rho_ice
+
+    IMPLICIT NONE
+
+    REAL :: mice      ! ice crystal mass [kg]
+    REAL :: temp      ! temperature [K]
+    REAL :: pres      ! air pressure [Pa]
+    REAL :: icecrystal_velo    ! fallspeed velocity of crystals [m/s]
+
+    !--Local
+    REAL :: c, g, d
+
+    c = (pres / 30000)**(-0.178) * (temp / 233.)**(-0.394)
+    IF ( mice .LT. 2.146e-13 ) THEN
+        g = 735.4
+        d = 0.42
+    ELSEIF ( mice .LT. 2.166e-9) THEN
+        g = 63292.4
+        d = 0.57
+    ELSEIF ( mice .LT. 4.264e-8 ) THEN
+        g = 329.8
+        d = 0.31
+    ELSE
+        g = 8.8
+        d = 0.096
+    ENDIF
+    icecrystal_velo = g * mice**d * c
+
+END FUNCTION ICECRYSTAL_VELO
 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++