Changeset 3094 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Timestamp:

Oct 22, 2023, 1:42:30 PM (14 months ago)

Author:

emillour

Message:

Mars PCM:
Fix issues with dates when in parallel with OpenMP (missing copyin
statement when opening parallel section in iniphysiq). Added some threadprivate
clauses in saved module variables in comcstfi_h.F90 and planete_h.F90.
Prettyfied solarlong.F and made it a module. Likewise for conf_phys.F
EM

File:

• trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/iniphysiq_mod.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/iniphysiq_mod.F90

@@ -10 +10 @@
                      prad,pg,pr,pcpp,iflag_phys)
 
+! Modules from the dynamics:
 use infotrac, only : nqtot, & ! number of advected tracers
                      tname, & ! tracer names
                      nqperes,nqfils! MVals: variables isotopes
+use temps_mod, only: day_ini, day_end, hour_ini
+
+! Modules from the dynamics-physics interface:
+use inigeomphy_mod, only: inigeomphy
+
+! Modules from the physics:
 use comgeomfi_h, only: ini_fillgeom
-use temps_mod, only: day_ini, day_end, hour_ini
-use phys_state_var_init_mod, only: phys_state_var_init
-use inigeomphy_mod, only: inigeomphy
 use geometry_mod, only: cell_area, & ! physics grid area (m2)
                         longitude, & ! longitudes (rad)
                         latitude ! latitudes (rad)
-! necessary to get klon_omp
-USE mod_phys_lmdz_para, ONLY: klon_omp ! number of columns (on local omp grid)
-USE dimphy, ONLY: init_dimphy
+use phys_state_var_init_mod, only: phys_state_var_init
+use mod_phys_lmdz_para, only: klon_omp ! number of columns (on local omp grid)
+use conf_phys_mod, only: conf_phys
 
 implicit none
@@ -62 +66 @@
   !------------------------------------------
 
-!$OMP PARALLEL
-  
 ! copy some fundamental parameters to physics 
 ! and do some initializations 
-
-! Initialize dimphy module => Now done in physics_distribution_mod
-!call init_dimphy(klon_omp,nlayer)
+! Note that for OpenMP threadprivate variables a COPYIN clause is required
+!$OMP PARALLEL COPYIN(day_ini,day_end,hour_ini)
 
 call phys_state_var_init(klon_omp,nlayer,nqtot,tname, &
@@ -80 +81 @@
 call conf_phys(klon_omp,nlayer,nqtot)
 
-! Initialize some "temporal and calendar" related variables
-!CALL init_time(day_ini,hour_ini,punjours,ptimestep)
-
 !$OMP END PARALLEL