Changeset 26714

Timestamp:

12/06/21 18:34:28 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: more change on Pig, still not working, but close

File:

• issm/trunk-jpl/src/jl/test501.jl (modified) (1 diff)

issm/trunk-jpl/src/jl/test501.jl

@@ -34 +34 @@
 
 #Boundary conditions
+pos = findall(vec(sum(Int64.(md.mask.ocean_levelset[md.mesh.elements].<0), dims=2)) .> 0.0)
+vertexonfloatingice=zeros(md.mesh.numberofvertices)
+vertexonfloatingice[md.mesh.elements[pos,:]] .= 1
+nodefront=(md.mesh.vertexonboundary .& (vertexonfloatingice.>0))
+md.mask.ice_levelset[findall(nodefront)] .= 0
+
 md.stressbalance.spcvx = NaN*ones(md.mesh.numberofvertices)
 md.stressbalance.spcvy = NaN*ones(md.mesh.numberofvertices)
-pos = findall(md.mesh.vertexonboundary)
+segmentsfront=md.mask.ice_levelset[md.mesh.segments[:,1:2]]==0
+segments = findall(vec(sum(Int64.(md.mask.ice_levelset[md.mesh.segments[:,1:2]].==0), dims=2)) .!=2)
+pos=md.mesh.segments[segments,1:2]
 md.stressbalance.spcvx[pos] .= 0.0
 md.stressbalance.spcvy[pos] .= 0.0