New:
---
Restart runs on SGI-ICE are working (mrun).
2d-decomposition is default on SGI-ICE systems. (init_pegrid)

Ocean-atmosphere coupling realized with MPI-1. mrun adjusted for this case
(-Y option). Adjustments in mrun, mbuild, and subjob for lcxt4.

DVRP arguments changed to single precision, mode pathlines added.

User can add additional routines in files user_additional_routines.f90


check_for_restart, check_parameters, data_output_dvrp, init_dvrp, init_pegrid, local_stop, modules, package_parin, palm, surface_coupler, timestep, user_additional_routines
Makefile, mrun, mbuild, subjob

New: init_coupling


Changed:
-------
User interface has been split into single files.

Makefile

New: user_3d_data_averaging user_actions user_advec_particles user_check_data_output user_check_data_output_pr user_data_output_2d user_data_output_3d user_data_output_dvrp user_define_netcdf_grid user_dvrp_coltab user_header user_init user_init_3d_model user_init_grid user_init_particles user_init_plant_canopy user_last_actions user_module user_parin user_particle_attributes user_read_restart_data user_spectra user_statistics

Removed: user_interface

Errors:
------

Size of pf3d adjusted to the required output size (1 gridpoint less, along
all three dimensions), because output of a subset of the data
(pf3d(nxa:nxe...) in the NF90_PUT_VAR statement caused segmentation fault
with the INTEL compiler. (combine_plot_fields)

Bugfix: error in zu index in case of section_xy = -1 (header)

Inconsistency removed: as the thermal stratification is not taken into account for the evaluation of the wall 
fluxes at vertical walls, the eddy viscosity km must not be used for the evaluation of the velocity gradients dudy, 
dwdy, dvdx and dwdx. (production_e)

Bugfix in calculating k index in case of oceans runs. (sort_particles)

advec_particles, combine_plot_fields, header, production_e