source: LMDZ6/branches/cirrus/libf/phylmd/evappot.F90 @ 5435

Last change on this file since 5435 was 4593, checked in by yann meurdesoif, 18 months ago

Replace #include (c preprocessor) by INCLUDE (fortran keyword)

in phylmd (except rrtm and ecrad) filtrez, dy3dmem and dyn3dcommon

Other directories will follow
YM
  • 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)
3
4IMPLICIT NONE
5
6INCLUDE "YOMCST.h"
7INCLUDE "YOETHF.h"
8INCLUDE "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.