Luckily, the data is sorted the right way for this to work…
dRhoM = reshape(dRho[:,4], length(Long), length(Lat), length(Zi));
Other than that, the print statement is a major source of slowing down. The loops are in the wrong order too. Remember, Julia lays out arrays in column-major-order.