| 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 | 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 |
|---|
| 33 | |
|---|
| 34 | New: inflow_turbulence |
|---|
| 35 | |
|---|
| 36 | |
|---|
| 37 | Changed: |
|---|
| 38 | ------- |
|---|
| 39 | Sharp lateral edges of the plant canopy introduced. Before smaller |
|---|
| 40 | values of the leaf area density than those applied in the prognostic |
|---|
| 41 | equation for the vertical velocity were used in the prognostic equations |
|---|
| 42 | for u and v at the boundaries of the plant_canopy. (init_3d_model) |
|---|
| 43 | |
|---|
| 44 | -cross_ts_* (read_var_list, write_var_list) |
|---|
| 45 | -dopr_time_count, old profil parameters (cross_..., dopr_crossindex, |
|---|
| 46 | profile_***) removed (read_3d_binary, write_3d_binary) |
|---|
| 47 | |
|---|
| 48 | hom_sum moved, volume_flow_area, volume_flow_initial moved from ..._3d_binary |
|---|
| 49 | to ..._var_list, hom_sum is now allocated in parin instead of init_3d_model |
|---|
| 50 | |
|---|
| 51 | npex, npey moved from inipar- to d3par-namelist (parin) |
|---|
| 52 | |
|---|
| 53 | myid_char14 removed. This was used since on IBM machines the process rank may |
|---|
| 54 | be changed when the final communicator is defined, in order to save the |
|---|
| 55 | preliminary processor-id for opening the binary output file for restarts |
|---|
| 56 | (unit 14), because otherwise a mismatch occured when reading this file in a |
|---|
| 57 | restart job. This mechanism is not needed any more since the subdomain |
|---|
| 58 | informations are now contained in file _0000. (parin, check_open) |
|---|
| 59 | |
|---|
| 60 | check_open, init_3d_model, modules, parin, read_var_list, read_3d_binary, write_var_list, write_3d_binary |
|---|
| 61 | |
|---|
| 62 | |
|---|
| 63 | Errors: |
|---|
| 64 | ------ |
|---|
| 65 | Bugfix: add ')' in user_read_restart_data. |
|---|
| 66 | |
|---|
| 67 | The leaf area density is explicitly set to its surface value at k=0. (check_parameters) |
|---|
| 68 | |
|---|
| 69 | Introduce prefix_chr to ensure unique dvrp_file path. |
|---|
| 70 | small bugfixes for user_interface sample code (comments): |
|---|
| 71 | - initialize ustvst with 0.0 as it is now computed only until nxr and nyn |
|---|
| 72 | - two ALLOCATE statements moved from user_read_restart_data back to user_init |
|---|
| 73 | - remove 'READ (13) u2_av' statement in user_read_restart_data |
|---|
| 74 | Bugfix: remove IF statement in plant_canopy_model_ij (plant_canopy_model) |
|---|
| 75 | Bugfix: divide sums(k,8) (e) and sums(k,34) (e*) by ngp_2dh_s_inner(k,sr) |
|---|
| 76 | (like other scalars) (flow_statistics) |
|---|
| 77 | Bugfix: dopr_time_count was written on the binary file, which caused that |
|---|
| 78 | NetCDF files newly created by restart files (no append of existing files!) |
|---|
| 79 | contained uneccessary time levels. (read_3d_binary, write_3d_binary) |
|---|
| 80 | Bugfix: extra '*' removed in user_statistics sample code (user_interface) |
|---|
| 81 | Bugfix: a stop command was missing in some cases of the parallel branch (local_stop) |
|---|
| 82 | Bugfix in volume flow control for non-cyclic boundary conditions (pres) |
|---|
| 83 | |
|---|
| 84 | |
|---|
| 85 | check_parameters, flow_statistics, init_dvrp, local_stop, plant_canopy_model, pres, read_3d_binary, user_interface, write_3d_binary |
|---|
| 86 | |
|---|
