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 turn
md.transient.ismovingfront = 1`. The ice thickness, speed, geometry will evolve dynamically, you do not need to do anything else.
I hope this helps!
Mathieu