Closed duochanatharvard closed 3 years ago
Nevermind, this function seems to find the weight for depth instead of only in the mixing layer.
Right, the original code attempts to get the volume of the Indian Ocean. It knows the surface definition of the Indian, and then needs to identify all subsurface gridcells below that surface.
c_ind = mixit(d_ind,it,jt,kt,ones(Nfield,1)); Seems the last input should be "inmixlyr" for the function to return an indicator of wether in the mix-layer.