NEW: function to calculate steady state ice shelf temperature as a function of surface and basal temperature and basal melting rate. Accounts for vertical advection and conduction of heat within the ice column, does not account for horizontal advection