Specification of traction boundary conditions over lateral surfaces.

pf4d

Dear users/developers :

I have not found any documentation regarding the manual specification of Neumann-type boundary conditions over lateral surfaces. From investigating the source code for md.SetMarineIceSheetBC(), it appears that I can do something like:

# set up element-wise multiplicative identities :

# rank-zero tensor vertex ones vector :
v_ones = np.ones(md.mesh.numberofvertices)

# rank-zero tensor element ones vector :
e_ones = np.ones(md.mesh.numberofelements)

# rank-two tensor ones vector :
A_ones = np.ones((md.mesh.numberofvertices, 6))

# rank-one tensor ones vector :
b_ones = np.ones((md.mesh.numberofvertices, 3))

md.stressbalance.spcvx       = np.nan * v_ones
md.stressbalance.spcvy       = np.nan * v_ones
md.stressbalance.spcvz       = np.nan * v_ones

md.basalforcings.groundedice_melting_rate = 0.0 * v_ones
md.basalforcings.floatingice_melting_rate = 0.0 * v_ones
md.stressbalance.referential              = np.nan * A_ones
md.stressbalance.loadingforce             = np.nan * b_ones

However, my simulation output (MISMIP+ geometry) has no-slip lateral boundaries:

Simulation output 1
Simulation output 2
Simulation output 3

Can you help?

Thanks!

pf4d

anyone?

mathieumorlighem

Hi Evan,
Sorry again for the delay. For a no slip boundary condition, don't you just need to set spcvx and spcvy to 0?
Best,
Mathieu