fixed "aquaplanet case" so that initializations (creation of
files startphy.nc and limit.nc) now also works in parallel (mpi,omp,mixed).
call to "iniaqua" is now done from within "iniphysiq" ; also added
some tests to check consistency between essential variables shared by
dynamics and physics (planetary radius, gravity, Cp, ...)
simillarily adapted "phydev" routines, and added necessary routines to also
be able to generate startphy.nc/restartphy.nc files there. Also removed
common file "comcstphy.h" and replaced it with a module "comcstphy.F90"