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

Last change on this file since 5225 was 5153, checked in by abarral, 4 months ago

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