My Project
 All Classes Files Functions Variables Macros
evappot.F90
Go to the documentation of this file.
1 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 REAL, DIMENSION(klon,nbsrf) :: qsat_ftsol
17 REAL, DIMENSION(klon) :: rhos, norme_u
18 REAL :: t_coup
19 
20  t_coup=234. ! Quelle horreur !!!!!
21 
22 DO nsrf = 1, nbsrf
23  DO i = 1, klon
24  IF (ftsol(i,nsrf).LT.t_coup) THEN
25  qsat_ftsol(i,nsrf) = qsats(ftsol(i,nsrf))/pplay(i)
26  ELSE
27  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 rhos(:) = pplay(:)/(rd*t_seri(:))
36 norme_u(:)=1.+sqrt(u_seri(:)*u_seri(:)+v_seri(:)*v_seri(:))
37 DO nsrf = 1, nbsrf
38  evap_pot(:,nsrf)=rhos(:)*cdragh(:)*norme_u(:)*(qsat_ftsol(:,nsrf)-q_seri(:))
39 ENDDO
40 RETURN
41 
42 END