|
Last change
on this file since 3493 was
1508,
checked in by emillour, 9 years ago
|
|
Common dynamics:
Updates in the dynamics (seq and ) to keep up with updates
in LMDZ5 (up to LMDZ5 trunk, rev 2325):
IMPORTANT: Modifications for isotopes are only done in dyn3d, not in dyn3dpar
as in LMDZ5 these modifications were done in dyn3dmem.
Related LMDZ5 revisions are r2270 and r2281
- in dynlonlat_phylonlat:
- add module "grid_atob_m.F90" (a regridding utility
so far only used by phylmd/ce0l.F90, used to be dyn3d_common/grid_atob.F)
- in misc:
- follow up updates on wxios.F (add missing_val module variable)
- in dyn3d_common:
- pression.F => pression.F90
- misc_mod.F90: moved from misc to dyn3d_common
- added new iso_verif_dyn.F
- covcont.F => covcont.F90
- infotrac.F90 : add handling of isotopes (reading of corresponding
traceur.def for planets not implemented)
- dynetat0.F => dynetat0.F90 with some code factorization
- dynredem.F => dynredem.F90 with some code factorization
- added dynredem_mod.F90: routines used by dynredem
- iniacademic.F90 : added isotopes-related initialization for Earth case
- in dyn3d:
- added check_isotopes.F
- modified (isotopes) advtrac.F90, caladvtrac.F
- guide_mod.F90: ported updates
- leapfrog.F : (isotopes) updates (NB: call integrd with nqtot tracers)
- qminimium.F : adaptations for isotopes (copied over, except that #include
comvert.h is not needed).
- vlsplt.F: adaptations for isotopes (copied over, except than #include
logic.h, comvert.h not needed, and replace "include comconst.h" with
use comconst_mod, ONLY: pi)
- vlspltqs.F : same as vlsplt.F, but also keeping added modification for CP(T)
- in dyn3dpar:
- leapfrog_p.F: remove unecessary #ifdef CPP_EARTH cpp flag. and call
integrd_p with nqtot tracers (only important for Earth)
- dynredem_p.F => dynredem_p.F90 and some code factorization
- and no isotopes-relates changes in dyn3dpar (since these changes
have been made in LMDZ5 dyn3dmem).
EM
|
|
File size:
1.2 KB
|
| Line | |
|---|
| 1 | SUBROUTINE pression( ngrid, ap, bp, ps, p ) |
|---|
| 2 | ! |
|---|
| 3 | !------------------------------------------------------------------------------- |
|---|
| 4 | ! Authors: P. Le Van , Fr.Hourdin |
|---|
| 5 | !------------------------------------------------------------------------------- |
|---|
| 6 | ! Purpose: Compute pressure p(l) at different levels from l = 1 (ground level) |
|---|
| 7 | ! to l = llm +1. Those levels correspond to the llm layers interfaces, |
|---|
| 8 | ! with p(ij,llm+1) = 0. and p(ij,1) = ps(ij) . |
|---|
| 9 | !------------------------------------------------------------------------------- |
|---|
| 10 | IMPLICIT NONE |
|---|
| 11 | include "dimensions.h" |
|---|
| 12 | include "paramet.h" |
|---|
| 13 | !=============================================================================== |
|---|
| 14 | ! Arguments: |
|---|
| 15 | INTEGER, INTENT(IN) :: ngrid !--- NUMBER OF GRID POINTS |
|---|
| 16 | REAL, INTENT(IN) :: ap(llmp1), bp(llmp1) !--- HYBRID COEFFICIENTS |
|---|
| 17 | REAL, INTENT(IN) :: ps(ngrid) !--- SURFACE PRESSURE |
|---|
| 18 | REAL, INTENT(OUT) :: p(ngrid,llmp1) !--- 3D PRESSURE FIELD |
|---|
| 19 | !=============================================================================== |
|---|
| 20 | ! Local variables: |
|---|
| 21 | INTEGER :: l |
|---|
| 22 | !=============================================================================== |
|---|
| 23 | DO l=1,llmp1; p(:,l) = ap(l) + bp(l) * ps(:); END DO |
|---|
| 24 | |
|---|
| 25 | END SUBROUTINE pression |
|---|
| 26 | |
|---|
| 27 | |
|---|
Note: See
TracBrowser
for help on using the repository browser.
