Update src/Setup_geometry.jl

Co-authored-by: Albert de Montserrat <58044444+albert-de-montserrat@users.noreply.github.com>

Author: boriskaus

src/Setup_geometry.jl

@@ -601,7 +601,7 @@ function add_cylinder!(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