/[lmdze]/trunk/phylmd/transp_lay.f90
ViewVC logotype

Contents of /trunk/phylmd/transp_lay.f90

Parent Directory Parent Directory | Revision Log Revision Log


Revision 328 - (show annotations)
Thu Jun 13 14:40:06 2019 UTC (4 years, 11 months ago) by guez
File size: 1610 byte(s)
Change all `.f` suffixes to `.f90`. (The opposite was done in revision
82.)  Because of change of philosopy in GNUmakefile: we already had a
rewritten rule for `.f`, so it does not make the makefile longer to
replace it by a rule for `.f90`. And it spares us options of
makedepf90 and of the compiler. Also we prepare the way for a simpler
`CMakeLists.txt`.

1 module transp_lay_m
2
3 IMPLICIT NONE
4
5 contains
6
7 SUBROUTINE transp_lay(paprs, t, q, u, v, geom, vtran_e, vtran_q, utran_e, &
8 utran_q)
9
10 ! Author: Z. X. Li (LMD/CNRS)
11 ! Date: April, 25th 1994
12 ! Objet : calculer le transport de l'\'energie et de la vapeur d'eau
13
14 USE dimphy, only: klon, klev
15 USE suphec_m, only: rcpd, rg, rlvtt
16
17 REAL, INTENT(IN):: paprs(klon, klev+1)
18 REAL, INTENT(IN):: t(klon, klev)
19 REAL, INTENT(IN):: q(klon, klev), u(klon, klev), v(klon, klev)
20 REAL, INTENT(IN):: geom(klon, klev)
21 REAL, INTENT(out):: vtran_e(klon, klev), vtran_q(klon, klev)
22 REAL, INTENT(out):: utran_e(klon, klev), utran_q(klon, klev)
23
24 ! Local:
25 INTEGER i, l
26 real esh
27
28 !------------------------------------------------------------------
29
30 DO l = 1, klev
31 DO i = 1, klon
32 utran_e(i, l) = 0.
33 utran_q(i, l) = 0.
34 vtran_e(i, l) = 0.
35 vtran_q(i, l) = 0.
36 END DO
37 END DO
38
39 DO l = 1, klev
40 DO i = 1, klon
41 esh = rcpd * t(i, l) + rlvtt * q(i, l) + geom(i, l)
42 utran_e(i, l) = utran_e(i, l) + u(i, l) * esh &
43 * (paprs(i, l) - paprs(i, l+1)) / rg
44 utran_q(i, l) = utran_q(i, l) + u(i, l) * q(i, l) &
45 * (paprs(i, l) - paprs(i, l+1)) / rg
46 vtran_e(i, l) = vtran_e(i, l) + v(i, l) * esh &
47 * (paprs(i, l) - paprs(i, l+1)) / rg
48 vtran_q(i, l) = vtran_q(i, l) + v(i, l) * q(i, l) &
49 * (paprs(i, l) - paprs(i, l+1)) / rg
50 END DO
51 END DO
52
53 END SUBROUTINE transp_lay
54
55 end module transp_lay_m

Properties

Name Value
svn:executable

  ViewVC Help
Powered by ViewVC 1.1.21