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_loc.F

    r1907 r1987  
    8080!      include "mpif.h"
    8181     
    82       INTEGER         longcles
    83       PARAMETER     ( longcles = 20 )
    84       REAL  clesphy0( longcles )
     82      INTEGER,PARAMETER :: longcles = 20
     83      REAL,INTENT(IN) :: clesphy0( longcles ) ! not used
     84      REAL,INTENT(IN) :: time_0 ! not used
     85
     86c   dynamical variables:
     87      REAL,INTENT(IN) :: ucov0(ijb_u:ije_u,llm)    ! zonal covariant wind
     88      REAL,INTENT(IN) :: vcov0(ijb_v:ije_v,llm)    ! meridional covariant wind
     89      REAL,INTENT(IN) :: teta0(ijb_u:ije_u,llm)    ! potential temperature
     90      REAL,INTENT(IN) :: q0(ijb_u:ije_u,llm,nqtot) ! advected tracers
     91      REAL,INTENT(IN) :: ps0(ijb_u:ije_u)          ! surface pressure (Pa)
     92      REAL,INTENT(IN) :: masse0(ijb_u:ije_u,llm)   ! air mass
     93      REAL,INTENT(IN) :: phis0(ijb_u:ije_u)        ! geopotentiat at the surface
    8594
    8695      real zqmin,zqmax
    87 
    88 c   variables dynamiques
    89       REAL :: vcov0(ijb_v:ije_v,llm),ucov0(ijb_u:ije_u,llm) ! vents covariants
    90       REAL :: teta0(ijb_u:ije_u,llm)                 ! temperature potentielle
    91       REAL :: q0(ijb_u:ije_u,llm,nqtot)              ! champs advectes
    92       REAL :: ps0(ijb_u:ije_u)                       ! pression  au sol
    93       REAL :: masse0(ijb_u:ije_u,llm)                ! masse d'air
    94       REAL :: phis0(ijb_u:ije_u)                     ! geopotentiel au sol
    9596
    9697!      REAL,SAVE,ALLOCATABLE :: p (:,:  )               ! pression aux interfac.des couches
     
    133134      REAL       time
    134135
    135       REAL  SSUM
    136       REAL time_0
     136      REAL  SSUM
    137137!      REAL,SAVE,ALLOCATABLE :: finvmaold(:,:)
    138138
     
    670670     
    671671      IF( forward. OR . leapf )  THEN
    672  
    673  
     672! Ehouarn: NB: fields sent to advtrac are those at the beginning of the time step
    674673         CALL caladvtrac_loc(q,pbaru,pbarv,
    675674     *        p, masse, dq,  teta,
     
    680679!      enddo
    681680
    682 c
     681! Ehouarn: Storage of mass flux for off-line tracers... not implemented...
     682
    683683      ENDIF ! of IF( forward. OR . leapf )
    684684
     
    11401140        endif
    11411141c$OMP BARRIER
     1142        CALL massdair_loc(p,masse)
     1143c$OMP BARRIER
    11421144
    11431145cc$OMP END PARALLEL
Note: See TracChangeset for help on using the changeset viewer.