source: LMDZ6/branches/Amaury_dev/libf/phylmd/Ocean_skin/phiw_m.F90 @ 5119

Last change on this file since 5119 was 5119, checked in by abarral, 4 months ago

enforce PRIVATE by default in several modules, expose PUBLIC as needed
move eigen.f90 to obsolete/
(lint) aslong the way

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.