I'd like to propose the following change to IOJpeg. It gets rid of annoying error messages when one of the dimensions is one.
if (CCTK_EQUALS(gridpoints, "hyperslab")) {
/* now do the actual I/O looping over all directions */ for (dir = 0; dir < maxdir; dir++) { /* get the directions to span the hyperslab */ if (dir == 0) { dir_i = 0; dir_j = 1; /* xy */ if(GH->cctk_gsh[0]==1 || GH->cctk_gsh[1]==1) continue; } else if (dir == 1) { dir_i = 0; dir_j = 2; /* xz */ if(GH->cctk_gsh[0]==1 || GH->cctk_gsh[2]==1) continue; } else { dir_i = 1; dir_j = 2; /* yz */ if(GH->cctk_gsh[1]==1 || GH->cctk_gsh[2]==1) continue; }
Cheers, Steve