Ignore:
Timestamp:
Jul 16, 2018 11:08:55 AM (6 years ago)
Author:
gronemeier
Message:

calculate km according to MOST within the surface layer

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

    r3045 r3130  
    2626! -----------------
    2727! $Id$
     28! move phi_m from turbulence_closure_mod
     29!
     30! 3045 2018-05-28 07:55:41Z Giersch
    2831! Error message revised
    2932!
     
    287290    PRIVATE
    288291
    289     PUBLIC init_surface_layer_fluxes, surface_layer_fluxes
     292    PUBLIC init_surface_layer_fluxes, phi_m, surface_layer_fluxes
    290293
    291294    INTERFACE init_surface_layer_fluxes
    292295       MODULE PROCEDURE init_surface_layer_fluxes
    293296    END INTERFACE init_surface_layer_fluxes
     297
     298    INTERFACE phi_m
     299       MODULE PROCEDURE phi_m
     300    END INTERFACE phi_m
    294301
    295302    INTERFACE surface_layer_fluxes
     
    22932300    END FUNCTION psi_h
    22942301
     2302
     2303!------------------------------------------------------------------------------!
     2304! Description:
     2305! ------------
     2306!> Calculates stability function for momentum
     2307!>
     2308!> @author Hauke Wurps
     2309!------------------------------------------------------------------------------!
     2310    FUNCTION phi_m( zeta )
     2311   
     2312       IMPLICIT NONE
     2313   
     2314       REAL(wp)            :: phi_m         !< Value of the function
     2315       REAL(wp)            :: zeta          !< Stability parameter z/L
     2316   
     2317       REAL(wp), PARAMETER :: a = 16.0_wp   !< constant
     2318       REAL(wp), PARAMETER :: c = 5.0_wp    !< constant
     2319   
     2320       IF ( zeta < 0.0_wp )  THEN
     2321          phi_m = 1.0_wp / SQRT( SQRT( 1.0_wp - a * zeta ) )
     2322       ELSE
     2323          phi_m = 1.0_wp + c * zeta
     2324       ENDIF
     2325   
     2326    END FUNCTION phi_m
     2327
    22952328 END MODULE surface_layer_fluxes_mod
Note: See TracChangeset for help on using the changeset viewer.