source: flexpart.git/src/psim.f90

10.4.1_pesei
Last change on this file was 92fab65, checked in by Ignacio Pisso <ip@…>, 9 months ago

add SPDX-License-Identifier to all .f90 files

  • Property mode set to 100644
File size: 948 bytes
Line 
1! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
2! SPDX-License-Identifier: GPL-3.0-or-later
3
4real function psim(z,al)
5
6  !**********************************************************************
7  !                                                                     *
8  ! DESCRIPTION: CALCULATION OF THE STABILITY CORRECTION FUNCTION FOR   *
9  !              MOMENTUM AS FUNCTION OF HEIGHT Z AND OBUKHOV SCALE     *
10  !              HEIGHT L                                               *
11  !                                                                     *
12  !**********************************************************************
13
14  use par_mod
15
16  implicit none
17
18  real :: z,al,zeta,x,a1,a2
19
20  zeta=z/al
21  if(zeta.le.0.) then
22  ! UNSTABLE CASE
23    x=(1.-15.*zeta)**0.25
24    a1=((1.+x)/2.)**2
25    a2=(1.+x**2)/2.
26    psim=log(a1*a2)-2.*atan(x)+pi/2.
27  else
28  ! STABLE CASE
29    psim=-4.7*zeta
30  endif
31
32end function psim
Note: See TracBrowser for help on using the repository browser.
hosted by ZAMG