Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

'ncar' ocean forcing fault

dpath2o

Daniel Atwater
New Member
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
 
Top