シミュレーションがどうもうまくいかない

今はRempel(2005)の太陽の差動回転と子午面環流を実現するシミュレーションの再現をやっているのだけどどうにもうまくいかない。

今の問題をとりあえず、メモしておく。
低緯度のtachoclineからそれより下のところでなぜか負の速度が発生してそれによって、そこらへんに正のエントロピーができてしまいエネルギーがうまくバランスしないというところ。

うーん