! pas reutilise ds nvelle formulation...
!!! tv pas reutiliser ds nvelle formulation... !! rho pas reutiliser ds nvelle formulation....
!!!! le pas reutiliser ds nvelle formulation....
!! u = vitesse du vent 1 er niveau !!!!! u_N utilisée pour le calcul des premier calcul...
!!! initialisation des parametres et premier calcul de z_0 cd et ch
!!! pour l'initialisation on fixe z_0 a 10-4 !!! on calcul ensuite cd a 10 m d'après smith !!! calcul de cd pour core pur
!!!! !!!!! on calcul le premier z_0 !!!!! si regime stable calcul du ch
!!!!! si regime instable calcul du ch
!!!! calcul ensuite de ce
! def ! def
! def ! def ! def
! def
!!!!!!! début des itérations....
!! calcul du zeta pour connaitre le signe et donc en deduire le regime...
! si regime stable alors valeurs pr chi et psi
!!! si regime instable alors valeurs pr psi et chi !!!
!!! calcul de cd pour core pur
!!! c'est sur cette ligne la qu'il y a un bug... ou sur la ligne du u_N
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real, | intent(in) | :: | du | |||
| real, | intent(in) | :: | dt | |||
| real, | intent(in) | :: | dq | |||
| real, | intent(in) | :: | t1 | |||
| real, | intent(in) | :: | q1 | |||
| real, | intent(in) | :: | zu | |||
| real, | intent(in) | :: | zt | |||
| real, | intent(in) | :: | zq | |||
| real, | intent(in) | :: | P | |||
| integer, | intent(in) | :: | n_it | |||
| logical, | intent(in) | :: | mixte | |||
| real, | intent(out), | DIMENSION(3) | :: | coeffs | ||
| real, | intent(out) | :: | rugosm | |||
| real, | intent(out) | :: | rugosh |