New: --- A heatflux can be prescribed at the top with new inipar parameters top_heatflux and use_top_fluxes. New 2d-arrays qswst, qswst_m, tswst and tswst_m are used to store this flux. Use of fluxes are controlled with new index variable nzt_diff. A Neumann boundary condition for temperature can be applied under these conditions. Additionally, a Dirichlet condition for temperature can be used at the top. check_parameters, diffusion_s, flow_statistics, header, init_grid, init_3d_model, modules, parin, production_e, prognostic_equations, read_var_list, read_3d_binary, swap_timelevel, write_var_list, write_3d_binary Changed: ------- Prognostic equations for all scalars are now solved up to gridpoint nzt (formerly nzt-1). Boundary conditions for scalars at top adjusted respectively (now applied only at nzt+1). The default top boundary condition for temperature has been renamed to "initial gradient". dvrp_output_local, which was commented out for a long time, is now activated for all streams. advec_s_pw, boundary_conds, calc_precipitation, check_parameters, diffusion_e, diffusion_s, impact_of_latent_heat, init_dvrp, init_pt_anomaly, modules, production_e, prognostic_equations, spline_z Errors: ------ Bugfix: 3d-array p is not a pointer any more. Bugfix in init_particles: MPI_REAL in MPI_ALLREDUCE replaced by MPI_INTEGER (caused error on NEC only) Bugfix: ddzw dimensioned 1:nzt"+1" diffusion_e, diffusion_s, diffusion_u, diffusion_v, diffusion_w, init_particles, modules