Parent Directory | Revision Log
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 |