Changes between Version 27 and Version 28 of doc/app/initialization_parameters

Timestamp:

Sep 13, 2010 2:55:54 PM (14 years ago)

Author:

kanani

Comment:

--

doc/app/initialization_parameters

@@ -258 +258 @@
 |----------------
 {{{#!td style="vertical-align:top"
-[=#<insert_parameter_name> '''<insert_parameter_name>''']
-}}}
-{{{#!td style="vertical-align:top"
-<insert type>
-}}}
-{{{#!td style="vertical-align:top"
-<insert value>
-}}}
-{{{#!td
-<insert explanation>
-}}}
-|----------------
-{{{#!td style="vertical-align:top"
-[=#<insert_parameter_name> '''<insert_parameter_name>''']
-}}}
-{{{#!td style="vertical-align:top"
-<insert type>
-}}}
-{{{#!td style="vertical-align:top"
-<insert value>
-}}}
-{{{#!td
-<insert explanation>
+[=#bc_q_b '''bc_q_b''']
+}}}
+{{{#!td style="vertical-align:top"
+C*20
+}}}
+{{{#!td style="vertical-align:top"
+'dirichlet'
+}}}
+{{{#!td
+Bottom boundary condition of the specific humidity / total water content.\\\\
+Allowed values are '' 'dirichlet' '' (q(k=0) = const. = [#q_surface q_surface] + [#q_surface_initial_change q_surface_initial_change]; the user may change this value during the run using [#user-defined] code) and '' 'neumann' '' (q(k=0)=q(k=1)). 
+When a constant surface latent heat flux is used ([#surface_waterflux surface_waterflux]), '''bc_q_b''' = '' 'neumann' '' must be used, because otherwise the resolved scale may contribute to the surface flux so that a constant value cannot be guaranteed.
+}}}
+|----------------
+{{{#!td style="vertical-align:top"
+[=#bc_q_t '''bc_q_t''']
+}}}
+{{{#!td style="vertical-align:top"
+C*20
+}}}
+{{{#!td style="vertical-align:top"
+'neumann'
+}}}
+{{{#!td
+Top boundary condition of the specific humidity / total water content.\\\\
+Allowed are the values '' 'dirichlet' '' (q(k=nz) and q(k=nz+1) do not change during the run) and '' 'neumann' ''. With the Neumann boundary condition the value of the humidity gradient at the top is calculated from the initial humidity profile (see [#q_surface q_surface], [#q_vertical_gradient q_vertical_gradient]) by: bc_q_t_val = ( q_init(k=nz) - q_init(k=nz-1)) / dzu(nz).
+Using this value (assumed constant during the run) the humidity boundary values are calculated as 
+
+      q(k=nz+1) =q(k=nz) + bc_q_t_val * dzu(nz+1)
+
+(up tp k=nz the prognostic equation for q is solved). 
 }}}
 |----------------