Changes between Version 7 and Version 8 of doc/app/examples/cbl

Timestamp:

Sep 16, 2010 11:56:36 AM (14 years ago)

Author:

suehring

Comment:

--

doc/app/examples/cbl

@@ -15 +15 @@
   [../../inipar#surface_heatflux surface_heatflux]{{{ = 0.1, }}} [../../inipar#bc_pt_b bc_pt_b]{{{ = 'neumann', / }}} \\
 
-[../../d3par &d3par]   [../../d3par#end_time end_time]{{{ = 3600.0,}}} \\
- create_disturbances = .T.,
- dt_disturb = 150.0, disturbance_energy_limit = 0.01,
- dt_run_control = 0.0,
- data_output = 'w_xy', 'w_xz', 'w_xz_av', 'pt_xy', 'pt_xz',
- dt_data_output = 900.0,
- dt_data_output_av = 1800.0,
- averaging_interval = 900.0,
- dt_averaging_input = 10.0,
- section_xy = 2, 10, section_xz = 20,
- data_output_2d_on_each_pe = .F.,
- dt_dopr = 900.0, averaging_interval_pr = 600.0,
- dt_averaging_input_pr = 10.0,
- data_output_pr = '#pt', 'w”pt”', 'w*pt*', 'wpt', 'w*2', 'pt*2',
- cross_profiles = ' pt ', ' w"pt" w*pt* wpt ', ' w*2 ', ' pt*2 ',
- cross_xtext = 'pot. temperature in K',
- 'heat flux in K ms>->1',
- 'velocity variance in m>2s>->2',
- 'temperature variance in K>2',
- z_max_do1d = 1500.0, /
+[../../d3par &d3par] 
+  [../../d3par#end_time end_time]{{{ = 3600.0,}}} \\
+  [../../d3par#create_disturbances create_disturbances]{{{ = .T., }}} \\
+  [../../d3par#dt_disturb dt_disturb]{{{ = 150.0, }}}  [../../d3par#disturbance_energy_limit disturbance_energy_limit]{{{ = 0.01, }}} \\
+  [../../d3par#dt_run_control dt_run_control]{{{ = 0.0,}}} \\
+  [../../d3par#data_output data_output]{{{ = 'w_xy', 'w_xz', 'w_xz_av', 'pt_xy', 'pt_xz', }}} \\
+  [../../d3par#dt_data_output dt_data_output]{{{ = 900.0, }}} \\
+  [../../d3par#dt_data_output_av  dt_data_output_av]{{{ = 1800.0, }}} \\
+  [../../d3par#averaging_interval averaging_interval]{{{ = 900.0, }}} \\
+  [../../d3par#dt_averaging_input dt_averaging_input]{{{ = 10.0, }}} \\
+  [../../d3par#section_xy section_xy]{{{ = 2, 10, }}}  [../../d3par#section_xz section_xz]{{{ = 20, }}} \\
+  [../../d3par#data_output_2d_on_each_pe data_output_2d_on_each_pe]{{{ = .F., }}} \\
+  [../../d3par#dt_dopr dt_dopr]{{{ = 900.0, [../../d3par#averaging_interval_pr averaging_interval_pr]{{{ = 600.0,}}} \\
+  [../../d3par#dt_averaging_input_pr dt_averaging_input_pr]{{{ = 10.0, }}} \\
+  [../../d3par#data_output_pr data_output_pr]{{{ = '#pt', 'w”pt”', 'w*pt*', 'wpt', 'w*2', 'pt*2', }}} \\
+  [../../d3par#z_max_do1d z_max_do1d]{{{ = 1500.0, / }}} \\
 
-The initialization parameters (&inipar) are located at the beginning of the file. For analysis of a convective boundary layer of approx. 1000 m thickness the horizontal size of the model domain should amount to at least 2 km x 2 km. In order to resolve the convective structures a grid spacing of dx = dy = dz = 50 m is enough, since the typical diameter of convective plumes is more than 100 m. Thereby the upper array index in the two horizontal directions needs to be nx = ny = 39. Since in each case the lower array index has the value 0, 40 grid points are used along both horizontal directions. In the vertical direction the domain must be high enough to include the entrainment processes at the top of the boundary layer as well as the propagation of gravity waves, which were stimulated by the convection. However, in the stably stratified region the grid resolution has not necessarily to be as high as within the boundary layer. This can be obtained by a vertical stretching of the grid starting from 1200 m via dz_stretch_level = 1200.0 m. This saves grid points and computing time. The upper boundary of the model is located at (see dz_stretch_factor) … m (computed by the model).
+The initialization parameters ([../../inipar &inipar]) are located at the beginning of the file. For analysis of a convective boundary layer of approx. 1000 m thickness the horizontal size of the model domain should amount to at least 2 km x 2 km. In order to resolve the convective structures a grid spacing of [../../inipar#dx dx]{{{ = }}} [../../inipar#dy dy]{{{ = }}}[../../inipar#dz dz]{{{ = 50 m}}} is enough, since the typical diameter of convective plumes is more than 100 m. Thereby the upper array index in the two horizontal directions needs to be [../../inipar#nx nx]{{{ = }}}[../../inipar#ny ny]{{{ = 39}}}. Since in each case the lower array index has the value 0, 40 grid points are used along both horizontal directions. In the vertical direction the domain must be high enough to include the entrainment processes at the top of the boundary layer as well as the propagation of gravity waves, which were stimulated by the convection. However, in the stably stratified region the grid resolution has not necessarily to be as high as within the boundary layer. This can be obtained by a vertical stretching of the grid starting from 1200 m via dz_stretch_level = 1200.0 m. This saves grid points and computing time. The upper boundary of the model is located at (see dz_stretch_factor) … m (computed by the model).
 
 Fast Fourier transformations are calculated using the Temperton-algorithm, which -on the IBM Regatta- is faster than the default system-specific algorithm (from IBM essl library).