Re: use of -C and -i option with psxy
Paul Wessel <[email protected]>
| Newsgroups | gmane.comp.gis.gmt.user |
|---|---|
| Message-ID | <[email protected]> |
Hi Cyril- Yes, this is a "feature" [http://gmt.soest.hawaii.edu/issues/519] in that you cannot ask for the same column more than once on input. We will be fixing this in the next minor release. For now the workaround is to use gmtconvert to OUTPUT the right columns and pipe to psxy, i.e. gmtconvert $sism -hi -i4,3,8,8s0.025 | psxy $region $projection -Sci -hi -Wdefault -C$depthscale -O -K -V >> $carte -p On Apr 3, 2014, at 7:08 AM, Langlois <[email protected]> wrote: > Dear GMT-users, > > I’m searching to plot french seismicity catalog. First, I had to convert the available text datafile into a .dat file with comma separator between columns (attached file below). Quake depth appear in the 6th column, magnitude in the last one (the 9th). > > So, I mapped the quakes with symbol sizes function of the magnitude and symbol color related to the depth of the quake, using : > > gmt psxy $sism $region $projection -i4,3,5,8s0.025 -Sci -hi -Wdefault -C$depthscale -O -K -V >> $carte > > where $sism call the data file and $depthscale correspond to my depth cpt file. The -i option count the column from zero, so longitude and latitude are the 4th and third columns, depth is the 5th and magnitude the 8th one. Symbols are scaled with magnitude using -i8s0.025 and colored by the -Cdepth.cpt option. > > Then, I tried to set both symbol color and symbol size according to magnitude. I believed that the same command with the option -i4,3,8,8s0.025 and a magnitude-related cpt file -Cmagnitude.cpt could do the job. But it doesn’t work and psxy say : « mismatch between actual (3) and expected (4) columns in the file » > > I don’t understand where the problem. > Thanks for help, > Cyril > > > To unsubscribe, send the message "signoff gmt-help" to [email protected] <France-seismicity.zip> > > 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]