Changeset 4502 for palm/trunk/SOURCE/basic_constants_and_equations_mod.f90

Timestamp:

Apr 17, 2020 4:14:16 PM (4 years ago)

Author:

schwenkel

Message:

Implementation of ice microphysics

File:

• palm/trunk/SOURCE/basic_constants_and_equations_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/basic_constants_and_equations_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implementation of ice microphysics
+! 
+! 4400 2020-02-10 20:32:41Z suehring
 ! Move routine to transform coordinates from netcdf_interface_mod to 
 ! basic_constants_and_equations_mod
@@ -92 +95 @@
     REAL(wp), PARAMETER ::  g_d_cp  = g   / c_p   !< precomputed g / c_p
     REAL(wp), PARAMETER ::  lv_d_cp = l_v / c_p   !< precomputed l_v / c_p
+    REAL(wp), PARAMETER ::  ls_d_cp = l_s / c_p   !< precomputed l_s / c_p
     REAL(wp), PARAMETER ::  lv_d_rd = l_v / r_d   !< precomputed l_v / r_d
     REAL(wp), PARAMETER ::  rd_d_rv = r_d / r_v   !< precomputed r_d / r_v
@@ -106 +110 @@
     PRIVATE magnus_0d, &
             magnus_1d, &
+            magnus_tl_0d, &
+            magnus_tl_1d, &
+            magnus_0d_ice, &
+            magnus_1d_ice, &
             ideal_gas_law_rho_0d, &
             ideal_gas_law_rho_1d, &
@@ -126 +134 @@
        MODULE PROCEDURE magnus_1d
     END INTERFACE magnus
+
+    INTERFACE magnus_tl
+       MODULE PROCEDURE magnus_tl_0d
+       MODULE PROCEDURE magnus_tl_1d
+    END INTERFACE magnus_tl
+
+    INTERFACE magnus_ice
+       MODULE PROCEDURE magnus_0d_ice
+       MODULE PROCEDURE magnus_1d_ice
+    END INTERFACE magnus_ice
 
     INTERFACE ideal_gas_law_rho
@@ -337 +355 @@
 ! Description:
 ! ------------
+!> This function computes the magnus formula (Press et al., 1992) using the
+!> (ice-) liquid water potential temperature.
+!> The magnus formula is needed to calculate the saturation vapor pressure over
+!> a plane liquid water surface
+!------------------------------------------------------------------------------!
+    FUNCTION magnus_tl_0d( t_l )
+
+       IMPLICIT NONE
+
+       REAL(wp), INTENT(IN) ::  t_l  !< liquid water temperature (K)
+
+       REAL(wp) ::  magnus_tl_0d
+!
+!--    Saturation vapor pressure for a specific temperature:
+       magnus_tl_0d =  610.78_wp * EXP( 17.269_wp * ( t_l - 273.16_wp ) /             &
+                                                    ( t_l - 35.86_wp  ) )
+
+    END FUNCTION magnus_tl_0d
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> This function computes the magnus formula (Press et al., 1992) using the
+!> (ice-) liquid water potential temperature.
+!> The magnus formula is needed to calculate the saturation vapor pressure over
+!> a plane liquid water surface
+!------------------------------------------------------------------------------!
+    FUNCTION magnus_tl_1d( t_l )
+
+       IMPLICIT NONE
+
+       REAL(wp), INTENT(IN), DIMENSION(:) ::  t_l  !< liquid water temperature (K)
+
+       REAL(wp), DIMENSION(size(t_l)) ::  magnus_tl_1d
+!
+!--    Saturation vapor pressure for a specific temperature:
+       magnus_tl_1d =  610.78_wp * EXP( 17.269_wp * ( t_l - 273.16_wp ) /      &
+                                                    ( t_l - 35.86_wp  ) )
+
+    END FUNCTION magnus_tl_1d
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> This function computes the magnus formula (Press et al., 1992).
+!> The magnus formula is needed to calculate the saturation vapor pressure over
+!> a plane ice surface
+!------------------------------------------------------------------------------!
+    FUNCTION magnus_0d_ice( t )
+
+       IMPLICIT NONE
+
+       REAL(wp), INTENT(IN) ::  t  !< temperature (K)
+
+       REAL(wp) ::  magnus_0d_ice
+!
+!--    Saturation vapor pressure for a specific temperature:
+       magnus_0d_ice =  611.2_wp * EXP( 22.46_wp * ( t - degc_to_k ) /         &
+                                                   ( t - 0.53_wp  ) )
+
+    END FUNCTION magnus_0d_ice
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> This function computes the magnus formula (Press et al., 1992).
+!> The magnus formula is needed to calculate the saturation vapor pressure over
+!> a plane ice surface
+!------------------------------------------------------------------------------!
+    FUNCTION magnus_1d_ice( t )
+
+       IMPLICIT NONE
+
+       REAL(wp), INTENT(IN), DIMENSION(:) ::  t  !< temperature (K)
+
+       REAL(wp), DIMENSION(size(t)) ::  magnus_1d_ice
+!
+!--    Saturation vapor pressure for a specific temperature:
+       magnus_1d_ice =  611.2_wp * EXP( 22.46_wp * ( t - degc_to_k ) /          &
+                                                   ( t - 0.53_wp  ) )
+
+    END FUNCTION magnus_1d_ice
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
 !> Compute the ideal gas law for scalar arguments.
 !------------------------------------------------------------------------------!