Changeset 1682 for trunk/LMDZ.COMMON/libf/dyn3dpar

Timestamp:

Mar 31, 2017, 11:31:38 AM (8 years ago)

Author:

emillour

Message:

All GCMs: set things up to enable pluging physics with dynamico

• dyn3d
• gcm.F90 : move I/O initialization (dates) to be done before physics

initialization

• dyn3dpar
• gcm.F : move I/O initialization (dates) to be done before physics

initialization

• dynphy_lonlat:
• inigeomphy_mod.F90 : add ind_cell_glo computation and transfer

to init_geometry

• phy_common:
• geometry_mod.F90 : add ind_cell_glo module variable to store global

column index

• print_control_mod.F90 : make initialization occur via init_print_control_mod

to avoid circular module dependencies

• init_print_control_mod.F90 : added to initialize print_control_mod module

variables

• mod_phys_lmdz_mpi_data.F90 : use print_control_mod (rather than iniprint.h)
• mod_phys_lmdz_para.F90 : use print_control_mod (rather than iniprint.h)
• mod_phys_lmdz_omp_data.F90 : add is_omp_master (alias of is_omp_root) module

variable and use print_control_mod (rather than
iniprint.h)

• physics_distribution_mod.F90 : add call to init_dimphy in

init_physics_distribution

• xios_writefield.F90 : generic routine to output field with XIOS (for debug)

• misc:
• handle_err_m.F90 : call abort_physic, rather than abort_gcm
• wxios.F90 : updates to enable unstructured grids

set module variable g_ctx_name to "LMDZ"
wxios_init(): remove call to wxios_context_init
wxios_context_init(): call xios_context_initialize with COMM_LMDZ_PHY
add routine wxios_set_context() to get handle and set context to XIOS
wxios_domain_param(): change arguments and generate the domain in-place
add wxios_domain_param_unstructured(): generate domain for unstructured case

NB: access is via "domain group" (whereas it is via "domain" in

wxios_domain_param)

• dynphy_lonlat/phy[std|mars|venus|titan]:
• iniphysiq_mod.F90 : Remove call to init_dimphy (which is now done in

phy_common/physics_distribution_mod.F90)

EM

File:

• trunk/LMDZ.COMMON/libf/dyn3dpar/gcm.F (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/dyn3dpar/gcm.F

@@ -441 +441 @@
 
 c-----------------------------------------------------------------------
-c   Initialisation de la physique :
-c   -------------------------------
-
-      IF ((iflag_phys==1).or.(iflag_phys>=100)) THEN
-!         latfi(1)=rlatu(1)
-!         lonfi(1)=0.
-!         zcufi(1) = cu(1)
-!         zcvfi(1) = cv(1)
-!         DO j=2,jjm
-!            DO i=1,iim
-!               latfi((j-2)*iim+1+i)= rlatu(j)
-!               lonfi((j-2)*iim+1+i)= rlonv(i)
-!               zcufi((j-2)*iim+1+i) = cu((j-1)*iip1+i)
-!               zcvfi((j-2)*iim+1+i) = cv((j-1)*iip1+i)
-!            ENDDO
-!         ENDDO
-!         latfi(ngridmx)= rlatu(jjp1)
-!         lonfi(ngridmx)= 0.
-!         zcufi(ngridmx) = cu(ip1jm+1)
-!         zcvfi(ngridmx) = cv(ip1jm-iim)
-
-         ! build airefi(), mesh area on physics grid
-!         CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,aire,airefi)
-         ! Poles are single points on physics grid
-!         airefi(1)=airefi(1)*iim
-!         airefi(ngridmx)=airefi(ngridmx)*iim
-
-! Physics
-#ifdef CPP_PHYS
-!         CALL iniphysiq(ngridmx,llm,daysec,day_ini,dtphys/nsplit_phys,
-!     &                latfi,lonfi,airefi,zcufi,zcvfi,rad,g,r,cpp,
-!     &                iflag_phys)
-         CALL iniphysiq(iim,jjm,llm,
-     &                distrib_phys(mpi_rank),comm_lmdz,
-     &                daysec,day_ini,dtphys/nsplit_phys,
-     &                rlatu,rlatv,rlonu,rlonv,aire,cu,cv,rad,g,r,cpp,
-     &                iflag_phys)
-#endif
-!         call_iniphys=.false.
-      ENDIF ! of IF (call_iniphys.and.(iflag_phys==1.or.iflag_phys>=100))
-
-
-c-----------------------------------------------------------------------
 c   Initialisation des I/O :
 c   ------------------------
@@ -519 +476 @@
 #endif
 
+
+c-----------------------------------------------------------------------
+c   Initialisation de la physique :
+c   -------------------------------
+
+      IF ((iflag_phys==1).or.(iflag_phys>=100)) THEN
+! Physics
+#ifdef CPP_PHYS
+         CALL iniphysiq(iim,jjm,llm,
+     &                distrib_phys(mpi_rank),comm_lmdz,
+     &                daysec,day_ini,dtphys/nsplit_phys,
+     &                rlatu,rlatv,rlonu,rlonv,aire,cu,cv,rad,g,r,cpp,
+     &                iflag_phys)
+#endif
+!         call_iniphys=.false.
+      ENDIF ! of IF (call_iniphys.and.(iflag_phys==1.or.iflag_phys>=100))
+
+
       if (planet_type=="mars") then
          ! For Mars we transmit day_ini