I believe it will be useful to smooth the initial condition fields before starting the model. I have had improved results with this choice in MOM6 (for instance) and will be hardcode it in pyroms_tools/initial/ic_file_build.py for now. I'll be using a gaussian filter in order to implement it.