Hi, Am I right in thinking that shelfice and GMREDI don;t work together? because gmredi always assumes that the top of the ocean is at k=1 (which can be fixed with a couple of mask appropriately located). Or am I missing something ? cheers, david