get rid of global sponge array, pass as argument workc instead. this can be used for FSI and passive scalar as well. also, flusi shows how much memory she allocated.