Running CLM single point from Restart file

[Esther_Bender]

Hei,

I am trying to run a CLM single point simulation from a restart file, I did a spinup run and now try to use the last restart file as the initial data file for another run, by using
finidat ="restart_file.nc"
in the user_nl_clm file.

for some reason, only some variables are read from the file, the soil temperature is used from the file, but the Carbon values, like TSOM, total ecosystem carbon etc. are starting from some arbitrary values and not the values of the end of the spin up run.

is there any reason for that?
I am using the script below for the run:
"

cd ~/CTSM_ExcessIce/cime/scripts


echo "create case" $1 $2 $3


./create_newcase --machine fram --case ~/cases/$1 --res CLM_USRDAT --compset 1850_DATM%GSWP3v1_CLM50%BGC_SICE_SOCN_SROF_SGLC_SWAV
--run-unsupported --user-mods-dirs /cluster/projects/nn2806k/estherbe/Input_Data/1x1Iskoras --compiler intel --handle-preexisting-dirs r --walltime 01:00:00 -q normal
#--queue devel


cd ~/cases/$1

cp user_mods/user_nl_datm_streams ~/cases/$1

./case.setup

./xmlchange STOP_N=1
./xmlchange STOP_OPTION=nyears
./xmlchange PTS_LON=24.250
./xmlchange PTS_LAT=69.383
./xmlchange RESUBMIT=3
./xmlchange NTASKS=1
./xmlchange RUN_TYPE=startup
./xmlchange JOB_WALLCLOCK_TIME="04:00:00"

./xmlchange RUN_STARTDATE="1901-01-01"
./xmlchange DATM_YR_ALIGN="1901"
./xmlchange DATM_YR_START="1901"
./xmlchange DATM_YR_END="1910"
./xmlchange CLM_ACCELERATED_SPINUP="off"
./xmlchange CLM_FORCE_COLDSTART="off"

./xmlchange DEBUG=FALSE
./xmlchange LND_DOMAIN_PATH=/1x1Iskoras
./xmlchange LND_DOMAIN_FILE="domain.lnd.fv0.9x1.25_gx1v7_Iskoras_c240325.nc"

./xmlchange DIN_LOC_ROOT_CLMFORC=/cluster/projects/nn2806k/estherbe/Input_data/1x1Iskoras


cat > user_nl_clm << EOF
use_excess_ice = .true.
use_excess_ice_tiles = .true.
soil_layerstruct_predefined = '49SL_10m'
fsurdat = '/cluster/home/estherbe/Input_data/Iskoras_geo_ORG_neu_cdf5.nc'
finidat='/cluster/work/users/XX/run/0711_Iskoras_spin_t7.clm2.r.2101-01-01-00000.nc'
hist_dov2xy=.false.
hist_nhtfrq=-0
hist_mfilt=12
hist_fincl1="H2OSOI", "TSOI", "SOILICE", "EXCESS_ICE", "SUBSIDENCE", "QDRAI_PERCH","FROST_TABLE","ZWT_PERCH", \
"ZWT", "H2OSNO", "SNOWDP", "SNOW_5D", "SNOW_ICE", "H2OSFC", "QFLX_LIQ_GRND", "QFLX_SNOW_GRND", "QFLX_EVAP_VEG",\
"QSOIL", "QOVER", "QRGWL", "QSNWCPICE", "FH2OSFC_NOSNOW", "CONC_CH4_SAT", "CONC_CH4_UNSAT",\
"TSOI_10CM", "TG", "QDRAI", "QDRAI_XS", "QICE", "QICE_MELT", "TBOT","TH2OSFC","THBOT","TOTEXICE_VOL",\
"TOTSOILICE","TOTSOILLIQ", "SNOW_DEPTH", "SNOW", "SNOWLIQ", "FSAT", "QH2OSFC", "QOVER", \
"SOILLIQ" , "TWS", "H2OSFC","FGEV" ,"FH2OSFC_NOSNOW", "QFLX_LIQ_GRND", "QFLX_SNOW_GRND", "QFLX_EVAP_VEG",\
"QSOIL", "QOVER", "QRGWL", "QSNWCPICE", "SNOWLIQ", "FSNO", "FSNO_EFF", "FH2OSFC", "FSAT", "H2OSFC", \
"QH2OSFC", "EXCESS_ICE", "TSA","SUBSACC","RAIN", "SNOW","QFLOOD","QIRRIG_FROM_SURFACE",'RAIN_FROM_ATM',\
"QFLX_EVAP_TOT", "NEE", "NEP", "NPP", "GPP", "TWS", "TLAI", "TOTCOLC", "TOTECOSYSC", "TOTSOMC", "TOTVEGC",\
"TOTCOLCH4", "CH4PROD", "CH4_SURF_DIFF_SAT", "CH4_SURF_DIFF_UNSAT", "CONC_CH4_SAT", "CONC_CH4_SAT", "CONC_CH4_UNSAT",\
"FCH4TOCO2", "PCH4", "PCO2", "ZWT_CH4_UNSAT"
EOF

"

 

[oleson]

That does seem odd. But if your initial file is from the end of a spinup run (where CLM_ACCELERATED_SPINUP="on"), and you transition to a normal run, the carbon and nitrogen pools are "de-accelerated". The technical note describes this as: When leaving the accelerated decomposition mode, the concentration of C and N in pools that had been accelerated are multiplied by the same term

, to bring the model into approximate equilibrium. See

2.21. Decomposition — ctsm release-clm5.0 documentation

So I would expect the carbon and nitrogen pools to look quite different, although they shouldn't be "arbitrarily" different.