@ronjareese @talbrecht
In my first prognostic runs, I get quickly spreading shelves throughout Antarctica. Hunting this phenomenon, I realized that we have some issues with units that become important for prognostic simulations.
With Ronja's latest commit, melt rates look good in the diagnostic variable basalmeltrate_shelf (0-8 m/year). However, if I print out what arrives in iMgeometry.cc around line 710
, I get at max 3e-8 m/s, which rougly corresponds to 1 m/year melt. This may explain well the fast gain of ice volume and fast spreading of shelves that I find.
Any ideas on that?
Related: the basalmeltrate_shelf is somehow not transferred to the bmelt variable. I find it equal to zero in shelfy areas.
@ronjareese @talbrecht In my first prognostic runs, I get quickly spreading shelves throughout Antarctica. Hunting this phenomenon, I realized that we have some issues with units that become important for prognostic simulations.
With Ronja's latest commit, melt rates look good in the diagnostic variable basalmeltrate_shelf (0-8 m/year). However, if I print out what arrives in iMgeometry.cc around line 710
, I get at max 3e-8 m/s, which rougly corresponds to 1 m/year melt. This may explain well the fast gain of ice volume and fast spreading of shelves that I find. Any ideas on that?
Related: the basalmeltrate_shelf is somehow not transferred to the bmelt variable. I find it equal to zero in shelfy areas.