use of -C and -i option with psxy
Langlois <[email protected]>
| Newsgroups | gmane.comp.gis.gmt.user |
|---|---|
| Message-ID | <[email protected]> |
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 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]