source: palm/trunk/SOURCE/CURRENT_MODIFICATIONS @ 388

Last change on this file since 388 was 388, checked in by raasch, 12 years ago

in-situ AND potential density are calculated and used in the ocean version

  • Property svn:keywords set to Id
File size: 9.2 KB
Line 
1New:
2---
3
4Adapted for machine lck
5(mrun, mbuild, subjob)
6
7bc_lr/bc_ns  in most subroutines replaced by LOGICAL variables bc_lr_cyc,
8bc_ns_cyc for speed optimization
9(check_parameters, diffusion_u, diffusion_v, diffusion_w, modules)
10
11Additional timestep criterion in case of simulations with plant canopy (timestep)
12
13Check for illegal entries in section_xy|xz|yz that exceed nz+1|ny+1|nx+1
14(check_parameters)
15
16Clipping of dvrp output implemented. Default colourtable for particles
17implemented, particle attributes (color, dvrp_size) can be set with new
18parameters particle_color, particle_dvrpsize, color_interval,
19dvrpsize_interval (init_dvrp, data_output_dvrp, modules, user_data_output_dvrp).
20Slicer attributes (dvrp) are set with new routine set_slicer_attributes_dvrp
21and are controlled with existing parameters slicer_range_limits.
22(set_slicer_attributes_dvrp)
23
24Ocean atmosphere coupling allows to use independent precursor runs in order
25to account for different spin-up times. The time when coupling has to be
26started is given by new inipar parameter coupling_start_time. The precursor
27ocean run has to be started using new mrun option "-y" in order to add
28appendix "_O" to all output files.
29(check_for_restart, check_parameters, data_output_2d, data_output_3d,
30data_output_profiles, data_output_ptseries, data_output_spectra,
31data_output_tseries, header, init_coupling, modules, mrun,
32parin, read_var_list, surface_coupler, time_integration, write_var_list)
33
34Polygon reduction for topography and ground plate isosurface. Reduction level
35for buildings can be chosen with parameter cluster_size. (init_dvrp)
36
37External pressure gradient (check_parameters, header, init_3d_model, modules,
38parin, prognostic_equations, read_var_list, write_var_list)
39
40New topography case 'single_street_canyon' (header, init_grid, modules, parin,
41read_var_list, user_check_parameters, user_header, user_init_grid, write_var_list)
42
43Option to predefine a target bulk velocity for conserve_volume_flow
44(check_parameters, header, init_3d_model, modules, parin, read_var_list,
45write_var_list)
46
47Option for user defined 2D data output in xy cross sections at z=nzb+1
48(data_output_2d, user_data_output_2d)
49
50xy cross section output of surface heatfluxes (latent, sensible)
51(average_3d_data, check_parameters, data_output_2d, modules, read_3d_binary,
52sum_up_3d_data, write_3d_binary)
53
54
55average_3d_data, check_for_restart, check_parameters, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, init_coupling, init_dvrp, init_grid, init_3d_model, header, mbuild, modules, mrun, package_parin, parin, prognostic_equations, read_3d_binary, read_var_list, subjob, surface_coupler, timestep, time_integration, user_check_parameters, user_data_output_2d, user_data_output_dvrp, user_header, user_init_grid, write_3d_binary, write_var_list
56
57New: set_particle_attributes, set_slicer_attributes_dvrp
58
59
60Changed:
61-------
62__lcmuk changed to __lc to avoid problems with Intel compiler on sgi-ice
63(poisfft)
64
65For extended NetCDF files, the updated title attribute includes an update of
66time_average_text where appropriate. (netcdf)
67
68In case of restart runs without extension, initial profiles are not written
69to NetCDF-file anymore. (data_output_profiles, modules, read_var_list, write_var_list)
70
71Small change in formatting of the message handling routine concering the output in the
72job protocoll. (message)
73
74initializing_actions='read_data_for_recycling' renamed to 'cyclic_fill', now
75independent of turbulent_inflow (check_parameters, header, init_3d_model)
76
772 NetCDF error numbers changed. (data_output_3d)
78
79A Link to the website appendix_a.html is printed for further information
80about the possible errors. (message)
81
82Temperature gradient criterion for estimating the boundary layer height
83replaced by the gradient criterion of Sullivan et al. (1998). (flow_statistics)
84
85NetCDF unit attribute in timeseries output in case of statistic regions added
86(netcdf)
87
88Output of NetCDF messages with aid of message handling routine.
89(check_open, close_file, data_output_2d, data_output_3d,
90data_output_profiles, data_output_ptseries, data_output_spectra,
91data_output_tseries, netcdf, output_particles_netcdf)
92
93Output of messages replaced by message handling routine.
94(advec_particles, advec_s_bc, buoyancy, calc_spectra, check_for_restart,
95check_open, coriolis, cpu_log, data_output_2d, data_output_3d, data_output_dvrp,
96data_output_profiles, data_output_spectra, fft_xy, flow_statistics, header,
97init_1d_model, init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid,
98netcdf, parin, plant_canopy_model, poisfft_hybrid, poismg, read_3d_binary,
99read_var_list, surface_coupler, temperton_fft, timestep, user_actions,
100user_data_output_dvrp, user_dvrp_coltab, user_init_grid, user_init_plant_canopy,
101user_parin, user_read_restart_data, user_spectra )
102
103Maximum number of tails is calculated from maximum number of particles and
104skip_particles_for_tail (init_particles)
105
106Value of vertical_particle_advection may differ for each particle group
107(advec_particles, header, modules)
108
109First constant in array den also defined as type double. (eqn_state_seawater)
110
111Parameter dvrp_psize moved from particles_par to dvrp_graphics_par. (package_parin)
112
113topography_grid_convention moved from userpar to inipar (check_parameters,
114header, parin, read_var_list, user_check_parameters, user_header,
115user_init_grid, user_parin, write_var_list)
116
117Default value of grid_matching changed to strict.
118
119Adjustments for runs on lcxt4 (necessary due to an software update on CRAY) and
120for coupled runs on ibmy (mrun, subjob)
121
122advec_particles, advec_s_bc, buoyancy, calc_spectra, check_for_restart, check_open, check_parameters, close_file, coriolis, cpu_log, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, eqn_state_seawater, fft_xy, flow_statistics, header, init_1d_model, init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid, message, mrun, netcdf, output_particles_netcdf, package_parin, parin, plant_canopy_model, poisfft, poisfft_hybrid, poismg, read_3d_binary, read_var_list, sort_particles, subjob, user_check_parameters, user_header, user_init_grid, user_parin, surface_coupler, temperton_fft, timestep, user_actions, user_data_output_dvrp, user_dvrp_coltab, user_init_grid, user_init_plant_canopy, user_parin, user_read_restart_data, user_spectra, write_var_list
123
124
125Errors:
126------
127Bugfix: Initial hydrostatic pressure profile in case of ocean runs is now
128calculated in 5 iteration steps. (init_ocean)
129
130Bugfix: wrong sign in buoyancy production of ocean part in case of not using
131the reference density (only in 3D routine production_e) (production_e)
132
133Bugfix: output of averaged 2d/3d quantities requires that an avaraging
134interval has been set, respective error message is included (check_parameters)
135
136Bugfix: Output on unit 14 only if requested by write_binary.
137(user_last_actions)
138
139Bugfix to avoid zero division by km_neutral (production_e)
140
141Bugfix for extended NetCDF files: In order to avoid 'data mode' errors if
142updated attributes are larger than their original size, NF90_PUT_ATT is called
143in 'define mode' enclosed by NF90_REDEF and NF90_ENDDEF calls. This implies a
144possible performance loss; an alternative strategy would be to ensure equal
145attribute size in a job chain. (netcdf)
146
147Bugfix: correction of initial volume flow for non-flat topography (init_3d_model)
148Bugfix: zero initialization of arrays within buildings for 'cyclic_fill' (init_3d_model)
149
150Bugfix: to_be_resorted => s_av for time-averaged scalars (data_output_2d, data_output_3d)
151
152Bugfix: error in formatting the output (message)
153
154Bugfix: avoid that ngp_2dh_s_inner becomes zero (init_3_model)
155
156Typographical error: unit of wpt in dots_unit (modules)
157
158Bugfix: error in check, if particles moved further than one subdomain length.
159This check must not be applied for newly released particles. (advec_particles)
160
161Bugfix: several tail counters are initialized, particle_tail_coordinates is
162only written to file if its third index is > 0, arrays for tails are allocated
163with a minimum size of 10 tails if there is no tail initially (init_particles,
164advec_particles)
165
166Bugfix: pressure included for profile output (check_parameters)
167
168Bugfix: Type of count and count_rate changed to default INTEGER on NEC machines
169(cpu_log)
170
171Bugfix: output if particle time series only if particle advection is switched
172on. (time_integration)
173
174Bugfix: qsws was calculated in case of constant heatflux = .FALSE. (prandtl_fluxes)
175
176Bugfix: averaging along z is not allowed for 2d quantities (e.g. u* and z0) (data_output_2d)
177
178Typographical errors (netcdf)
179
180If the inversion height calculated by the prerun is zero, inflow_damping_height
181must be explicitly specified (init_3d_model)
182
183Small bugfix concerning 3d 64bit netcdf output format (header)
184
185Bugfix: dt_fixed removed from the restart file, because otherwise, no change
186from a fixed to a variable timestep would be possible in restart runs.
187(read_var_list, write_var_list)
188
189Bugfix: initial setting of time_coupling in coupled restart runs (time_integration)
190
191advec_particles, check_parameters, cpu_log, data_output_2d, data_output_3d, header, init_3d_model, init_particles, init_ocean, modules, netcdf, prandtl_fluxes, production_e, read_var_list, time_integration, user_last_actions, write_var_list
Note: See TracBrowser for help on using the repository browser.