Bug #1140

grdimage using orthographic does not render full dataset

Added by Ryan about 1 year ago. Updated 8 months ago.

Status:ClosedStart date:2017-08-01
Priority:NormalDue date:
Assignee:Paul% Done:

100%

Category:-
Target version:Candidate for next bugfix release
Affected version:6.x-svn Platform:Linux

Description

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.

SSTinst_daily.0000158184.r2160x1080.gTop.pdf - Sample buggy output pdf (378 KB) Ryan, 2017-08-01 07:10

SSTinst_daily.0000158184.r2160x1080.nc - Sample input file (9 MB) Ryan, 2017-08-01 07:11

sstPlot.sh Magnifier - Sample runscript (552 Bytes) Ryan, 2017-08-01 07:12

History

#1 Updated by Paul about 1 year ago

  • Status changed from New to Feedback

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.

#2 Updated by Ryan about 1 year ago

Paul wrote:

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 10 months 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.

#4 Updated by Joaquim 8 months ago

  • Status changed from Feedback to Closed

Fixed

Also available in: Atom PDF