1 |
|
288 |
SUBROUTINE evappot(klon,nbsrf,ftsol,pplay,cdragh, & |
2 |
|
|
& t_seri,q_seri,u_seri,v_seri,evap_pot) |
3 |
|
|
|
4 |
|
|
IMPLICIT NONE |
5 |
|
|
|
6 |
|
|
INCLUDE "YOMCST.h" |
7 |
|
|
INCLUDE "YOETHF.h" |
8 |
|
|
INCLUDE "FCTTRE.h" |
9 |
|
|
|
10 |
|
|
|
11 |
|
|
INTEGER :: klon, nbsrf |
12 |
|
|
REAL, DIMENSION(klon,nbsrf) :: ftsol,evap_pot |
13 |
|
|
REAL, DIMENSION(klon) :: pplay,t_seri,wind,q_seri,u_seri,v_seri,cdragh |
14 |
|
|
|
15 |
|
|
INTEGER :: nsrf,i |
16 |
|
576 |
REAL, DIMENSION(klon,nbsrf) :: qsat_ftsol |
17 |
|
288 |
REAL, DIMENSION(klon) :: rhos, norme_u |
18 |
|
|
REAL :: t_coup |
19 |
|
|
|
20 |
|
|
t_coup=234. ! Quelle horreur !!!!! |
21 |
|
|
|
22 |
✓✓ |
1440 |
DO nsrf = 1, nbsrf |
23 |
✓✓ |
1146528 |
DO i = 1, klon |
24 |
✓✓ |
1146240 |
IF (ftsol(i,nsrf).LT.t_coup) THEN |
25 |
|
56480 |
qsat_ftsol(i,nsrf) = qsats(ftsol(i,nsrf))/pplay(i) |
26 |
|
|
ELSE |
27 |
|
1088608 |
qsat_ftsol(i,nsrf) = qsatl(ftsol(i,nsrf))/pplay(i) |
28 |
|
|
ENDIF |
29 |
|
|
ENDDO |
30 |
|
|
ENDDO |
31 |
|
|
! ========================================================== c |
32 |
|
|
! Calcul de l'evaporation Potentielle |
33 |
|
|
|
34 |
|
|
|
35 |
✓✓ |
286560 |
rhos(:) = pplay(:)/(RD*t_seri(:)) |
36 |
✓✓ |
286560 |
norme_u(:)=1.+sqrt(u_seri(:)*u_seri(:)+v_seri(:)*v_seri(:)) |
37 |
✓✓ |
1440 |
DO nsrf = 1, nbsrf |
38 |
✓✓ |
1146528 |
evap_pot(:,nsrf)=rhos(:)*cdragh(:)*norme_u(:)*(qsat_ftsol(:,nsrf)-q_seri(:)) |
39 |
|
|
ENDDO |
40 |
|
288 |
RETURN |
41 |
|
|
|
42 |
|
|
END |