Re: grdgradient and -R specifications

Perhaps the $grid file was not created as a geographic grid? Does grdinfo say degrees_east and degrees_east as units? If not you might need to add -fg to the command.

-p

On Dec 3, 2013, at 11:18 AM, Cyril Langlois <[email protected]> wrote:

> Hello,
> 
> I want to use the global Bouguer anomalies model WGM2012 available at BGI site (http://bgi.obs-mip.fr/fr/data-products/Grids-and-models/wgm2012) with GMT5.1.0 for the Hawaii-Emperor islands. Then, I want to extract the region -R150/220/0/60 from the grid file. So I write : 
> 
> gmt grdgradient $grid -R150/220/0/60 -A270 -GBouguer_intens.grd -Ne0.6
> 
> where $grid calls WGM2012_Bouguer_ponc_2min.grd
> 
> But the output image shows data only between 150°W and 180°, as if grdgradient could not pass from -180/180 to 0/360 longitudes data. Without -R specification, the map is correct (but the running time is delayed) ! I don’t understand...
> 
> My script is attached below. 
> 
> Thanks for help,
> Cyril
> 
> To unsubscribe, send the message "signoff gmt-help" to [email protected] <GMT5_Hawaii-Empereur_Bouguer_WGM2012.sh>
> 
> Dr. Cyril LANGLOIS
> Professeur agrégé - enseignant de Sciences de la Terre 
> [email protected]
> Tél. portable : 06 89 05 68 22
> Tél. bureau : 04 72 72 85 28
> § --------------------------- §
> Préparation à l'Agrégation SV-STU
> Ecole Normale Supérieure de Lyon -  Site Monod
> 46, allée d'Italie
> 69364 Lyon cedex 07
> § --------------------------- §
> 
> To unsubscribe, send the message "signoff gmt-help" to [email protected]


To unsubscribe, send the message "signoff gmt-help" to [email protected]