write DoY for ice on/off:
mod to getZmean function to handle empty