def __init__(self,x,y,z,cells,nfaces,mask,method='nearest',grdfile=None): self.method=method # Initialise the trisearch array GridSearch.__init__(self,x,y,cells,nfaces=nfaces,force_inside=True) if self.method == 'linear': Grid.__init__(self,grdfile) self.datatmp = np.zeros(mask.shape,dtype=np.double) self.z = np.sort(z) self.z[-1]=10.0 # Set the surface layer to large self.Nkmax = z.size-1 self.mask3d = mask self.maskindex = -1*np.ones(self.mask3d.shape,dtype=np.int32) rr=0 for ii in range(self.mask3d.shape[0]): for jj in range(self.mask3d.shape[1]): if self.mask3d[ii,jj]: self.maskindex[ii,jj]=rr rr+=1
def __init__(self,x,y,cells,nfaces,method='nearest',grdfile=None): self.method=method # Initialise the trisearch array GridSearch.__init__(self,x,y,cells,nfaces=nfaces,force_inside=True)