Changes between Version 2 and Version 3 of doc/tec/changelog_2009-2012
- Timestamp:
- Feb 23, 2021 3:47:07 PM (4 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
doc/tec/changelog_2009-2012
v2 v3 142 142 {{{#!td style="vertical-align:top" 143 143 '''Bugfixes:'''\\ 144 It is not allowed to use [wiki:doc/app/ini par#cloud_scheme cloud_scheme] = 'seifert_beheng' without [wiki:doc/app/inipar#precipitation precipitation] = .TRUE. (check_parameters).145 146 Boundaries of ql-cross-sections in case of [wiki:doc/app/ini par#cloud_scheme cloud_scheme] = 'seifert_beheng'] adjusted (data_output_2d).147 148 '''Changes in the two-moments cloud physics scheme ([wiki:doc/app/ini par#cloud_scheme cloud_scheme] = 'seifert_beheng':'''\\149 Sedimentation of rain according to Stevens and Seifert (2008) has been implemented. Monotonicity of the sedimentation process is insured by limiting the slopes of qr and nr. The limiter can be steered with the parameter [wiki:doc/app/ini par#limiter_sedimentation limiter_sedimentation]. The default is [wiki:doc/app/inipar#limiter_sedimentation limiter_sedimentation] = .TRUE.. Therefore, the Courant number of the sedimentation process can be greater than one. The Courant number can be controlled by [wiki:doc/app/inipar#c_sedimentation c_sedimentation] (default: [wiki:doc/app/inipar#c_sedimentation c_sedimentation] = 2.0) (microphysics, init_cloud_physics, modules, parin, read_var_list, write_var_list).144 It is not allowed to use [wiki:doc/app/initialization_parameters#cloud_scheme cloud_scheme] = 'seifert_beheng' without [wiki:doc/app/initialization_parameters#precipitation precipitation] = .TRUE. (check_parameters). 145 146 Boundaries of ql-cross-sections in case of [wiki:doc/app/initialization_parameters#cloud_scheme cloud_scheme] = 'seifert_beheng'] adjusted (data_output_2d). 147 148 '''Changes in the two-moments cloud physics scheme ([wiki:doc/app/initialization_parameters#cloud_scheme cloud_scheme] = 'seifert_beheng':'''\\ 149 Sedimentation of rain according to Stevens and Seifert (2008) has been implemented. Monotonicity of the sedimentation process is insured by limiting the slopes of qr and nr. The limiter can be steered with the parameter [wiki:doc/app/initialization_parameters#limiter_sedimentation limiter_sedimentation]. The default is [wiki:doc/app/initialization_parameters#limiter_sedimentation limiter_sedimentation] = .TRUE.. Therefore, the Courant number of the sedimentation process can be greater than one. The Courant number can be controlled by [wiki:doc/app/initialization_parameters#c_sedimentation c_sedimentation] (default: [wiki:doc/app/initialization_parameters#c_sedimentation c_sedimentation] = 2.0) (microphysics, init_cloud_physics, modules, parin, read_var_list, write_var_list). 150 150 151 151 Turbulence effects on autoconversion and accretion according to Seifert, Nuijens and Stevens (2010) have been implemented and can be activated with [wiki:doc/app/d3par#turbulence turbulence] = .TRUE. (microphysics, diffusion_e, init_3d_model, modules, parin, read_var_list, time_integration, write_var_list). 152 152 153 The number density of cloud droplets can be adjusted with [wiki:doc/app/ini par#nc nc]. The default ([wiki:doc/app/inipar#nc nc] = 70.0E6 ) is applicable for marine conditions (parin, read_var_list, write_var_list).153 The number density of cloud droplets can be adjusted with [wiki:doc/app/initialization_parameters#nc nc]. The default ([wiki:doc/app/initialization_parameters#nc nc] = 70.0E6 ) is applicable for marine conditions (parin, read_var_list, write_var_list). 154 154 155 155 '''Attention:'''\\ … … 196 196 Two-moment cloud physics (Seifert and Beheng, 2006) have been implemented. Therefore, two new prognostic equations (nr: rain drop number density; qr: rain water content) and associated basic output quantities are introduced (advec_ws, average_3d_data, boundary_conds, check_parameter, data_output_2d, data_output_3d, flow_statistics, init_3d_model, init_cloud_physics, modules, netcdf, parin, prognostic_equations, read_3d_binary, read_var_list, sum_up_3d_data, swap_timelevel, time_integration, user_actions, write_3d_binary, write_var_list). Microphysical tendencies affecting the prognostic equations of liquid water potential temperature, specific humidity, rain drop number density and rain water content are computed in the new module microphysics (additional changes: calc_liquid_water_content, Makefile). According to the CFL criterion, the timestep has to be reduced in case of precipitation (timestep). 197 197 198 The new cloud physics scheme can be used with [wiki:doc/app/ini par#humidty humidity] = .TRUE. and [wiki:doc/app/inipar#cloud_physics cloud_physics] = .TRUE.. It is the the default value of the new steering parameter [wiki:doc/app/inipar#cloud_scheme cloud_scheme] = 'seifert_beheng' ([wiki:doc/app/inipar#cloud_scheme cloud_scheme] = 'kessler' allows the use of the former cloud physics scheme). It is not allowed to use the new cloud physics without [wiki:doc/app/inipar#precipitation precipitation] = .TRUE. ([wiki:doc/app/errmsg#PA0363 PA0363]).199 200 The ventilation effect on evaporation of rain drops according to Seifert (2008) has been implemented and can be activated with [wiki:doc/app/ini par#ventilation_effect ventilation_effect].201 202 Drizzle of raindrops according to Heus et al. (2010) has been implemented and can be activated with [wiki:doc/app/ini par#drizzle drizzle].198 The new cloud physics scheme can be used with [wiki:doc/app/initialization_parameters#humidty humidity] = .TRUE. and [wiki:doc/app/initialization_parameters#cloud_physics cloud_physics] = .TRUE.. It is the the default value of the new steering parameter [wiki:doc/app/initialization_parameters#cloud_scheme cloud_scheme] = 'seifert_beheng' ([wiki:doc/app/initialization_parameters#cloud_scheme cloud_scheme] = 'kessler' allows the use of the former cloud physics scheme). It is not allowed to use the new cloud physics without [wiki:doc/app/initialization_parameters#precipitation precipitation] = .TRUE. ([wiki:doc/app/errmsg#PA0363 PA0363]). 199 200 The ventilation effect on evaporation of rain drops according to Seifert (2008) has been implemented and can be activated with [wiki:doc/app/initialization_parameters#ventilation_effect ventilation_effect]. 201 202 Drizzle of raindrops according to Heus et al. (2010) has been implemented and can be activated with [wiki:doc/app/initialization_parameters#drizzle drizzle]. 203 203 204 204 '''Attention:'''\\ … … 283 283 mask locations determined based on scalar positions (init_masks) 284 284 285 save memory by not allocating pt_2 in case of [wiki:doc/app/ini par#neutral neutral] = .T. (init_3d_model, swap_timelevel; see also r940)285 save memory by not allocating pt_2 in case of [wiki:doc/app/initialization_parameters#neutral neutral] = .T. (init_3d_model, swap_timelevel; see also r940) 286 286 287 287 minor reformatting (check_for_restart) … … 652 652 {{{#!td style="vertical-align:top" 653 653 '''New:''' \\ 654 Roughness length for scalar quantities z0h can be set by the new inipar-parameter [wiki:doc/app/ini par#z0h_factor z0h_factor].654 Roughness length for scalar quantities z0h can be set by the new inipar-parameter [wiki:doc/app/initialization_parameters#z0h_factor z0h_factor]. 655 655 (average_3d_data, check_parameters, data_output_2d, header, init_1d_model, init_3d_model, init_masks, modules, parin, prandtl_fluxes, read_3d_binary, read_var_list, sum_up_3d_data, write_3d_binary, write_var_list) 656 656 657 657 An inflow damping zone for the potential temperature has been implemented for the case of non-cyclic lateral boundaries. 658 It is steered by new inipar-parameter [wiki:doc/app/ini par#pt_damping_width pt_damping_width] and [wiki:doc/app/inipar#pt_damping_factor pt_damping_factor].658 It is steered by new inipar-parameter [wiki:doc/app/initialization_parameters#pt_damping_width pt_damping_width] and [wiki:doc/app/initialization_parameters#pt_damping_factor pt_damping_factor]. 659 659 (check_parameters, header, init_3d_model, modules, parin, prognostic_equations, read_var_list, write_var_list) 660 660 … … 813 813 }}} 814 814 {{{#!td style="vertical-align:top" 815 Prognostic equation for temperature can be switched off with new inipar parameter [wiki:doc/app/ini par#neutral neutral] in order to spare cpu time in case of simulating pure neutrally stratified flows. This will also switch off calculation of all buoyancy related terms.815 Prognostic equation for temperature can be switched off with new inipar parameter [wiki:doc/app/initialization_parameters#neutral neutral] in order to spare cpu time in case of simulating pure neutrally stratified flows. This will also switch off calculation of all buoyancy related terms. 816 816 (check_parameters, modules, header, parin, production_e, prognostic_equations, read_var_list, write_var_list) 817 817 }}} … … 873 873 }}} 874 874 {{{#!td style="vertical-align:top" 875 Masking method can be switched on for multigrid-solver using inipar parameter [../../app/ini par/#masking_method masking_method]. This feature has been implemented because some topography data have caused job aborts which are probably due to wrong boundary condition settings in the multigrid-solver. This problem will be adressed in a future revision. (header.f90, init_grid.f90, modules.f90, read_var_list.f90, write_var_list.f90).875 Masking method can be switched on for multigrid-solver using inipar parameter [../../app/initialization_parameters#masking_method masking_method]. This feature has been implemented because some topography data have caused job aborts which are probably due to wrong boundary condition settings in the multigrid-solver. This problem will be adressed in a future revision. (header.f90, init_grid.f90, modules.f90, read_var_list.f90, write_var_list.f90). 876 876 }}} 877 877 |---------------- … … 1430 1430 }}} 1431 1431 {{{#!td style="vertical-align:top" 1432 Droplet growth by condensation may include curvature and solution effects. These are steered by the new inipar-parameter [../../app/ini par/#curvature_solution_effects curvature_solution_effects]. The droplet growth equation by condensation is solved with the Rosenbrock method for stiff o.d.e's with monitoring local truncation error to adjust stepsize (see Numerical recipes in FORTRAN, 2nd edition, p. 731). (advec_particles, check_parameters, header, init_cloud_physics, init_particles, modules, parin, read_var_list, write_var_list)1432 Droplet growth by condensation may include curvature and solution effects. These are steered by the new inipar-parameter [../../app/initialization_parameters#curvature_solution_effects curvature_solution_effects]. The droplet growth equation by condensation is solved with the Rosenbrock method for stiff o.d.e's with monitoring local truncation error to adjust stepsize (see Numerical recipes in FORTRAN, 2nd edition, p. 731). (advec_particles, check_parameters, header, init_cloud_physics, init_particles, modules, parin, read_var_list, write_var_list) 1433 1433 1434 1434 Mean, minimum, and maximum particle radius added as output quantity to the particle time series. (data_output_ptseries, modules) … … 1591 1591 }}} 1592 1592 {{{#!td style="vertical-align:top" 1593 New: the utility program '''check_namelist_files''' has been implemented. It is used by mrun to check the parameter files (e.g. *_p3d, *_p3df) before the batch job is submitted/the interactive job is started. The source code of palm has thus been modified by means of the '''pre-processor directive !__check''', which is only used for the compilation of the check_namelist_files program. Basically the directive is used to exclude the MPI-functions for compilation with a serial compiler (on the local host). Note: the namelist file check is currently not possible for coupled runs and for [../../app/ini par/#initializing_actions initializing_actions='cyclic_fill']! A new mrun option '-z' has been implemented to switch off check_namelist_files program (PALM will of course still check the namelist files when the simulation starts) (Makefile_check, mbuild, mrun, check_namelist_files, init_pegrid, check_parameters, init_masks, poisfft, poisfft_hybrid, local_stop, init_grid, check_open, modules).1593 New: the utility program '''check_namelist_files''' has been implemented. It is used by mrun to check the parameter files (e.g. *_p3d, *_p3df) before the batch job is submitted/the interactive job is started. The source code of palm has thus been modified by means of the '''pre-processor directive !__check''', which is only used for the compilation of the check_namelist_files program. Basically the directive is used to exclude the MPI-functions for compilation with a serial compiler (on the local host). Note: the namelist file check is currently not possible for coupled runs and for [../../app/initialization_parameters#initializing_actions initializing_actions='cyclic_fill']! A new mrun option '-z' has been implemented to switch off check_namelist_files program (PALM will of course still check the namelist files when the simulation starts) (Makefile_check, mbuild, mrun, check_namelist_files, init_pegrid, check_parameters, init_masks, poisfft, poisfft_hybrid, local_stop, init_grid, check_open, modules). 1594 1594 }}} 1595 1595 |---------------- … … 1774 1774 }}} 1775 1775 {{{#!td style="vertical-align:top" 1776 Scalar quantities ({{{pt}}}, {{{q}}}, {{{sa}}}) can be excluded from Rayleigh damping by setting new parameter [../../app/ini par/#scalar_rayleigh_damping scalar_rayleigh_damping] = .FALSE.1776 Scalar quantities ({{{pt}}}, {{{q}}}, {{{sa}}}) can be excluded from Rayleigh damping by setting new parameter [../../app/initialization_parameters#scalar_rayleigh_damping scalar_rayleigh_damping] = .FALSE. 1777 1777 (init_3d_model, modules, parin, prognostic_equations) 1778 1778 }}} … … 1947 1947 }}} 1948 1948 {{{#!td style="vertical-align:top" 1949 Output of liquid water potential temperature lpt as cross section, volume and masked data in case of [../../app/ini par/#cloud_physics cloud_physics]=.T. enabled. See also [../../app/d3par/#data_output data_output] and [../../app/d3par/#data_output_masks data_output_masks]. (average_3d_data, check_parameters, data_output_2d, data_output_3d, data_output_mask, init_masks, modules, netcdf, read_3d_binary, sum_up_3d_data, write_3d_binary)1949 Output of liquid water potential temperature lpt as cross section, volume and masked data in case of [../../app/initialization_parameters#cloud_physics cloud_physics]=.T. enabled. See also [../../app/d3par/#data_output data_output] and [../../app/d3par/#data_output_masks data_output_masks]. (average_3d_data, check_parameters, data_output_2d, data_output_3d, data_output_mask, init_masks, modules, netcdf, read_3d_binary, sum_up_3d_data, write_3d_binary) 1950 1950 }}} 1951 1951 |---------------- … … 1987 1987 }}} 1988 1988 {{{#!td style="vertical-align:top" 1989 Flow field initialization with given (e.g. measured) profiles implemented. Profile data for u-,v-velocity components + respective heights are given with new inipar-parameters [../../app/ini par/#u_profile u_profile], [../../app/inipar/#v_profile v_profile], and [../../app/inipar/#uv_heights uv_heights]. Final profiles are calculated from these given profiles by linear interpolation. (check_parameters, header, init_3d_model, modules, parin)1989 Flow field initialization with given (e.g. measured) profiles implemented. Profile data for u-,v-velocity components + respective heights are given with new inipar-parameters [../../app/initialization_parameters#u_profile u_profile], [../../app/initialization_parameters#v_profile v_profile], and [../../app/initialization_parameters#uv_heights uv_heights]. Final profiles are calculated from these given profiles by linear interpolation. (check_parameters, header, init_3d_model, modules, parin) 1990 1990 }}} 1991 1991 |---------------- … … 2944 2944 '''Most important changes:''' New advection scheme added. Different number of processors and different horizontal resolution in ocean and atmosphere is now allowed. Bottom boundary conditions for u and v changed from mirror to dirichlet boundary conditions. Inflow turbulence is now defined by fluctuations around spanwise mean. 2945 2945 2946 Advection scheme ( see [../../app/ini par/#momentum_advec momentum_advec] and [../../app/inipar/#scalar_advec scalar_advec] ) :2946 Advection scheme ( see [../../app/initialization_parameters#momentum_advec momentum_advec] and [../../app/initialization_parameters#scalar_advec scalar_advec] ) : 2947 2947 2948 2948 New advection scheme added ( Wicker and Skamarock 5th order ). Turbulent fluxes are computed directly inside advection routines, removed from flow_statistics. (advec_ws) 2949 2949 2950 Check for illegal combination of ws-scheme and [../../app/ini par/#timestep_scheme timestep_scheme].2951 Check for [../../app/ini par/#topography topography] and ws-scheme.2952 Check for not cyclic boundary conditions in combination with ws-scheme and [../../app/ini par/#loop_optimization loop_optimization] = 'vector'.2953 Check for [../../app/ini par#call_psolver_at_all_substeps call_psolver_at_all_substeps] and [../../app/inipar#momentum_advec momentum_advec] = 'ws-scheme'. (check_parameters)2950 Check for illegal combination of ws-scheme and [../../app/initialization_parameters#timestep_scheme timestep_scheme]. 2951 Check for [../../app/initialization_parameters#topography topography] and ws-scheme. 2952 Check for not cyclic boundary conditions in combination with ws-scheme and [../../app/initialization_parameters#loop_optimization loop_optimization] = 'vector'. 2953 Check for [../../app/initialization_parameters#call_psolver_at_all_substeps call_psolver_at_all_substeps] and [../../app/initialization_parameters#momentum_advec momentum_advec] = 'ws-scheme'. (check_parameters) 2954 2954 2955 2955 When advection is computed with ws-scheme, turbulent fluxes are already … … 3075 3075 init_pt_anomaly, init_rankine, poismg, sor, time_integration) 3076 3076 3077 Default of [../../app/ini par/#momentum_advec momentum_advec] and [../../app/inipar/#scalar_advec scalar_advec] changed to 'ws-scheme'. (modules)3077 Default of [../../app/initialization_parameters#momentum_advec momentum_advec] and [../../app/initialization_parameters#scalar_advec scalar_advec] changed to 'ws-scheme'. (modules) 3078 3078 3079 3079 nxl-1, nxr+1, nys-1, nyn+1 replaced by nxlg, nxrg, nysg, nyng … … 3131 3131 }}} 3132 3132 {{{#!td style="vertical-align:top" 3133 Bugfix in checking for [../../app/ini par#conserve_volume_flow_mode conserve_volume_flow_mode]. (check_parameters)3133 Bugfix in checking for [../../app/initialization_parameters#conserve_volume_flow_mode conserve_volume_flow_mode]. (check_parameters) 3134 3134 3135 3135 Changed length of string run_description_header. (modules) … … 3275 3275 }}} 3276 3276 {{{#!td style="vertical-align:top" 3277 Optional barriers included in order to speed up collective operations {{{MPI_ALLTOALL}}} and {{{MPI_ALLREDUCE}}}. This feature is controlled with new initial parameter [[../../app/ini par#collective_wait|collective_wait]]. Default is .FALSE, but .TRUE. on SGI-type systems. (advec_particles, advec_s_bc, buoyancy, check_for_restart, cpu_statistics, data_output_2d, data_output_ptseries, flow_statistics,3277 Optional barriers included in order to speed up collective operations {{{MPI_ALLTOALL}}} and {{{MPI_ALLREDUCE}}}. This feature is controlled with new initial parameter [[../../app/initialization_parameters#collective_wait|collective_wait]]. Default is .FALSE, but .TRUE. on SGI-type systems. (advec_particles, advec_s_bc, buoyancy, check_for_restart, cpu_statistics, data_output_2d, data_output_ptseries, flow_statistics, 3278 3278 global_min_max, inflow_turbulence, init_3d_model, init_particles, init_pegrid, init_slope, parin, pres, poismg, set_particle_attributes, timestep, read_var_list, user_statistics, write_compressed, write_var_list) 3279 3279