Hi all,
I'm using ISSM to estimate Balance Velocities, and I just wanted to understand a bit better what the model is doing and make sure that my results make sense. I typed "md = solve(md, 'Balancevelocity');" and I got some results. But I noticed that in order to solve for BalanceVelocity, I had to feed the variables md.initialization.vel and md.balancethickness.thickening_rate. I just set md.balancethickness.thickening_rate to zeros, but I'm not sure if this is the right approach here. Should I derive that from the mass accumulation field instead?
I wonder which method ISSM uses to estimate Balance Velocities, and how the model uses md.initialization.vel and md.balancethickness.thickening_rate in that calculation. I wrote a function myself to calculate balance velocities based on Joughin et al. (1997) approach, so I just wanted to compare ISSM function with mine.
Thank you in advance.
Best,
-Ana