Error with mppncscatter when splitting high dimension variables

[philenso]

Hi all,

I am trying to scatter the SIS2 restart file ice_model.res.nc with mppncscatter:

#mppncscatter -X lonh -Y lath -x 80 -y 32 -i 4 -j 4 ice_model.res.nc

I got an error message like this:

Error. Failed to get var enth_snow double data from input file.

But it works well when splitting the MOM6 restart file MOM.res.nc.

It seems that mppncscatter fails to split variables with dimensions higher than 4.

Following are some information of the files I try to splitted:

#ls -lh MOM.res.nc ice_model.res.nc
-rw-r--r-- 1 philenso climate 24G Mar 21 01:56 MOM.res.nc
-rw-r--r-- 1 philenso climate 33G Mar 21 02:09 ice_model.res.nc

#ncdump -h ice_model.res.nc
netcdf ice_model.res {
dimensions:
xaxis_1 = 7200 ;
yaxis_1 = 5400 ;
zaxis_1 = 6 ;
zaxis_2 = 5 ;
aaxis_1 = 1 ;
aaxis_2 = 4 ;
Time = UNLIMITED ; // (1 currently)
variables:
double xaxis_1(xaxis_1) ;
xaxis_1:long_name = "xaxis_1" ;
xaxis_1:units = "none" ;
xaxis_1:cartesian_axis = "X" ;
double yaxis_1(yaxis_1) ;
yaxis_1:long_name = "yaxis_1" ;
yaxis_1:units = "none" ;
yaxis_1:cartesian_axis = "Y" ;
double zaxis_1(zaxis_1) ;
zaxis_1:long_name = "zaxis_1" ;
zaxis_1:units = "none" ;
zaxis_1:cartesian_axis = "Z" ;
double zaxis_2(zaxis_2) ;
zaxis_2:long_name = "zaxis_2" ;
zaxis_2:units = "none" ;
zaxis_2:cartesian_axis = "Z" ;
double aaxis_1(aaxis_1) ;
aaxis_1:long_name = "aaxis_1" ;
aaxis_1:units = "none" ;
aaxis_1:cartesian_axis = "N" ;
double aaxis_2(aaxis_2) ;
aaxis_2:long_name = "aaxis_2" ;
aaxis_2:units = "none" ;
aaxis_2:cartesian_axis = "N" ;
double Time(Time) ;
Time:long_name = "Time" ;
Time:units = "time level" ;
Time:cartesian_axis = "T" ;
double flux_u(Time, yaxis_1, xaxis_1) ;
flux_u:long_name = "flux_u" ;
flux_u:units = "none" ;
flux_u:checksum = "F41D2749BC340B41" ;
double flux_v(Time, yaxis_1, xaxis_1) ;
flux_v:long_name = "flux_v" ;
flux_v:units = "none" ;
flux_v:checksum = "A455A8D45ACCA49D" ;
double flux_t(Time, yaxis_1, xaxis_1) ;
flux_t:long_name = "flux_t" ;
flux_t:units = "none" ;
flux_t:checksum = "C0AE999C38ABBF83" ;
double flux_q(Time, yaxis_1, xaxis_1) ;
flux_q:long_name = "flux_q" ;
flux_q:units = "none" ;
flux_q:checksum = "400D166E70141F76" ;
double flux_salt(Time, yaxis_1, xaxis_1) ;
flux_salt:long_name = "flux_salt" ;
flux_salt:units = "none" ;
flux_salt:checksum = "49428C23F69D7C37" ;
double flux_lw(Time, yaxis_1, xaxis_1) ;
flux_lw:long_name = "flux_lw" ;
flux_lw:units = "none" ;
flux_lw:checksum = "8591002E7D4EC0C9" ;
double lprec(Time, yaxis_1, xaxis_1) ;
lprec:long_name = "lprec" ;
lprec:units = "none" ;
lprec:checksum = "B5BD8DBCCFA61593" ;
double fprec(Time, yaxis_1, xaxis_1) ;
fprec:long_name = "fprec" ;
fprec:units = "none" ;
fprec:checksum = "CA7D4B12A7DC11CE" ;
double runoff(Time, yaxis_1, xaxis_1) ;
runoff:long_name = "runoff" ;
runoff:units = "none" ;
runoff:checksum = "5E6D2BAE809BA463" ;
double calving(Time, yaxis_1, xaxis_1) ;
calving:long_name = "calving" ;
calving:units = "none" ;
calving:checksum = " 0" ;
double runoff_hflx(Time, yaxis_1, xaxis_1) ;
runoff_hflx:long_name = "runoff_hflx" ;
runoff_hflx:units = "none" ;
runoff_hflx:checksum = " 0" ;
double calving_hflx(Time, yaxis_1, xaxis_1) ;
calving_hflx:long_name = "calving_hflx" ;
calving_hflx:units = "none" ;
calving_hflx:checksum = " 0" ;
double p_surf(Time, yaxis_1, xaxis_1) ;
p_surf:long_name = "p_surf" ;
p_surf:units = "none" ;
p_surf:checksum = "6A5D5D1025784467" ;
double flux_sw_vis_dir(Time, yaxis_1, xaxis_1) ;
flux_sw_vis_dir:long_name = "flux_sw_vis_dir" ;
flux_sw_vis_dir:units = "none" ;
flux_sw_vis_dir:checksum = "89FFDE96CAF191BB" ;
double flux_sw_vis_dif(Time, yaxis_1, xaxis_1) ;
flux_sw_vis_dif:long_name = "flux_sw_vis_dif" ;
flux_sw_vis_dif:units = "none" ;
flux_sw_vis_dif:checksum = "A2352A2A0A556701" ;
double flux_sw_nir_dir(Time, yaxis_1, xaxis_1) ;
flux_sw_nir_dir:long_name = "flux_sw_nir_dir" ;
flux_sw_nir_dir:units = "none" ;
flux_sw_nir_dir:checksum = "95718612960C588E" ;
double flux_sw_nir_dif(Time, yaxis_1, xaxis_1) ;
flux_sw_nir_dif:long_name = "flux_sw_nir_dif" ;
flux_sw_nir_dif:units = "none" ;
flux_sw_nir_dif:checksum = "95718612960C588E" ;
double part_size(Time, zaxis_1, yaxis_1, xaxis_1) ;
part_size:long_name = "part_size" ;
part_size:units = "none" ;
part_size:checksum = "BA314C406557F6AD" ;
double h_pond(Time, zaxis_2, yaxis_1, xaxis_1) ;
h_pond:long_name = "h_pond" ;
h_pond:units = "H_to_kg_m2 kg m-2" ;
h_pond:checksum = " 0" ;
double h_snow(Time, zaxis_2, yaxis_1, xaxis_1) ;
h_snow:long_name = "h_snow" ;
h_snow:units = "H_to_kg_m2 kg m-2" ;
h_snow:checksum = "CD455862CC768865" ;
double enth_snow(Time, aaxis_1, zaxis_2, yaxis_1, xaxis_1) ;
enth_snow:long_name = "enth_snow" ;
enth_snow:units = "none" ;
enth_snow:checksum = "358802D65F223A58" ;
double h_ice(Time, zaxis_2, yaxis_1, xaxis_1) ;
h_ice:long_name = "h_ice" ;
h_ice:units = "H_to_kg_m2 kg m-2" ;
h_ice:checksum = "6851D63E214905D1" ;
double H_to_kg_m2(Time) ;
H_to_kg_m2:long_name = "The conversion factor from SIS2 mass-thickness units to kg m-2." ;
H_to_kg_m2:units = "none" ;
H_to_kg_m2:checksum = "3FF0000000000000" ;
double enth_ice(Time, aaxis_2, zaxis_2, yaxis_1, xaxis_1) ;
enth_ice:long_name = "enth_ice" ;
enth_ice:units = "J kg-1" ;
enth_ice:checksum = "F0BFA808B6ADE462" ;
double sal_ice(Time, aaxis_2, zaxis_2, yaxis_1, xaxis_1) ;
sal_ice:long_name = "sal_ice" ;
sal_ice:units = "kg/kg" ;
sal_ice:checksum = "FEC75016E224617B" ;
double u_ice_C(Time, yaxis_1, xaxis_1) ;
u_ice_C:long_name = "u_ice_C" ;
u_ice_C:units = "none" ;
u_ice_C:checksum = "F83D0D857BFFF60C" ;
double v_ice_C(Time, yaxis_1, xaxis_1) ;
v_ice_C:long_name = "v_ice_C" ;
v_ice_C:units = "none" ;
v_ice_C:checksum = "48771EAB3292BD49" ;
double m_to_Z(Time) ;
m_to_Z:long_name = "The conversion factor from m to SIS2 height units." ;
m_to_Z:units = "Z meter-1" ;
m_to_Z:checksum = "3FF0000000000000" ;
double m_to_L(Time) ;
m_to_L:long_name = "The conversion factor from m to SIS2 length units." ;
m_to_L:units = "L meter-1" ;
m_to_L:checksum = "3FF0000000000000" ;
double s_to_T(Time) ;
s_to_T:long_name = "The conversion factor from s to SIS2 time units." ;
s_to_T:units = "T second-1" ;
s_to_T:checksum = "3FF0000000000000" ;
double kg_m3_to_R(Time) ;
kg_m3_to_R:long_name = "The conversion factor from kg m-3 to SIS2 density units." ;
kg_m3_to_R:units = "R m3 kg-1" ;
kg_m3_to_R:checksum = "3FF0000000000000" ;
double str_d(Time, yaxis_1, xaxis_1) ;
str_d:long_name = "str_d" ;
str_d:units = "none" ;
str_d:checksum = "C0D2E0428DDBF3B1" ;
double str_t(Time, yaxis_1, xaxis_1) ;
str_t:long_name = "str_t" ;
str_t:units = "none" ;
str_t:checksum = "CADE369FB2602A04" ;
double str_s(Time, yaxis_1, xaxis_1) ;
str_s:long_name = "str_s" ;
str_s:units = "none" ;
str_s:checksum = "DAFE63A92135418D" ;
double rough_mom(Time, zaxis_1, yaxis_1, xaxis_1) ;
rough_mom:long_name = "rough_mom" ;
rough_mom:units = "none" ;
rough_mom:checksum = "E028F13F27773323" ;
double rough_heat(Time, zaxis_1, yaxis_1, xaxis_1) ;
rough_heat:long_name = "rough_heat" ;
rough_heat:units = "none" ;
rough_heat:checksum = "36D80D215EEC1E1F" ;
double rough_moist(Time, zaxis_1, yaxis_1, xaxis_1) ;
rough_moist:long_name = "rough_moist" ;
rough_moist:units = "none" ;
rough_moist:checksum = "53577FBA322A06C6" ;
double coszen(Time, yaxis_1, xaxis_1) ;
coszen:long_name = "coszen" ;
coszen:units = "none" ;
coszen:checksum = " 0" ;
double T_skin(Time, zaxis_2, yaxis_1, xaxis_1) ;
T_skin:long_name = "T_skin" ;
T_skin:units = "none" ;
T_skin:checksum = "757577E352E88C5B" ;

#ncdump -h MOM.res.nc
netcdf MOM.res {
dimensions:
lath = 5400 ;
lonh = 7200 ;
Layer = 75 ;
Time = UNLIMITED ; // (1 currently)
variables:
double lath(lath) ;
lath:long_name = "Latitude" ;
lath:units = "degrees_north" ;
lath:cartesian_axis = "Y" ;
double lonh(lonh) ;
lonh:long_name = "Longitude" ;
lonh:units = "degrees_east" ;
lonh:cartesian_axis = "X" ;
double Layer(Layer) ;
Layer:long_name = "Layer z-rho" ;
Layer:units = "meter" ;
Layer:cartesian_axis = "Z" ;
Layer:positive = "up" ;
double Time(Time) ;
Time:long_name = "Time" ;
Time:units = "days" ;
Time:cartesian_axis = "T" ;
double Temp(Time, Layer, lath, lonh) ;
Temp:long_name = "Potential Temperature" ;
Temp:units = "degC" ;
Temp:checksum = "2FC28A5AA6EF48CF" ;