Jinlun, we have been using 3rd order direct space time with flux limiter http://mitgcm.org/sealion/online_documents/node78.html for the cube runs. Advection schemes are a big black box for me. Do you know if central differencing is less diffusive than above scheme? Dimitris