'ncar' ocean forcing fault

[dpath2o]

Hello,

I am attempting to force CICE6 on a global 1/10-degree tripolar grid, with a global atmosphere and ocean datasets re-gridded to that same grid.

Amongst other machine specific inputs in the ice_in file, critical to this post/help request, I believe, is that I am using atm_data_type set to 'JRA55_tx1' and ocn_data_type set to 'ncar'. This crashes with the following error (truncated):

Atmospheric data files:
/scratch/jk72/da1339/cice-dirs/input/AFIM/forcing/0p1/JRA55/8XDAILY/JRA55_gx3_0
3hr_forcing_2005.nc
Set current forcing data year = 2005
(JRA55_data) reading forcing file 1st ts =
/scratch/jk72/da1339/cice-dirs/input/AFIM/forcing/0p1/JRA55/8XDAILY/JRA55_gx3_0
3hr_forcing_2005.nc
(JRA55_data)fdbg JRA55_bulk_data
(JRA55_data)fdbg fsw 0.000000000000000E+000 1116.57788085938
(JRA55_data)fdbg flw 109.821899414062 449.184173583984
(JRA55_data)fdbg fsnow 0.000000000000000E+000 6.004014518111944E-003
(JRA55_data)fdbg Tair 229.463770548503 305.760060628255
(JRA55_data)fdbg uatm -21.5435422261556 21.5619373321533
(JRA55_data)fdbg vatm -22.8203296661377 21.7584559122721
(JRA55_data)fdbg Qa 5.210744348005392E-005 2.182209864258766E-002
ocn_data_ncar
Tf -2.95014794332433 -1.079747287104882E-003
sst -2.17914343392977 32.8922670836090
sss 1.995333172230711E-002 51.7001094766842
hmix 20.0000000000000 2942.82056189096
uocn -1.83191154240280 1.67568895268825
vocn -1.64283971905067 1.71251867599385
ss_tltx -5.810538027816875E-004 1.654350355434713E-004
ss_tlty -6.879669175987391E-004 1.655585753416640E-004
qdp -332.760181796166 329.432021561489

Finished writing ./history/iceh_ic.2005-01-01-03600.nc
(ice_memusage_print) memory use (MB) = -0.0010 1877.2145 (cice_init):end
[gadi-mmem-clx-0004:772730:0:773275] Caught signal 11 (Segmentation fault: Sent by the kernel at address (nil))
.
.

.
.

--------------------------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
forrtl: error (78): process killed (SIGTERM)
Image PC Routine Line Source
libpthread-2.28.s 000014F7B96C4CE0 Unknown Unknown Unknown
libuct.so.0.0.0 000014F7B296D230 Unknown Unknown Unknown
libucp.so.0.0.0 000014F7B2BD6E9A ucp_worker_progre Unknown Unknown
libopen-pal.so.40 000014F7B306AB83 opal_progress Unknown Unknown
libmpi.so.40.30.3 000014F7BA2D8CB5 ompi_request_defa Unknown Unknown
libmpi.so.40.30.3 000014F7BA28D1F7 PMPI_Waitall Unknown Unknown
libmpi_mpifh.so 000014F7BA6108AD Unknown Unknown Unknown
cice 00000000004420F1 ice_boundary_mp_i 1386 ice_boundary.F90
cice 00000000008A3585 ice_transport_rem 576 ice_transport_remap.F90
cice 000000000088F64E ice_transport_dri 545 ice_transport_driver.F90
cice 00000000004142AC cice_runmod_mp_ci 297 CICE_RunMod.F90
cice 000000000040E806 MAIN__ 49 CICE.F90
cice 000000000040E79D Unknown Unknown Unknown
libc-2.28.so 000014F7B9327CF3 __libc_start_main Unknown Unknown
cice 000000000040E6BE Unknown Unknown Unknown

However, when I run using atm_data_type set to 'JRA55_tx1' and ocn_data_type set to 'default', the run completes without any errors.

My first assumption, is that there must be something wrong with my ocean forcing file, but I have gone over it, and can find no faults with it. Here is what the header of that file looks like:

netcdf acom2_ocn_frcg_cice6_2005 {
dimensions:
ni = 3600 ;
nj = 2700 ;
time = UNLIMITED ; // (31 currently)

variables:
float S(time,nj,ni) ;
S:_FillValue = -2.e+16f ;
S:units = "psu" ;
S:long_name = "sea surface salinity" ;
S:coordinates = "xc yc" ;

float T(time,nj,ni) ;
T:_FillValue = -2.e+16f ;
T:units = "C" ;
T:long_name = "sea surface temperature" ;
T:coordinates = "xc yc" ;

float U(time,nj,ni) ;
U:_FillValue = -2.e+16f ;
U:units = "m/s" ;
U:long_name = "meridional surface current" ;
U:coordinates = "xc yc" ;

float V(time,nj,ni) ;
V:_FillValue = -2.e+16f ;
V:units = "m/s" ;
V:long_name = "zonal surface current" ;
V:coordinates = "xc yc" ;

double dhdx(time,nj,ni) ;
dhdx:_FillValue = -2.e+16 ;
dhdx:units = "m" ;
dhdx:long_name = "zonal sea surface slope" ;
dhdx:coordinates = "xc yc" ;

double dhdy(time,nj,ni) ;
dhdy:_FillValue = -2.e+16 ;
dhdy:units = "m" ;
dhdy:long_name = "meridional sea surface slope" ;
dhdy:coordinates = "xc yc" ;

float hblt(time,nj,ni) ;
hblt:_FillValue = -2.e+16f ;
hblt:units = "m" ;
hblt:long_name = "mixed layer depth" ;
hblt:coordinates = "xc yc" ;

float qdp(time,nj,ni) ;
qdp:_FillValue = -2.e+16f ;
qdp:units = "W/m^2" ;
qdp:long_name = "deep ocean heat flux" ;
qdp:coordinates = "xc yc" ;

int64 time(time) ;
time:long_name = "time" ;
time:cartesian_axis = "T" ;
time:calendar_type = "GREGORIAN" ;
time:units = "days since 2005-01-01 12:00:00" ;
time:calendar = "proleptic_gregorian" ;

double xc(nj,ni) ;
xc:_FillValue = NaN ;
xc:units = "degrees_east" ;

double yc(nj,ni) ;
yc:_FillValue = NaN ;
yc:units = "degrees_north" ;
} // group /

Does anyone have an idea where the problem might lie?

Thanks,
Dan

 

[dbailey]

Oh, is your ocean data file on the tx1 grid? That is 360x240 grid points. CICE standalone does not interpolate.

 

[dpath2o]

Rats! It looks like my reply did not stick. Hopefully this one does as I would really love some help here nutting this problem out. (I'm presently in Antarctica, Casey Station and the bandwidth is pretty poor).

Thanks David.

The ocean *and* atmosphere datasets are *not* on the Tx1 grid. They are both on a 3600x2700 global tripole grid (along with bathymetry and coastal masking) which is being supplied to CICE. Again, CICE is happily working when *just* supplied the atmosphere, it is only when I attempt to 'turn on' the 'ncar' ocean dataset.