source: issm/trunk-jpl/src/jl/md/parameterization.jl@ 26713

# setmask {{{
"""
SETMASK - establish boundaries between grounded and floating ice.

   By default, ice is considered grounded. The contour floatingicename defines nodes 
   for which ice is floating. The contour groundedicename defines nodes inside an floatingice, 
   that are grounded (ie: ice rises, islands, etc ...)
   All input files are in the Argus format (extension .exp).

   Usage:
      md=setmask(md,floatingicename,groundedicename)

   Examples:
      md=setmask(md,'all','');
      md=setmask(md,'Iceshelves.exp','Islands.exp');
"""
function setmask(md::model,floatingicename::String,groundedicename::String)

        elementonfloatingice = FlagElements( md, floatingicename)
        elementongroundedice = FlagElements( md, groundedicename)

        elementonfloatingice = convert( Array{Float64}, (elementonfloatingice.>0) .&  (elementongroundedice.==0.))
        elementongroundedice = convert( Array{Float64}, elementonfloatingice.==0.)

        vertexonfloatingice=zeros(md.mesh.numberofvertices)
        vertexongroundedice=zeros(md.mesh.numberofvertices)

        vertexongroundedice[md.mesh.elements[findall(elementongroundedice.>0),:]] .= 1.
        vertexonfloatingice[findall(vertexongroundedice.==0.)] .= 1.

        #define levelsets
        md.mask.ocean_levelset = vertexongroundedice
        md.mask.ocean_levelset[findall(vertexongroundedice .==0.)] .= -1.
        md.mask.ice_levelset = -1*ones(md.mesh.numberofvertices)

        return md
end
#}}}
# FlagElements {{{
function FlagElements(md::model,region::String)

        if isempty(region)
                flags = zeros(md.mesh.numberofelements)
        elseif region == "all"
                flags = ones(md.mesh.numberofelements)
        else
                xcenter = md.mesh.x[md.mesh.elements]*[1;1;1]/3
                ycenter = md.mesh.y[md.mesh.elements]*[1;1;1]/3
                flags = ContourToNodes(xcenter, ycenter, region, 2.)
        end

        return flags
end
#}}}