Changes between Version 94 and Version 95 of doc/app/chemistry_parameters
- Timestamp:
- Jan 8, 2019 7:31:59 AM (5 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
doc/app/chemistry_parameters
v94 v95 78 78 Bottom boundary condition of the chemical species (`cs`) concentration. 79 79 80 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.80 Allowed values are 'dirichlet' (cs(k=0) = const. = [#cs_surface cs_surface] + cs_surface_initial_change; When a constant surface concentration flux is used ([#surface_csflux surface_csflux]) or emissions are applied ([#do_emis do_emis] = .T.), '''bc_cs_b''' = 'neumann' must be used. 81 81 }}} 82 82 |---------------- … … 95 95 Allowed are the values 'dirichlet' (cs(k=nz+1) does not change during the run), 'neumann' (cs(k=nz+1) = cs(k=nz)), 96 96 and 'initial_gradient' . 97 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 as97 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_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 98 98 cs(k=nz+1) = cs(k=nz) + bc_cs_t_val * dzu(nz+1) 99 99 (up to k=nz the prognostic equation for the chemical species concentration is solved).