Hi, I am trying to implement a model in High Order approximation and, in this particular study I wish to force the glacier terminus movement with the levelset method. However, it is not clear to me what is the proper way of specifying levelset after the extrusion. The terminus data only provides the surface ice-ocean boundary, but all subsequent layers need their spclevelset value and the glacier front thickness / geometry evolves dynamically, thus below the surface the boundary levelset cannot be prescribed. What is the right approach here?

Thank you!
Donglai

Hi @donglaiyang

I am not quite sure to understand your problem. The level set should be uniform vertically. I highly recommend that you first start doing your simulation in 2D using SSA, and once you are happy with it, you can extrude the model and run exactly the same code after setting the stress balance equations to HO.

In order to "force" the motion of the calving front you need to make md.levelset.spclevelset a time series (i.e. with an additional row with the time assigned to each specified levelsetand turnmd.transient.ismovingfront = 1`. The ice thickness, speed, geometry will evolve dynamically, you do not need to do anything else.

I hope this helps!
Mathieu

Ah sorry some other error led me to think that spclevelset needs to be prescribed at every level, but thanks