Confusion over units in SIMPLE

[matthiasmengel]

@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

double basal_mass_balance   = use_basal_melt_rate ? -basal_melt_rate(i, j) : 0.0;

, 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.