Preparation soil map data using SoilGrid datasets

[Mohammedhhu]

Dear Developers,

I am working with the SoilGrid datasets for soil mapping and have downloaded six layers of bulk density, sand, and clay content using Google Earth Engine. I processed these datasets in ArcMap to fill missing pixels, clip, resample, and mask them with the river basin boundary, then converted them to text files.

I updated the file paths and names in _main_prepare_lut.f90 and main_group_soil.f90 and attempted to compile the code using the following steps:
Compilation Steps

gfortran -c mo_kind.f90
gfortran _main_prepare_lut.f90 -o prepare_lut.exe
prepare_lut.exe

These steps ran successfully, generating soil_id.asc and soil_lut.txt. However, I encountered errors in the next steps:

gfortran -c mo_moment.f90
gfortran main_group_soil.f90 -o main_group_soil.exe
main_group_soil.exe

The compilation of main_group_soil.f90 fails with an undefined reference error related to functions from mo_moment.f90. I would like to confirm if I am following the correct procedure for preparing inputs and compiling the code.

Any guidance on resolving this issue would be greatly appreciated.

Thank you.

Best regards,
Mohammed Abdallah

[rohinikumar82]

Possibly, you have a soil type in your dataset that occurs only once, and therefore, the standard deviation routine can not be rightly executed. One way would be to comment out - lines in the code for the SD estimation.

If you answered, do not forget to close the question. Thanks.

Best,
Rohini