source: LMDZ6/trunk/libf/phylmd/evappot.f90 @ 5301

Last change on this file since 5301 was 5285, checked in by abarral, 4 days ago

As discussed internally, remove generic ONLY: ... for new _mod_h 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: 1023 bytes
Line 
1SUBROUTINE evappot(klon,nbsrf,ftsol,pplay,cdragh,  &
2       &    t_seri,q_seri,u_seri,v_seri,evap_pot)
3
4USE yoethf_mod_h
5USE yomcst_mod_h
6IMPLICIT NONE
7
8
9INCLUDE "FCTTRE.h"
10
11
12INTEGER :: klon, nbsrf
13REAL, DIMENSION(klon,nbsrf) :: ftsol,evap_pot
14REAL, DIMENSION(klon) :: pplay,t_seri,wind,q_seri,u_seri,v_seri,cdragh
15
16INTEGER :: nsrf,i
17REAL, DIMENSION(klon,nbsrf) :: qsat_ftsol
18REAL, DIMENSION(klon) :: rhos, norme_u
19REAL :: t_coup
20
21      t_coup=234.   ! Quelle horreur !!!!!
22
23DO nsrf = 1, nbsrf
24   DO i = 1, klon
25      IF (ftsol(i,nsrf).LT.t_coup) THEN
26         qsat_ftsol(i,nsrf) = qsats(ftsol(i,nsrf))/pplay(i)
27      ELSE
28         qsat_ftsol(i,nsrf) = qsatl(ftsol(i,nsrf))/pplay(i)
29      ENDIF
30   ENDDO
31ENDDO
32! ========================================================== c
33! Calcul de l'evaporation Potentielle
34
35
36rhos(:) = pplay(:)/(RD*t_seri(:))
37norme_u(:)=1.+sqrt(u_seri(:)*u_seri(:)+v_seri(:)*v_seri(:))
38DO nsrf = 1, nbsrf
39  evap_pot(:,nsrf)=rhos(:)*cdragh(:)*norme_u(:)*(qsat_ftsol(:,nsrf)-q_seri(:))
40ENDDO
41RETURN
42
43END
Note: See TracBrowser for help on using the repository browser.