Changeset 2454 for trunk/MESOSCALE

Timestamp:

Feb 16, 2021, 12:50:58 AM (4 years ago)

Author:

aslmd

Message:

MESOSCALE Mars: getin_p introduced in physics in commit r2304 so have to apply to MESOSCALE Mars the same fix as MESOSCALE generic in r2267 (fixing ticket #37) that is, a fix that initializes MPI communicator so that getin_p works, otherwise parameters from .def are not loaded

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_mars_lmd_new/iniphysiq_mod.F (modified) (2 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_mars_lmd_new/iniphysiq_mod.F

@@ -14 +14 @@
    !real*8,intent(in) :: ptimestep !physics time step (s) [dtphys]
 
+! this is necessary to initialize getin_p
+USE mod_phys_lmdz_para, ONLY: Init_phys_lmdz_para
+
 implicit none
+
+INCLUDE 'mpif.h'
 
 real,intent(in) :: prad ! radius of the planet (m)
@@ -44 +49 @@
 
 !! read callphys.def
+CALL Init_phys_lmdz_para(1,1,1,MPI_COMM_WORLD) ! initializes MPI communicator
+                                               ! ... necessary for getin_p in conf_phys
 call conf_phys(ngrid,nlayer,nq)