|
Last change
on this file since 4095 was
1987,
checked in by Ehouarn Millour, 11 years ago
|
|
Add updating pressure, mass and Exner function (ie: all variables which depend on surface pressure) after adding physics tendencies (which include a surface pressure tendency).
Note that this change induces slight changes in GCM results with respect to previous svn version of the code, even if surface pressure tendency is zero (because of recomputation of polar values as an average over polar points on the dynamics grid).
EM
|
-
Property copyright set to
Name of program: LMDZ
Creation date: 1984
Version: LMDZ5
License: CeCILL version 2
Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
See the license file in the root directory
|
|
File size:
1.2 KB
|
| Rev | Line | |
|---|
| [1632] | 1 | SUBROUTINE pression_loc( ngrid, ap, bp, ps, p ) |
|---|
| [1987] | 2 | USE parallel_lmdz, ONLY: ij_begin, ij_end, ijb_u, ije_u, |
|---|
| 3 | & pole_nord, pole_sud, omp_chunk |
|---|
| [1632] | 4 | c |
|---|
| 5 | |
|---|
| 6 | c Auteurs : P. Le Van , Fr.Hourdin . |
|---|
| 7 | |
|---|
| 8 | c ************************************************************************ |
|---|
| 9 | c Calcule la pression p(l) aux differents niveaux l = 1 ( niveau du |
|---|
| 10 | c sol) a l = llm +1 ,ces niveaux correspondant aux interfaces des (llm) |
|---|
| 11 | c couches , avec p(ij,llm +1) = 0. et p(ij,1) = ps(ij) . |
|---|
| 12 | c ************************************************************************ |
|---|
| 13 | c |
|---|
| 14 | IMPLICIT NONE |
|---|
| 15 | c |
|---|
| 16 | #include "dimensions.h" |
|---|
| 17 | #include "paramet.h" |
|---|
| 18 | c |
|---|
| [1987] | 19 | INTEGER,INTENT(IN) :: ngrid ! not used |
|---|
| [1632] | 20 | INTEGER l,ij |
|---|
| 21 | |
|---|
| [1987] | 22 | REAL,INTENT(IN) :: ap( llmp1 ), bp( llmp1 ), ps( ijb_u:ije_u ) |
|---|
| 23 | REAL,INTENT(OUT) :: p( ijb_u:ije_u,llmp1 ) |
|---|
| [1632] | 24 | |
|---|
| 25 | INTEGER ijb,ije |
|---|
| 26 | |
|---|
| 27 | |
|---|
| 28 | ijb=ij_begin-iip1 |
|---|
| 29 | ije=ij_end+2*iip1 |
|---|
| 30 | |
|---|
| 31 | if (pole_nord) ijb=ij_begin |
|---|
| 32 | if (pole_sud) ije=ij_end |
|---|
| 33 | |
|---|
| 34 | c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) |
|---|
| 35 | DO l = 1, llmp1 |
|---|
| 36 | DO ij = ijb, ije |
|---|
| 37 | p(ij,l) = ap(l) + bp(l) * ps(ij) |
|---|
| 38 | ENDDO |
|---|
| 39 | ENDDO |
|---|
| 40 | c$OMP END DO NOWAIT |
|---|
| 41 | RETURN |
|---|
| 42 | END |
|---|
Note: See
TracBrowser
for help on using the repository browser.
