Changeset 1987 for LMDZ5/trunk/libf/dyn3d/advtrac.F90

Timestamp:

Feb 24, 2014, 4:05:47 PM (10 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:

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

LMDZ5/trunk/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