Changeset 26738 for issm/trunk-jpl/src/jl/core/elements.jl

Timestamp:

12/17/21 12:21:15 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: account for passive elements now if the ice levelset crosses the domain

File:

• issm/trunk-jpl/src/jl/core/elements.jl (modified) (3 diffs)

issm/trunk-jpl/src/jl/core/elements.jl

@@ -46 +46 @@
         end
 end #}}}
+function InputUpdateFromVector(element::Tria, vector::Vector{Float64}, enum::IssmEnum, layout::IssmEnum) #{{{
+
+        lidlist = element.vertexids
+        data = Vector{Float64}(undef, 3)
+
+        if(layout==VertexSIdEnum)
+                for i in 1:3
+                        data[i] = vector[element.vertices[i].sid]
+                        @assert isfinite(data[i])
+                end
+                SetTriaInput(element.inputs, enum, P1Enum, lidlist, data)
+        else
+                error("layout ", layout, " not supported yet");
+        end
+end #}}}
 function Update(element::Tria, inputs::Inputs, index::Int64, md::model, finiteelement::IssmEnum) #{{{
 
@@ -108 +123 @@
         InputServe!(element,input)
         return input
+
+end # }}}
+function GetInputListOnNodes!(element::Tria, vector::Vector{Float64}, enum::IssmEnum) # {{{
+
+        #Get Input first 
+        input = GetInput(element, enum)
+
+        #Get value at each vertex (i.e. P1 Nodes)
+        gauss=GaussTria(P1Enum)
+        for i in 1:gauss.numgauss
+                vector[i] = GetInputValue(input, gauss, i)
+        end
 
 end # }}}
@@ -197 +224 @@
         end
 end#}}}
+function IsIceInElement(element::Tria) #{{{
+        #We consider that an element has ice if at least one of its nodes has a negative level set
+
+        input=GetInput(element, MaskIceLevelsetEnum)
+
+        if GetInputMin(input)<0
+                return true
+        else
+                return false
+        end
+
+end#}}}
 function GetIcefrontCoordinates!(element::Tria, xyz_front::Matrix{Float64}, xyz_list::Matrix{Float64}, levelsetenum::IssmEnum) #{{{