Changeset 106 for palm/trunk/SOURCE/CURRENT_MODIFICATIONS

Timestamp:

Aug 16, 2007 2:30:26 PM (17 years ago)

Author:

raasch

Message:

preliminary update of bugfixes and extensions for non-cyclic BCs

File:

• palm/trunk/SOURCE/CURRENT_MODIFICATIONS (modified) (2 diffs)

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -13 +13 @@
 inipar parameter top_momentumflux_u|v.
 
-boundary_conds, check_open, check_parameters, diffusion_u, diffusion_v, flow_statistics, header, init_pegrid, init_3d_model, modules, palm, parin, prognostic_equations, read_var_list, read_3d_binary, swap_timelevel, time_integration, write_var_list, write_3d_binary
+Non-cyclic boundary conditions can be used along all horizontal directions.
+
+Quantities w*p* and w"e can be output as vertical profiles.
+
+boundary_conds, check_open, check_parameters, diffusion_u, diffusion_v, flow_statistics, header, init_pegrid, init_3d_model, modules, palm, parin, pres, prognostic_equations, read_var_list, read_3d_binary, swap_timelevel, time_integration, write_var_list, write_3d_binary
 
 New:
@@ -21 +25 @@
 Changed:
 -------
+Remaining 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.
+
+advec_particles, init_particles, time_integration
 
 
 Errors:
 ------
-+dots_num_palm in module user, +module netcdf_control in user_init (both in user_interface.f90)
+Bugs 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).
 
-check_parameters, user_interface
+advec_u_pw, advec_u_up, advec_v_pw, advec_v_up, boundary_conds, buoyancy, check_parameters, coriolis, diffusion_u, diffusion_v, init_pegrid, init_3d_model, modules, production_e, prognostic_equations, user_interface