Hi, The model crash and the unrealistic SSH have both been fixed by setting: #undef INCLUDE_CD_CODE in CPP_OPTIONS.h Thanks to Jake Gebbie for pointing that out. It is no longer necessary to use implicit viscosity/diffusion which is good because I still can't get either to work properly. Cheers Dan