Hi Mathieu,
thank you for your answers and the documentation in the email 😃
mathieumorlighem remember that we are comparing Dv/dt to gravity, with speeds of 1000s m/year at most, inertia will always be negligible.
Makes sense, thanks
mathieumorlighem We tried some tests some time ago and it seemed like it was not worth implementing it. We can revisit this of course.
I thought it might if some equations are expensive, but it does of course introduce the additional error of skipping steps or averaging. I don't have a good feeling yet for what changes at which scale (in particular where seasons might even come in), so I'll just trust you on that, especially since you've already tested it out and decided against it.
mathieumorlighem you can output anything that is an "Input" in the C++ core. Let me know if there is something you want to see and I can tell you if it exists (or you can familiarize yourself with the code since it sounds like you have a good background in numerical modeling).
I was mostly thinking about sea level rise, but otherwise just generally curious about what's happening in the model and what people are interested in understanding. In the core files I found sealevelchange_core.cpp, and now noticed that there's the isslc
variable in the transient parameters - so that answers that question.
I've browsed through some of the c++ files, and I've found a list of outputs in /trunk/src/c/classes/FemModel
, line 2395 onwards (I tested some of the case names without the suffix "Enum" as "requested outputs" and that worked). This is for future reference if someone else is looking.
The part about the ice volume calculation and loadmodel
makes sense, thanks. For the latter it just took me a while to figure out why the same model would run only sometimes, so I wanted to mention it.
Thank you a lot for your help!
Nicole