an initial implementation; not correct yet for distributed because then the wall melt should be a cavity opening term; even in routing, if you add the wall melt then the bwat field blows up, more-or-less as expected; power alpha mostly controls this; why the problem occurs at the edge of null strip in thinish ice is not clear to me