So far, only the initial hydrostatic pressure is entered into this equation.

Salinity (sa) and potential density (rho) are included as new 2d/3d output quantities. Vertical profiles of salinity (sa), salinity fluxes (w"sa", w*sa*, wsa), and potential density (rho) can also be output.advec_s_bc, average_3d_data,

boundary_conds, buoyancy,

check_parameters, data_output_2d,

data_output_3d, diffusion_e,

diffusivities, flow_statistics, header,

init_grid, init_3d_model, modules,

netcdf, parin, production_e,

prognostic_equations,

read_3d_binary, read_var_list,

sum_up_3d_data, swap_timelevel,

time_integration, user_interface,

write_var_list, write_3d_binary

new:

eqn_state_seawater, init_ocean CInipar-parameter use_pt_reference renamed use_reference.

Internal variable hydro_press renamed hyp, routine calc_mean_pt_profile renamed calc_mean_profile.

The format of the RUN_CONTROL file has been adjusted for ocean runs.

Defaults of .._vertical_gradient_levels changed from -1.0 to -9999999.9.advec_particles, buoyancy,

calc_liquid_water_content,

check_parameters, diffusion_e,

diffusivities, header,

init_cloud_physics, modules,

production_e, prognostic_equations,

run_control EBugfix: height above topography instead of height above level k=0 is used for calculating the mixing length (diffusion_e and diffusivities).

Bugfix: error in boundary condition for TKE removed (advec_s_bc).advec_s_bc, diffusion_e, diffusivities +Atmos. Oceanic Technol., 23, 1709-1728).

So far, only the initial hydrostatic pressure is entered into this equation.

Salinity (sa) and potential density (rho) are included as new 2d/3d output quantities. Vertical profiles of salinity (sa), salinity fluxes (w"sa", w*sa*, wsa), and potential density (rho) can also be output.advec_s_bc, average_3d_data,

boundary_conds, buoyancy,

check_parameters, data_output_2d,

data_output_3d, diffusion_e,

diffusivities, flow_statistics, header,

init_grid, init_3d_model, modules,

netcdf, parin, production_e,

prognostic_equations,

read_3d_binary, read_var_list,

sum_up_3d_data, swap_timelevel,

time_integration, user_interface,

write_var_list, write_3d_binary

new:

eqn_state_seawater, init_ocean CInipar-parameter use_pt_reference renamed use_reference.

Internal variable hydro_press renamed hyp, routine calc_mean_pt_profile renamed calc_mean_profile.

The format of the RUN_CONTROL file has been adjusted for ocean runs.

Defaults of .._vertical_gradient_levels changed from -1.0 to -9999999.9.advec_particles, buoyancy,

calc_liquid_water_content,

check_parameters, diffusion_e,

diffusivities, header,

init_cloud_physics, modules,

production_e, prognostic_equations,

run_control EBugfix: height above topography instead of height above level k=0 is used for calculating the mixing length (diffusion_e and diffusivities).

Bugfix: error in boundary condition for TKE removed (advec_s_bc).advec_s_bc, diffusion_e, diffusivities05/10/07SR3.4NThis version allows runs for a coupled atmosphere-ocean LES.

Coupling frequency is controlled by new d3par-parameter dt_coupling, the coupling mode ('atmosphere_to_ocean' or 'ocean_to_atmosphere') for the respective processes is read from environment variable coupling_mode, which is set by the mpiexec-command in mrun. Communication between the two models is handled by the intercommunicator comm_inter. Local files opened by the ocean model get the additional suffic "_O". +Saturation at k=nzb is assumed for the atmosphere coupled to ocean. +Usage of the coupled model is described in the new section 3.8.

A momentum flux can be set as top boundary condition using the new inipar parameters top_momentumflux_u|v.

Non-cyclic boundary conditions can be used along all horizontal directions.

Quantities w*p* and w"e can be output as vertical profiles.

Initial profiles are reset to constant profiles in case that initializing_actions /= 'set_constant_profiles'. (init_rankine)

Optionally calculate km and kh from initial TKE e_init.

Attention:

+The format of the testresult-file (run control output) has changed +due to this modification.boundary_conds, +check_open, check_parameters, diffusion_u, diffusion_v, +flow_statistics, header, init_pegrid, init_rankine, init_3d_model, +local_stop, modules, palm, parin, prandtl_fluxes, pres, +prognostic_equations, read_var_list, read_3d_binary, swap_timelevel, +timestep, time_integration, write_var_list, write_3d_binary

new:

surface_coupler CRemaining variables iran changed to iran_part (advec_particles, init_particles).

In case that the presure solver is not called for every Runge-Kutta substep (call_psolver_at_all_substeps = .F.), +it is called after the first substep instead of the last. In that case, +random perturbations are also added to the velocity field after the +first substep.

Initialization of km, kh = 0.00001 for ocean = .T. (for ocean = .F. it remains 0.01).

Allow data_output_pr = 'q', 'wq', 'w"q"', 'w*q*' for humidity = .T. (instead of cloud_physics = .T.).advec_particles, check_for_restart, check_parameters, init_particles, init_3d_model, time_integration

EBugs from code parts for non-cyclic boundary conditions are removed:

loops for u and v are starting from index nxlu, nysv, +respectively. The radiation boundary condition is used for every +Runge-Kutta substep. Velocity phase speeds for the radiation boundary +conditions are calculated for the first Runge-Kutta substep only and +reused for the further substeps. New arrays c_u, c_v, and c_w +are defined for this purpose. Several index errors are removed from the +radiation boundary condition code parts. Upper bounds for calculating u_0 and v_0 (in production_e) are nxr+1 and nyn+1 because otherwise these values are not available in case of non-cyclic boundary conditions.

+dots_num_palm in module user, +module netcdf_control in user_init (both in user_interface)

Bugfix: wrong sign removed from the buoyancy production term in the case use_reference = .T. (production_e)

Bugfix: Error message concerning output of particle concentration (pc) modified (check_parameters).

Bugfix: Rayleigh damping for ocean fixed.advec_u_pw, +advec_u_up, advec_v_pw, advec_v_up, boundary_conds, buoyancy, +check_parameters, coriolis, diffusion_u, diffusion_v, header, +init_pegrid, init_3d_model, modules, production_e, +prognostic_equations, user_interface