GCC Code Coverage Report
Directory: ./ Exec Total Coverage
File: phylmd/Ocean_skin/phiw_m.F90 Lines: 0 5 0.0 %
Date: 2023-06-30 12:51:15 Branches: 0 2 0.0 %

Line Branch Exec Source
1
module Phiw_m
2
3
  implicit none
4
5
contains
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
24
end module Phiw_m