Hi, I am experimenting with a simple model setup to simulate a moving ice front over land. As an initial condition I chose a round "blob" of ice in the centre of the domain which terminates at some radius. I would like the rest of the domain to be empty (land), but I see that the value of md.masstransport.min_thickness is applied everywhere where I defined a md.geometry.thickness of zero.
Is there a way to have parts of the domain empty and perhaps use the md.transient.ismovingfront = 1 module to simulate a propagating ice front into these regions?
I tried setting the md.mask.ice_levelset fields according to the initial geometry where -1 is the blob, 0 is at the front and 1 everywhere else. I get the following error prompt during the first iteration, and I am not sure what that means.
[0] ??? Error using ==> ./classes/matrix/ElementMatrix.cpp:233
[0] CheckConsistency error message: NaN found in Element Matrix
Also, if I raise my initial condition by a meter, then after the first iteration the md.results.TransientSolution(1).Thickness field just equals 1 everywhere, even though my initial condition in the centre is way above 1.
Thanks,
Jeremie
