Changes between Version 51 and Version 52 of doc/tec/changelog_2018

Timestamp:

Jan 3, 2011 3:29:53 PM (14 years ago)

Author:

suehring

Comment:

--

doc/tec/changelog_2018

@@ -19 +19 @@
 }}}
 {{{#!td style="vertical-align:top;width: 75px"
-r667
+r670
 }}}
 {{{#!td style="vertical-align:top"
@@ -28 +28 @@
 }}}
 {{{#!td style="vertical-align:top"
-New advection scheme added ( Wicker and Skamarock 5th order ) (advec_ws.f90)
-Please note the further information in [../../app/inipar/#momentum_advec inipar] .
+'''Most important changes:''' New advection scheme added ( Wicker and Skamarock 5th order ). For further information
+see [../../app/inipar/#momentum_advec inipar]. (advec_ws)
 
 Changes in coupling: Different number of processors and different horizontal resolution in ocean and atmosphere is now allowed.
@@ -36 +36 @@
 Bottom boundary conditions for u and v  changed from mirror to dirichlet boundary conditions.
 
-Detailed information will follow at beginning of January 2011.
-
-
+Inflow turbulence is now defined by fluctuations around spanwise mean.
+
+Removed mirror boundary conditions for u and v at the bottom in case of
+ibc_uv_b == 0. Instead, dirichlet boundary conditions (u=v=0) are set
+in init_3d_model. (boundary_conds)
+
+Exchange of parameters between ocean and atmosphere via PE0
+Check for illegal combination of ws-scheme and [../../app/inipar/#timestep_scheme timestep_scheme].
+Check for [../../app/inipar/#topography topography] and ws-scheme.
+Check for not cyclic boundary conditions in combination with ws-scheme and [../../app/inipar/#loop_optimization loop_optimization] = 'vector'.
+Check for [../../app/inipar#call_psolver_at_all_substeps call_psolver_at_all_substeps] and [../../app/inipar#momentum_advec momentum_advec] = 'ws-scheme'.
+Different processor/grid topology in atmosphere and ocean is now allowed! (check_parameters)
+
+Dynamic exchange of ghost points with nbgp_local to ensure that no useless 
+ghost points exchanged in case of multigrid. type_yz(0) and type_xz(0) used for
+normal grid, the remaining types used for the several grid levels. 
+Exchange is done via MPI-Vectors with a dynamic value of ghost points which 
+depend on the advection scheme. Exchange of left and right PEs is 10% faster
+with MPI-Vectors than without. (exchange_horiz, exchange_horiz_2d)
+
+When advection is computed with ws-scheme, turbulent fluxes are already 
+computed in the respective advection routines and buffered in arrays
+sums_xxxx_ws_l(). This is due to a consistent treatment of statistics 
+with the numerics and to avoid unphysical kinks near the surface. Furthermore the computation of z_i is only done if the heat flux 
+exceeds a minimum value. This affects only simulations of a neutral boundary 
+layer and is due to reasons of computations in the advection scheme. 
+(flow_statistics)
+
+Using nbgp recycling planes for a better resolution of the turbulent flow near
+the inflow. (inflow_turbulence)
+
+Definition of new array bounds nxlg, nxrg, nylg, nyng on each PE. 
+Furthermore the allocation of arrays and steering of loops is done with these
+parameters. In case of dirichlet boundary condition at the bottom zu(0)=0.0
+dzu_mg has to be set explicitly for a equally spaced grid near bottom.
+ddzu_pres added to use a equally spaced grid near bottom. (init_grid)
+
+Moved determination of target_id's from init_coupling.
+Determination of parameters needed for coupling (coupling_topology, ngp_a, ngp_o)
+with different grid/processor-topology in ocean and atmosphere
+Adaption of ngp_xy, ngp_y to a dynamic number of ghost points.
+The maximum_grid_level changed from 1 to 0. 0 is the normal grid, 1 to 
+maximum_grid_level the grids for multigrid, in which 0 and 1 are normal grids.
+This distinction is due to reasons of data exchange and performance for the 
+normal grid and grids in poismg.
+The definition of MPI-Vectors adapted to a dynamic number of ghost points.
+New MPI-Vectors for data exchange between left and right boundaries added. 
+This is due to reasons of performance (10% faster).
+ATTENTION: nnz_x undefined problem still has to be solved!!!!!!!!
+TEST OUTPUT (TO BE REMOVED) logging mpi2 ierr values. (init_pegrid)
+
+Call ws_init to initialize arrays needed for statistical evaluation and
+optimization when ws-scheme is used. (init_3d_model)
+
+Module advec_ws added. (Makefile)
+
+Removed u_nzb_p1_for_vfc and v_nzb_p1_for_vfc
+For coupling with different resolution in ocean and atmosphere:
++nx_a, +nx_o, +ny_a, +ny_o, +ngp_a, +ngp_o, +total_2d_o, +total_2d_a,
++coupling_topology 
+Buffer arrays for the left sided advective fluxes added in arrays_3d.
++flux_s_u, +flux_s_v, +flux_s_w, +diss_s_u, +diss_s_v, +diss_s_w,
++flux_s_pt, +diss_s_pt, +flux_s_e, +diss_s_e, +flux_s_q, +diss_s_q,
++flux_s_sa, +diss_s_sa.
++exchange_mg added in control_parameters to steer the data exchange.
+Parameters +nbgp, +nxlg, +nxrg, +nysg, +nyng added in indices.
+flag array +boundary_flags added in indices to steer the degradation of order 
+of the advective fluxes when non-cyclic boundaries are used.
+MPI-datatypes +type_y, +type_y_int and +type_yz for data_exchange added in 
+pegrid.
++sums_wsus_ws_l, +sums_wsvs_ws_l, +sums_us2_ws_l, +sums_vs2_ws_l, 
++sums_ws2_ws_l, +sums_wspts_ws_l, +sums_wssas_ws_l, +sums_wsqs_ws_l 
+and +weight_substep added in statistics to steer the statistical evaluation 
+of turbulent fluxes in the advection routines.
+LOGICALS +ws_scheme_sca and +ws_scheme_mom added to get a better performance
+in prognostic_equations. (modules)
+
+New allocation of tend when ws-scheme and multigrid is used. This is due to 
+reasons of performance of the data_exchange. The same is done with p after
+poismg is called.
+nxl-1, nxr+1, nys-1, nyn+1 replaced by nxlg, nxrg, nysg, nyng when no 
+multigrid is used.
+Call of SOR routine is referenced with ddzu_pres. (pres)
+
+Calls of the advection routines with WS5 added.
+Calls of ws_statistics added to set the statistical arrays to zero after each 
+time step. (prognostic_equations)
+
+Additional case for nonequivalent processor and grid topology in ocean and 
+atmosphere added (coupling_topology = 1).
+Added exchange of u and v from Ocean to Atmosphere. (surface_coupler)
+
+
+
+
+}}}
+|----------------
+{{{#!td style="vertical-align:top;width: 50px"
+
+}}}
+{{{#!td style="vertical-align:top;width: 50px"
+
+}}}
+{{{#!td style="vertical-align:top;width: 75px"
+
+}}}
+{{{#!td style="vertical-align:top"
+
+}}}
+{{{#!td style="vertical-align:top"
+C
+}}}
+{{{#!td style="vertical-align:top"
+Default of [../../app/inipar/#momentum_advec momentum_advec] and [../../app/inipar/#scalar_advec scalar_advec] changed to 'ws-scheme'. (modules)
+
+nxl-1, nxr+1, nys-1, nyn+1 replaced by nxlg, nxrg, nysg, nyng 
+(advec_particles, asselin_filter.f90, average_3d_data.f90, boundary_conds.f90, 
+calc_liquid_water_content.f90, calc_spectra, data_output_2d, data_output_3d,
+diffusion_e.f90, diffusion_u.f90, diffusion_v.f90, diffusion_w.f90, 
+diffusivities, init_3d_model, init_rankine, init_slope, prandtl_fluxes, pres,
+production_e, sor, subsidence, sum_up_3d_data, surface_coupler, timestep, 
+user_..., wall_fluxes)
+
+Calls of exchange_horiz modified. (advec_particles, data_output_2d, 
+data_output_3d, data_output_mask, diffusivities, init_3d_model, 
+init_pt_anomaly, init_rankine, poismg, sor, time_integration)
+
+Allocation of arrays, output of array bounds and skip values adapted to a 
+dynamic number of ghost points with nbgp. (check_open, data_output_3d,
+global_min_max, read_3d_binary, user_read_restart_data, write_compressed)
+
+Changed surface boundary conditions for u and v from mirror bc to dirichelt bc,
+therefore u(uzb,:,:) and v(nzb,:,:) is now representative for the height z0. 
+(prandtl_fluxes)
+
+determination of target_id's moved to init_pegrid. (init_coupling)
+
+Initial volume flow is now calculated by using the variable hom_sum.
+Therefore the correction of initial volume flow for non-flat topography
+removed (removed u_nzb_p1_for_vfc and v_nzb_p1_for_vfc)
+Changed surface boundary conditions for u and v in case of ibc_uv_b == 0 from
+mirror bc to dirichlet boundary conditions (u=v=0), so that k=nzb is
+representative for the height z0. (init_3d_model)
+}}}
+|----------------
+{{{#!td style="vertical-align:top;width: 50px"
+
+}}}
+{{{#!td style="vertical-align:top;width: 50px"
+
+}}}
+{{{#!td style="vertical-align:top;width: 75px"
+
+}}}
+{{{#!td style="vertical-align:top"
+
+}}}
+{{{#!td style="vertical-align:top"
+B
+}}}
+{{{#!td style="vertical-align:top"
+Bugfix in checking for [../../app/inipar#conserve_volume_flow_mode conserve_volume_flow_mode]. (check_parameters)
+
+Changed length of string run_description_header. (modules)
+
+Bugfix: After pressure correction no volume flow correction in case of 
+non-cyclic boundary conditions (has to be done only before 
+pressure correction). (pres)
+
+Bug removed in declaration of ddzw(), nz replaced by nzt+1. (sor)
 }}}
 |----------------