| 1 | New: |
|---|
| 2 | --- |
|---|
| 3 | The plant canopy model has been extended to the prognostic equation for |
|---|
| 4 | the potential temperature and the prognostic equation for the concentration |
|---|
| 5 | of a passive scalar, as proposed by Watanabe (2004) and Shaw and Schumann |
|---|
| 6 | (1992). (prognostic_equations, plant_canopy_model) |
|---|
| 7 | |
|---|
| 8 | Within a job chain, the number of processors or the virtual processor grid |
|---|
| 9 | is now allowed to be changed. This modifies the domain decomposition and |
|---|
| 10 | requires that the PEs within a restart job have to know from which file(s) they |
|---|
| 11 | have to read their restart data (up to now there is a one to one relation |
|---|
| 12 | between the restart files and the PEs, i.e. processor 4 only reads data |
|---|
| 13 | from file "_0004"). Informations about the subdomain index limits of all |
|---|
| 14 | restart files are now written by PE0 (to file "_0000"). Also, all variables |
|---|
| 15 | containing global values which do not depend on the subdomain are written |
|---|
| 16 | to file "_0000" only. In a restart run, this file is opened by all PEs to |
|---|
| 17 | read these global values and the informations about the subdomain index |
|---|
| 18 | limits. Depending on these index limit information it is decided which |
|---|
| 19 | restart file(s) has/have to be opened in order to read the complete data |
|---|
| 20 | necessary for the current subdomain. |
|---|
| 21 | |
|---|
| 22 | The turbulence recycling method of Kataoka and Mizuno (2002) has been |
|---|
| 23 | implemented. A pre-run with smaller model domain can be used to initialize |
|---|
| 24 | the real run, in order to spare the CPU time needed to get the turbulence |
|---|
| 25 | into a quasi-stationary state (see initializing_actions = |
|---|
| 26 | 'read_data_for_recycling'). |
|---|
| 27 | |
|---|
| 28 | Particle advection adjusted for ocean runs. |
|---|
| 29 | |
|---|
| 30 | User-defined spectra. |
|---|
| 31 | |
|---|
| 32 | Argument -c introduced to most of the subjob calls, which allows the user to |
|---|
| 33 | choose his own job catalog by setting job_catalog in the configuration file |
|---|
| 34 | (default is ~/job_queue). Workaround for mpixec with -env option. |
|---|
| 35 | Adjustments for lcxt4 (Bergen Center for Computational Science) (mrun, subjob) |
|---|
| 36 | |
|---|
| 37 | advec_particles, calc_spectra, check_open, check_parameters, data_output_spectra, header, init_particles, init_pegrid, init_3d_model, modules, netcdf, parin, particle_boundary_conds, plant_canopy_model, prognostic_equations, read_var_list, read_3d_binary, time_integration, user_interface, write_var_list, write_3d_binary |
|---|
| 38 | |
|---|
| 39 | New: inflow_turbulence |
|---|
| 40 | |
|---|
| 41 | |
|---|
| 42 | Changed: |
|---|
| 43 | ------- |
|---|
| 44 | Strict grid matching along z is not needed for mg-solver. (check_parameters) |
|---|
| 45 | |
|---|
| 46 | Sharp lateral edges of the plant canopy introduced. Before smaller |
|---|
| 47 | values of the leaf area density than those applied in the prognostic |
|---|
| 48 | equation for the vertical velocity were used in the prognostic equations |
|---|
| 49 | for u and v at the boundaries of the plant_canopy. (init_3d_model) |
|---|
| 50 | |
|---|
| 51 | -cross_ts_* (read_var_list, write_var_list) |
|---|
| 52 | -dopr_time_count, old profil parameters (cross_..., dopr_crossindex, |
|---|
| 53 | profile_***) removed (read_3d_binary, write_3d_binary) |
|---|
| 54 | |
|---|
| 55 | hom_sum moved, volume_flow_area, volume_flow_initial moved from ..._3d_binary |
|---|
| 56 | to ..._var_list, hom_sum is now allocated in parin instead of init_3d_model |
|---|
| 57 | |
|---|
| 58 | npex, npey moved from inipar- to d3par-namelist (parin) |
|---|
| 59 | |
|---|
| 60 | myid_char14 removed. This was used since on IBM machines the process rank may |
|---|
| 61 | be changed when the final communicator is defined, in order to save the |
|---|
| 62 | preliminary processor-id for opening the binary output file for restarts |
|---|
| 63 | (unit 14), because otherwise a mismatch occured when reading this file in a |
|---|
| 64 | restart job. This mechanism is not needed any more since the subdomain |
|---|
| 65 | informations are now contained in file _0000. (parin, check_open) |
|---|
| 66 | |
|---|
| 67 | Transpositions for the 2D domain decomposition have been optimized by using |
|---|
| 68 | f_inv as an automatic array instead of providing the memory by a dummy argument. |
|---|
| 69 | This spares one copy loop per transposition. Order of indices in the 3D loops |
|---|
| 70 | in some of the transpose routines have been rearranged for better cache utilization. |
|---|
| 71 | Both have been suggested by Roland Richter (SGI) as part of the |
|---|
| 72 | HLRN-II benchmark process. (transpose) |
|---|
| 73 | |
|---|
| 74 | Workaround for getting information about the coupling mode. (palm) |
|---|
| 75 | |
|---|
| 76 | Format adjustments in order to allow CPU# > 999. (cpu_statistics) |
|---|
| 77 | |
|---|
| 78 | advec_s_ups, advec_u_ups, advec_v_ups, advec_w_ups, calc_spectra, check_open, check_parameters, cpu_statistics, init_3d_model, modules, palm, parin, poisfft, read_var_list, read_3d_binary, transpose, write_var_list, write_3d_binary |
|---|
| 79 | |
|---|
| 80 | |
|---|
| 81 | Errors: |
|---|
| 82 | ------ |
|---|
| 83 | Bugfix: calculation of horizontal fluxes at vertical walls (diffusion_s) |
|---|
| 84 | |
|---|
| 85 | Bugfix: zero assignments to tendency arrays in case of restarts (init_3d_model) |
|---|
| 86 | |
|---|
| 87 | Bugfix: wrong grid level used in multigrid solver, routine restrict (poismg) |
|---|
| 88 | |
|---|
| 89 | Bugfix: Construction of vertical profiles when 10 gradients have been |
|---|
| 90 | specified in the parameter list (ug, vg, pt, q, sa, lad). (check_parameters) |
|---|
| 91 | |
|---|
| 92 | Bugfix: Determination of the number of vertical levels for which spectra |
|---|
| 93 | shall be output - 100 replaced by 10. (netcdf) |
|---|
| 94 | |
|---|
| 95 | Bugfix: Lengths of those characters that are needed for the output of the |
|---|
| 96 | characteristic levels of potential temperature, the geostrophic wind, |
|---|
| 97 | scalar concentration, humidity and leaf area density to the header file |
|---|
| 98 | has been increased. (header) |
|---|
| 99 | |
|---|
| 100 | Bugfix: add ')' in user_read_restart_data. |
|---|
| 101 | |
|---|
| 102 | The leaf area density is explicitly set to its surface value at k=0. (check_parameters) |
|---|
| 103 | |
|---|
| 104 | Introduce prefix_chr to ensure unique dvrp_file path. |
|---|
| 105 | small bugfixes for user_interface sample code (comments): |
|---|
| 106 | - initialize ustvst with 0.0 as it is now computed only until nxr and nyn |
|---|
| 107 | - two ALLOCATE statements moved from user_read_restart_data back to user_init |
|---|
| 108 | - remove 'READ (13) u2_av' statement in user_read_restart_data |
|---|
| 109 | Bugfix: remove IF statement in plant_canopy_model_ij (plant_canopy_model) |
|---|
| 110 | Bugfix: divide sums(k,8) (e) and sums(k,34) (e*) by ngp_2dh_s_inner(k,sr) |
|---|
| 111 | (like other scalars) (flow_statistics) |
|---|
| 112 | Bugfix: dopr_time_count was written on the binary file, which caused that |
|---|
| 113 | NetCDF files newly created by restart files (no append of existing files!) |
|---|
| 114 | contained uneccessary time levels. (read_3d_binary, write_3d_binary) |
|---|
| 115 | Bugfix: extra '*' removed in user_statistics sample code (user_interface) |
|---|
| 116 | Bugfix: a stop command was missing in some cases of the parallel branch (local_stop) |
|---|
| 117 | Bugfix in volume flow control for non-cyclic boundary conditions (pres) |
|---|
| 118 | Bugfix: misplaced #endif directives (combine_plot_fields) |
|---|
| 119 | |
|---|
| 120 | |
|---|
| 121 | check_parameters, diffusion_s, flow_statistics, init_dvrp, init_3d_model, local_stop, plant_canopy_model, poismg, pres, read_3d_binary, user_interface, write_3d_binary |
|---|
| 122 | |
|---|
