| 1 | New: |
|---|
| 2 | --- |
|---|
| 3 | |
|---|
| 4 | Allow new case bc_uv_t = 'dirichlet_0' for channel flow. |
|---|
| 5 | |
|---|
| 6 | For unknown variables (CASE DEFAULT) call new subroutine user_data_output_dvrp |
|---|
| 7 | |
|---|
| 8 | Pressure boundary conditions for vertical walls added to the multigrid solver. |
|---|
| 9 | They are applied using new wall flag arrays (wall_flags_..) which are defined |
|---|
| 10 | for each grid level. New argument gls added to routine user_init_grid |
|---|
| 11 | (user_interface). |
|---|
| 12 | |
|---|
| 13 | Frequence of sorting particles can be controlled with new particles_par |
|---|
| 14 | parameter dt_sort_particles. Sorting is moved from the SGS timestep loop in |
|---|
| 15 | advec_particles after the end of this loop. |
|---|
| 16 | |
|---|
| 17 | advec_particles, check_parameters, data_output_dvrp, header, init_3d_model, init_grid, init_particles, init_pegrid, modules, package_parin, read_var_list, user_interface, write_var_list |
|---|
| 18 | |
|---|
| 19 | |
|---|
| 20 | Changed: |
|---|
| 21 | ------- |
|---|
| 22 | |
|---|
| 23 | Redefine initial nzb_local as the actual total size of topography (later the |
|---|
| 24 | extent of topography in nzb_local is reduced by 1dx at the E topography walls |
|---|
| 25 | and by 1dy at the N topography walls to form the basis for nzb_s_inner); |
|---|
| 26 | for consistency redefine 'single_building' case. |
|---|
| 27 | |
|---|
| 28 | Vertical profiles now based on nzb_s_inner; they are divided by |
|---|
| 29 | ngp_2dh_s_inner (scalars, procucts of scalars) and ngp_2dh (staggered velocity |
|---|
| 30 | components and their products, procucts of scalars and velocity components), |
|---|
| 31 | respectively. |
|---|
| 32 | |
|---|
| 33 | Allow two instead of one digit to specify isosurface and slicer variables. |
|---|
| 34 | |
|---|
| 35 | Status of 3D-volume NetCDF data file only depends on switch netcdf_64bit_3d (check_open) |
|---|
| 36 | |
|---|
| 37 | prognostic_equations include the respective wall_*flux in the parameter list of |
|---|
| 38 | calls of diffusion_s. Same as before, only the values of wall_heatflux(0:4) |
|---|
| 39 | can be assigned. At present, wall_humidityflux, wall_qflux, wall_salinityflux, |
|---|
| 40 | and wall_scalarflux are kept zero. diffusion_s uses the respective wall_*flux |
|---|
| 41 | instead of wall_heatflux. This update serves two purposes: |
|---|
| 42 | - it avoids errors in calculations with humidity/scalar/salinity and prescribed |
|---|
| 43 | non-zero wall_heatflux, |
|---|
| 44 | - it prepares PALM for a possible assignment of wall fluxes of |
|---|
| 45 | humidity/scalar/salinity in a future release. |
|---|
| 46 | |
|---|
| 47 | buoyancy, check_open, data_output_dvrp, diffusion_s, diffusivities, flow_statistics, header, init_3d_model, init_dvrp, init_grid, modules, prognostic_equations |
|---|
| 48 | |
|---|
| 49 | |
|---|
| 50 | Errors: |
|---|
| 51 | ------ |
|---|
| 52 | |
|---|
| 53 | Several bugfixes in the ocean part: Initial density rho is calculated |
|---|
| 54 | (init_ocean). Error in initializing u_init and v_init removed |
|---|
| 55 | (check_parameters). Calculation of density flux now starts from |
|---|
| 56 | nzb+1 (production_e). |
|---|
| 57 | |
|---|
| 58 | Bugfix: pleft/pright changed to pnorth/psouth in sendrecv of particle tail |
|---|
| 59 | numbers along y, small bugfixes in the SGS part (advec_particles) |
|---|
| 60 | |
|---|
| 61 | Bugfix: model_string needed a default value (combine_plot_fields) |
|---|
| 62 | |
|---|
| 63 | Bugfix: wavenumber calculation for even nx in routines maketri (poisfft) |
|---|
| 64 | |
|---|
| 65 | Bugfix: assignment of fluxes at walls |
|---|
| 66 | |
|---|
| 67 | Bugfix: absolute value of f must be used when calculating the Blackadar mixing length (init_1d_model) |
|---|
| 68 | |
|---|
| 69 | advec_particles, check_parameters, combine_plot_fields, diffusion_s, init_ocean, init_1d_model, poisfft, production_e |
|---|
