Here are the results so far, in iterations before blowing up: deltaT implicit advection component both none mom temp 900 211 217 164 190 700 1459 1530 1452 190 600 2310 1070 * 3660 * hasn't blow up so far (~20000)