Update src/Setup_geometry.jl

Co-authored-by: Albert de Montserrat <[email protected]>

Author: boriskaus

src/Setup_geometry.jl

@@ -430,7 +430,7 @@ function add_sphere!(Phase, Temp, Grid::AbstractGeneralGrid;    # required input
     end
 
     # Set the phase. Different routines are available for that - see below.
-    Phase[ind_flat] = compute_phase(Phase[ind_flat], Temp[ind_flat], X[ind], Y[ind], Z[ind], phase)
+    @views Phase[ind_flat] .= compute_phase(Phase[ind_flat], Temp[ind_flat], X[ind], Y[ind], Z[ind], phase)
 
     return nothing
 end