grdgradient and -R specifications

Cyril Langlois <[email protected]>
Newsgroups gmane.comp.gis.gmt.user
Message-ID <[email protected]>
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



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]
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.