Changeset 1987 for LMDZ5/trunk/libf/dyn3dmem/pression_loc.F

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:

• LMDZ5/trunk/libf/dyn3dmem/pression_loc.F (modified) (2 diffs)

LMDZ5/trunk/libf/dyn3dmem/pression_loc.F

@@ -1 +1 @@
       SUBROUTINE pression_loc( ngrid, ap, bp, ps, p )
-      USE parallel_lmdz
+      USE parallel_lmdz, ONLY: ij_begin, ij_end, ijb_u, ije_u,
+     &                         pole_nord, pole_sud, omp_chunk
 c
 
@@ -16 +17 @@
 #include "paramet.h"
 c
-      INTEGER ngrid
+      INTEGER,INTENT(IN) :: ngrid ! not used
       INTEGER l,ij
  
-      REAL ap( llmp1 ), bp( llmp1 ), ps( ijb_u:ije_u )
-      REAL p( ijb_u:ije_u,llmp1 ) 
+      REAL,INTENT(IN) :: ap( llmp1 ), bp( llmp1 ), ps( ijb_u:ije_u )
+      REAL,INTENT(OUT) :: p( ijb_u:ije_u,llmp1 ) 
       
       INTEGER ijb,ije