Questions about Soil Layers and Interpolation in CLM Model Output

[liu845968376]

Hello, I am a new user of the CLM model. The following questions might seem a bit naive, but I hope you can help me with them:

In my CLM model output files, the H2OSOI (soil moisture) variable has 15 layers, but I checked the user manual and found that it should have 20 layers. Also, the values for the soil layers seem strange; they are irrational numbers, including:
0.00710063
0.02792499959
0.06225857511
0.1188650653
0.2121933997
0.3660658002
0.6197584867
1.038027048
1.727635264
2.864607096
4.739156723
7.829766273
12.92532063
21.32646942
35.17761993

My questions are:

• Is this output correct?
• Are the units of the soil layers in the output in real physical layers, and are the units in meters?
• If I need to interpolate the data to other layers, do I simply need to use the command cdo intlevel,0.05,0.25,0.7,1.5 X.nc X1.nc?

Finally, why is a non-uniform grid like this used for the soil layers?

 

[slevis]

I don't, yet, have answers to your questions, but I looked at one of my simulations and found the following:

1) In the lnd_in I see
soil_layerstruct_predefined = '20SL_8.5m'
Your simulation should say the same. If not, then you should add it to your case's user_nl_clm.

2) In my simulation's lnd.log I see
dzsoi: 2.000000000000000E-002 4.000000000000000E-002 6.000000000000000E-002
8.000000000000000E-002 0.120000000000000 0.160000000000000
0.200000000000000 0.240000000000000 0.280000000000000
0.320000000000000 0.360000000000000 0.400000000000000
0.440000000000000 0.540000000000000 0.640000000000000
0.740000000000000 0.840000000000000 0.940000000000000
1.04000000000000 1.14000000000000 2.39000000000000
4.67553390593274 7.63519052838329 11.1400000000000
15.1154248593737
and this looks different from yours.

I don't know where your simulation got 15 layers from, but you are able to set your own soil layers through the namelist variables
soil_layerstruct_userdefined
soil_layerstruct_userdefined_nlevsoi