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 | module massbarxy_m |
2 | |
3 | IMPLICIT NONE |
4 | |
5 | contains |
6 | |
7 | SUBROUTINE massbarxy(masse, massebxy) |
8 | |
9 | ! From LMDZ4/libf/dyn3d/massbarxy.F, version 1.1.1.1 2004/05/19 12:53:07 |
10 | |
11 | ! 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 |