how to force air with a SST forcing with excess 1year?

[Yaqi Jin]

Hello everyone, I am preforming a atomosphere experiment forced by SST from 2020-2021. And want to force the model with 2y SST cyclically, ie. the model first year is forced by 2020 SST, and the model second year is forced by 2021. And the model third year is forced by 2020, the fourth year is forced by 2021...... and so forth.

My try:
First, I put a forcing field with sized lonxlatx31year, and its resolution in time is monthly.

Second, I modified the SSTICE_YEAR_START、SSTICE_YEAR_END、SSTICE_YEAR_ALIGN and RUN_STARTDATE as lower.

./xmlchange RUN_STARTDATE=0000-01-01
./xmlchange SSTICE_YEAR_ALIGN=3
./xmlchange SSTICE_YEAR_START=0
./xmlchange SSTICE_YEAR_END=1

Finally , I build it successfully but I cannot run successfully. At the same time the error is showed as lower description.

ocn.log.
(docn_comp_init) ocn_in = docn_ocn_in
(docn_comp_init) decomp = 1d
(docn_comp_init) restfilm = undefined
(docn_comp_init) restfils = undefined
(shr_strdata_readnml) reading input namelist file: docn_ocn_in
(shr_stream_init) Reading file docn.streams.txt.prescribed
(shr_file_freeUnit) unit 94 was not in use
(docn_comp_init) ocn mode = SSTDATA
(shr_strdata_init) calling shr_dmodel_mapSet for fill
(shr_strdata_init) calling shr_dmodel_mapSet for remap
(shr_strdata_print) ----------------------------------------------------------
(shr_strdata_print) name = SDOCN data
(shr_strdata_print) dataMode = SSTDATA
(shr_strdata_print) domainFile = /public/home/yq/cesm1_2_2/inputdata/share/domains/domain.ocn.gx1v6.090206.nc
(shr_strdata_print) nxg = 320
(shr_strdata_print) nyg = 384
(shr_strdata_print) calendar = NO_LEAP
(shr_strdata_print) io_type = 6
(shr_strdata_print) eccen = 1.000000E+36
(shr_strdata_print) mvelpp = 1.000000E+36
(shr_strdata_print) lambm0 = 1.000000E+36
(shr_strdata_print) obliqr = 1.000000E+36
(shr_strdata_print) nstreams = 1
(shr_strdata_print) pio_iotype = 6
(shr_strdata_print) streams ( 1) = docn.streams.txt.prescribed 3 0 1
(shr_strdata_print) taxMode ( 1) = cycle
(shr_strdata_print) dtlimit ( 1) = 1.500000E+00
(shr_strdata_print) strnxg ( 1) = 360
(shr_strdata_print) strnyg ( 1) = 180
(shr_strdata_print) dofill ( 1) = T
(shr_strdata_print) fillalgo( 1) = nn
(shr_strdata_print) fillmask( 1) = nomask
(shr_strdata_print) fillread( 1) = NOT_SET
(shr_strdata_print) fillwrit( 1) = NOT_SET
(shr_strdata_print) domaps ( 1) = T
(shr_strdata_print) mapalgo ( 1) = bilinear
(shr_strdata_print) mapmask ( 1) = dstmask
(shr_strdata_print) mapread ( 1) = NOT_SET
(shr_strdata_print) mapwrit ( 1) = NOT_SET
(shr_strdata_print) tintalgo( 1) = linear
(shr_strdata_print)
(shr_strdata_print) nvectors = 0
(shr_strdata_print) ----------------------------------------------------------
(docn_comp_init) initialize gsmaps
(docn_comp_init) copy domains
(docn_comp_init) allocate AVs
(shr_stream_findBounds) ERROR: LVD not found, all data is after yearLast
(shr_sys_abort) ERROR: (shr_stream_findBounds) ERROR: LVD not found, all data is after yearLast
(shr_sys_abort) WARNING: calling shr_mpi_abort() and stopping


My question:
I want to know "how to run with a forcing field with excess 12 months"?
And how to make it to read cyclically the forcing field, whether need to divided multiple files or a forcing field with a 2y-time dimension is ok?
And how to solve the above error?

Thanks for your help and any advice!!!!!!!!!!!!!!!!!!!!!

 

[dbailey]

I believe you want:

./xmlchange RUN_STARTDATE=2020-01-01
./xmlchange SSTICE_YEAR_ALIGN=2020
./xmlchange SSTICE_YEAR_START=2020
./xmlchange SSTICE_YEAR_END=2021

Then make sure it that taxmode = "cycle" in docn_in.

 

[Yaqi Jin]

I believe you want:

./xmlchange RUN_STARTDATE=2020-01-01
./xmlchange SSTICE_YEAR_ALIGN=2020
./xmlchange SSTICE_YEAR_START=2020
./xmlchange SSTICE_YEAR_END=2021

Then make sure it that taxmode = "cycle" in docn_in.

Thanks for your reply. I am sure that the taxmode had been set into cycle. But I don't know where I had seen it.

But this is not showed in docn_in.
&docn_nml
decomp = '1d'
ocn_in = 'docn_ocn_in'
/


This is information in docn.streams.txt.prescribed.
<dataSource>
GENERIC
</dataSource>
<domainInfo>
<variableNames>
time time
xc lon
yc lat
area area
mask mask
</variableNames>
<filePath>
/public/home/yqjin/cesm1_2_2/inputdata/ocn/docn7
</filePath>
<fileNames>
domain.ocn.1x1.111007.nc
</fileNames>
</domainInfo>
<fieldInfo>
<variableNames>
SST_cpl t
</variableNames>
<filePath>
/public/home/yqjin/cesm1_2_2/inputdata/atm/cam/sst
</filePath>
<fileNames>
sst_HadOIBl_bc_2020_2021_sensitive_EXP.nc
</fileNames>
<offset>
0
</offset>
</fieldInfo>

There is any way to query or set?

 

[Yaqi Jin]

I believe you want:

./xmlchange RUN_STARTDATE=2020-01-01
./xmlchange SSTICE_YEAR_ALIGN=2020
./xmlchange SSTICE_YEAR_START=2020
./xmlchange SSTICE_YEAR_END=2021

Then make sure it that taxmode = "cycle" in docn_in.

Hi , dbailey. I have another question. Should I change the SSTICE_YEAR_START=2020 but not SSTICE_YEAR_START=0 in spite that the time of SST forcing file's date is from year 0.

 

[dbailey]

If you are using a variant of the HadOIBC this should be real years. What if you do:

ncdump -v time sst_HadOIBl_bc_2020_2021_sensitive_EXP.nc

This should give you days since 0001-01-01. I suspect it is not 0.

 

[Yaqi Jin]

If you are using a variant of the HadOIBC this should be real years. What if you do:

ncdump -v time sst_HadOIBl_bc_2020_2021_sensitive_EXP.nc

This should give you days since 0001-01-01. I suspect it is not 0.

Hi dbailey,
Thanks for your reply!
I had solved this problem with your help. For a startup exp, I can set the RUN_STARTDATE and align the model year with it. And give atomosphere a specified SSTA forcing by setting SSTICE_YEAR_START and SSTICE_YEAR_END.
Thanks for your help again.