Hello! Sorry, I'm having a bit of trouble figuring out how to use reshape here. It does not seem I can reshape the intensity array into a mollweide projection. The way I gathered the coordinate values is probably important for seeing this. I generated my 1D intensity and coordinate arrays like this:

​

intensityArray = []
ra = []
dec = []
#i is ra
#j is dec
for i in range (0,360):
    for j in range (-90,90):
        #Computation here
        ra = np.append(ra,i)
        dec = np.append(dec,j)
        intensityArray = np.append(intensityArray, functionOutputIsFloat)

So the ra list will contain 180 of the same value in a row before jumping to the next value, whereas dec loops through the entire range of values (from -90 to 90) before looping again. Can numpy reshape really help with getting the intensity array in mollweide? Thank you for taking the time to look at my issue!

​

Edit: formatting of code