source: LMDZ.3.3/trunk/libf/phylmd/transp.F @ 556

Last change on this file since 556 was 2, checked in by lmdz, 25 years ago

Initial revision

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