| 51 | {{{#!td style="vertical-align:top; width: 150px" |
| 52 | [=#cs_surface '''bc_cs_b'''] |
| 53 | }}} |
| 54 | {{{#!td style="vertical-align:top; width: 50px" |
| 55 | R |
| 56 | }}} |
| 57 | {{{#!td style="vertical-align:top; width: 75px" |
| 58 | 0.0 |
| 59 | }}} |
| 60 | {{{#!td |
| 61 | Bottom boundary condition of the chemical species ('cs')concentration. |
| 62 | |
| 63 | Allowed values are 'dirichlet' (cs(k=0) = const. = cs_surface + cs_surface_initial_change; When a constant surface concentration flux is used (surface_cs_flux) or emissions are applied ('do_emis' = .T.), bc_cs_b = 'neumann' must be used. }}} |
| 64 | |
| 65 | {{{#!td style="vertical-align:top; width: 150px" |
| 66 | [=#cs_surface '''bc_cs_t'''] |
| 67 | }}} |
| 68 | {{{#!td style="vertical-align:top; width: 50px" |
| 69 | R |
| 70 | }}} |
| 71 | {{{#!td style="vertical-align:top; width: 75px" |
| 72 | 0.0 |
| 73 | }}} |
| 74 | {{{#!td |
| 75 | Top boundary condition of the scalar concentration. |
| 76 | |
| 77 | Allowed are the values 'dirichlet' (cs(k=nz+1) does not change during the run), 'neumann' (cs(k=nz+1) = cs(k=nz)), |
| 78 | |
| 79 | ### RFo No idea whther this also works für cs |
| 80 | and 'initial_gradient' . |
| 81 | With the 'initial_gradient' boundary condition the value of the scalar concentration gradient at the top is calculated from the initial scalar concentration profile (see cs_surface, cs_vertical_gradient) by: bc_cs_t_val = (cs_init(k=nz) - ss_init(k=nz-1)) / dzu(nz). Using this value (assumed constant during the run) the concentration boundary values are calculated as |
| 82 | |
| 83 | |
| 84 | cs(k=nz+1) = cs(k=nz) + bc_cs_t_val * dzu(nz+1) |
| 85 | |
| 86 | (up to k=nz the prognostic equation for the chemical species concentration is solved). |
| 87 | ### |
| 88 | |
| 89 | When a constant cs flux is used at the top boundary (top_csflux), bc_cs_t = 'neumann' must be used, because otherwise the resolved scale may contribute to the top flux so that a constant value cannot be guaranteed. |
| 90 | }}} |
| 91 | |