background density for empty cells

cpinte · Nov 21, 2024 · 702f99c · 702f99c
1 parent d12427b
commit 702f99c
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/density.f90 b/src/density.f90
@@ -1821,7 +1821,7 @@ subroutine normalize_dust_density(disk_dust_mass)
   do l=1,n_grains_tot
      somme=0.0_dp
      do icell=1,n_cells
-        if (densite_pouss(l,icell) <= 0.0_dp) densite_pouss(l,icell) = 0.0_dp
+        if (densite_pouss(l,icell) <= 0.0_dp) densite_pouss(l,icell) = tiny_real *  nbre_grains(l)
         somme=somme+densite_pouss(l,icell)*volume(icell)
      enddo !icell
      if (somme > tiny_dp) densite_pouss(l,:) = densite_pouss(l,:) / somme * nbre_grains(l) ! nbre_grains pour avoir Sum densite_pouss = 1  dans le disque

diff --git a/src/read_idefix.f90 b/src/read_idefix.f90
@@ -200,7 +200,7 @@ subroutine read_idefix_model()
              icell = cell_map(i,j,phik)
 
              densite_gaz(icell) = rho(i,i2,i3) * udens
-             densite_pouss(:,icell) = rho(i,i2,i3) * udens
+             densite_pouss(:,icell) = rho(i,i2,i3) * udens  ! flat dust distribution
 
              ! todo : check in cyl
              vfield3d(icell,1)  = vx1(i,i2,i3) * uvelocity ! vr