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 |
25 |
|
|
|