Create the laplacian operator. Lap = divand_laplacian(mask,pmn,nu,iscyclic) Form a Laplacian using finite differences assumes that gradient is zero at "coastline" Input: mask: binary mask delimiting the domain. 1 is inside and 0 outside. For oceanographic application, this is the land-sea mask. pmn: scale factor of the grid. nu: diffusion coefficient of the Laplacian field of the size mask or cell arrays of fields Output: Lap: sparce matrix represeting a Laplaciant
Package: divand