SpinupStability_v10.ncl error

[tenghao]

I tried to use ncl SpinupStability_v10.ncl to judge whether the operation has reached a balanced state, but it failed. Prompt:



fatal:conform_dims: The array to be conformed must have the same number of dimensions as indicated by the length of the last argument

fatal:["Execute.c":8637]:Execute: Error occurred at or near line 229 in file SpinupStability_v10.ncl



I have used the atmospheric driving data for three years (2010-2012) and saved it as a result every month. How can I solve this error and make SpinupStability_v10.ncl work normally?




In addition, I saw in the forum that it must be saved as a file every month or every year before it can be run. How to change this situation, for example, save it in an h0 file for 12 months, read the monthly results in this file and run SpinupStability_v10.ncl.

 

[oleson]

How many times have your cycled over the atmospheric data? As the comment in the script says:

You need a run consisting of at least three cycles of atmospheric data.

Is this a single-point simulation?

The script supports either monthly history files, one file per month, or annual average history files. If you want to add a third option, you'lll have to modify the script yourself.

 

[oleson]

It looks like the error is associated with this line:

landareaC = conform_dims(dimsizes(totecosysc),landarea,(/1/)) ; conforming dimensions of landarea to totecosysc

For the single point simulation option you've selected, the variables landfrac and area (which are used to calculate landarea) should have a single dimension named lndgrid, and the variable TOTECOSYSC should have the dimensions (time,lndgrid).

 

[tenghao]

I'm running regional simulation, but I set region=”SPT ".Is this the cause of the problem? For regional simulation, should region be set to Arctic?

My atmospheric forcing data is from 2010 to 2012. I set subper=3, and output the result file for 108 months, which is saved as a file every month, but it has not been run yet.The simulation time should be at least three times as long as the atmospheric forcing data, right?

 

[oleson]

Yes, that's probably the problem. You can try Arctic and set your lat/lons for your region appropriately, although I don't think that option has been tested for awhile, See the discussion beginning Dec 7, 2023 in this post for some possibly helpful information regarding modifying the script for regional runs:

Questions about BGC spinup

Dear all, Hi! I am trying to run a I2000Bgc compset spinup following the section 1.5.5 Spinup of CLM5.0-BGC of Clm5.0 User's guide.I want to spinup using the forcing data from 1979 to 1983.and try to simulate the model year from 1984 to 2014. Here are some setting: ./create_newcase -case...

bb.cgd.ucar.edu

Yes, the simulation time should be at least three times as long as the atmospheric forcing you are looping over.

 

[tenghao]

Thank you very much for your answer. I will try it.

 

[AquaTech_007]

Hi @oleson I have some confusion regarding the subper years. I am running single-point simulations for three years of atmospheric data (2021–2023), and I’ve tried spinning it up for 300 years. However, the variables are still not in an equilibrium state and remain noisy. I’m using a subper as 3 years, should I adjust this? Any insights would be greatly appreciated. Thank you in advance!

Following is the case built, user_nl_clm, SpinupStability.ncl and the graph of variables.

Case Built
./xmlchange STOP_OPTION=nyears,STOP_N=300
./xmlchange CLM_BLDNML_OPTS='-bgc bgc -crop -clm_start_type cold -clm_accelerated_spinup on'
./xmlchange CLM_USRDAT_NAME=1x1_Data_Busso_Bosco
./xmlchangeATM_DOMAIN_FILE=domain.lnd.1x1_Busso_Bosco_Woodland_1x1_Busso_Bosco_Woodland.240227.nc,LND_DOMAIN_FILE=domain.lnd.1x1_Busso_Bosco_Woodland_1x1_Busso_Bosco_Woodland.240227.nc
./xmlchange DATM_CLMNCEP_YR_START=2021
./xmlchange DATM_CLMNCEP_YR_ALIGN=2021
./xmlchange DATM_CLMNCEP_YR_END=2023
./xmlchange RUN_STARTDATE=2021-01-01

user_nl_clm
fsurdat='/home/hydrology/cesm-inputdata/lnd/clm2/surfdata_map/Busso_Bosco_Woodland_hist_78pfts_CMIP6_simyr2000_c240227.nc'

!hist_empty_htapes = .true.
hist_nhtfrq = 0
hist_mfilt = 3600

SpinupStability.ncl
caseid = "BussoTEST_2021_2023_ACC"
username = "hydrology"
annual_hist = True
region = "SPT" ; Global, Arctic, or SPT (single point)
subper = 3 ; Subsampling period in years

 

[AquaTech_007]

The Equilibrium summary follow as:

(0) EQUILIBRIUM SUMMARY
(0) =======================================================================
(0) TOTECOSYSC is in equilibrium. Eq. Yr. = 2021
(0) TOTSOMC is in equilibrium. Eq. Yr. = 2021
(0) TOTVEGC is in equilibrium. Eq. Yr. = 2021
(0) FATAL: TLAI is NOT in equilibrium
(0) GPP is in equilibrium. Eq. Yr. = 2021
(0) WARNING: TWS is NOT in equilibrium or is missing
(0) H2OSNO is in equilibrium. Eq. Yr. = 2021
(0) =========================================
(0) Finish Time: gio 24 lug 2025, 14.49.59, CEST
(0) =========================================
(0) Successfully ran the script

 

[oleson]

You've set annual_hist = True, which requires annual mean (not monthly) history files. However, you've set hist_nhtfrq = 0, which outputs monthly history files, this should be hist_nhtfrq = -8760.