Bug fix: checking density at initialisation

cpinte · May 30, 2024 · 23b5cae · 23b5cae
1 parent 13e61a6
commit 23b5cae
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 4 deletions.
diff --git a/src/density.f90 b/src/density.f90
@@ -965,7 +965,6 @@ subroutine is_density_file_Voronoi()
      write(*,*) "Found 1D density structure, using a Voronoi mesh"
      lVoronoi = .true.
      l3D = .true.
-
      n_cells = naxes(1)
   else
      write(*,*) "Found a structured mesh"

diff --git a/src/dust_transfer.f90 b/src/dust_transfer.f90
@@ -109,8 +109,6 @@ subroutine transfert_poussiere()
 
   laffichage=.true.
 
-  if (ldensity_file) call is_density_file_Voronoi()
-
   if (lVoronoi) then
      if (lmhd_voronoi) then
         call setup_mhd_to_mcfost() !uses sph_to_voronoi

diff --git a/src/init_mcfost.f90 b/src/init_mcfost.f90
@@ -3,7 +3,7 @@ module init_mcfost
   use parametres
   use naleat
   use grains, only : aggregate_file, mueller_aggregate_file
-  use density, only : species_removed, T_rm
+  use density, only : species_removed, T_rm, is_density_file_Voronoi
   use molecular_emission
   !$ use omp_lib
   use benchmarks
@@ -1477,6 +1477,8 @@ subroutine initialisation_mcfost()
      call read_para(para)
   endif
 
+  if (ldensity_file) call is_density_file_Voronoi()
+
   if (lfargo3d) then
      l3D = .true.
      if (n_zones > 1) call error("fargo3d mode only work with 1 zone")