Page tree

Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.
Comment: Minor updates
Excerpt
The NRSOILLAYEREADFLAG control flag controls the number of layers to read from MESH_parameters_CLASS.ini or MESH_parameters.txt when reading soil-related parameters and initial values.


Page contents:

Table of Contents

Available options

OptionDescriptionSupported since
0No effect, the flag is disabled. Three values are expected for all soil-related parameters and initial values when reading from MESH_parameters_CLASS.ini or MESH_parameters.txt.r747
1Read as many values as the total number of soil layers (NSL) defined in MESH_input_soil_levels.txt for all soil-related parameters and initial values from MESH_parameters_CLASS.ini or MESH_parameters.txt.r747
<N>

Read the specified number of <N> values for all soil-related parameters and initial values from MESH_parameters_CLASS.ini and MESH_parameters.txt. The value must be greater than 3 and will be automatically limited to the total number of soil layers (NSL) defined in MESH_input_soil_levels.txt:

Math
3 < N <= "NSL"
r1027

Examples

Behaviour without this option (NRSOILAYEREADFLAG 0)

Without activating NRSOILAYEREADFLAG, three values for each soil-related parameter and initial value field are required in MESH_parameters_CLASS.ini or MESH_parameters.txt, regardless of the number of soil layers defined in MESH_input_soil_levels.txt:

...

Info
The mapping of these fields from the lowest parameterized layer to lower layers happens automatically inside the MESH code. This is exactly the same process that occurred in MESH prior to the addition of the NRSOILAYEREADFLAG control flag.

NRSOILAYEREADFLAG 1

When NRSOILAYEREADFLAG is activated with option 1, no mapping is done and the same number of values must exist for the soil-related parameters and initial values in MESH_parameters_CLASS.ini as the total number of soil layers (NSL) defined in MESH_input_soil_levels.txt.

...

Example

NRSOILAYEREADFLAG is activated with option 1 in MESH_input_run_options.ini:

...

No Format
languagetext
titleMESH_parameters_CLASS.ini (excerpt)
      70.0      70.0      70.0      70.0  !> SAND(1:4)
      30.0      30.0      30.0      30.0  !> CLAY(1:4)
       0.0       0.0       0.0       0.0  !> ORGM(1:4)
       5.0       5.0       5.0       5.0       5.0       0.0       5.0  !> TBAR(1:4)/TCAN/TSNO/TPND
       0.2       0.2       0.2       0.2       0.0       0.0       0.0       0.0       0.0  !> THLQ(1:4)/THIC(1:4)/ZPND

NRSOILAYEREADFLAG <n>

When NRSOILAYEREADFLAG is activated with an option > 3 but less than the total number of soil layers (NSL), a similar mapping is done as when NRSOILAYEREADFLAG is disabled, but the lowest parameterized layer used to assign values to lower layers is the layer specified as <n>.

Info
If <n> is set equal to NSL, the the option is equivalent to "NRSOILAYEREADFLAG 1".

...

Example

NRSOILAYEREADFLAG is activated with option <n> in MESH_input_run_options.ini, where <n> is set to 4 to read values for 4 soil layers:

...

Code Block
languagefortran
do j = 5, NSL
  sand(:, j) = sand(:, 4)
  clay(:, j) = clay(:, 4)
  orgm(:, j) = orgm(:, 4)
  tbar(:, j) = tbar(:, 4)
  thlq(:, j) = thlq(:, 4)
  thic(:, j) = thic(:, 4)
end do

Verifying the state of NRSOILAYEREADFLAG

Info
The following only applies if CLASS is activated as the HLSS in the model. The following information is not printed for runs using SVS as the HLSS.

The state of NRSOILAYEREADFLAG can be inferred from the CLASS prognostic state summary printed to MESH_output_echo_print.txt at the end of a normal simulation.

If NRSOILAYEREADFLAG is disabled or if only 3 layers are defined in MESH_input_soil_layers.txt, only three values for the respective fields are printed to the file:

No Format
languagetext
titleMESH_output_echo_print.txt (excerpt)
    1.770     1.913    2.122     0.000   -0.011     0.723  !> TBAR(1:3)/TCAN/TSNO/TPND
    0.444     0.447    0.447     0.028    0.000     0.392     0.095  !>THLQ(1:3)/THIC(1:3)/ZPND
    0.000     0.000  151.256     0.840  893.673     0.000  !>RCAN/SNCAN/SNO/ALBS/RHOS/GRO

If NRSOILAYEREADFLAG is activated with option 1, and 4 layers are defined in MESH_input_soil_levels.txt, four values for the respective fields are printed to the file:

No Format
languagetext
titleMESH_output_echo_print.txt (excerpt)
    2.585     4.023    4.165     5.566    0.000     0.000     0.000  !>TBAR(1:4)/TCAN/TSNO/TPND
    0.117     0.139    0.042     0.342    0.000     0.000     0.000     0.000     0.000  !>THLQ(1:4)/THIC(1:4)/ZPND
    0.000     0.000    0.192     0.685  286.956     1.000  !> RCAN/SNCAN/SNO/ALBS/RHOS/GRO

If NRSOILAYEREADFLAG is activated with option <n> set to 4, and 6 layers are defined in MESH_input_soil_levels.txt, only four values for the respective fields are printed to the file:

No Format
languagetext
titleMESH_output_echo_print.txt (excerpt)
    2.585     4.023    4.165     5.566    0.000     0.000     0.000  !>TBAR(1:4)/TCAN/TSNO/TPND
    0.117     0.139    0.042     0.342    0.000     0.000     0.000     0.000     0.000  !>THLQ(1:4)/THIC(1:4)/ZPND
    0.000     0.000    0.192     0.685  286.956     1.000  !> RCAN/SNCAN/SNO/ALBS/RHOS/GRO

Considerations

Decreasing the number of soil layers in an existing setup with "NRSOILAYEREADFLAG 1"

...

This cascade would continue until no more numbers could be read from the file, when the model would again stop with a 'read'-related error.

Verifying the state of NRSOILAYEREADFLAG

Info
The following only applies if CLASS is activated as the HLSS in the model. The following information is not printed for runs using SVS as the HLSS.

The state of NRSOILAYEREADFLAG can be inferred from the CLASS prognostic state summary printed to MESH_output_echo_print.txt at the end of a normal simulation.

If NRSOILAYEREADFLAG is disabled or if only 3 layers are defined in MESH_input_soil_layers.txt, only three values for the respective fields are printed to the file:

No Format
languagetext
titleMESH_output_echo_print.txt (excerpt)
    1.770     1.913    2.122     0.000   -0.011     0.723  !> TBAR(1:3)/TCAN/TSNO/TPND
    0.444     0.447    0.447     0.028    0.000     0.392     0.095  !>THLQ(1:3)/THIC(1:3)/ZPND
    0.000     0.000  151.256     0.840  893.673     0.000  !>RCAN/SNCAN/SNO/ALBS/RHOS/GRO

If NRSOILAYEREADFLAG is activated with option 1, and 4 layers are defined in MESH_input_soil_levels.txt, four values for the respective fields are printed to the file:

No Format
languagetext
titleMESH_output_echo_print.txt (excerpt)
    2.585     4.023    4.165     5.566    0.000     0.000     0.000  !>TBAR(1:4)/TCAN/TSNO/TPND
    0.117     0.139    0.042     0.342    0.000     0.000     0.000     0.000     0.000  !>THLQ(1:4)/THIC(1:4)/ZPND
    0.000     0.000    0.192     0.685  286.956     1.000  !> RCAN/SNCAN/SNO/ALBS/RHOS/GRO

If NRSOILAYEREADFLAG is activated with option <n> set to 4, and 6 layers are defined in MESH_input_soil_levels.txt, only four values for the respective fields are printed to the file:

No Format
languagetext
titleMESH_output_echo_print.txt (excerpt)
    2.585     4.023    4.165     5.566    0.000     0.000     0.000  !>TBAR(1:4)/TCAN/TSNO/TPND
    0.117     0.139    0.042     0.342    0.000     0.000     0.000     0.000     0.000  !>THLQ(1:4)/THIC(1:4)/ZPND
    0.000     0.000    0.192     0.685  286.956     1.000  !> RCAN/SNCAN/SNO/ALBS/RHOS/GRO

Related pages

Content by Label
showLabelsfalse
showSpacefalse
excerptTypesimple
cqllabel in ("control-flag","mesh-input","= "soil-levels" ) and space = "MESH"