About Soil various problem.

[FuYingWang]

Hello,
After I completed a simple simulation, I worked on improving the simulation accuracy. I improved the 'HGT_M' and 'LANDUSEF' variables in the 'geo_em.d01.nc' file using more accurate data. Next, I'm going to use more precise soil variables instead of the defaults, but I'm running into problems.

1. there seem to be only a few variables I've found regarding soil (e.g.CLAYFRAC, SOILCBOT,SOILCL1).
   I'd like to know if noahmp is using the 'CLAYFRAC' and 'SOILCBOT' variables??.
   
   2.I have more precise data for soil variables (such as sand%, clay%, soil classification), how should I use them to replace the default values.
   Thanks.

[tslin2]

There are options in the namelist.hrldas to change:
SOIL_DATA_OPTION = 1, ! Noah-MP options for defining soil properties [default = 1]
! geogrid must have been run with GEOGRID.TBL.ARW.noahmp, use with caution
! **1 -> use input dominant soil texture
! 2 -> use input soil texture that varies with depth
! 3 -> use soil composition (sand, clay, orgm) and pedotransfer functions (OPT_PEDO)
! 4 -> use input soil properties (BEXP_3D, SMCMAX_3D, etc.) (not valid in WRF)

Default option 1 use top layer soil texture for all soil layers in the Noah-MP
Option 2 use SOILC1-4, so soil texture change with soil layer.
Option 3 use clay and sand % from SOILCOMP, it has 8 values.
1-4 is sand % for each 4 soil layer, 5-8 is clay % for each 4 soil layer

You could use option 3 to change the SOILCOMP from default values to your soil data

[FuYingWang]

Thanks a lot.
I have add soil variable in geo_em.d01.nc about ' SOILCOMP'. And ncview can read it.

After creating GLDAS_forcing and change SOIL_DATA_OPTION = 3, I run the model.
But, When I run ./hrldas.exe namelist.hrldas, an error apper.

Then, I ncview HRLDAS_setup_2023010100_d1

It don't have SOILCOMP.
What details did I miss? And what should I do

[FuYingWang]

Thanks a lot.
I have add soil variable in geo_em.d01.nc about ' SOILCOMP'. And ncview can read it.
image
After creating GLDAS_forcing and change SOIL_DATA_OPTION = 3, I run the model.
But, When I run ./hrldas.exe namelist.hrldas, an error apper.

Then, I ncview HRLDAS_setup_2023010100_d1

It don't have SOILCOMP.

What details did I miss? And what should I do

[cenlinhe]

please directly add SOILCOMP to HRLDAS_setup file. Our current HRLDAS_forcing preprocessor does not have this SOILCOMP variable, which will not pass it from geo_em file to HRLDAS_setup file. We will add it to HRLDAS_forcing package later. Thanks!

[FuYingWang]

Sorry to bother you again, I followed your instructions and added 'SOILCOMP' to the 'HRLDA_SETUP' file.
SOILCOMP, it has 8 values. 1-4 is sand % for each 4 soil layer, 5-8 is clay % for each 4 soil layer.
But, the 'PedoTransferSR2006Mod.F90' seem also need 'Orgm', I only have sand and clay.
I can run ./rrldas.exe, but I don't seem to have 'Orgm' data. what would happen? What will be the value of Orgm?
Thanks.

[tslin2]

This assumes zero in the current model, it is not necessary to have it.

noahmp/drivers/hrldas/WaterVarInTransferMod.F90

Line 230 in fa2d209

SoilOrg = 0.0

if you like to replace it, you will revise this value. Or you can revise the code to read other values