Hi Chris, You were right: the values of theta in the corner part of the halo regions affect slightly the solution but should not. It did the test in global_ocean.cs32x15 and it seems to be related to GMRedi (it's OK when I turn off GMRedi). Do you remember if the problem was in an atmospheric set-up or an oceanic one ? See you, Jean-Michel