Prevent from using specific sea ice boundary conditions in zdftke and zdfgls if nn_ice=0 (i.e. no ice model but key_si3 || key_cice)
With key_si3 or key_cice and nn_ice==0, one should not be allowed to use specific sea-ice boundary conditions in zdftke and zdfgls (nn_mxlice >1). Sea ice thickness array, even if being undefined, could indeed be used.
I recommend adding at initialization stage something like:
IF ((nn_mxlice>1).AND.(nn_ice==0)) CALL ctl_stop(...)