/[lmdze]/trunk/dyn3d/vitvert.f
ViewVC logotype

Contents of /trunk/dyn3d/vitvert.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 91 - (show annotations)
Wed Mar 26 17:18:58 2014 UTC (10 years, 1 month ago) by guez
File size: 887 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 module vitvert_m
2
3 IMPLICIT NONE
4
5 contains
6
7 SUBROUTINE vitvert(convm, w)
8
9 ! From libf/dyn3d/vitvert.F, version 1.1.1.1, 2004/05/19 12:53:05
10 ! Authors: P. Le Van, F. Hourdin
11
12 ! Objet : calcul de la vitesse verticale aux niveaux sigma
13
14 ! La vitesse verticale est orientée de haut en bas. Au sol, au
15 ! niveau sigma(1), w(i, j, 1) = 0. Au sommet, au niveau
16 ! sigma(llm+1), la vitesse verticale est aussi égale à 0 et n'est
17 ! pas stockée dans le tableau w.
18
19 USE dimens_m, ONLY : llm
20 USE disvert_m, ONLY : bp
21 USE paramet_m, ONLY : ip1jmp1
22
23 real, intent(in):: convm(ip1jmp1, llm)
24 REAL, intent(out):: w(ip1jmp1, llm)
25
26 ! Local:
27 INTEGER l
28
29 !------------------------------------------------------
30
31 forall (l = 2: llm) w(:, l) = convm(:, l) - bp(l) * convm(:, 1)
32 w(:, 1) = 0.
33
34 END SUBROUTINE vitvert
35
36 end module vitvert_m

  ViewVC Help
Powered by ViewVC 1.1.21