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

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/dyn3dmem/calfis_loc.F (modified) (4 diffs)

LMDZ5/trunk/libf/dyn3dmem/calfis_loc.F

@@ -41 +41 @@
       Use Write_field_p
       USE Times
-      USE infotrac
-      USE control_mod
+      USE infotrac, ONLY: nqtot, niadv, tname
+      USE control_mod, ONLY: planet_type, nsplit_phys
 
       IMPLICIT NONE
@@ -115 +115 @@
 c    Arguments :
 c    -----------
-      LOGICAL  lafin
-!      REAL heure
-      REAL, intent(in):: jD_cur, jH_cur
-      REAL pvcov(iip1,jjb_v:jje_v,llm)
-      REAL pucov(iip1,jjb_u:jje_u,llm)
-      REAL pteta(iip1,jjb_u:jje_u,llm)
-      REAL pmasse(iip1,jjb_u:jje_u,llm)
-      REAL pq(iip1,jjb_u:jje_u,llm,nqtot)
-      REAL pphis(iip1,jjb_u:jje_u)
-      REAL pphi(iip1,jjb_u:jje_u,llm)
-c
-      REAL pdvcov(iip1,jjb_v:jje_v,llm)
-      REAL pducov(iip1,jjb_u:jje_u,llm)
-      REAL pdteta(iip1,jjb_u:jje_u,llm)
-      REAL pdq(iip1,jjb_u:jje_u,llm,nqtot)
-      REAL flxw(iip1,jjb_u:jje_u,llm)  ! Flux de masse verticale sur la grille dynamique
-c
-      REAL pps(iip1,jjb_u:jje_u)
-      REAL pp(iip1,jjb_u:jje_u,llmp1)
-      REAL ppk(iip1,jjb_u:jje_u,llm)
-c
-      REAL pdvfi(iip1,jjb_v:jje_v,llm)
-      REAL pdufi(iip1,jjb_u:jje_u,llm)
-      REAL pdhfi(iip1,jjb_u:jje_u,llm)
-      REAL pdqfi(iip1,jjb_u:jje_u,llm,nqtot)
-      REAL pdpsfi(iip1,jjb_u:jje_u)
-
-      INTEGER        longcles
-      PARAMETER    ( longcles = 20 )
-      REAL clesphy0( longcles )
+      LOGICAL,INTENT(IN) ::  lafin ! .true. for the very last call to physics
+      REAL,INTENT(IN):: jD_cur, jH_cur
+      REAL,INTENT(IN):: pvcov(iip1,jjb_v:jje_v,llm) ! covariant meridional velocity
+      REAL,INTENT(IN):: pucov(iip1,jjb_u:jje_u,llm) ! covariant zonal velocity
+      REAL,INTENT(IN):: pteta(iip1,jjb_u:jje_u,llm) ! potential temperature
+      REAL,INTENT(IN):: pmasse(iip1,jjb_u:jje_u,llm) ! mass in each cell ! not used
+      REAL,INTENT(IN):: pq(iip1,jjb_u:jje_u,llm,nqtot) ! tracers
+      REAL,INTENT(IN):: pphis(iip1,jjb_u:jje_u) ! surface geopotential
+      REAL,INTENT(IN):: pphi(iip1,jjb_u:jje_u,llm) ! geopotential
+
+      REAL,INTENT(IN) :: pdvcov(iip1,jjb_v:jje_v,llm) ! dynamical tendency on vcov ! not used
+      REAL,INTENT(IN) :: pducov(iip1,jjb_u:jje_u,llm) ! dynamical tendency on ucov
+      REAL,INTENT(IN) :: pdteta(iip1,jjb_u:jje_u,llm) ! dynamical tendency on teta ! not used
+      REAL,INTENT(IN) :: pdq(iip1,jjb_u:jje_u,llm,nqtot) ! dynamical tendency on tracers ! not used
+
+      REAL,INTENT(IN) :: pps(iip1,jjb_u:jje_u) ! surface pressure (Pa)
+      REAL,INTENT(IN) :: pp(iip1,jjb_u:jje_u,llmp1) ! pressure at mesh interfaces (Pa)
+      REAL,INTENT(IN) :: ppk(iip1,jjb_u:jje_u,llm) ! Exner at mid-layer
+      REAL,INTENT(IN) :: flxw(iip1,jjb_u:jje_u,llm)  ! Vertical mass flux on dynamics grid
+
+      ! tendencies (in */s) from the physics
+      REAL,INTENT(OUT) :: pdvfi(iip1,jjb_v:jje_v,llm) ! tendency on covariant meridional wind
+      REAL,INTENT(OUT) :: pdufi(iip1,jjb_u:jje_u,llm) ! tendency on covariant zonal wind
+      REAL,INTENT(OUT) :: pdhfi(iip1,jjb_u:jje_u,llm) ! tendency on potential temperature (K/s)
+      REAL,INTENT(OUT) :: pdqfi(iip1,jjb_u:jje_u,llm,nqtot) ! tendency on tracers
+      REAL,INTENT(OUT) :: pdpsfi(iip1,jjb_u:jje_u) ! tendency on surface pressure (Pa/s)
+
+      INTEGER,PARAMETER :: longcles = 20
+      REAL,INTENT(IN) :: clesphy0( longcles ) ! unused
 
 
@@ -221 +220 @@
 c
 cIM diagnostique PVteta, Amip2
-      INTEGER ntetaSTD
-      PARAMETER(ntetaSTD=3)
-      REAL rtetaSTD(ntetaSTD)
-      DATA rtetaSTD/350., 380., 405./ ! Earth-specific values, beware !!
+      INTEGER,PARAMETER :: ntetaSTD=3
+      REAL,SAVE :: rtetaSTD(ntetaSTD)=(/350.,380.,405./) ! Earth-specific, beware !!
       REAL PVteta(klon,ntetaSTD)
       
@@ -230 +227 @@
       REAL SSUM
 
-      LOGICAL firstcal, debut
-      DATA firstcal/.true./
-      SAVE firstcal,debut
+      LOGICAL,SAVE :: firstcal=.true., debut=.true.
 c$OMP THREADPRIVATE(firstcal,debut)