294 | | [=#<insert_parameter_name> '''<insert_parameter_name>'''] |
295 | | }}} |
296 | | {{{#!td style="vertical-align:top" |
297 | | <insert type> |
298 | | }}} |
299 | | {{{#!td style="vertical-align:top" |
300 | | <insert value> |
301 | | }}} |
302 | | {{{#!td |
303 | | <insert explanation> |
304 | | }}} |
305 | | |---------------- |
306 | | {{{#!td style="vertical-align:top" |
307 | | [=#<insert_parameter_name> '''<insert_parameter_name>'''] |
308 | | }}} |
309 | | {{{#!td style="vertical-align:top" |
310 | | <insert type> |
311 | | }}} |
312 | | {{{#!td style="vertical-align:top" |
313 | | <insert value> |
314 | | }}} |
315 | | {{{#!td |
316 | | <insert explanation> |
| 294 | [=#bc_s_b '''bc_s_b'''] |
| 295 | }}} |
| 296 | {{{#!td style="vertical-align:top" |
| 297 | C*20 |
| 298 | }}} |
| 299 | {{{#!td style="vertical-align:top" |
| 300 | 'dirichlet' |
| 301 | }}} |
| 302 | {{{#!td |
| 303 | Bottom boundary condition of the scalar concentration.\\\\ |
| 304 | Allowed values are '' 'dirichlet' '' (s(k=0) = const. = [#s_surface s_surface] + [#s_surface_initial_change s_surface_initial_change]; the user may change this value during the run using [#user-defined] code) and '' 'neumann' '' (s(k=0) = s(k=1)). |
| 305 | When a constant surface concentration flux is used ([#surface_scalarflux surface_scalarflux]), '''bc_s_b''' = '' 'neumann' '' must be used, because otherwise the resolved scale may contribute to the surface flux so that a constant value cannot be guaranteed. |
| 306 | }}} |
| 307 | |---------------- |
| 308 | {{{#!td style="vertical-align:top" |
| 309 | [=#bc_s_t '''bc_s_t'''] |
| 310 | }}} |
| 311 | {{{#!td style="vertical-align:top" |
| 312 | C*20 |
| 313 | }}} |
| 314 | {{{#!td style="vertical-align:top" |
| 315 | 'neumann' |
| 316 | }}} |
| 317 | {{{#!td |
| 318 | Top boundary condition of the scalar concentration.\\\\ |
| 319 | Allowed are the values '' 'dirichlet' '' (s(k=nz) and s(k=nz+1) do not change during the run) and '' 'neumann' ''. With the Neumann boundary condition the value of the scalar concentration gradient at the top is calculated from the initial scalar concentration profile (see [#s_surface s_surface], [#s_vertical_gradient s_vertical_gradient]) by: bc_s_t_val = (s_init(k=nz) - s_init(k=nz-1)) / dzu(nz). |
| 320 | Using this value (assumed constant during the run) the concentration boundary values are calculated as |
| 321 | |
| 322 | s(k=nz+1) = s(k=nz) + bc_s_t_val * dzu(nz+1) |
| 323 | |
| 324 | (up to k=nz the prognostic equation for the scalar concentration is solved). |