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

Last change on this file since 5146 was 5119, checked in by abarral, 11 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.