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