/[lmdze]/trunk/Sources/phylmd/transp_lay.f
ViewVC logotype

Contents of /trunk/Sources/phylmd/transp_lay.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 38 - (show annotations)
Thu Jan 6 17:52:19 2011 UTC (13 years, 4 months ago) by guez
Original Path: trunk/libf/phylmd/transp_lay.f
File size: 1650 byte(s)
Extracted ASCII art from "inigeom" into a separate text file in the
documentation.

"test_disvert" now creates a separate file for layer thicknesses.

Moved variables from module "yomcst" to module "suphec_m" because this
is where those variables are defined. Kept in "yomcst" only parameters
of Earth orbit. Gave the attribute "parameter" to some variables of
module "suphec_m".

Variables of module "yoethf" were defined in procedure "suphec". Moved
these definitions to a new procedure "yoethf" in module "yoethf_m".

1 SUBROUTINE transp_lay (paprs,tsol,
2 e t, q, u, v, geom,
3 s vtran_e, vtran_q, utran_e, utran_q)
4 c
5 use dimens_m
6 use dimphy
7 use SUPHEC_M
8 IMPLICIT none
9 c======================================================================
10 c Auteur(s): Z.X.Li (LMD/CNRS)
11 c Date: le 25 avril 1994
12 c Objet: Calculer le transport de l'energie et de la vapeur d'eau
13 c======================================================================
14 c
15 c
16 REAL, intent(in):: paprs(klon,klev+1)
17 real 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 c
22 INTEGER i, l
23 c ------------------------------------------------------------------
24 REAL geom(klon,klev), esh
25 c ------------------------------------------------------------------
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
34 c
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
48 c
49 RETURN
50 END

Properties

Name Value
svn:executable

  ViewVC Help
Powered by ViewVC 1.1.21