So, this is still sort of clunky, and the reason is I don't want it to create a new instance of Geometry with every loop, so it is taking in r_max as an array and still uses the counter in the loop to change geometry. Same with joining the clustering realizations, I only want to make two of those, not four, so it's not going to loop over every element of x_interp_list. Let me know if a more elegant way is jumping out at you.
So, this is still sort of clunky, and the reason is I don't want it to create a new instance of Geometry with every loop, so it is taking in r_max as an array and still uses the counter in the loop to change geometry. Same with joining the clustering realizations, I only want to make two of those, not four, so it's not going to loop over every element of x_interp_list. Let me know if a more elegant way is jumping out at you.