source: LMDZ6/branches/blowing_snow/libf/phylmd/Ocean_skin/phiw_m.F90 @ 5450

Last change on this file since 5450 was 3834, checked in by Laurent Fairhead, 4 years ago

Renaming module files so that makelmdz and create_make_gcm finds them

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.