source: LMDZ4/branches/V3_test/libf/phylmd/transp_lay.F @ 1215

Last change on this file since 1215 was 704, checked in by Laurent Fairhead, 18 years ago

Inclusion des modifs de Y. Meurdesoif pour la version V3
LF

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.6 KB
Line 
1      SUBROUTINE transp_lay (paprs,tsol,
2     e                   t, q, u, v, geom,
3     s                   vtran_e, vtran_q, utran_e, utran_q)
4c
5      USE dimphy
6      IMPLICIT none
7c======================================================================
8c Auteur(s): Z.X.Li (LMD/CNRS)
9c Date: le 25 avril 1994
10c Objet: Calculer le transport de l'energie et de la vapeur d'eau
11c======================================================================
12c
13cym#include "dimensions.h"
14cym#include "dimphy.h"
15#include "YOMCST.h"
16c
17      REAL paprs(klon,klev+1), tsol(klon)
18      REAL t(klon,klev), q(klon,klev), u(klon,klev), v(klon,klev)
19      REAL utran_e(klon,klev), utran_q(klon,klev)
20      REAL vtran_e(klon,klev), vtran_q(klon,klev)
21c
22      INTEGER i, l
23c     ------------------------------------------------------------------
24      REAL geom(klon,klev), esh
25c     ------------------------------------------------------------------
26      DO l = 1, klev
27      DO i = 1, klon
28         utran_e(i,l) = 0.0
29         utran_q(i,l) = 0.0
30         vtran_e(i,l) = 0.0
31         vtran_q(i,l) = 0.0
32      ENDDO
33      ENDDO
34c
35      DO l = 1, klev
36      DO i = 1, klon
37         esh = RCPD*t(i,l) + RLVTT*q(i,l) + geom(i,l)
38         utran_e(i,l)=utran_e(i,l)+ u(i,l)*esh*
39     .                (paprs(i,l)-paprs(i,l+1))/RG
40         utran_q(i,l)=utran_q(i,l)+ u(i,l)*q(i,l)
41     .                *(paprs(i,l)-paprs(i,l+1))/RG
42         vtran_e(i,l)=vtran_e(i,l)+ v(i,l)*esh*
43     .                (paprs(i,l)-paprs(i,l+1))/RG
44         vtran_q(i,l)=vtran_q(i,l)+ v(i,l)*q(i,l)
45     .                *(paprs(i,l)-paprs(i,l+1))/RG
46      ENDDO
47      ENDDO
48c
49      RETURN
50      END
Note: See TracBrowser for help on using the repository browser.