The output problem of spin up

[Lian Gong]

Hi,

I spin up a single-point in clm5.0, using

./create_newcase --case eight.fiftten --res CLM_USRDAT --compset 2000_DATM%1PT_CLM50%BGC_SICE_SOCN_SROF_SGLC_SWAV --run-unsupported --mach syau

and customized the user_nl_clm:
cat << EOF >> user_nl_${hlm}
fsurdat= '${HLM_SURDAT_PATH}/surfdata_1x1_qianyan_16pfts_Irrig_CMIP6_simyr1850_c230623.nc'
spinup_state = 2
hist_mfilt = 240
hist_nhtfrq = 0
hist_empty_htapes = .true.
hist_fincl1 = 'TOTECOSYSC', 'TOTECOSYSN', 'TOTSOMC', 'TOTSOMN', 'TOTVEGC', 'TOTVEGN', 'TLAI', 'GPP', 'CPOOL', 'NPP', 'TWS', 'H2OSNO','H2OSOI'
EOF

After 100 spin-up, the GPP seems attain an equilibrium state. However, when I check the TLA, all the values are 3-4 instead of the observed range of 5 and 6。Besides, the gpp_24 in the r file, values only show half, some are designated as NAN (Not a Number). Why?

I really don’t know what happened, and any assistance would be greatly appreciated！

Here is the simulated result that I have attached:

 

[oleson]

There could be any number of reasons why TLAI doesn't match up with observations. There could be problems with the input data, the forcing, or the model parameterizations, that's the research project.
The gpp24 on the restart file is generally dimensioned by pft (plant functional type). Maybe your analysis program is interpreting missing values as NaNs. Missing values should be 1.e36, i.e.,

gpp24:_FillValue = 1.e+36 ;
gpp24:missing_value = 1.e+36 ;