source: LMDZ6/branches/Amaury_dev/libf/phylmd/evappot.F90 @ 5143

Last change on this file since 5143 was 5143, checked in by abarral, 3 months ago

Put YOEGWD.h, FCTTRE.h into modules
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 1.0 KB
Line 
1SUBROUTINE evappot(klon,nbsrf,ftsol,pplay,cdragh,  &
2      t_seri,q_seri,u_seri,v_seri,evap_pot)
3USE lmdz_YOETHF
4USE lmdz_fcttre, ONLY: foeew, foede, qsats, qsatl, dqsats, dqsatl, thermcep
5
6IMPLICIT NONE
7
8INCLUDE "YOMCST.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)<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.