Changes between Version 29 and Version 30 of doc/app/examples/cbl
- Timestamp:
- Feb 24, 2021 3:10:24 PM (4 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
doc/app/examples/cbl
v29 v30 8 8 Fast Fourier transformations ([../../initialization_parameters#fft_method fft_method]) are calculated using the Temperton-algorithm. The initial profiles for wind and temperature can be assigned via [../../initialization_parameters#initializing_actions initializing_actions]'' = 'set_constant_profiles'.'' The wind speed, constant with height, amounts to [../../initialization_parameters#ug_surface ug_surface]'' = ''[../../initialization_parameters#vg_surface vg_surface]'' = 0.0 m/s.'' In order to allow for a fast onset of convection, a neutral stratified layer up to ''z = 800 m'' capped by an inversion with ''dtheta/dz = 1K/100 m'' is given: [../../initialization_parameters#pt_vertical_gradient pt_vertical_gradient]'' = 0.0, 1.0,'' [../../initialization_parameters#pt_vertical_gradient_level pt_vertical_gradient_level]'' = 0.0, 800.0.'' The surface temperature, which by default amounts to ''300 K'', provides the fixed point for the temperature profile (see [../../initialization_parameters#pt_surface pt_surface]). Convection is driven by a given, near-surface sensible heat flux via [../../initialization_parameters#surface_heatflux surface_heatflux]'' = 0.1 K m/s''. A given surface sensible heta flux requires the bottom boundary condition for potential temperature to be [../../initialization_parameters#bc_pt_b bc_pt_b]'' = 'neumann'.'' Thus all initialization parameters are determined. These can not be changed during the run (also not for restart runs). 9 9 10 Now the run parameters ([../../ d3par &runtime_parameters] ) must be specified. To produce a quasi stationary boundary layer the simulated time should be at least one hour, i.e. [../../d3par#end_time end_time]'' = 3600 s.'' To stimulate convection, the initially homogeneous (zero) wind field must be disturbed ([../../d3par#create_disturbances create_disturbances]'' = .T.''). These perturbations should be repeated in a temporal interval of [../../d3par#dt_disturb dt_disturb]'' = 150.0 s'' until the energy of the perturbations exceeds the value [../../d3par#disturbance_energy_limit disturbance_energy_limit]'' = 0.01 m^2^/s^2^''. After each time step run time information (e.g. size of the timestep, maximum velocities, etc.) are to be written to the local file [../../iofiles#RUN_CONTROL RUN_CONTROL] ([../../d3par#dt_run_control dt_run_control]'' = 0.0 s'').10 Now the run parameters ([../../runtime_parameters &runtime_parameters] ) must be specified. To produce a quasi stationary boundary layer the simulated time should be at least one hour, i.e. [../../runtime_parameters#end_time end_time]'' = 3600 s.'' To stimulate convection, the initially homogeneous (zero) wind field must be disturbed ([../../runtime_parameters#create_disturbances create_disturbances]'' = .T.''). These perturbations should be repeated in a temporal interval of [../../runtime_parameters#dt_disturb dt_disturb]'' = 150.0 s'' until the energy of the perturbations exceeds the value [../../runtime_parameters#disturbance_energy_limit disturbance_energy_limit]'' = 0.01 m^2^/s^2^''. After each time step run time information (e.g. size of the timestep, maximum velocities, etc.) are to be written to the local file [../../iofiles#RUN_CONTROL RUN_CONTROL] ([../../runtime_parameters#dt_run_control dt_run_control]'' = 0.0 s''). 11 11 12 Instantaneous cross section data of vertical velocity (''w'') and potential temperature (''pt'') are to be output for horizontal (''xy'') and vertical (''xz'') cross sections, and additionally, time averaged (''av'') vertical cross section data are to be output for the vertical velocity: [../../ d3par#data_output data_output]'' = 'w_xy', 'w_xz', 'w_xz_av', 'theta_xy', 'theta_xz'.'' Output of instantaneous (time averaged) data is done after each ''900 (1800)s'': [../../d3par#dt_data_output dt_data_output]'' = 900.0,'' [../../d3par#dt_data_output_av dt_data_output_av]'' = 1800.0.'' The averaged data are time averaged over the last ''900.0 s'', where the temporal interval of data entering the average is ''10 s'': [../../d3par#averaging_interval averaging_interval]'' = 900.0,'' [../../d3par#dt_averaging_input dt_averaging_input]'' = 10.0.'' Horizontal cross sections are output for vertical levels with grid index ''k=2'' and ''k=10,'' vertical cross sections are output for index ''j=20'': [../../d3par#section_xy section_xy]'' = 2, 10, '' [../../d3par#section_xz section_xz]'' = 20.'' For runs on more than one processor, cross section data are collected and output on PE0: [../../d3par#data_output_2d_on_each_pe data_output_2d_on_each_pe]'' = .F.''.12 Instantaneous cross section data of vertical velocity (''w'') and potential temperature (''pt'') are to be output for horizontal (''xy'') and vertical (''xz'') cross sections, and additionally, time averaged (''av'') vertical cross section data are to be output for the vertical velocity: [../../runtime_parameters#data_output data_output]'' = 'w_xy', 'w_xz', 'w_xz_av', 'theta_xy', 'theta_xz'.'' Output of instantaneous (time averaged) data is done after each ''900 (1800)s'': [../../runtime_parameters#dt_data_output dt_data_output]'' = 900.0,'' [../../runtime_parameters#dt_data_output_av dt_data_output_av]'' = 1800.0.'' The averaged data are time averaged over the last ''900.0 s'', where the temporal interval of data entering the average is ''10 s'': [../../runtime_parameters#averaging_interval averaging_interval]'' = 900.0,'' [../../runtime_parameters#dt_averaging_input dt_averaging_input]'' = 10.0.'' Horizontal cross sections are output for vertical levels with grid index ''k=2'' and ''k=10,'' vertical cross sections are output for index ''j=20'': [../../runtime_parameters#section_xy section_xy]'' = 2, 10, '' [../../runtime_parameters#section_xz section_xz]'' = 20.'' For runs on more than one processor, cross section data are collected and output on PE0: [../../runtime_parameters#data_output_2d_on_each_pe data_output_2d_on_each_pe]'' = .F.''. 13 13 14 Output of vertical profiles is to be done after each ''900 s''. The profiles shall be temporally averaged over the last ''600 s'', whereby the temporal interval of the profiles entering the average has to be ''10 s'': [../../ d3par#dt_dopr dt_dopr]'' = 900.0 s,'' [../../d3par#averaging_interval_pr averaging_interval_pr]'' = 600.0 s'', [../../d3par#dt_averaging_input_pr dt_averaging_input_pr]'' = 10.0 s.'' The temperature profile including the initial temperature profile (therefore'' '#theta' ''), the subgrid scale, resolved and total vertical sensible heat flux as well as the variances of the vertical velocity and the potential temperature are to be output: [../../d3par#data_output_pr data_output_pr]'' = '#theta', 'w"theta”', 'w*theta*', 'wtheta', 'w*2', 'theta*2'.''14 Output of vertical profiles is to be done after each ''900 s''. The profiles shall be temporally averaged over the last ''600 s'', whereby the temporal interval of the profiles entering the average has to be ''10 s'': [../../runtime_parameters#dt_dopr dt_dopr]'' = 900.0 s,'' [../../runtime_parameters#averaging_interval_pr averaging_interval_pr]'' = 600.0 s'', [../../runtime_parameters#dt_averaging_input_pr dt_averaging_input_pr]'' = 10.0 s.'' The temperature profile including the initial temperature profile (therefore'' '#theta' ''), the subgrid scale, resolved and total vertical sensible heat flux as well as the variances of the vertical velocity and the potential temperature are to be output: [../../runtime_parameters#data_output_pr data_output_pr]'' = '#theta', 'w"theta”', 'w*theta*', 'wtheta', 'w*2', 'theta*2'.''