| 1 | == Canopy Parameters == |
| 2 | [[TracNav(doc/app/partoc|nocollapse)]] ----noch ändern? |
| 3 | \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\ |
| 4 | '''NAMELIST group name: canopy_par''' \\ |
| 5 | |
| 6 | ||='''Parameter Name''' =||='''[../fortrantypes FORTRAN Type]''' =||='''Default Value''' =||='''Explanation''' =|| |
| 7 | |---------------- |
| 8 | {{{#!td style="vertical-align:top;width: 150px" |
| 9 | [=#canopy_mode '''canopy_mode'''] |
| 10 | }}} |
| 11 | {{{#!td style="vertical-align:top;width: 50px" |
| 12 | C*20 |
| 13 | }}} |
| 14 | {{{#!td style="vertical-align:top;width: 75px" |
| 15 | 'block' |
| 16 | }}} |
| 17 | {{{#!td |
| 18 | Canopy mode.\\\\ |
| 19 | Besides using the default value, that will create a horizontally homogeneous plant canopy that extends over the total horizontal extension of the model domain, the user may add code to the user interface (see [#3.5.1 3.5.1]) subroutine {{{user_init_plant_canopy}}} to allow further canopy modes.\\\\ |
| 20 | The setting of '''canopy_mode''' becomes only active, if [#plant_canopy plant_canopy] has been set ''.T.'' and a non-zero [#drag_coefficient drag_coefficient] has been defined. |
| 21 | }}} |
| 22 | |---------------- |
| 23 | {{{#!td style="vertical-align:top" |
| 24 | [=#cthf '''cthf'''] |
| 25 | }}} |
| 26 | {{{#!td style="vertical-align:top" |
| 27 | R |
| 28 | }}} |
| 29 | {{{#!td style="vertical-align:top" |
| 30 | 0.0 |
| 31 | }}} |
| 32 | {{{#!td |
| 33 | Average heat flux that is prescribed at the top of the plant canopy.\\\\ |
| 34 | If [#plant_canopy plant_canopy] is set ''.T.'', the user can prescribe a heat flux at the top of the plant canopy. |
| 35 | It is assumed that solar radiation penetrates the canopy and warms the foliage which, in turn, warms the air in contact with it.\\\\ |
| 36 | '''Note:'''\\ |
| 37 | Instead of using the value prescribed by [#surface_heatflux surface_heatflux], the near surface heat flux is determined from an exponential function that is dependent on the cumulative leaf_area_index (Shaw and Schumann (1992, Boundary Layer Meteorol., 61, 47-64)). |
| 38 | }}} |
| 39 | |---------------- |
| 40 | {{{#!td style="vertical-align:top" |
| 41 | [=#drag_coefficient '''drag_coefficient'''] |
| 42 | }}} |
| 43 | {{{#!td style="vertical-align:top" |
| 44 | R |
| 45 | }}} |
| 46 | {{{#!td style="vertical-align:top" |
| 47 | 0.0 |
| 48 | }}} |
| 49 | {{{#!td |
| 50 | Drag coefficient used in the {{{plant_canopy_model}}}.\\\\ |
| 51 | This parameter has to be non-zero, if the parameter [#plant_canopy plant_canopy] is set ''.T.''. |
| 52 | }}} |
| 53 | |---------------- |
| 54 | {{{#!td style="vertical-align:top" |
| 55 | [=#lad_surface '''lad_surface'''] |
| 56 | }}} |
| 57 | {{{#!td style="vertical-align:top" |
| 58 | R |
| 59 | }}} |
| 60 | {{{#!td style="vertical-align:top" |
| 61 | 0.0 |
| 62 | }}} |
| 63 | {{{#!td |
| 64 | Surface value of the leaf area density (in m^2^/m^3^).\\\\ |
| 65 | This parameter assigns the value of the leaf area density '''lad''' at the surface (k=0). Starting from this value, the leaf area density profile is constructed with [#lad_vertical_gradient lad_vertical_gradient] and [#lad_vertical_gradient_level lad_vertical_gradient_level]. |
| 66 | }}} |
| 67 | |---------------- |
| 68 | {{{#!td style="vertical-align:top" |
| 69 | [=#lad_vertical_gradient '''lad_vertical_gradient'''] |
| 70 | }}} |
| 71 | {{{#!td style="vertical-align:top" |
| 72 | R(10) |
| 73 | }}} |
| 74 | {{{#!td style="vertical-align:top" |
| 75 | 10 * 0.0 |
| 76 | }}} |
| 77 | {{{#!td |
| 78 | Gradient(s) of the leaf area density (in m^2^/m^4^).\\\\ |
| 79 | This leaf area density gradient holds starting from the height level defined by [#lad_vertical_gradient_level lad_vertical_gradient_level] (precisely: for all uv levels k where zu(k) > lad_vertical_gradient_level, lad(k) is set: lad(k) = lad(k-1) + dzu(k) * '''lad_vertical_gradient''') up to the level defined by [#pch_index pch_index]. Above that level lad(k) will automatically be set to 0.0. A total of 10 different gradients for 11 height intervals (10 intervals if lad_vertical_gradient_level(1) = 0.0) can be assigned. The leaf area density at the surface is assigned via [#lad_surface lad_surface]. |
| 80 | }}} |
| 81 | |---------------- |
| 82 | {{{#!td style="vertical-align:top" |
| 83 | [=#lad_vertical_gradient_level '''lad_vertical_gradient_level'''] |
| 84 | }}} |
| 85 | {{{#!td style="vertical-align:top" |
| 86 | R(10) |
| 87 | }}} |
| 88 | {{{#!td style="vertical-align:top" |
| 89 | 10 * 0.0 |
| 90 | }}} |
| 91 | {{{#!td |
| 92 | Height level from which on the gradient of the leaf area density defined by [#lad_vertical_gradient lad_vertical_gradient] is effective (in m).\\\\ |
| 93 | The height levels have to be assigned in ascending order. The default values result in a leaf area density that is constant with height up to the top of the plant canopy layer defined by [#pch_index pch_index]. For the piecewise construction of temperature profiles see [#lad_vertical_gradient lad_vertical_gradient]. |
| 94 | }}} |
| 95 | |---------------- |
| 96 | {{{#!td style="vertical-align:top" |
| 97 | [=#leaf_surface_concentration '''leaf_surface_concentration'''] |
| 98 | }}} |
| 99 | {{{#!td style="vertical-align:top" |
| 100 | R |
| 101 | }}} |
| 102 | {{{#!td style="vertical-align:top" |
| 103 | 0.0 |
| 104 | }}} |
| 105 | {{{#!td |
| 106 | Concentration of a passive scalar at the surface of a leaf (in K m/s).\\\\ |
| 107 | This parameter is only of importance in cases in that both, [#plant_canopy plant_canopy] and [#passive_scalar passive_scalar], are set ''.T..'' The value of the concentration of a passive scalar at the surface of a leaf is required for the parametrisation of the sources and sinks of scalar concentration due to the canopy. |
| 108 | }}} |
| 109 | |---------------- |
| 110 | {{{#!td style="vertical-align:top" |
| 111 | [=#pch_index '''pch_index'''] |
| 112 | }}} |
| 113 | {{{#!td style="vertical-align:top" |
| 114 | I |
| 115 | }}} |
| 116 | {{{#!td style="vertical-align:top" |
| 117 | 0 |
| 118 | }}} |
| 119 | {{{#!td |
| 120 | Grid point index (scalar) of the upper boundary of the plant canopy layer.\\\\ |
| 121 | Above '''pch_index''' the arrays of leaf area density and [#drag_coeffient drag_coeffient] are automatically set to zero in case of [#plant_canopy plant_canopy] = ''.T.''. Up to '''pch_index''' a leaf area density profile can be prescribed by using the parameters [#lad_surface lad_surface], [#lad_vertical_gradient lad_vertical_gradient] and [#lad_vertical_gradient_level lad_vertical_gradient_level]. |
| 122 | }}} |
| 123 | |---------------- |
| 124 | {{{#!td style="vertical-align:top" |
| 125 | [=#plant_canopy '''plant_canopy'''] |
| 126 | }}} |
| 127 | {{{#!td style="vertical-align:top" |
| 128 | L |
| 129 | }}} |
| 130 | {{{#!td style="vertical-align:top" |
| 131 | .F. |
| 132 | }}} |
| 133 | {{{#!td |
| 134 | Switch for the plant canopy model.\\\\ |
| 135 | If '''plant_canopy''' is set ''.T.'', the plant canopy model of Watanabe (2004, BLM 112, 307-341) is used.\\ |
| 136 | The impact of a plant canopy on a turbulent flow is considered by an additional drag term in the momentum equations and an additional sink term in the prognostic equation for the subgrid-scale TKE. These additional terms depend on the leaf drag coefficient (see [#drag_coefficient drag_coefficient]), and the leaf area density (see [#lad_surface lad_surface], [#lad_vertical_gradient lad_vertical_gradient], [#lad_vertical_gradient_level lad_vertical_gradient_level]). The top boundary of the plant canopy is determined by the parameter [#pch_index pch_index]. For all heights equal or larger than zw(k=pch_index), the leaf area density is 0 (i.e. there is no canopy at these heights!).\\ |
| 137 | By default, a horizontally homogeneous plant canopy is prescribed, if '''plant_canopy''' is set ''.T.''. However, the user can define other types of plant canopies (see [#canopy_mode canopy_mode]).\\\\ |
| 138 | If '''plant_canopy''' and [#passive_scalar passive_scalar] are set ''.T.'', the canopy acts as an additional scalar source or sink, respectively. The source/sink strength depends on the scalar concentration at the leaf surface, which generally does not vary with time in PALM, and which can be specified with parameter [#leaf_surface_concentration leaf_surface_concentration].\\\\ |
| 139 | Additional heating by the plant canopy is taken into account, if the default value of parameter [#cthf cthf] is altered in the parameter file. In that case, the value of [#surface_heatflux surface_heatflux] specified in the parameter file is not used in the model. Instead, the near-surface heat flux is derived from an exponential function that depends on the cumulative leaf area index.\\\\ |
| 140 | '''plant_canopy''' = ''.T.'' is only allowed with a non-zero [#drag_coefficient drag_coefficient]. |
| 141 | }}} |
| 142 | |---------------- |
| 143 | {{{#!td style="vertical-align:top" |
| 144 | [=#scalar_exchange_coefficient '''scalar_exchange_coefficient'''] |
| 145 | }}} |
| 146 | {{{#!td style="vertical-align:top" |
| 147 | R |
| 148 | }}} |
| 149 | {{{#!td style="vertical-align:top" |
| 150 | 0.0 |
| 151 | }}} |
| 152 | {{{#!td |
| 153 | Scalar exchange coefficient for a leaf (dimensionless).\\\\ |
| 154 | This parameter is only of importance in cases in that both, [#plant_canopy plant_canopy] and [#passive_scalar passive_scalar], are set ''.T.''. The value of the scalar exchange coefficient is required for the parametrisation of the sources and sinks of scalar concentration due to the canopy. |
| 155 | }}} |
| 156 | |---------------- |
| 157 | |
| 158 | |