* Update landmask handling to have land = 1, ocean = 0, to be consistent with the idea that if L is the landmask, then X[L] are the land pixels. * Update apply landmask to match the new method * Find locations where landmask is multiplied instead of using apply landmask and correct those