source: LMDZ5/branches/testing/libf/phylmd/transp_lay.F90 @ 4025

Last change on this file since 4025 was 2408, checked in by Laurent Fairhead, 9 years ago

Merged trunk changes r2298:2396 into testing branch
  • 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
  • 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 
[1992]1
[776]2! $Header$
[1992]3
4SUBROUTINE transp_lay(paprs, tsol, t, q, u, v, geom, vtran_e, vtran_q, &
5    utran_e, utran_q)
6
7  USE dimphy
8  IMPLICIT NONE
9  ! ======================================================================
10  ! Auteur(s): Z.X.Li (LMD/CNRS)
11  ! Date: le 25 avril 1994
12  ! Objet: Calculer le transport de l'energie et de la vapeur d'eau
13  ! ======================================================================
14
15  include "YOMCST.h"
16
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)
21
22  INTEGER i, l
23  ! ------------------------------------------------------------------
24  REAL geom(klon, klev), esh
25  ! ------------------------------------------------------------------
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    END DO
33  END DO
34
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*(paprs(i,l)-paprs(i,l+1))/ &
39        rg
40      utran_q(i, l) = utran_q(i, l) + u(i, l)*q(i, l)*(paprs(i,l)-paprs(i,l+1 &
41        ))/rg
42      vtran_e(i, l) = vtran_e(i, l) + v(i, l)*esh*(paprs(i,l)-paprs(i,l+1))/ &
43        rg
44      vtran_q(i, l) = vtran_q(i, l) + v(i, l)*q(i, l)*(paprs(i,l)-paprs(i,l+1 &
45        ))/rg
46    END DO
47  END DO
48
49  RETURN
50END SUBROUTINE transp_lay
Note: See TracBrowser for help on using the repository browser.