question: SOILC_vr

[tenghao]

I used clm 5.0 to create the following case：

./create_newcase --case test02_I2000Clm50BgcCrop --compset 2000_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_SGLC_SWAV --res CLM_USRDAT --machine myintel --compiler intel --run-unsupported



Defined in user_nl_clm:

fsurdat = '/**/**/****.nc'

flanduse_timeseries = '/**/**/****.nc'

spinup_state = 2



hist_empty_htapes=.true.

hist_fincl1 = 'TOTSOMC', 'SOILC_vr'

hist_nhtfrq = -24

hist_mfilt = 10



Prompt error:

htapes_fieldlist ERROR: SOILC_vr in fincl( 2 ) for history tape

1 not found

ENDRUN:

ERROR: ERROR in histFileMod.F90 at line 681



Can't this compsets generate the variable SOILC_vr? Which compsets can generate SOILC_vr? Does SOILC_vr refer to organic carbon in soil?

Additionally, I am unable to access the CLM namelist definitions website. Is the website being maintained?

Thank you for your response.

Best wishes.

 

[tenghao]

Q1: When I add the output variables SOIL1C, SOIL2C and SOIL3C to the user_nl_clm, the model can run normally and output the results, but SOIL1C, SOIL2C and SOIL3C are two-dimensional variables, which are not in line with my expectations. What should I do if I still want to output the variable SOILC_vr? Is this caused by a version problem? Can I output SOILC_vr variable using clm5.0.35?


Q2: How to customize the soil thickness in clm? It seems that it needs to be defined in user_nl_clm by soil_layerstruct_userdefined and soil_layerstruct_userdefined_nlevsoi, but I don't know how to set these two variables to meet my requirements (for example, I want to define the soil layer size as 20 layers, with soil layer thicknesses of 0-0.2m, 0.2-0.4m, 0.4-0.6m, 0.6-0.8m, 0.8-1.0m, 1-1.5m, 1.5-2m, 2-2.5m, 2.5-3m, 3-3.5m, 3.5-4m, 4-4.5m, 4.5-5m, 5-5.5m, 5.5-6m, 6-6.5m, 6.5-7m, 7-7.5m, 7.5-8m, 8-8.6m) Do I need to change the F90 file related to soil layers or only set it in user_nl_clm? What should I do if I need to make changes to the F90 file?

 

[oleson]

From your first post. SOILC_vr should be available in release-clm5.0.35 with that compset:

call hist_addfld2d (fname='SOILC_vr', units='gC/m^3', type2d='levsoi', &
avgflag='A', long_name='SOIL C (vertically resolved)', &
ptr_col=this%decomp_soilc_vr_col)

It looks like it is conditional on this in the lnd_in so you might check to see if you have this:

use_vertsoilc = .true.

It should actually be output by default, you shouldn't have to specify it.
I've tried it here and it outputs by default fine. Maybe you have a typo in your user_nl_clm? Try not specifying any changes to the history tapes in your user_nl_clm.

The web site for the CESM2.1.5/CLM5 namelist definitions seems to be up:

CESM2.1.5 CLM Namelist Definitions

docs.cesm.ucar.edu

This post may be the most useful for instructions on how to customize the soil thicknesses:

Lower boundary conditions (differences between CLM4.5 and CLM5)

Hello, Pardon my insufficient knowledge of the Fortran coding but I have to ask the following questions here as I only have this forum to check my inquiries with: I have found the "Modified CLM4.5 files to change the bottom heat flux and add 24 bedrock layers of thickness 12.5 m at the bottom...

bb.cgd.ucar.edu

The user there has successfully customize the soil layer structure in CLM5 so you may want to ping them directly.