02/06/2009

• Move continuity calculation from implicitsurface class to domain class.
• Modify calcContVel() to use u, v and w directly from domain object.
• We need the continuity update for vertical velocity w after nonhydrostatic update (continuity must be applied after exchanging of boundary conditions).

04/06/2009

• Theta-method is introduced in implicit surface calculation.
• Error in implicitDiffusionZ has corrected to take account boundary conditions properly.

05/06/2009

• Theta-method is introduced in implicit calculation of vertical diffusion.

29/06/2009

• Pressure correction method have improved to have 2nd-order accuracy in time following Armfield 2000.