+model version: @@ -2254,21 +2253,50 @@ modules, parin, production_e, prognostic_equations, read_var_list, read_3d_binary, swap_timelevel, write_var_list, write_3d_binary +modules, production_e, prognostic_equations, spline_z

Current
-model version: | 3.1c | ||||

3.2 | |||||

Last
change of this document: | $Id$ | ||||

C | Prognostic equations for all
-scalars are now solved up to gridpoint nzt (formerly nzt-1).
+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 +applied only at nzt+1). The default top boundary condition for temperature has been renamed to 'initial_gradient'. Calls -of dvrp_output_local, +of dvrp_output_local, which were commented out for a long time, are 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 | |||

E | Bugfix: 3d-array p is not a
-pointer any more. Bugfix in init_particles: -MPI_REAL -argument in MPI_ALLREDUCE -replaced by MPI_INTEGER -(caused error on NEC only). Bugfix: ddzw now -dimensioned 1:nzt"+1". | diffusion_e, diffusion_s, -diffusion_u, diffusion_v, diffusion_w, init_particles, modules | |||

E | Bugfix: 3d-array p is not a
+pointer any more. Bugfix in init_particles: +MPI_REAL +argument in MPI_ALLREDUCE +replaced by MPI_INTEGER +(caused error on NEC only). Bugfix: ddzw now +dimensioned 1:nzt"+1". | diffusion_e, diffusion_s, +diffusion_u, diffusion_v, diffusion_w, init_particles, modules | |||

29/03/07 | SR | 3.2 | N | The code, including scripts, utility programs and documentation is now under full control of subversion. This revision includes a quite large number of new features, changes and bugfixes. Particle reflection from vertical walls implemented. The particle SGS model is also adjusted to these vertical walls. Wall functions for vertical walls now include diabatic conditions. New subroutines wall_fluxes, wall_fluxes_e. New 4D-array rif_wall. Calculation/output of precipitation amount, precipitation rate and z0 (by setting "pra*", "prr*", "z0*" with data_output). The time interval on which the precipitation amount is defined is set by new d3par-parameter precipitation_amount_interval. New inipar-parameter loop_optimization to control the loop optimization method. New +inipar-parameter pt_reference. If given, this value is used as the +reference that is used in the buoyancy terms (otherwise, the instantaneous +horizontally averaged temperature is used). New d3par-parameter netcdf_64bit_3d +to switch on and especially to switch off 64bit offset for 3D +files only (because still some pd-software is unable to handle 64bit +offset). New d3par-parameter dt_max to set the maximum value allowed for the timestep. New user interfaces user_advec_particles and user_init_3d_model. The new initializing action "by_user" calls user_init_3d_model and allows the initial setting of all 3d arrays under complete user control. Samples added to the user interface which show how to add user-define time series quantities. Topography height informations are stored on arrays zu_s_inner and zw_w_inner and output to the 2d/3d NetCDF files. Unit 9 is generally opened for output of debug informations (file DEBUG_<pe#>). | Makefile advec_particles, +average_3d_data, buoyancy, calc_precipitation, check_open, +check_parameters, data_output_2d, diffusion_e, diffusion_u, +diffusion_v, diffusion_w, diffusivities, header, impact_of_latent_heat, +init_particles, init_3d_model, modules, netcdf, parin, production_e, +read_var_list, read_3d_binary, sum_up_3d_data, user_interface, +write_var_list, write_3d_binary new: wall_fluxes |

C | General revision of non-cyclic horizontal boundary conditions: Radiation boundary conditions are now used instead of Neumann conditions at the outflow (calculation needs velocity values for t-dt, which are stored on new arrays u_m_l, u_m_r, etc.). The calculation of a mean outflow is not needed any more. Volume flow control is added for the outflow boundary (currently only for the north boundary!!). Additional gridpoints along x and y (uxrp, vynp) are not needed any more and removed from the code. Routine "boundary_conds" now operates on timelevel t+dt and is not split in two parts (main, uvw_outflow) any more. Neumann boundary conditions are used at inflow/outflow in case of non-cyclic boundary conditions for all 2d-arrays that are handled by exchange_horiz_2d. Attention: the +non-cyclic boundary conditions of this revision still contain some bugs +which will be removed within the next revisions! The Bott-Chlond-scheme can now be used in the vectorization-optimized part (NEC, etc.). Equations are solved up to gridpoint nzt. The +FFT-method for solving the Poisson-equation is now working with Neumann +boundary conditions both at the bottom and the top. This requires +adjustments of the tridiagonal coefficients and subtracting the +horizontally averaged mean from the vertical velocity field. Particles-package is now part of the default code ("-p particles" is not needed as an mrun-option any more). New attribute age_m in particle_type. Move call of user_actions( 'after_integration' ) below the increments of times and counters. user_actions is now called for each statistic region and has as an argument the number of the respective region (sr). d3par-parameter data_output_ts is removed. Timeseries output for profil-software removed. Timeseries are now switched on by setting dt_dots. Timeseries data is collected in flow_statistics. Initial velocities at nzb+1 +are regarded for volume flow control in case they have been set zero +before (to avoid small timesteps); see new internal parameters u/v_nzb_p1_for_vfc. q is not allowed any more to become negative (prognostic_equations). poisfft_init is only called if the fft-solver is switched on (init_pegrid). d3par-parameter moisture renamed to humidity. Subversion global revision number is read from mrun and added to the run description header and to the run control (_rc) file. __vtk directives removed from main program. The uitility routine interpret_config now reads PALM environment variables from NAMELIST instead using the system call GETENV. | advec_s_bc, advec_u_pw, +advec_u_up, advec_v_pw, advec_v_up, asselin_filter, boundary_conds, check_parameters, +coriolis, data_output_dvrp, data_output_ptseries, data_output_ts, +data_output_2d, data_output_3d, diffusion_u, diffusion_v, +disturb_field, exchange_horiz, exchange_horiz_2d, flow_statistics, header, init_grid, +init_particles, init_pegrid, init_rankine, init_pt_anomaly, +init_1d_model, init_3d_model, modules, palm, package_parin, parin, +poisfft, poismg, prandtl_fluxes, pres, production_e, +prognostic_equations, read_var_list, read_3d_binary, sor, +swap_timelevel, time_integration, write_var_list, write_3d_binary | |||

E | Bugfix: preset of tendencies te_em, te_um, and te_vm (init_1d_model). Bugfix in sample for reading user defined data from restart file (user_init). Bugfix in setting diffusivities for cases with the outflow damping layer extending over more than one subdomain (init_3d_model). Check for possible negative humidities in the initial humidity profile. In Makefile, default suffixes are removed from the suffix list to avoid calling of m2c in cases of existing .mod files (problem on Linux systems). | Makefile check_parameters, init_1d_model, init_3d_model, user_interface |