source: palm/trunk/SOURCE/CURRENT_MODIFICATIONS @ 108

Last change on this file since 108 was 108, checked in by letzel, 17 years ago
  • Improved coupler: evaporation - salinity-flux coupling for humidity = .T.,

avoid MPI hangs when coupled runs terminate, add DOC/app/chapter_3.8;

  • Optional calculation of km and kh from initial TKE e_init;
  • Default initialization of km,kh = 0.00001 for ocean = .T.;
  • Allow data_output_pr= q, wq, w"q", w*q* for humidity = .T.;
  • Bugfix: Rayleigh damping for ocean fixed.
  • Property svn:keywords set to Id
File size: 3.2 KB
Line 
1New:
2---
3Allows runs for a coupled atmosphere-ocean LES,
4coupling frequency is controlled by new d3par-parameter dt_coupling,
5the coupling mode (atmosphere_to_ocean or ocean_to_atmosphere) for the
6respective processes is read from environment variable coupling_mode,
7which is set by the mpiexec-command,
8communication between the two models is done using the intercommunicator
9comm_inter,
10local files opened by the ocean model get the additional suffic "_O".
11Assume saturation at k=nzb_s_inner(j,i) for atmosphere coupled to ocean.
12
13A momentum flux can be set as top boundary condition using the new
14inipar parameter top_momentumflux_u|v.
15
16Non-cyclic boundary conditions can be used along all horizontal directions.
17
18Quantities w*p* and w"e can be output as vertical profiles.
19
20Initial profiles are reset to constant profiles in case that initializing_actions /= 'set_constant_profiles'. (init_rankine)
21
22Optionally calculate km and kh from initial TKE e_init.
23
24boundary_conds, check_open, check_parameters, diffusion_u, diffusion_v, flow_statistics, header, init_pegrid, init_rankine, init_3d_model, modules, palm, parin, pres, prandtl_fluxes, prognostic_equations, read_var_list, read_3d_binary, swap_timelevel, time_integration, write_var_list, write_3d_binary
25
26New:
27surface_coupler
28
29
30Changed:
31-------
32Remaining variables iran changed to iran_part (advec_particles, init_particles).
33
34In case that the presure solver is not called for every Runge-Kutta substep
35(call_psolver_at_all_substeps = .F.), it is called after the first substep
36instead of the last. In that case, random perturbations are also added to the
37velocity field after the first substep.
38
39Initialization of km,kh = 0.00001 for ocean = .T. (for ocean = .F. it remains 0.01).
40
41Allow data_output_pr= q, wq, w"q", w*q* for humidity = .T. (instead of cloud_physics = .T.).
42
43Modifications to terminate coupled runs.
44
45advec_particles, check_for_restart, check_parameters, init_particles, init_3d_model, local_stop, timestep, time_integration
46
47
48Errors:
49------
50Bugs from code parts for non-cyclic boundary conditions are removed: loops for
51u and v are starting from index nxlu, nysv, respectively. The radiation boundary
52condition is used for every Runge-Kutta substep. Velocity phase speeds for
53the radiation boundary conditions are calculated for the first Runge-Kutta
54substep only and reused for the further substeps. New arrays c_u, c_v, and c_w
55are defined for this purpose. Several index errors are removed from the
56radiation boundary condition code parts. Upper bounds for calculating
57u_0 and v_0 (in production_e) are nxr+1 and nyn+1 because otherwise these
58values are not available in case of non-cyclic boundary conditions.
59
60+dots_num_palm in module user, +module netcdf_control in user_init (both in user_interface)
61
62Bugfix: wrong sign removed from the buoyancy production term in the case use_reference = .T. (production_e)
63
64Bugfix: Error message concerning output of particle concentration (pc) modified (check_parameters).
65
66Bugfix: Rayleigh damping for ocean fixed.
67
68advec_u_pw, advec_u_up, advec_v_pw, advec_v_up, boundary_conds, buoyancy, check_parameters, coriolis, diffusion_u, diffusion_v, header, init_pegrid, init_3d_model, modules, production_e, prognostic_equations, user_interface
Note: See TracBrowser for help on using the repository browser.