Changeset 2418 for LMDZ5/trunk/libf/dynphy_lonlat/calfis_loc.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_loc.F (modified) (4 diffs)

LMDZ5/trunk/libf/dynphy_lonlat/calfis_loc.F

@@ -47 +47 @@
       USE infotrac, ONLY: nqtot, niadv, tname
       USE control_mod, ONLY: planet_type, nsplit_phys
+#ifdef CPP_PHYS
+      USE callphysiq_mod, ONLY: call_physiq
+#endif 
 
 #ifdef CPP_PARA
@@ -138 +141 @@
       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
+      REAL,INTENT(IN) :: flxw(iip1,jjb_u:jje_u,llm) ! Vertical mass flux on lower mesh interfaces (kg/s) (on llm because flxw(:,:,llm+1)=0)
 
       ! tendencies (in */s) from the physics
@@ -240 +243 @@
       INTEGER :: jjb,jje
 
-      LOGICAL tracerdyn  
 c
 c-----------------------------------------------------------------------
@@ -709 +711 @@
          lafin_split=lafin.and.isplit==nsplit_phys
 
-      if (planet_type=="earth") then
-
-      CALL physiq (klon,
-     .             llm,
-     .             debut_split,
-     .             lafin_split,
-     .             jD_cur,
-     .             jH_cur_split,
-     .             zdt_split,
-     .             zplev_omp,
-     .             zplay_omp,
-     .             zphi_omp,
-     .             zphis_omp,
-     .             presnivs_omp,
-     .             zufi_omp,
-     .             zvfi_omp,
-     .             zrfi_omp,
-     .             ztfi_omp,
-     .             zqfi_omp,
-     .             flxwfi_omp,
-     .             zdufi_omp,
-     .             zdvfi_omp,
-     .             zdtfi_omp,
-     .             zdqfi_omp,
-     .             zdpsrf_omp,
-     .             pducov)
-
-      else if ( planet_type=="generic" ) then
-
-      CALL physiq (klon,     !! 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_p 
-     .             jH_cur_split,   !! ptime "fraction of day"
-     .             zdt_split,      !! ptimestep
-     .             zplev_omp,  !! pplev
-     .             zplay_omp,  !! pplay
-     .             zphi_omp,   !! pphi
-     .             zufi_omp,   !! pu
-     .             zvfi_omp,   !! pv
-     .             ztfi_omp,   !! pt
-     .             zqfi_omp,   !! pq
-     .             flxwfi_omp, !! pw !! or 0. anyway this is for diagnostic. not used in physiq.
-     .             zdufi_omp,  !! pdu
-     .             zdvfi_omp,  !! pdv
-     .             zdtfi_omp,  !! pdt
-     .             zdqfi_omp,  !! pdq
-     .             zdpsrf_omp, !! pdpsrf
-     .             tracerdyn)      !! tracerdyn <-- utilite ???
-
-      endif ! of if (planet_type=="earth")
+        CALL call_physiq(klon,llm,nqtot,tname,
+     &                   debut_split,lafin_split,
+     &                   jD_cur,jH_cur_split,zdt_split,
+     &                   zplev_omp,zplay_omp,
+     &                   zphi_omp,zphis_omp,
+     &                   presnivs_omp,
+     &                   zufi_omp,zvfi_omp,zrfi_omp,ztfi_omp,zqfi_omp,
+     &                   flxwfi_omp,pducov,
+     &                   zdufi_omp,zdvfi_omp,zdtfi_omp,zdqfi_omp,
+     &                   zdpsrf_omp)