Parent Directory | Revision Log
Moved procedure inigeom into module comgeom. In disvert, renamed s_sampling to vert_sampling, following LMDZ. Removed choice strato1. In case read, read ap and bp instead of s (following LMDZ). Added argument phis to start_init_orog and start_init_dyn, and removed variable phis of module start_init_orog_m. In etat0 and start_init_orog, renamed relief to zmea_2d. In start_init_dyn, renamed psol to ps. In start_init_orog, renamed relief_hi to relief. No need to set phis(iim + 1, :) = phis(1, :), already done in grid_noro. Documentation for massbar out of SVN, in massbar.txt. Documentation was duplicated in massdair, but not relevant in massdair. In conflx, no need to initialize pen_[ud] and pde_[ud]. In flxasc, used intermediary variable fact (following LMDZ). In grid_noro, added local variable zmea0 for zmea not smoothed and computed zphi from zmea instead of zmea0 (following LMDZ). This changes the results of ce0l. Removed arguments pen_u and pde_d of phytrac and nflxtr, which were not used.
1 | guez | 78 | module massbarxy_m |
2 | guez | 3 | |
3 | guez | 78 | IMPLICIT NONE |
4 | guez | 3 | |
5 | guez | 78 | contains |
6 | guez | 3 | |
7 | guez | 78 | SUBROUTINE massbarxy(masse, massebxy) |
8 | guez | 3 | |
9 | guez | 78 | ! From LMDZ4/libf/dyn3d/massbarxy.F, version 1.1.1.1 2004/05/19 12:53:07 |
10 | guez | 3 | |
11 | guez | 78 | ! Calcule les moyennes en x et y de la masse d'air dans chaque maille. |
12 | ! Auteurs : P. Le Van, F. Hourdin. | ||
13 | |||
14 | USE dimens_m, ONLY: iim, llm | ||
15 | USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmp1 | ||
16 | USE comgeom, ONLY: alpha1, alpha2, alpha3, alpha4 | ||
17 | |||
18 | REAL, intent(in):: masse(ip1jmp1, llm) | ||
19 | real, intent(out):: massebxy(ip1jm, llm) | ||
20 | |||
21 | ! Local: | ||
22 | integer l, ij | ||
23 | |||
24 | !------------------------------------------------------------------- | ||
25 | |||
26 | DO l = 1, llm | ||
27 | DO ij = 1, ip1jm - 1 | ||
28 | massebxy(ij, l) = masse(ij, l) * alpha2(ij) & | ||
29 | + masse(ij + 1, l) * alpha3(ij + 1) & | ||
30 | + masse(ij + iip1, l) * alpha1(ij + iip1) & | ||
31 | + masse(ij + iip2, l) * alpha4(ij + iip2) | ||
32 | end DO | ||
33 | |||
34 | ! correction pour massebxy(iip1, j) | ||
35 | DO ij = iip1, ip1jm, iip1 | ||
36 | massebxy(ij, l) = massebxy(ij - iim, l) | ||
37 | end DO | ||
38 | end DO | ||
39 | |||
40 | END SUBROUTINE massbarxy | ||
41 | |||
42 | end module massbarxy_m |
ViewVC Help | |
Powered by ViewVC 1.1.21 |