Changeset 2418 for LMDZ5/trunk/libf/dynphy_lonlat/calfis.F

Timestamp:

Jan 3, 2016, 11:16:34 AM (8 years ago)

Author:

Ehouarn Millour

Message:

Improving the physics/dynamics interface:

• added module callphysiq_mod.F90 in dynphy_lonlat/phy* which contains the routine "call_physiq" which is called by calfis* and calls the physics. This way different "physiq" routine from different physics packages may be called: The calfis* routines now exposes all available fields that might be transmitted to physiq but which is actually send (ie: expected/needed by physiq) is decided in call_physiq.
• turned "physiq.F90" into module "physiq_mod.F90" for better control of "physiq" arguments. Extracted embeded "gr_fi_ecrit" as self-standing routine (but note that this routine actually only works in serial mode).

EM

File:

• LMDZ5/trunk/libf/dynphy_lonlat/calfis.F (modified) (4 diffs)

LMDZ5/trunk/libf/dynphy_lonlat/calfis.F

@@ -31 +31 @@
       USE infotrac, ONLY: nqtot, niadv, tname
       USE control_mod, ONLY: planet_type, nsplit_phys
- 
+#ifdef CPP_PHYS
+      USE callphysiq_mod, ONLY: call_physiq
+#endif 
 
       IMPLICIT NONE
@@ -121 +123 @@
       REAL,INTENT(IN) :: pp(iip1,jjp1,llmp1) ! pressure at mesh interfaces (Pa)
       REAL,INTENT(IN) :: ppk(iip1,jjp1,llm) ! Exner at mid-layer
-      REAL,INTENT(IN) :: flxw(iip1,jjp1,llm)  ! Vertical mass flux on dynamics grid
+      REAL,INTENT(IN) :: flxw(iip1,jjp1,llm) ! Vertical mass flux on lower mesh interfaces (kg/s) (on llm because flxw(:,:,llm+1)=0)
 
       ! tendencies (in */s) from the physics
@@ -167 +169 @@
       LOGICAL,SAVE :: firstcal=.true., debut=.true.
 !      REAL rdayvrai
-
-      LOGICAL tracerdyn
 
 c
@@ -468 +468 @@
          lafin_split=lafin.and.isplit==nsplit_phys
 
-      if (planet_type=="earth") then
-
-         CALL physiq (ngridmx,
-     .             llm,
-     .             debut_split,
-     .             lafin_split,
-     .             jD_cur,
-     .             jH_cur_split,
-     .             zdt_split,
-     .             zplev,
-     .             zplay,
-     .             zphi,
-     .             zphis,
-     .             presnivs,
-     .             zufi,
-     .             zvfi, zrfi,
-     .             ztfi,
-     .             zqfi,
-     .             flxwfi,
-     .             zdufi,
-     .             zdvfi,
-     .             zdtfi,
-     .             zdqfi,
-     .             zdpsrf,
-     .             pducov)
-
-      else if ( planet_type=="generic" ) then
-
-         CALL physiq (ngridmx,     !! ngrid
-     .             llm,            !! nlayer
-     .             nqtot,          !! nq
-     .             tname,          !! tracer names from dynamical core (given in infotrac)
-     .             debut_split,    !! firstcall 
-     .             lafin_split,    !! lastcall
-     .             jD_cur,         !! pday. see leapfrog
-     .             jH_cur_split,   !! ptime "fraction of day"
-     .             zdt_split,      !! ptimestep
-     .             zplev,          !! pplev
-     .             zplay,          !! pplay
-     .             zphi,           !! pphi
-     .             zufi,           !! pu
-     .             zvfi,           !! pv
-     .             ztfi,           !! pt
-     .             zqfi,           !! pq
-     .             flxwfi,         !! pw !! or 0. anyway this is for diagnostic. not used in physiq.
-     .             zdufi,          !! pdu
-     .             zdvfi,          !! pdv
-     .             zdtfi,          !! pdt
-     .             zdqfi,          !! pdq
-     .             zdpsrf,         !! pdpsrf
-     .             tracerdyn)      !! tracerdyn <-- utilite ???
-
-      endif ! of if (planet_type=="earth")
+        CALL call_physiq(ngridmx,llm,nqtot,tname,
+     &                   debut_split,lafin_split,
+     &                   jD_cur,jH_cur_split,zdt_split,
+     &                   zplev,zplay,
+     &                   zphi,zphis,
+     &                   presnivs,
+     &                   zufi,zvfi,zrfi,ztfi,zqfi,
+     &                   flxwfi,pducov,
+     &                   zdufi,zdvfi,zdtfi,zdqfi,zdpsrf)
+                             
+!      if (planet_type=="earth") then
+!
+!         CALL physiq (ngridmx,
+!     .             llm,
+!     .             debut_split,
+!     .             lafin_split,
+!     .             jD_cur,
+!     .             jH_cur_split,
+!     .             zdt_split,
+!     .             zplev,
+!     .             zplay,
+!     .             zphi,
+!     .             zphis,
+!     .             presnivs,
+!     .             zufi,
+!     .             zvfi, zrfi,
+!     .             ztfi,
+!     .             zqfi,
+!     .             flxwfi,
+!     .             zdufi,
+!     .             zdvfi,
+!     .             zdtfi,
+!     .             zdqfi,
+!     .             zdpsrf,
+!     .             pducov)
+!
+!      else if ( planet_type=="generic" ) then
+!
+!         CALL physiq (ngridmx,     !! ngrid
+!     .             llm,            !! nlayer
+!     .             nqtot,          !! nq
+!     .             tname,          !! tracer names from dynamical core (given in infotrac)
+!     .             debut_split,    !! firstcall 
+!     .             lafin_split,    !! lastcall
+!     .             jD_cur,         !! pday. see leapfrog
+!     .             jH_cur_split,   !! ptime "fraction of day"
+!     .             zdt_split,      !! ptimestep
+!     .             zplev,          !! pplev
+!     .             zplay,          !! pplay
+!     .             zphi,           !! pphi
+!     .             zufi,           !! pu
+!     .             zvfi,           !! pv
+!     .             ztfi,           !! pt
+!     .             zqfi,           !! pq
+!     .             flxwfi,         !! pw !! or 0. anyway this is for diagnostic. not used in physiq.
+!     .             zdufi,          !! pdu
+!     .             zdvfi,          !! pdv
+!     .             zdtfi,          !! pdt
+!     .             zdqfi,          !! pdq
+!     .             zdpsrf,         !! pdpsrf
+!     .             tracerdyn)      !! tracerdyn <-- utilite ???
+!
+!      endif ! of if (planet_type=="earth")
 
          zufi(:,:)=zufi(:,:)+zdufi(:,:)*zdt_split