source: LMDZ6/branches/SETHET_DECOUPLE/libf/phylmd/evappot.F90 @ 5182

Last change on this file since 5182 was 1907, checked in by lguez, 11 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
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

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.
  • 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
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.