26 |
|
|
27 |
! Local: |
! Local: |
28 |
INTEGER i |
INTEGER i |
29 |
REAL zxmod, taux, tauy |
REAL wind |
30 |
|
|
31 |
!--------------------------------------------------------------- |
!--------------------------------------------------------------- |
32 |
|
|
33 |
DO i = 1, size(u) |
DO i = 1, size(u) |
34 |
zxmod = 1. + SQRT(u(i)**2 + v(i)**2) |
wind = SQRT(u(i)**2 + v(i)**2) |
35 |
taux = u(i) * zxmod * cd_m(i) |
ustarhb(i) = wind * (1 + wind) * cd_m(i) |
|
tauy = v(i) * zxmod * cd_m(i) |
|
|
ustarhb(i) = SQRT(taux**2 + tauy**2) |
|
36 |
ENDDO |
ENDDO |
37 |
|
|
38 |
end function ustarhb |
end function ustarhb |