grdimage using orthographic does not render full dataset
|Target version:||Candidate for next bugfix release|
I'm trying to use -JG (orthographic) and instead of getting the full dataset rendered I'm only getting a small slice; the slice changes depending upon parameters set with -J.
The dataset ought to be fine as it renders in other software such as ncview and panoply.
pscoast and psbasemap also seem to work fine. Script is attached along with sample data and current output.
Using GMT5.2.1 (on a remote machine, outside control to run 5.4.x or later) but also present in 4.5.7.
We've had a look at this but not finished. Seems like your grid is not fully COARDS compliant as it has no dimension called s, just a common dimension called Z for all the grids.
Not sure if it'll help, but I used NCO to strip the s variable of the (in this case) useless dimensions of time and Z (placeholder for elevation) and still get the same result when plotting.
ncwa -a Z,time SSTinst_daily.0000158184.r2160x1080.nc test.nc
It also renders fine in other software (eg. Panoply and ncview).
#3 Updated by Paul about 1 month ago
- Assignee set to Paul
- Target version set to Candidate for next bugfix release
- % Done changed from 0 to 100
- Affected version changed from 5.2.x to 6.x-svn
Sorry so long time. Your grid was too sloppy for our checks that determine if we have a global grid. I have allowed for less precision now and if you append =gd to have it read via GDAL then the plotting works. It was basically a failure to detect a global grid. The grid needs to be read via GDAL since GMT insists on a COARDS netcdf convention that your grid is not following.