315 | | The main principle are: |
316 | | -- The T scheme is always as depth or deeper than the W scheme |
317 | | -- W and T share the same common discretization for the common part: the layers are the same; For the nodes: W has the firt node at the surface while T has it at the middle of the first layer; Else they have the same node position |
318 | | -- The upper part of the scheme follow a geometric discretization of reason 2 (like the current W scheme) |
319 | | -- Below a certain depth, depth_lin, we have for the W and T a linear discretization, with the depth of each layers equal to that of the lowest layer of the "geometrical part" of the scheme. Note that depth_lin may be equal to depth max with thus no linear discretization section. |
320 | | -- Below another depth, depth_geom, we have again the geometrical discretization with a reason xx to be defined. This allows to model a deep soil for permafrost modelling. Note that Depth_geom may also be equal to depth_max with thus no bottom geometrical discretization. |
| 315 | The main principles are: |
| 316 | - The T scheme is always as depth or deeper than the W scheme |
| 317 | - W and T share the same common discretization for the common part: the layers are the same; For the nodes: W has the firt node at the surface while T has it at the middle of the first layer; Else they have the same node position |
| 318 | - The upper part of the scheme follow a geometric discretization of reason 2 (like the current W scheme) |
| 319 | - Below a certain depth, depth_lin, we have for the W and T a linear discretization, with the depth of each layers equal to that of the lowest layer of the "geometrical part" of the scheme. Note that depth_lin may be equal to depth max with thus no linear discretization section. |
| 320 | - Below another depth, depth_geom, we have again the geometrical discretization with a reason xx to be defined. This allows to model a deep soil for permafrost modelling. Note that Depth_geom may also be equal to depth_max with thus no bottom geometrical discretization. |
337 | | -- Heat capacity of each layer (C) should be a function of the layer water content: thus a function of the layer total water content ("W" variable of CWRR divided by the layer depth) |
338 | | -- Conductivity between layers (Lambda) should be a function of the "Theta" of layer i and i+1 (weighted arythmetic mean using the depth of layers) |
339 | | -- For the part of the thermal scheme that is below the hydrology scheme: we should take a constant in time value of "Theta" and "water content", to avoid changing heat capacity over time and thus not conserving the energy (without any corresponding imput/output of energy by water flow); the values to choose is open for discussion/tests ? |
| 337 | - Heat capacity of each layer (C) should be a function of the layer water content: thus a function of the layer total water content ("W" variable of CWRR divided by the layer depth) |
| 338 | - Conductivity between layers (Lambda) should be a function of the "Theta" of layer i and i+1 (weighted arythmetic mean using the depth of layers) |
| 339 | - For the part of the thermal scheme that is below the hydrology scheme: we should take a constant in time value of "Theta" and "water content", to avoid changing heat capacity over time and thus not conserving the energy (without any corresponding imput/output of energy by water flow); the values to choose is open for discussion/tests ? |