Implements ocean model BMR (-ocean bmr)
which reads forcing "shelfbmassflux" from file,
option "-adjust_bmr" enables adjustment of
shelfbmassflux wrt reference ice geometry
(which is variable "draft" from RTopo-1)