source: lmdz_wrf/trunk/WRFV3/lmdz/evappot.F90 @ 1465

Last change on this file since 1465 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

File size: 1.0 KB
Line 
1SUBROUTINE evappot(klon,nbsrf,ftsol,pplay,cdragh,  &
2       &    t_seri,q_seri,u_seri,v_seri,evap_pot)
3
4IMPLICIT NONE
5
6#include "YOMCST.h"
7#include "YOETHF.h"
8#include "FCTTRE.h"
9
10
11INTEGER :: klon, nbsrf
12REAL, DIMENSION(klon,nbsrf) :: ftsol,evap_pot
13REAL, DIMENSION(klon) :: pplay,t_seri,wind,q_seri,u_seri,v_seri,cdragh
14
15INTEGER :: nsrf,i
16REAL, DIMENSION(klon,nbsrf) :: qsat_ftsol
17REAL, DIMENSION(klon) :: rhos, norme_u
18REAL :: t_coup
19
20      t_coup=234.   ! Quelle horreur !!!!!
21
22DO 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
30ENDDO
31! ========================================================== c
32! Calcul de l'evaporation Potentielle
33
34
35rhos(:) = pplay(:)/(RD*t_seri(:))
36norme_u(:)=1.+sqrt(u_seri(:)*u_seri(:)+v_seri(:)*v_seri(:))
37DO nsrf = 1, nbsrf
38  evap_pot(:,nsrf)=rhos(:)*cdragh(:)*norme_u(:)*(qsat_ftsol(:,nsrf)-q_seri(:))
39ENDDO
40RETURN
41
42END
Note: See TracBrowser for help on using the repository browser.