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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 91 - (hide annotations)
Wed Mar 26 17:18:58 2014 UTC (10 years, 2 months ago) by guez
Original Path: trunk/phylmd/transp_lay.f
File size: 1555 byte(s)
Removed unused variables lock_startdate and time_stamp of module
calendar.

Noticed that physiq does not change the surface pressure. So removed
arguments ps and dpfi of subroutine addfi. dpfi was always 0. The
computation of ps in addfi included some averaging at the poles. In
principle, this does not change ps but in practice it does because of
finite numerical precision. So the results of the simulation are
changed. Removed arguments ps and dpfi of calfis. Removed argument
d_ps of physiq.

du at the poles is not computed by dudv1, so declare only the
corresponding latitudes in dudv1. caldyn passes only a section of the
array dudyn as argument.

Removed variable niadv of module iniadvtrac_m.

Declared arguments of exner_hyb as assumed-shape arrays and made all
other horizontal sizes in exner_hyb dynamic. This allows the external
program test_disvert to use exner_hyb at a single horizontal position.

1 guez 81 SUBROUTINE transp_lay(paprs, tsol, t, q, u, v, geom, vtran_e, vtran_q, &
2     utran_e, utran_q)
3    
4     USE dimens_m
5     USE dimphy
6     USE suphec_m
7     IMPLICIT NONE
8     ! ======================================================================
9     ! Auteur(s): Z.X.Li (LMD/CNRS)
10     ! Date: le 25 avril 1994
11     ! Objet: Calculer le transport de l'energie et de la vapeur d'eau
12     ! ======================================================================
13    
14    
15 guez 91 REAL, INTENT (IN):: paprs(klon, klev+1)
16 guez 81 REAL tsol(klon)
17 guez 91 REAL, INTENT (IN):: t(klon, klev)
18     REAL, INTENT (IN):: q(klon, klev), u(klon, klev), v(klon, klev)
19 guez 81 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
50     END SUBROUTINE transp_lay

Properties

Name Value
svn:executable

  ViewVC Help
Powered by ViewVC 1.1.21