Replace sum(ge) - sum(ge(1, :)) by sum(ge(:iim, :)) in procedure sortvarc. Remove intermediary variables vor and ztotl.
trunk/dyn3d/sortvarc.f