Ignore:
Timestamp:
Feb 24, 2014, 4:05:47 PM (11 years ago)
Author:
Ehouarn Millour
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/dyn3dmem/leapfrog_mod.F90

    r1907 r1987  
    11MODULE leapfrog_mod
    22
    3   REAL,POINTER,SAVE :: ucov(:,:)
    4   REAL,POINTER,SAVE :: vcov(:,:)
    5   REAL,POINTER,SAVE :: teta(:,:)
    6   REAL,POINTER,SAVE :: ps(:)
    7   REAL,POINTER,SAVE :: masse(:,:)
    8   REAL,POINTER,SAVE :: phis(:)
    9   REAL,POINTER,SAVE :: q(:,:,:)
    10   REAL,POINTER,SAVE :: p(:,:)
    11   REAL,POINTER,SAVE :: pks(:)
    12   REAL,POINTER,SAVE :: pk(:,:)
    13   REAL,POINTER,SAVE :: pkf(:,:)
    14   REAL,POINTER,SAVE :: phi(:,:)
    15   REAL,POINTER,SAVE :: w(:,:)
     3  REAL,POINTER,SAVE :: ucov(:,:) ! zonal covariant wind
     4  REAL,POINTER,SAVE :: vcov(:,:) ! meridional covariant wind
     5  REAL,POINTER,SAVE :: teta(:,:) ! potential temperature
     6  REAL,POINTER,SAVE :: ps(:) ! surface pressure
     7  REAL,POINTER,SAVE :: masse(:,:) ! air mass
     8  REAL,POINTER,SAVE :: phis(:) ! geopotential at the surface
     9  REAL,POINTER,SAVE :: q(:,:,:) ! advected tracers
     10  REAL,POINTER,SAVE :: p(:,:) ! interlayer pressure
     11  REAL,POINTER,SAVE :: pks(:) ! Exner at the surface
     12  REAL,POINTER,SAVE :: pk(:,:) ! Exner at mid-layer
     13  REAL,POINTER,SAVE :: pkf(:,:) ! filtered Exner
     14  REAL,POINTER,SAVE :: phi(:,:) ! geopotential
     15  REAL,POINTER,SAVE :: w(:,:) ! vertical velocity
    1616  REAL,POINTER,SAVE :: pbaru(:,:)
    1717  REAL,POINTER,SAVE :: pbarv(:,:)
Note: See TracChangeset for help on using the changeset viewer.