Error when running F_1850_WACCM compset using sstice_ts.nc dataset from B compset

[andrew_kren@colorado_edu]

Hi,I ran into an issue when running my F_1850_WACCM compset. I first followed the instructions to create an sstice.nc file from a B compset using the directions provided in the CESM user guide. My resulting dataset is a 200-yr run of monthly mean ice and sst data. I believe I put the data in the right format as it runs for the first 15 days. On the 16th day it gives me this error in the ice log file:arwt salt flx in dt(kg)=    0.00000000000000000E+00   0.00000000000000000E+00
arwt salt flx error    =   -2.23838047746190892E-06  -1.35363960735993708E-05
 ----------------------------
 
istep1:       744    idate:     10116    sec:     43200
                                             Arctic                 Antarctic
total ice area  (km^2) =    1.45002526842200644E+07   8.81659065380526148E+06
total ice extent(km^2) =    1.54457210554687399E+07   1.43049566412500646E+07
total ice volume (m^3) =    2.90005053684401289E+13   8.81659065380526172E+12
total snw volume (m^3) =    2.46290928533233984E+12   1.05295144771957654E+12
tot kinetic energy (J) =    0.00000000000000000E+00   0.00000000000000000E+00
rms ice speed    (m/s) =        0.00000000000000000       0.00000000000000000
average albedo         =        0.82451510319631771       0.78852941649317887
max ice volume     (m) =        2.00000000000000000       1.00000000000000000
max ice speed    (m/s) =        0.00000000000000000       0.00000000000000000
max strength    (kN/m) =        0.00000000000000000       0.00000000000000000
 ----------------------------
arwt rain h2o kg in dt =    6.92986077532349091E+10   7.07287948044940186E+10
arwt snow h2o kg in dt =    2.44500870968654663E+11   3.16646818289773743E+11
arwt evap h2o kg in dt =    7.68122293026825309E+08  -6.33337438266671143E+10
arwt frzl h2o kg in dt =    0.00000000000000000E+00   0.00000000000000000E+00
arwt frsh h2o kg in dt =    2.61076582570898781E+10   2.14736284309985077E+11
arwt ice mass (kg)     =    2.65934634228596000E+16   8.08481362953942500E+15
arwt snw mass (kg)     =    8.12760064159672125E+14   3.47473977747460250E+14
arwt tot mass (kg)     =    2.74062234870192720E+16   8.43228760728688500E+15
arwt tot mass chng(kg) =    1.83956426256000000E+11  -4.25921078380000000E+10
arwt water flux        =    2.88459942757826538E+11   1.09305584957615570E+11
 (=rain+snow+evap+frzl-fresh) 
water flux error       =    3.81313085881109280E-06   1.80138178238074970E-05
 ----------------------------
arwt atm heat flux (W) =   -3.80268436595113812E+14   6.66410829579165625E+13
arwt ocn heat flux (W) =    0.00000000000000000E+00   0.00000000000000000E+00
arwt frzl heat flux(W) =    0.00000000000000000E+00   0.00000000000000000E+00
arwt tot energy    (J) =   -8.88504311304809559E+21  -2.48431294458485801E+21
arwt net heat      (J) =   -6.84483185871204864E+17   1.19953949324249808E+17
arwt tot energy chng(J)=   -6.57350359909400576E+17   6.27190052909219840E+16
arwt heat error        =    3.05376413108884984E-06  -2.30385403570370614E-05
 ----------------------------
arwt salt mass (kg)    =    1.06373853691438406E+14   3.23392545181576992E+13
arwt salt mass chng(kg)=   -1.35145014048000008E+08  -4.35730259048000038E+08
arwt salt flx in dt(kg)=    0.00000000000000000E+00   0.00000000000000000E+00
arwt salt flx error    =   -1.27047210717793143E-06  -1.34737261430480631E-05
 ----------------------------
(shr_dmodel_readLBUB) reading file: /projects/ankr3992/inputdata//atm/cam/sst/sstice_ts.nc       2
 (shr_strdata_advance) ERROR: dt limit1    275.500000000000    
   29.5000000000000        1.50000000000000    
 (shr_strdata_advance) ERROR: dt limit2        10116       43200       10215
           0
(shr_sys_abort) ERROR: (shr_strdata_advance) ERROR dt limit
(shr_sys_abort) WARNING: calling shr_mpi_abort() and stopping

Does this have something to do with my time array in the sstice_ts.nc file, meaning is 10215 incorrect for the mid-month time? Should it be 10214? Or is there something else wrong?

 

[santos]

There is not an officially supported way of doing this (despite the user's guide entry), and I'm not sure that I know enough about this error to be much help in any case. I can tell you to check the time coordinate on the file, versus the SSTICE_YEAR_START, SSTICE_YEAR_ALIGN, and start date.I'm sure that the bcgen tool is relevant here; here's a recent post about that:http://bb.cgd.ucar.edu/problem-recipe-using-b-compset-output-create-sstice-forcing-files-f-compset

 

[andrew_kren@colorado_edu]

Yes I am aware of that post and had some back and forth with Sean Davis. What should the year start and year align be?

 

[santos]

START should be the first data set year you want to use (probably the first year on the file), and ALIGN should be the model year that corresponds to START (probably the year in RUN_STARTDATE).

 

[andrew_kren@colorado_edu]

Ok I fixed my run start date and this time the model ran till April 16 and then stopped with the same dtlimit error. By any chance do you know what may be going on? In the model error, it appears that when it got to April 16, it was looking for the data in 10116, meaning the next year, which should not be right. Any suggestions for how to debug would be most appreciated.

 

[santos]

It's not clear to me, especially as I'm not familiar with this code. When you say that 01 01 16 is the next year, does that mean that you are starting the model at year 0? Or do you mean that you believe that this is the date it is searching for in the file, and that the file starts at year 0? (Of course, either way you'd expect it to be looking for xx0515 instead.)