New:
---

A heatflux can be prescribed at the top with new inipar parameters
top_heatflux and use_top_fluxes. New 2d-arrays tsws and tswst are used to
store this flux. A Neumann boundary condition for
temperature can be applied under these conditions.

modules.f90 


Changed:
-------

Prognostic equations for all scalars are now solved up to gridpoint nzt
(formerly nzt-1).

The default top boundary condition for temperature has been renamed to
"initial gradient".

modules.f90


Errors:
------

Bugfix: 3d-array p is not a pointer any more.