Changeset 1189 for trunk


Ignore:
Timestamp:
Feb 25, 2014, 10:59:48 AM (11 years ago)
Author:
emillour
Message:

Common dynamics: a couple of bug fixes:

  • Correctly account for the change in pressure, mass, etc. after modifying surface pressure following a call to the physics.
  • Corrected tracer advection, which is computed using values at the beginning of the time step, so it is done at Matsuno forward step.

EM

Location:
trunk
Files:
15 edited

Legend:

Unmodified
Added
Removed

  • trunk/DOC/chantiers/commit_importants.log

    r1107 r1189  
    13181318  Some Mars-specific variables (ecritphy,ecritstart,timestart) added in
    13191319  control_mod.F and (hour_ini) in temps.h
    1320  
     1320
     1321**********************
     1322**** commit_v1189 ****
     1323**********************
     1324Ehouarn : updates in the dynamics:
     1325- Correctly account for the change in pressure, masse, etc. after modifying surface
     1326  pressure following a call to the physics
     1327- Corrected tracer advection, which is computed using values at the beginning of the
     1328  time step, so it is done at Matsuno forward step.
     1329

  • trunk/LMDZ.COMMON/libf/dyn3d/addfi.F

    r7 r1189  
    5555c    -----------
    5656c
    57       REAL pdt
     57      REAL,INTENT(IN) :: pdt ! time step for the integration (s)
    5858c
    59       REAL pvcov(ip1jm,llm),pucov(ip1jmp1,llm)
    60       REAL pteta(ip1jmp1,llm),pq(ip1jmp1,llm,nqtot),pps(ip1jmp1)
     59      REAL,INTENT(INOUT) :: pvcov(ip1jm,llm) ! covariant meridional wind
     60      REAL,INTENT(INOUT) :: pucov(ip1jmp1,llm) ! covariant zonal wind
     61      REAL,INTENT(INOUT) :: pteta(ip1jmp1,llm) ! potential temperature
     62      REAL,INTENT(INOUT) :: pq(ip1jmp1,llm,nqtot) ! tracers
     63      REAL,INTENT(INOUT) :: pps(ip1jmp1) ! surface pressure (Pa)
     64c respective tendencies (.../s) to add
     65      REAL,INTENT(IN) :: pdvfi(ip1jm,llm)
     66      REAL,INTENT(IN) :: pdufi(ip1jmp1,llm)
     67      REAL,INTENT(IN) :: pdqfi(ip1jmp1,llm,nqtot)
     68      REAL,INTENT(IN) :: pdhfi(ip1jmp1,llm)
     69      REAL,INTENT(IN) :: pdpfi(ip1jmp1)
    6170c
    62       REAL pdvfi(ip1jm,llm),pdufi(ip1jmp1,llm)
    63       REAL pdqfi(ip1jmp1,llm,nqtot),pdhfi(ip1jmp1,llm),pdpfi(ip1jmp1)
    64 c
    65       LOGICAL leapf,forward
     71      LOGICAL,INTENT(IN) :: leapf,forward ! not used
    6672c
    6773c
    … …  
    7177      REAL xpn(iim),xps(iim),tpn,tps
    7278      INTEGER j,k,iq,ij
    73       REAL qtestw, qtestt
    74       PARAMETER ( qtestw = 1.0e-15 )
    75       PARAMETER ( qtestt = 1.0e-40 )
     79      REAL,PARAMETER :: qtestw = 1.0e-15
     80      REAL,PARAMETER :: qtestt = 1.0e-40
    7681
    7782      REAL SSUM

  • trunk/LMDZ.COMMON/libf/dyn3d/advtrac.F90

    r1019 r1189  
    99  !            M.A Filiberti (04/2002)
    1010  !
    11   USE infotrac
    12   USE control_mod
     11  USE infotrac, ONLY: nqtot, iadv
     12  USE control_mod, ONLY: iapp_tracvl, day_step
    1313
    1414
    … …  
    3030  !     Arguments
    3131  !-------------------------------------------------------------------
     32  INTEGER,INTENT(OUT) :: iapptrac
     33  REAL,INTENT(IN) :: pbaru(ip1jmp1,llm)
     34  REAL,INTENT(IN) :: pbarv(ip1jm,llm)
     35  REAL,INTENT(INOUT) :: q(ip1jmp1,llm,nqtot)
     36  REAL,INTENT(IN) :: masse(ip1jmp1,llm)
     37  REAL,INTENT(IN) :: p( ip1jmp1,llmp1 )
     38  REAL,INTENT(IN) :: teta(ip1jmp1,llm)
     39  REAL,INTENT(IN) :: pk(ip1jmp1,llm)
     40  REAL,INTENT(OUT) :: flxw(ip1jmp1,llm)
     41  !-------------------------------------------------------------------
    3242  !     Ajout PPM
    3343  !--------------------------------------------------------
    3444  REAL massebx(ip1jmp1,llm),masseby(ip1jm,llm)
    35   !--------------------------------------------------------
    36   INTEGER iapptrac
    37   REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
    38   REAL q(ip1jmp1,llm,nqtot),masse(ip1jmp1,llm)
    39   REAL p( ip1jmp1,llmp1 ),teta(ip1jmp1,llm)
    40   REAL pk(ip1jmp1,llm)
    41   REAL flxw(ip1jmp1,llm)
    42 
    4345  !-------------------------------------------------------------
    4446  !     Variables locales

  • trunk/LMDZ.COMMON/libf/dyn3d/caldyn.F

    r5 r1189  
    11!
    2 ! $Header$
     2! $Id: $
    33!
    4 c
    5 c
    64      SUBROUTINE caldyn
    75     $ (itau,ucov,vcov,teta,ps,masse,pk,pkf,tsurpk,phis ,
    … …  
    108      IMPLICIT NONE
    119
    12 c=======================================================================
    13 c
    14 c  Auteur :  P. Le Van
    15 c
    16 c   Objet:
    17 c   ------
    18 c
    19 c   Calcul des tendances dynamiques.
    20 c
    21 c Modif 04/93 F.Forget
    22 c=======================================================================
     10!=======================================================================
     11!
     12!  Auteur :  P. Le Van
     13!
     14!   Objet:
     15!   ------
     16!
     17!   Calcul des tendances dynamiques.
     18!
     19! Modif 04/93 F.Forget
     20!=======================================================================
    2321
    24 c-----------------------------------------------------------------------
    25 c   0. Declarations:
    26 c   ----------------
     22!-----------------------------------------------------------------------
     23!   0. Declarations:
     24!   ----------------
    2725
    2826#include "dimensions.h"
    … …  
    3230#include "comgeom.h"
    3331
    34 c   Arguments:
    35 c   ----------
     32!   Arguments:
     33!   ----------
    3634
    37       LOGICAL conser
     35      LOGICAL,INTENT(IN) :: conser ! triggers printing some diagnostics
     36      INTEGER,INTENT(IN) :: itau ! time step index
     37      REAL,INTENT(IN) :: vcov(ip1jm,llm) ! covariant meridional wind
     38      REAL,INTENT(IN) :: ucov(ip1jmp1,llm) ! covariant zonal wind
     39      REAL,INTENT(IN) :: teta(ip1jmp1,llm) ! potential temperature
     40      REAL,INTENT(IN) :: ps(ip1jmp1) ! surface pressure
     41      REAL,INTENT(IN) :: phis(ip1jmp1) ! geopotential at the surface
     42      REAL,INTENT(IN) :: pk(ip1jmp1,llm) ! Exner at mid-layer
     43      REAL,INTENT(IN) :: pkf(ip1jmp1,llm) ! filtered Exner
     44      REAL,INTENT(IN) :: tsurpk(ip1jmp1,llm) ! cpp * temperature / pk
     45      REAL,INTENT(IN) :: phi(ip1jmp1,llm) ! geopotential
     46      REAL,INTENT(OUT) :: masse(ip1jmp1,llm) ! air mass
     47      REAL,INTENT(OUT) :: dv(ip1jm,llm) ! tendency on vcov
     48      REAL,INTENT(OUT) :: du(ip1jmp1,llm) ! tendency on ucov
     49      REAL,INTENT(OUT) :: dteta(ip1jmp1,llm) ! tenddency on teta
     50      REAL,INTENT(OUT) :: dp(ip1jmp1) ! tendency on ps
     51      REAL,INTENT(OUT) :: w(ip1jmp1,llm) ! vertical velocity
     52      REAL,INTENT(OUT) :: pbaru(ip1jmp1,llm) ! mass flux in the zonal direction
     53      REAL,INTENT(OUT) :: pbarv(ip1jm,llm) ! mass flux in the meridional direction
     54      REAL,INTENT(IN) :: time ! current time
    3855
    39       INTEGER itau
    40       REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm),teta(ip1jmp1,llm)
    41       REAL ps(ip1jmp1),phis(ip1jmp1)
    42       REAL pk(ip1jmp1,llm),pkf(ip1jmp1,llm)
    43       REAL tsurpk(ip1jmp1,llm)
     56!   Local:
     57!   ------
     58
    4459      REAL vcont(ip1jm,llm),ucont(ip1jmp1,llm)
    45       REAL phi(ip1jmp1,llm),masse(ip1jmp1,llm)
    46       REAL dv(ip1jm,llm),du(ip1jmp1,llm)
    47       REAL dteta(ip1jmp1,llm),dp(ip1jmp1)
    48       REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
    49       REAL time
    50 
    51 c   Local:
    52 c   ------
    53 
    5460      REAL ang(ip1jmp1,llm),p(ip1jmp1,llmp1)
    5561      REAL massebx(ip1jmp1,llm),masseby(ip1jm,llm),psexbarxy(ip1jm)
    5662      REAL vorpot(ip1jm,llm)
    57       REAL w(ip1jmp1,llm),ecin(ip1jmp1,llm),convm(ip1jmp1,llm)
     63      REAL ecin(ip1jmp1,llm),convm(ip1jmp1,llm)
    5864      REAL bern(ip1jmp1,llm)
    5965      REAL massebxy(ip1jm,llm)
    … …  
    6268      INTEGER   ij,l
    6369
    64 c-----------------------------------------------------------------------
    65 c   Calcul des tendances dynamiques:
    66 c   --------------------------------
     70!-----------------------------------------------------------------------
     71!   Compute dynamical tendencies:
     72!--------------------------------
    6773
     74      ! compute contravariant winds ucont() and vcont
    6875      CALL covcont  ( llm    , ucov    , vcov , ucont, vcont        )
     76      ! compute pressure p()
    6977      CALL pression ( ip1jmp1, ap      , bp   ,  ps  , p            )
     78      ! compute psexbarxy() XY-area weighted-averaged surface pressure (what for?)
    7079      CALL psextbar (   ps   , psexbarxy                            )
     80      ! compute mass in each atmospheric mesh: masse()
    7181      CALL massdair (    p   , masse                                )
     82      ! compute X and Y-averages of mass, massebx() and masseby()
    7283      CALL massbar  (   masse, massebx , masseby                    )
     84      ! compute XY-average of mass, massebxy()
    7385      call massbarxy(   masse, massebxy                             )
     86      ! compute mass fluxes pbaru() and pbarv()
    7487      CALL flumass  ( massebx, masseby , vcont, ucont ,pbaru, pbarv )
     88      ! compute dteta() , horizontal converging flux of theta
    7589      CALL dteta1   (   teta , pbaru   , pbarv, dteta               )
     90      ! compute convm(), horizontal converging flux of mass
    7691      CALL convmas  (   pbaru, pbarv   , convm                      )
    7792
     93      ! compute pressure variation due to mass convergence
    7894      DO ij =1, ip1jmp1
    7995         dp( ij ) = convm( ij,1 ) / airesurg( ij )
    8096      ENDDO
    8197
     98      ! compute vertical velocity w()
    8299      CALL vitvert ( convm  , w                                  )
     100      ! compute potential vorticity vorpot()
    83101      CALL tourpot ( vcov   , ucov  , massebxy  , vorpot         )
     102      ! compute rotation induced du() and dv()
    84103      CALL dudv1   ( vorpot , pbaru , pbarv     , du     , dv    )
     104      ! compute kinetic energy ecin()
    85105      CALL enercin ( vcov   , ucov  , vcont     , ucont  , ecin  )
     106      ! compute Bernouilli function bern()
    86107      CALL bernoui ( ip1jmp1, llm   , phi       , ecin   , bern  )
     108      ! compute and add du() and dv() contributions from Bernouilli and pressure
    87109      CALL dudv2   ( tsurpk , pkf   , bern      , du     , dv    )
    88110
    … …  
    94116      ENDDO
    95117
    96 
     118      ! compute vertical advection contributions to du(), dv() and dteta()
    97119      CALL advect( ang, vcov, teta, w, massebx, masseby, du, dv,dteta )
    98120
    99 C  WARNING probleme de peridocite de dv sur les PC/linux. Pb d'arrondi
    100 C          probablement. Observe sur le code compile avec pgf90 3.0-1
     121!  WARNING probleme de peridocite de dv sur les PC/linux. Pb d'arrondi
     122!          probablement. Observe sur le code compile avec pgf90 3.0-1
    101123
    102124      DO l = 1, llm
    103125         DO ij = 1, ip1jm, iip1
    104126           IF( dv(ij,l).NE.dv(ij+iim,l) )  THEN
    105 c         PRINT *,'!!!ATTENTION!!! probleme de periodicite sur vcov', 
    106 c    ,   ' dans caldyn'
    107 c         PRINT *,' l,  ij = ', l, ij, ij+iim,dv(ij+iim,l),dv(ij,l)
     127!         PRINT *,'!!!ATTENTION!!! probleme de periodicite sur vcov', 
     128!    ,   ' dans caldyn'
     129!         PRINT *,' l,  ij = ', l, ij, ij+iim,dv(ij+iim,l),dv(ij,l)
    108130          dv(ij+iim,l) = dv(ij,l)
    109131          endif
    110132         enddo
    111133      enddo
    112 c-----------------------------------------------------------------------
    113 c   Sorties eventuelles des variables de controle:
    114 c   ----------------------------------------------
     134
     135!-----------------------------------------------------------------------
     136!   Output some control variables:
     137!---------------------------------
    115138
    116139      IF( conser )  THEN
    … …  
    120143      ENDIF
    121144
    122       RETURN
    123145      END

  • trunk/LMDZ.COMMON/libf/dyn3d/calfis.F

    r1129 r1189  
    2929c    Auteur :  P. Le Van, F. Hourdin
    3030c   .........
    31       USE infotrac
    32       USE control_mod
     31      USE infotrac, ONLY: nqtot, niadv, tname
     32      USE control_mod, ONLY: planet_type, nsplit_phys
    3333      USE write_field
    3434      USE cpdet_mod, only: t2tpot,tpot2t
    … …  
    106106c    Arguments :
    107107c    -----------
    108       LOGICAL  lafin
    109 
    110 
    111       REAL pvcov(iip1,jjm,llm)
    112       REAL pucov(iip1,jjp1,llm)
    113       REAL pteta(iip1,jjp1,llm)
    114       REAL pmasse(iip1,jjp1,llm)
    115       REAL pq(iip1,jjp1,llm,nqtot)
    116       REAL pphis(iip1,jjp1)
    117       REAL pphi(iip1,jjp1,llm)
    118 
    119       REAL pdvcov(iip1,jjm,llm)
    120       REAL pducov(iip1,jjp1,llm)
    121       REAL pdteta(iip1,jjp1,llm)
     108      LOGICAL,INTENT(IN) ::  lafin ! .true. for the very last call to physics
     109      REAL,INTENT(IN) :: jD_cur, jH_cur
     110      REAL,INTENT(IN) :: pvcov(iip1,jjm,llm) ! covariant meridional velocity
     111      REAL,INTENT(IN) :: pucov(iip1,jjp1,llm) ! covariant zonal velocity
     112      REAL,INTENT(IN) :: pteta(iip1,jjp1,llm) ! potential temperature
     113      REAL,INTENT(IN) :: pmasse(iip1,jjp1,llm) ! mass in each cell ! not used
     114      REAL,INTENT(IN) :: pq(iip1,jjp1,llm,nqtot) ! tracers
     115      REAL,INTENT(IN) :: pphis(iip1,jjp1) ! surface geopotential
     116      REAL,INTENT(IN) :: pphi(iip1,jjp1,llm) ! geopotential
     117
     118      REAL,INTENT(IN) :: pdvcov(iip1,jjm,llm) ! dynamical tendency on vcov
     119      REAL,INTENT(IN) :: pducov(iip1,jjp1,llm) ! dynamical tendency on ucov
     120      REAL,INTENT(IN) :: pdteta(iip1,jjp1,llm) ! dynamical tendency on teta
    122121! commentaire SL: pdq ne sert que pour le calcul de pcvgq,
    123122! qui lui meme ne sert a rien dans la routine telle qu'elle est
    124123! ecrite, et que j'ai donc commente....
    125       REAL pdq(iip1,jjp1,llm,nqtot)
    126 
    127       REAL pps(iip1,jjp1)
    128       REAL pp(iip1,jjp1,llmp1)
    129       REAL ppk(iip1,jjp1,llm)
    130 
    131 c TENDENCIES in */s
    132       REAL pdvfi(iip1,jjm,llm)
    133       REAL pdufi(iip1,jjp1,llm)
    134       REAL pdhfi(iip1,jjp1,llm)
    135       REAL pdqfi(iip1,jjp1,llm,nqtot)
    136       REAL pdpsfi(iip1,jjp1)
     124      REAL,INTENT(IN) :: pdq(iip1,jjp1,llm,nqtot) ! dynamical tendency on tracers
     125      ! NB: pdq is only used to compute pcvgq which is in fact not used...
     126
     127      REAL,INTENT(IN) :: pps(iip1,jjp1) ! surface pressure (Pa)
     128      REAL,INTENT(IN) :: pp(iip1,jjp1,llmp1) ! pressure at mesh interfaces (Pa)
     129      REAL,INTENT(IN) :: ppk(iip1,jjp1,llm) ! Exner at mid-layer
     130      REAL,INTENT(IN) :: flxw(iip1,jjp1,llm)  ! Vertical mass flux on dynamics grid
     131
     132      ! tendencies (in */s) from the physics
     133      REAL,INTENT(OUT) :: pdvfi(iip1,jjm,llm) ! tendency on covariant meridional wind
     134      REAL,INTENT(OUT) :: pdufi(iip1,jjp1,llm) ! tendency on covariant zonal wind
     135      REAL,INTENT(OUT) :: pdhfi(iip1,jjp1,llm) ! tendency on potential temperature (K/s)
     136      REAL,INTENT(OUT) :: pdqfi(iip1,jjp1,llm,nqtot) ! tendency on tracers
     137      REAL,INTENT(OUT) :: pdpsfi(iip1,jjp1) ! tendency on surface pressure (Pa/s)
    137138
    138139
    … …  
    173174
    174175cIM diagnostique PVteta, Amip2
    175       INTEGER ntetaSTD
    176       PARAMETER(ntetaSTD=3)
    177       REAL rtetaSTD(ntetaSTD)
    178       DATA rtetaSTD/350., 380., 405./ ! Earth-specific values, beware !!
     176      INTEGER,PARAMETER :: ntetaSTD=3
     177      REAL,SAVE :: rtetaSTD(ntetaSTD)=(/350.,380.,405./) ! Earth-specific, beware !!
    179178      REAL PVteta(ngridmx,ntetaSTD)
    180179
    181       REAL flxw(iip1,jjp1,llm)  ! Flux de masse verticale sur la grille dynamique
    182180      REAL flxwfi(ngridmx,llm)  ! Flux de masse verticale sur la grille physiq
    183181     
    184182      REAL SSUM
    185183
    186       LOGICAL firstcal, debut
    187       DATA firstcal/.true./
    188       SAVE firstcal,debut
     184      LOGICAL,SAVE :: firstcal=.true., debut=.true.
    189185!      REAL rdayvrai
    190       REAL, intent(in):: jD_cur, jH_cur
    191186
    192187      LOGICAL tracerdyn ! for generic/mars physics call ; possibly to get rid of

  • trunk/LMDZ.COMMON/libf/dyn3d/conf_gcm.F

    r1056 r1189  
    174174      starttime = 0
    175175      CALL getin('starttime',starttime)
    176 
     176     
     177      ! Mars: time of start for run in "start.nc" (when there are multiple time
     178      !       steps stored in the file)
     179      timestart=-9999 ! default value; if <0, use last stored time
     180      call getin("timestart",timestart)
     181     
    177182!Config  Key  = less1day
    178183!Config  Desc = Possibilite d'integrer moins d'un jour

  • trunk/LMDZ.COMMON/libf/dyn3d/dissip.F

    r1 r1189  
    11!
    2 ! $Header$
     2! $Id: $
    33!
    44      SUBROUTINE dissip( vcov,ucov,teta,p, dv,du,dh )
    … …  
    3535c   ----------
    3636
    37       REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm),teta(ip1jmp1,llm)
    38       REAL  p( ip1jmp1,llmp1 )
    39       REAL dv(ip1jm,llm),du(ip1jmp1,llm),dh(ip1jmp1,llm)
     37      REAL,INTENT(IN) :: vcov(ip1jm,llm) ! covariant meridional wind
     38      REAL,INTENT(IN) :: ucov(ip1jmp1,llm) ! covariant zonal wind
     39      REAL,INTENT(IN) :: teta(ip1jmp1,llm) ! potential temperature
     40      REAL,INTENT(IN) :: p(ip1jmp1,llmp1) ! pressure
     41      ! tendencies (.../s) on covariant winds and potential temperature
     42      REAL,INTENT(OUT) :: dv(ip1jm,llm)
     43      REAL,INTENT(OUT) :: du(ip1jmp1,llm)
     44      REAL,INTENT(OUT) :: dh(ip1jmp1,llm)
    4045
    4146c   Local:

  • trunk/LMDZ.COMMON/libf/dyn3d/leapfrog.F

    r1107 r1189  
    1212      use IOIPSL
    1313#endif
    14       USE infotrac
     14      USE infotrac, ONLY: nqtot
    1515      USE guide_mod, ONLY : guide_main
    16       USE write_field
    17       USE control_mod, only: planet_type,nday,day_step,iperiod,iphysiq,
     16      USE write_field, ONLY: writefield
     17      USE control_mod, ONLY: planet_type,nday,day_step,iperiod,iphysiq,
    1818     &                       less1day,fractday,ndynstep,iconser,
    1919     &                       dissip_period,offline,ip_ebil_dyn,
    … …  
    7575! #include "clesphys.h"
    7676
    77       real zqmin,zqmax
    78 
    79 c   variables dynamiques
    80       REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm) ! vents covariants
    81       REAL teta(ip1jmp1,llm)                 ! temperature potentielle
    82       REAL q(ip1jmp1,llm,nqtot)               ! champs advectes
    83       REAL ps(ip1jmp1)                       ! pression  au sol
    84       REAL p (ip1jmp1,llmp1  )               ! pression aux interfac.des couches
    85       REAL pks(ip1jmp1)                      ! exner au  sol
    86       REAL pk(ip1jmp1,llm)                   ! exner au milieu des couches
    87       REAL pkf(ip1jmp1,llm)                  ! exner filt.au milieu des couches
    88       REAL masse(ip1jmp1,llm)                ! masse d'air
    89       REAL phis(ip1jmp1)                     ! geopotentiel au sol
    90       REAL phi(ip1jmp1,llm)                  ! geopotentiel
    91       REAL w(ip1jmp1,llm)                    ! vitesse verticale
     77      REAL,INTENT(IN) :: time_0 ! not used
     78
     79c   dynamical variables:
     80      REAL,INTENT(INOUT) :: ucov(ip1jmp1,llm)    ! zonal covariant wind
     81      REAL,INTENT(INOUT) :: vcov(ip1jm,llm)      ! meridional covariant wind
     82      REAL,INTENT(INOUT) :: teta(ip1jmp1,llm)    ! potential temperature
     83      REAL,INTENT(INOUT) :: ps(ip1jmp1)          ! surface pressure (Pa)
     84      REAL,INTENT(INOUT) :: masse(ip1jmp1,llm)   ! air mass
     85      REAL,INTENT(INOUT) :: phis(ip1jmp1)        ! geopotentiat at the surface
     86      REAL,INTENT(INOUT) :: q(ip1jmp1,llm,nqtot) ! advected tracers
     87
     88      REAL p (ip1jmp1,llmp1  )               ! interlayer pressure
     89      REAL pks(ip1jmp1)                      ! exner at the surface
     90      REAL pk(ip1jmp1,llm)                   ! exner at mid-layer
     91      REAL pkf(ip1jmp1,llm)                  ! filtered exner at mid-layer
     92      REAL phi(ip1jmp1,llm)                  ! geopotential
     93      REAL w(ip1jmp1,llm)                    ! vertical velocity
    9294! ADAPTATION GCM POUR CP(T)
    9395      REAL temp(ip1jmp1,llm)                 ! temperature 
    9496      REAL tsurpk(ip1jmp1,llm)               ! cpp*T/pk 
     97
     98      real zqmin,zqmax
    9599
    96100c variables dynamiques intermediaire pour le transport
    … …  
    133137
    134138      REAL  SSUM
    135       REAL time_0
    136139!     REAL finvmaold(ip1jmp1,llm)
    137140
    … …  
    410413c   -------------------------------------------------------------
    411414
    412 !      IF( forward. OR . leapf )  THEN
    413 ! Ehouarn: NB: at this point p with ps are not synchronized
    414 !              (whereas mass and ps are...)
    415       IF((.not.forward).OR. leapf )  THEN
    416         ! Ehouarn: gather mass fluxes during backward Matsuno or LF step
     415      IF( forward. OR . leapf )  THEN
     416! Ehouarn: NB: fields sent to advtrac are those at the beginning of the time step
    417417         CALL caladvtrac(q,pbaru,pbarv,
    418418     *        p, masse, dq,  teta,
    … …  
    547547     $                  ucov, vcov, teta , q   ,ps ,
    548548     $                 dufi, dvfi, dtetafi , dqfi ,dpfi  )
    549 
     549          ! since addfi updates ps(), also update p(), masse() and pk()
     550          CALL pression (ip1jmp1,ap,bp,ps,p)
     551          CALL massdair(p,masse)
     552          if (pressure_exner) then
     553            CALL exner_hyb(ip1jmp1,ps,p,alpha,beta,pks,pk,pkf)
     554          else
     555            CALL exner_milieu(ip1jmp1,ps,p,beta,pks,pk,pkf)
     556          endif
     557         
    550558c      Couche superieure :
    551559c      -------------------
    … …  
    600608          CALL exner_milieu( ip1jmp1, ps, p, beta, pks, pk, pkf )
    601609        endif
    602 
     610        CALL massdair(p,masse)
    603611
    604612c-----------------------------------------------------------------------

  • trunk/LMDZ.COMMON/libf/dyn3dpar/addfi_p.F

    r1019 r1189  
    5555c    -----------
    5656c
    57       REAL pdt
    58 c
    59       REAL pvcov(ip1jm,llm),pucov(ip1jmp1,llm)
    60       REAL pteta(ip1jmp1,llm),pq(ip1jmp1,llm,nqtot),pps(ip1jmp1)
    61 c
    62       REAL pdvfi(ip1jm,llm),pdufi(ip1jmp1,llm)
    63       REAL pdqfi(ip1jmp1,llm,nqtot),pdhfi(ip1jmp1,llm),pdpfi(ip1jmp1)
    64 c
    65       LOGICAL leapf,forward
     57      REAL,INTENT(IN) :: pdt ! time step for the integration (s)
     58c
     59      REAL,INTENT(INOUT) :: pvcov(ip1jm,llm) ! covariant meridional wind
     60      REAL,INTENT(INOUT) :: pucov(ip1jmp1,llm) ! covariant zonal wind
     61      REAL,INTENT(INOUT) :: pteta(ip1jmp1,llm) ! potential temperature
     62      REAL,INTENT(INOUT) :: pq(ip1jmp1,llm,nqtot) ! tracers
     63      REAL,INTENT(INOUT) :: pps(ip1jmp1) ! surface pressure (Pa)
     64c respective tendencies (.../s) to add
     65      REAL,INTENT(IN) :: pdvfi(ip1jm,llm)
     66      REAL,INTENT(IN) :: pdufi(ip1jmp1,llm)
     67      REAL,INTENT(IN) :: pdqfi(ip1jmp1,llm,nqtot)
     68      REAL,INTENT(IN) :: pdhfi(ip1jmp1,llm)
     69      REAL,INTENT(IN) :: pdpfi(ip1jmp1)
     70c
     71      LOGICAL,INTENT(IN) :: leapf,forward ! not used
    6672c
    6773c
    … …  
    7177      REAL xpn(iim),xps(iim),tpn,tps
    7278      INTEGER j,k,iq,ij
    73       REAL qtestw, qtestt
    74       PARAMETER ( qtestw = 1.0e-15 )
    75       PARAMETER ( qtestt = 1.0e-40 )
     79      REAL,PARAMETER :: qtestw = 1.0e-15
     80      REAL,PARAMETER :: qtestt = 1.0e-40
    7681
    7782      REAL SSUM

  • trunk/LMDZ.COMMON/libf/dyn3dpar/advtrac_p.F90

    r1019 r1189  
    1616  USE Vampir
    1717  USE times
    18   USE infotrac
    19   USE control_mod
     18  USE infotrac, ONLY: nqtot, iadv
     19  USE control_mod, ONLY: iapp_tracvl, day_step, planet_type
    2020  IMPLICIT NONE
    2121  !
    … …  
    3434  !     Arguments
    3535  !-------------------------------------------------------------------
     36  INTEGER,INTENT(OUT) :: iapptrac
     37  REAL,INTENT(IN) :: pbaru(ip1jmp1,llm)
     38  REAL,INTENT(IN) :: pbarv(ip1jm,llm)
     39  REAL,INTENT(INOUT) :: q(ip1jmp1,llm,nqtot)
     40  REAL,INTENT(IN) :: masse(ip1jmp1,llm)
     41  REAL,INTENT(IN) :: p( ip1jmp1,llmp1 )
     42  REAL,INTENT(IN) :: teta(ip1jmp1,llm)
     43  REAL,INTENT(IN) :: pk(ip1jmp1,llm)
     44  REAL,INTENT(OUT) :: flxw(ip1jmp1,llm)
     45  !-------------------------------------------------------------------
    3646  !     Ajout PPM
    3747  !--------------------------------------------------------
    3848  REAL massebx(ip1jmp1,llm),masseby(ip1jm,llm)
    39   !--------------------------------------------------------
    40   INTEGER iapptrac
    41   REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
    42   REAL q(ip1jmp1,llm,nqtot),masse(ip1jmp1,llm)
    43   REAL p( ip1jmp1,llmp1 ),teta(ip1jmp1,llm)
    44   REAL pk(ip1jmp1,llm)
    45   REAL               :: flxw(ip1jmp1,llm)
    46 
    4749  !-------------------------------------------------------------
    4850  !     Variables locales

  • trunk/LMDZ.COMMON/libf/dyn3dpar/caldyn_p.F

    r1019 r1189  
    11!
    2 ! $Header$
    3 !
    4 c
    5 c
     2! $Id: $
     3!
    64#undef DEBUG_IO
    75c#define DEBUG_IO
    … …  
    1513      IMPLICIT NONE
    1614
    17 c=======================================================================
    18 c
    19 c  Auteur :  P. Le Van
    20 c
    21 c   Objet:
    22 c   ------
    23 c
    24 c   Calcul des tendances dynamiques.
    25 c
    26 c Modif 04/93 F.Forget
    27 c=======================================================================
    28 
    29 c-----------------------------------------------------------------------
    30 c   0. Declarations:
    31 c   ----------------
     15!=======================================================================
     16!
     17!  Auteur :  P. Le Van
     18!
     19!   Objet:
     20!   ------
     21!
     22!   Calcul des tendances dynamiques.
     23!
     24! Modif 04/93 F.Forget
     25!=======================================================================
     26
     27!-----------------------------------------------------------------------
     28!   0. Declarations:
     29!   ----------------
    3230
    3331#include "dimensions.h"
    … …  
    3735#include "comgeom.h"
    3836
    39 c   Arguments:
    40 c   ----------
    41 
    42       LOGICAL conser
    43 
    44       INTEGER itau
    45       REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm),teta(ip1jmp1,llm)
    46       REAL ps(ip1jmp1),phis(ip1jmp1)
    47       REAL pk(ip1jmp1,llm),pkf(ip1jmp1,llm)
    48       REAL tsurpk(ip1jmp1,llm)
     37!   Arguments:
     38!   ----------
     39
     40      LOGICAL,INTENT(IN) :: conser ! triggers printing some diagnostics
     41      INTEGER,INTENT(IN) :: itau ! time step index
     42      REAL,INTENT(IN) :: vcov(ip1jm,llm) ! covariant meridional wind
     43      REAL,INTENT(IN) :: ucov(ip1jmp1,llm) ! covariant zonal wind
     44      REAL,INTENT(IN) :: teta(ip1jmp1,llm) ! potential temperature
     45      REAL,INTENT(IN) :: ps(ip1jmp1) ! surface pressure
     46      REAL,INTENT(IN) :: phis(ip1jmp1) ! geopotential at the surface
     47      REAL,INTENT(IN) :: pk(ip1jmp1,llm) ! Exner at mid-layer
     48      REAL,INTENT(IN) :: pkf(ip1jmp1,llm) ! filtered Exner
     49      REAL,INTENT(IN) :: tsurpk(ip1jmp1,llm) ! cpp * temperature / pk
     50      REAL,INTENT(IN) :: phi(ip1jmp1,llm) ! geopotential
     51      REAL,INTENT(OUT) :: masse(ip1jmp1,llm) ! air mass
     52      REAL,INTENT(OUT) :: dv(ip1jm,llm) ! tendency on vcov
     53      REAL,INTENT(OUT) :: du(ip1jmp1,llm) ! tendency on ucov
     54      REAL,INTENT(OUT) :: dteta(ip1jmp1,llm) ! tenddency on teta
     55      REAL,INTENT(OUT) :: dp(ip1jmp1) ! tendency on ps
     56      REAL,INTENT(OUT) :: w(ip1jmp1,llm) ! vertical velocity
     57      REAL,INTENT(OUT) :: pbaru(ip1jmp1,llm) ! mass flux in the zonal direction
     58      REAL,INTENT(OUT) :: pbarv(ip1jm,llm) ! mass flux in the meridional direction
     59      REAL,INTENT(IN) :: time ! current time
     60
     61!   Local:
     62!   ------
     63
    4964      REAL,SAVE :: vcont(ip1jm,llm),ucont(ip1jmp1,llm)
    50       REAL phi(ip1jmp1,llm),masse(ip1jmp1,llm)
    51       REAL dv(ip1jm,llm),du(ip1jmp1,llm)
    52       REAL dteta(ip1jmp1,llm),dp(ip1jmp1)
    53       REAL w(ip1jmp1,llm)
    54       REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
    55       REAL time
    56 
    57 c   Local:
    58 c   ------
    59 
    6065      REAL,SAVE :: ang(ip1jmp1,llm)
    6166      REAL,SAVE :: p(ip1jmp1,llmp1)
    … …  
    7075      INTEGER   ij,l,ijb,ije,ierr
    7176
    72 c-----------------------------------------------------------------------
    73 c   Calcul des tendances dynamiques:
    74 c   --------------------------------
     77!-----------------------------------------------------------------------
     78!   Compute dynamical tendencies:
     79!--------------------------------
     80
     81      ! compute contravariant winds ucont() and vcont
    7582      CALL covcont_p  ( llm    , ucov    , vcov , ucont, vcont        )
     83      ! compute pressure p()
    7684      CALL pression_p ( ip1jmp1, ap      , bp   ,  ps  , p            )
    7785cym      CALL psextbar (   ps   , psexbarxy                          )
    7886c$OMP BARRIER
     87      ! compute mass in each atmospheric mesh: masse()
    7988      CALL massdair_p (    p   , masse                                )
     89      ! compute X and Y-averages of mass, massebx() and masseby()
    8090      CALL massbar_p  (   masse, massebx , masseby                    )
     91      ! compute XY-average of mass, massebxy()
    8192      call massbarxy_p(   masse, massebxy                             )
     93      ! compute mass fluxes pbaru() and pbarv()
    8294      CALL flumass_p  ( massebx, masseby , vcont, ucont ,pbaru, pbarv )
     95      ! compute dteta() , horizontal converging flux of theta
    8396      CALL dteta1_p   (   teta , pbaru   , pbarv, dteta               )
     97      ! compute convm(), horizontal converging flux of mass
    8498      CALL convmas1_p  (   pbaru, pbarv   , convm                      )
    8599c$OMP BARRIER     
    … …  
    108122      ijb=ij_begin
    109123      ije=ij_end
    110            
     124      ! compute pressure variation due to mass convergence
    111125      DO ij =ijb, ije
    112126         dp( ij ) = convm( ij,1 ) / airesurg( ij )
    … …  
    115129c$OMP BARRIER
    116130c$OMP FLUSH
     131
     132      ! compute vertical velocity w()
    117133      CALL vitvert_p ( convm  , w                                  )
     134      ! compute potential vorticity vorpot()
    118135      CALL tourpot_p ( vcov   , ucov  , massebxy  , vorpot         )
     136      ! compute rotation induced du() and dv()
    119137      CALL dudv1_p   ( vorpot , pbaru , pbarv     , du     , dv    )
    120138
    … …  
    129147c$OMP BARRIER
    130148#endif     
     149     
     150      ! compute kinetic energy ecin()
    131151      CALL enercin_p ( vcov   , ucov  , vcont     , ucont  , ecin  )
     152      ! compute Bernouilli function bern()
    132153      CALL bernoui_p ( ip1jmp1, llm   , phi       , ecin   , bern  )
     154      ! compute and add du() and dv() contributions from Bernouilli and pressure
    133155      CALL dudv2_p   ( tsurpk , pkf   , bern      , du     , dv    )
    134156
    … …  
    159181c$OMP END DO
    160182
     183      ! compute vertical advection contributions to du(), dv() and dteta()
    161184      CALL advect_new_p(ang,vcov,teta,w,massebx,masseby,du,dv,dteta)
    162185
    … …  
    190213      ENDIF
    191214
    192       RETURN
    193215      END

  • trunk/LMDZ.COMMON/libf/dyn3dpar/calfis_p.F

    r1129 r1189  
    3737      Use Write_field_p
    3838      USE Times
    39       USE infotrac
    40       USE control_mod
     39      USE infotrac, ONLY: nqtot, niadv, tname
     40      USE control_mod, ONLY: planet_type, nsplit_phys
    4141      USE cpdet_mod, only: tpot2t_p, t2tpot_p
    4242
    … …  
    116116c    Arguments :
    117117c    -----------
    118       LOGICAL  lafin
    119 !      REAL heure
    120       REAL, intent(in):: jD_cur, jH_cur
    121       REAL pvcov(iip1,jjm,llm)
    122       REAL pucov(iip1,jjp1,llm)
    123       REAL pteta(iip1,jjp1,llm)
    124       REAL pmasse(iip1,jjp1,llm)
    125       REAL pq(iip1,jjp1,llm,nqtot)
    126       REAL pphis(iip1,jjp1)
    127       REAL pphi(iip1,jjp1,llm)
    128 
    129       REAL pdvcov(iip1,jjm,llm)
    130       REAL pducov(iip1,jjp1,llm)
    131       REAL pdteta(iip1,jjp1,llm)
     118      LOGICAL,INTENT(IN) ::  lafin ! .true. for the very last call to physics
     119      REAL,INTENT(IN) :: jD_cur, jH_cur
     120      REAL,INTENT(IN) :: pvcov(iip1,jjm,llm) ! covariant meridional velocity
     121      REAL,INTENT(IN) :: pucov(iip1,jjp1,llm) ! covariant zonal velocity
     122      REAL,INTENT(IN) :: pteta(iip1,jjp1,llm) ! potential temperature
     123      REAL,INTENT(IN) :: pmasse(iip1,jjp1,llm) ! mass in each cell ! not used
     124      REAL,INTENT(IN) :: pq(iip1,jjp1,llm,nqtot) ! tracers
     125      REAL,INTENT(IN) :: pphis(iip1,jjp1) ! surface geopotential
     126      REAL,INTENT(IN) :: pphi(iip1,jjp1,llm) ! geopotential
     127
     128      REAL,INTENT(IN) :: pdvcov(iip1,jjm,llm) ! dynamical tendency on vcov
     129      REAL,INTENT(IN) :: pducov(iip1,jjp1,llm) ! dynamical tendency on ucov
     130      REAL,INTENT(IN) :: pdteta(iip1,jjp1,llm) ! dynamical tendency on teta
    132131! commentaire SL: pdq ne sert que pour le calcul de pcvgq,
    133132! qui lui meme ne sert a rien dans la routine telle qu'elle est
    134133! ecrite, et que j'ai donc commente....
    135       REAL pdq(iip1,jjp1,llm,nqtot)
    136       REAL flxw(iip1,jjp1,llm)  ! Flux de masse verticale sur la grille dynamique
    137 
    138       REAL pps(iip1,jjp1)
    139       REAL pp(iip1,jjp1,llmp1)
    140       REAL ppk(iip1,jjp1,llm)
    141 
    142       REAL pdvfi(iip1,jjm,llm)
    143       REAL pdufi(iip1,jjp1,llm)
    144       REAL pdhfi(iip1,jjp1,llm)
    145       REAL pdqfi(iip1,jjp1,llm,nqtot)
    146       REAL pdpsfi(iip1,jjp1)
     134      REAL,INTENT(IN) :: pdq(iip1,jjp1,llm,nqtot) ! dynamical tendency on tracers
     135      ! NB: pdq is only used to compute pcvgq which is in fact not used...
     136
     137      REAL,INTENT(IN) :: pps(iip1,jjp1) ! surface pressure (Pa)
     138      REAL,INTENT(IN) :: pp(iip1,jjp1,llmp1) ! pressure at mesh interfaces (Pa)
     139      REAL,INTENT(IN) :: ppk(iip1,jjp1,llm) ! Exner at mid-layer
     140      REAL,INTENT(IN) :: flxw(iip1,jjp1,llm)  ! Vertical mass flux on dynamics grid
     141
     142      ! tendencies (in */s) from the physics
     143      REAL,INTENT(OUT) :: pdvfi(iip1,jjm,llm) ! tendency on covariant meridional wind
     144      REAL,INTENT(OUT) :: pdufi(iip1,jjp1,llm) ! tendency on covariant zonal wind
     145      REAL,INTENT(OUT) :: pdhfi(iip1,jjp1,llm) ! tendency on potential temperature (K/s)
     146      REAL,INTENT(OUT) :: pdqfi(iip1,jjp1,llm,nqtot) ! tendency on tracers
     147      REAL,INTENT(OUT) :: pdpsfi(iip1,jjp1) ! tendency on surface pressure (Pa/s)
    147148
    148149#ifdef CPP_PHYS
    … …  
    225226
    226227cIM diagnostique PVteta, Amip2
    227       INTEGER ntetaSTD
    228       PARAMETER(ntetaSTD=3)
    229       REAL rtetaSTD(ntetaSTD)
    230       DATA rtetaSTD/350., 380., 405./ ! Earth-specific values, beware !!
     228      INTEGER,PARAMETER :: ntetaSTD=3
     229      REAL,SAVE :: rtetaSTD(ntetaSTD)=(/350.,380.,405./) ! Earth-specific, beware !!
    231230      REAL PVteta(klon,ntetaSTD)
    232231     
    … …  
    234233      REAL SSUM
    235234
    236       LOGICAL firstcal, debut
    237       DATA firstcal/.true./
    238       SAVE firstcal,debut
     235      LOGICAL,SAVE :: firstcal=.true., debut=.true.
    239236c$OMP THREADPRIVATE(firstcal,debut)
    240237     

  • trunk/LMDZ.COMMON/libf/dyn3dpar/conf_gcm.F

    r1056 r1189  
    202202      CALL getin('starttime',starttime)
    203203
     204      ! Mars: time of start for run in "start.nc" (when there are multiple time
     205      !       steps stored in the file)
     206      timestart=-9999 ! default value; if <0, use last stored time
     207      call getin("timestart",timestart)
     208     
    204209!Config  Key  = less1day
    205210!Config  Desc = Possibilite d'integrer moins d'un jour

  • trunk/LMDZ.COMMON/libf/dyn3dpar/dissip_p.F

    r1019 r1189  
     1!
     2! $Id: $
     3!
    14      SUBROUTINE dissip_p( vcov,ucov,teta,p, dv,du,dh )
    25c
    … …  
    3437c   ----------
    3538
    36       REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm),teta(ip1jmp1,llm)
    37       REAL  p( ip1jmp1,llmp1 )
    38       REAL dv(ip1jm,llm),du(ip1jmp1,llm),dh(ip1jmp1,llm)
     39      REAL,INTENT(IN) :: vcov(ip1jm,llm) ! covariant meridional wind
     40      REAL,INTENT(IN) :: ucov(ip1jmp1,llm) ! covariant zonal wind
     41      REAL,INTENT(IN) :: teta(ip1jmp1,llm) ! potentail temperature
     42      REAL,INTENT(IN) :: p(ip1jmp1,llmp1) ! pressure
     43      ! tendencies (.../s) on covariant winds and potential temperature
     44      REAL,INTENT(OUT) :: dv(ip1jm,llm)
     45      REAL,INTENT(OUT) :: du(ip1jmp1,llm)
     46      REAL,INTENT(OUT) :: dh(ip1jmp1,llm)
    3947
    4048c   Local:

  • trunk/LMDZ.COMMON/libf/dyn3dpar/leapfrog_p.F

    r1107 r1189  
    8383#include "academic.h"
    8484     
    85       real zqmin,zqmax
    86 
    87 c   variables dynamiques
    88       REAL :: vcov(ip1jm,llm),ucov(ip1jmp1,llm) ! vents covariants
    89       REAL :: teta(ip1jmp1,llm)                 ! temperature potentielle
    90       REAL :: q(ip1jmp1,llm,nqtot)              ! champs advectes
    91       REAL :: ps(ip1jmp1)                       ! pression  au sol
    92       REAL,SAVE :: p (ip1jmp1,llmp1  )               ! pression aux interfac.des couches
    93       REAL,SAVE :: pks(ip1jmp1)                      ! exner au  sol
    94       REAL,SAVE :: pk(ip1jmp1,llm)                   ! exner au milieu des couches
    95       REAL,SAVE :: pkf(ip1jmp1,llm)                  ! exner filt.au milieu des couches
    96       REAL :: masse(ip1jmp1,llm)                ! masse d'air
    97       REAL :: phis(ip1jmp1)                     ! geopotentiel au sol
    98       REAL,SAVE :: phi(ip1jmp1,llm)                  ! geopotentiel
    99       REAL,SAVE :: w(ip1jmp1,llm)                    ! vitesse verticale
     85      REAL,INTENT(IN) :: time_0 ! not used
     86
     87c   dynamical variables:
     88      REAL,INTENT(INOUT) :: ucov(ip1jmp1,llm)    ! zonal covariant wind
     89      REAL,INTENT(INOUT) :: vcov(ip1jm,llm)      ! meridional covariant wind
     90      REAL,INTENT(INOUT) :: teta(ip1jmp1,llm)    ! potential temperature
     91      REAL,INTENT(INOUT) :: ps(ip1jmp1)          ! surface pressure (Pa)
     92      REAL,INTENT(INOUT) :: masse(ip1jmp1,llm)   ! air mass
     93      REAL,INTENT(INOUT) :: phis(ip1jmp1)        ! geopotentiat at the surface
     94      REAL,INTENT(INOUT) :: q(ip1jmp1,llm,nqtot) ! advected tracers
     95     
     96      REAL,SAVE :: p (ip1jmp1,llmp1  )       ! interlayer pressure
     97      REAL,SAVE :: pks(ip1jmp1)              ! exner at the surface
     98      REAL,SAVE :: pk(ip1jmp1,llm)           ! exner at mid-layer
     99      REAL,SAVE :: pkf(ip1jmp1,llm)          ! filtered exner at mid-layer
     100      REAL,SAVE :: phi(ip1jmp1,llm)          ! geopotential
     101      REAL,SAVE :: w(ip1jmp1,llm)            ! vertical velocity
    100102! ADAPTATION GCM POUR CP(T)
    101103      REAL,SAVE :: temp(ip1jmp1,llm)                 ! temperature 
    102104      REAL,SAVE :: tsurpk(ip1jmp1,llm)               ! cpp*T/pk 
     105
     106      real zqmin,zqmax
    103107
    104108c variables dynamiques intermediaire pour le transport
    … …  
    144148      REAL       time
    145149
    146       REAL  SSUM
    147       REAL time_0
     150      REAL  SSUM
    148151!      REAL,SAVE :: finvmaold(ip1jmp1,llm)
    149152
    … …  
    721724c   -------------------------------------------------------------
    722725
    723 !      IF( forward. OR . leapf )  THEN
    724       IF((.not.forward).OR. leapf )  THEN
    725         ! Ehouarn: gather mass fluxes during backward Matsuno or LF step
    726 cc$OMP PARALLEL DEFAULT(SHARED)
    727 c
     726      IF( forward. OR . leapf )  THEN
     727! Ehouarn: NB: fields sent to advtrac are those at the beginning of the time step
    728728         CALL advtrac_p( pbaru,pbarv,
    729729     *             p,  masse,q,iapptrac, teta,
    … …  
    737737
    738738      ENDIF ! of IF( forward. OR . leapf )
    739 cc$OMP END PARALLEL
    740739
    741740c-----------------------------------------------------------------------
    … …  
    10371036     $                  ucov, vcov, teta , q   ,ps ,
    10381037     $                 dufi, dvfi, dtetafi , dqfi ,dpfi  )
     1038          ! since addfi updates ps(), also update p(), masse() and pk()
     1039          CALL pression_p(ip1jmp1,ap,bp,ps,p)
     1040c$OMP BARRIER
     1041          CALL massdair_p(p,masse)
     1042c$OMP BARRIER
     1043          if (pressure_exner) then
     1044            CALL exner_hyb_p(ip1jmp1,ps,p,alpha,beta,pks,pk,pkf)
     1045          else
     1046            CALL exner_milieu_p(ip1jmp1,ps,p,beta,pks,pk,pkf)
     1047          endif
     1048c$OMP BARRIER
    10391049
    10401050c      Couche superieure :
    … …  
    11951205          CALL exner_milieu_p( ip1jmp1, ps, p, beta, pks, pk, pkf )
    11961206        endif
     1207c$OMP BARRIER
     1208        CALL massdair_p(p,masse)
    11971209c$OMP BARRIER
    11981210
Note: See TracChangeset for help on using the changeset viewer.