About SoilGrid routines

[fmanquehual]

Hello, I hope you are well.

I am using the routines “_main_prepare_lut.f90” and “main_group_soil.f90” located in the folder “mhm-v5.11.1/pre-proc/process_soilgrid” to process the ‘SoilGrid’ data. However, when I try to use it, I have a doubt about the LUT table generated, as it does not provide the ‘HORIZON’, 'UD' or 'LD' columns, and I think this is why I get an error when using the outputs in mHM. If you could clarify whether the outputs I have are correct or incorrect, either when generating it or how to represent it properly in ‘mhm.nml’, I would appreciate it. Below are some outputs that may clarify the situation:

After compiling and preparing the data, it gives me the following outputs:
#===================================================#
... jun 16 23:18 soil_lut.txt
... jun 16 23:18 soil_id.asc
... jun 16 23:20 soil_class.asc (<-- this is a copy from soil_id.asc)
... jun 16 23:20 soil_classdefinition.txt (<-- this is a copy from soil_lut.txt)
... jun 16 23:20 soil_class_horizon_depth_00600-01000mm.asc
... jun 16 23:20 soil_class_horizon_depth_00300-00600mm.asc
... jun 16 23:20 soil_class_horizon_depth_00150-00300mm.asc
... jun 16 23:20 soil_class_horizon_depth_00050-00150mm.asc
... jun 16 23:20 soil_class_horizon_depth_00000-00050mm.asc
... jun 16 23:20 soil_id_reduced_grid.asc
... jun 16 23:20 soil_frequency_reduced_grid.asc
... jun 16 23:20 bd_stdev_reduced_grid.asc
... jun 16 23:20 bd_mean_reduced_grid.asc
... jun 16 23:20 soil_classdefinition_reduced_grid.txt
#===================================================#

The content of ‘soil_classdefinition_reduced_grid.txt’ (which I use as “soil_classdefinition.txt” in “mhm.nml”) is:
#===================================================#
nSoil_Types 885
ID CLAY[%] SAND[%] Bd_mu[gcm-3] Bd_sd[gcm-3] no.Sample
1 5 7 0.098 0.001 1
2 6 1 0.080 0.001 3
3 6 2 0.084 0.004 3
4 6 6 0.083 0.001 1
5 6 7 0.103 0.001 1
6 6 8 0.091 0.011 6
7 6 9 0.086 0.008 7
8 6 10 0.087 0.009 5
9 6 11 0.091 0.014 4
10 6 12 0.080 0.001 1
#===================================================#

Error when I ran mHM:
#===================================================#
Given data directories:
#--------------
DOMAIN 1
#--------------
Morphological directory: nhuble_alto/input/morph (found)
Land cover directory: nhuble_alto/input/luse (found)
Level-2 header directory: nhuble_alto/input/meteo/pre (found)
Precipitation directory: nhuble_alto/input/meteo/pre (found)
Temperature directory: nhuble_alto/input/meteo/tavg (found)
PET directory: nhuble_alto/input/meteo/pet (found)
Output directory: nhuble_alto/output_b1 (found)

Read data ...
Reading data ...
Reading dem and lcover ...
Reading dem for domain: 1 ...
Reading lcover for domain: 1 ...
read_soil_LUT: ERROR occurred: Mixed horizon depths in soil type 1 and horizon no. 5
STOP 1
#===================================================#

This is a segment of 'mhm.nml' about soil depth:
#===================================================#
iFlag_soilDB = 0
!
!> [mm] soil depth down to which organic matter is possible
tillageDepth = 100 ! 200
!
!> No. of soil horizons to be modeled
nSoilHorizons_mHM = 5
!
! soil depth information
!> IF (iFlag_soilDB = 0)
!> Provide below the soil_Depth() for 1,2,..,n-1 soil horizons. Depth of the last layer(n) is taken from the soil LUT
!> IF (iFlag_soilDB = 1)
!> Provide below the soil_Depth() for every 1,2..n-1,n soil horizons. You must have soil_class-id gridded file for each layer
!> Also check your tillage layer depth. It should conform with one of the below specified soil_Depth.
!
!> Soil_Horizon Depth[mm] ! bottom depth of soil horizons w.r.t. ground surface (positive downwards)
soil_Depth(1) = 50 ! 200
soil_Depth(2) = 150
soil_Depth(3) = 300
soil_Depth(4) = 600
soil_Depth(5) = 1000
#===================================================#

[fmanquehual]

The answer is simple, just I need change "iFlag_soilDB = 0" by "iFlag_soilDB = 1". Thank you Sr. Rohini.

By the way, "tillageDepth" need the same depth of some soil_Depth(:).