Changeset 3130 for palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

Timestamp:

Jul 16, 2018 11:08:55 AM (7 years ago)

Author:

gronemeier

Message:

calculate km according to MOST within the surface layer

File:

• palm/trunk/SOURCE/surface_layer_fluxes_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! move phi_m from turbulence_closure_mod
+! 
+! 3045 2018-05-28 07:55:41Z Giersch
 ! Error message revised
 ! 
@@ -287 +290 @@
     PRIVATE
 
-    PUBLIC init_surface_layer_fluxes, surface_layer_fluxes
+    PUBLIC init_surface_layer_fluxes, phi_m, surface_layer_fluxes
 
     INTERFACE init_surface_layer_fluxes
        MODULE PROCEDURE init_surface_layer_fluxes
     END INTERFACE init_surface_layer_fluxes
+
+    INTERFACE phi_m
+       MODULE PROCEDURE phi_m
+    END INTERFACE phi_m
 
     INTERFACE surface_layer_fluxes
@@ -2293 +2300 @@
     END FUNCTION psi_h
 
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Calculates stability function for momentum
+!>
+!> @author Hauke Wurps
+!------------------------------------------------------------------------------!
+    FUNCTION phi_m( zeta ) 
+   
+       IMPLICIT NONE 
+   
+       REAL(wp)            :: phi_m         !< Value of the function
+       REAL(wp)            :: zeta          !< Stability parameter z/L
+   
+       REAL(wp), PARAMETER :: a = 16.0_wp   !< constant
+       REAL(wp), PARAMETER :: c = 5.0_wp    !< constant
+   
+       IF ( zeta < 0.0_wp )  THEN
+          phi_m = 1.0_wp / SQRT( SQRT( 1.0_wp - a * zeta ) )
+       ELSE
+          phi_m = 1.0_wp + c * zeta
+       ENDIF
+   
+    END FUNCTION phi_m
+
  END MODULE surface_layer_fluxes_mod