Changeset 26706

Timestamp:

12/03/21 11:26:44 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added IsIcefront function

Location:

issm/trunk-jpl/src/jl/solve

Files:

• analyses/stressbalanceanalysis.jl (modified) (1 diff)
• elements.jl (modified) (1 diff)

issm/trunk-jpl/src/jl/solve/analyses/stressbalanceanalysis.jl

@@ -174 +174 @@
                         pe.values[2*i  ] += -gauss.weights[ig]*Jdet*rho_ice*g*H*ds[2]*basis[i]
                 end
+        end
+
+        if(IsIcefront(element))
+                error("not supported")
         end
 

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

@@ -160 +160 @@
         return phi
 end#}}}
+function IsIcefront(element::Tria) #{{{
+
+        level = Vector{Float64}(undef,3)
+        GetInputListOnVertices!(element, level, MaskIceLevelsetEnum)
+
+        nbice = 0
+        for i in 1:3
+                if(level[i]<0.) nbice+=1 end
+        end
+
+        if(nbice==1)
+                return true
+        else
+                return false
+        end
+end#}}}
 
 #Finite Element stuff