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