BUGFIX - fixes #434

danieljprice · Sep 12, 2023 · a5ee42a · a5ee42a
1 parent 377b2c0
commit a5ee42a
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/src/setup/set_disc.F90 b/src/setup/set_disc.F90
@@ -323,6 +323,13 @@ subroutine set_disc(id,master,mixture,nparttot,npart,npart_start,rmin,rmax, &
  !
  sigma_normdust = 1.d0
  if (do_mixture) then
+    if (present(r_grid)) then
+       rad_tmp = r_grid
+    else
+       do i=1,maxbins
+          rad_tmp(i) = R_indust + (i-1) * (R_outdust-R_indust)/real(maxbins-1)
+       enddo
+    endif
     !--sigma_normdust set from dust disc mass
     call get_disc_mass(disc_mdust,enc_m_tmp,rad_tmp,Q_tmp,sigmaprofiledust, &
                        sigma_normdust,star_m,p_indexdust,q_inddust, &