1. Rewrite source code.
From now on diagnose.f90 include dry_test and whschubert1982 boussinesq
approx.

2. Add testing unit (under developing...)

3. New problem : integral check problem.