source: LMDZ6/trunk/libf/phylmd/Ocean_skin/phiw_m.f90 @ 5447

Last change on this file since 5447 was 5268, checked in by abarral, 2 months ago

.f90 <-> .F90 depending on cpp key use

File size: 470 bytes
Line 
1module Phiw_m
2
3  implicit none
4
5contains
6
7  elemental real function Phiw(zL)
8
9    ! Fonction de stabilit\'e (Monin Obukhov)
10
11    real, intent(in):: zL ! profondeur  / longueur de Monin-Obukhov
12
13    !----------------------------------------------
14
15    if (zL < 0.) then
16       Phiw = (1. - 16. * zL)**(- 0.5)
17    else
18       Phiw = 1. + (5. * zL + 4. * zL**2)/(1. + 3. * zL + 0.25 * zL**2)
19       ! (Takaya et al. 2010)
20    end if
21
22  end function Phiw
23
24end module Phiw_m
Note: See TracBrowser for help on using the repository browser.