Changeset 1189 for trunk/LMDZ.COMMON/libf/dyn3d/advtrac.F90

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

File:

• trunk/LMDZ.COMMON/libf/dyn3d/advtrac.F90 (modified) (2 diffs)

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

@@ -9 +9 @@
   !            M.A Filiberti (04/2002)
   !
-  USE infotrac
-  USE control_mod
+  USE infotrac, ONLY: nqtot, iadv
+  USE control_mod, ONLY: iapp_tracvl, day_step
 
 
@@ -30 +30 @@
   !     Arguments
   !-------------------------------------------------------------------
+  INTEGER,INTENT(OUT) :: iapptrac
+  REAL,INTENT(IN) :: pbaru(ip1jmp1,llm)
+  REAL,INTENT(IN) :: pbarv(ip1jm,llm)
+  REAL,INTENT(INOUT) :: q(ip1jmp1,llm,nqtot)
+  REAL,INTENT(IN) :: masse(ip1jmp1,llm)
+  REAL,INTENT(IN) :: p( ip1jmp1,llmp1 )
+  REAL,INTENT(IN) :: teta(ip1jmp1,llm)
+  REAL,INTENT(IN) :: pk(ip1jmp1,llm)
+  REAL,INTENT(OUT) :: flxw(ip1jmp1,llm)
+  !-------------------------------------------------------------------
   !     Ajout PPM
   !--------------------------------------------------------
   REAL massebx(ip1jmp1,llm),masseby(ip1jm,llm)
-  !--------------------------------------------------------
-  INTEGER iapptrac
-  REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
-  REAL q(ip1jmp1,llm,nqtot),masse(ip1jmp1,llm)
-  REAL p( ip1jmp1,llmp1 ),teta(ip1jmp1,llm)
-  REAL pk(ip1jmp1,llm)
-  REAL flxw(ip1jmp1,llm)
-
   !-------------------------------------------------------------
   !     Variables locales