Re: grdcut: Error cutting over +/-180° longitude
Paul Wessel <[email protected]>
| Newsgroups | gmane.comp.gis.gmt.user |
|---|---|
| Message-ID | <[email protected]> |
Usually, one has to ensure that west < east, even if picking corners, so the -124 is causing trouble. It works better with 236 since 236 > 5 while -124 < 5. -p On Mar 23, 2014, at 9:03 AM, Kristof Koch <[email protected]> wrote: > Dear Paul, > > I'd say so: > > $ grdinfo ETOPO1_Ice_g_gmt4.grdETOPO1_Ice_g_gmt4.grd: Title: ETOPO1_Ice_g_gmt4.grd > ETOPO1_Ice_g_gmt4.grd: Command: > ETOPO1_Ice_g_gmt4.grd: Remark: > ETOPO1_Ice_g_gmt4.grd: Gridline node registration used [Geographic grid] > ETOPO1_Ice_g_gmt4.grd: Grid file format: ni = GMT netCDF format (32-bit integer), COARDS, CF-1.5 > ETOPO1_Ice_g_gmt4.grd: x_min: -180 x_max: 180 x_inc: 0.0166666666667 name: Longitude [degrees_east] nx: 21601 > ETOPO1_Ice_g_gmt4.grd: y_min: -90 y_max: 90 y_inc: 0.0166666666667 name: Latitude [degrees_north] ny: 10801 > ETOPO1_Ice_g_gmt4.grd: z_min: -10898 z_max: 8271 name: z > ETOPO1_Ice_g_gmt4.grd: scale_factor: 1 add_offset: 0 > ETOPO1_Ice_g_gmt4.grd: format: classic > > The Version of GMT I'm using is 5.1.1_r12890 [64-bit] on Ubuntu 12.04.4 LTS (GNU/Linux 2.6.32-042stab084.17 x86_64) > > Am 23.03.2014 um 19:52 schrieb Paul Wessel <[email protected]>: > >> Does grdinfo say your etopo1 grid is geographic? If not you will have to add -fg. >> >> Sent from my iPad >> >>> On Mar 23, 2014, at 12:31 AM, Kristof Koch <[email protected]> wrote: >>> >>> Dear List, >>> >>> I'm trying to cut the ETOPO1 grid. My script works as desired as long as -R doesn't span over +/-180° longitude, however I get the following error trying to cut over +/-180° longitude: >>> >>> grdcut: Error: grid x range <= 0.0 >>> grdcut (GMT_grdcut): Unacceptable PLOT_DEGREE_FORMAT template. A not allowed [/tmp/apollodor-8GOUnDRc/tmp_elevation_data_raw.grd] >>> >>> the relevant snippet of code: >>> >>> TERRAIN_BOUNDARIES=75/17/-124/79 >>> >>> grdcut \ >>> ETOPO1_Ice_g_gmt4.grd \ >>> -G"${WORK_DIR}"/tmp_elevation_data_raw.grd \ >>> -R"${TERRAIN_BOUNDARIES}"r >>> >>> After studying the grdcut manpage, I'm not quite sure how to make it work. Any hints are welcome. >>> >>> All the best, >>> Kristof Koch >>> >>> -- >>> My PGP Public Key: http://kristofkoch.de/pgp/ >>> >>> >>> To unsubscribe, send the message "signoff gmt-help" to [email protected] >> >> To unsubscribe, send the message "signoff gmt-help" to [email protected] > > > -- > My PGP Public Key: http://kristofkoch.de/pgp/ > > > To unsubscribe, send the message "signoff gmt-help" to [email protected] To unsubscribe, send the message "signoff gmt-help" to [email protected]