# How to plot fault geometry with fault slip rate in GMT? **URL:** \ **Category:** Q&A\ **Tags:** gmt\ **Created:** [May 7, 2024, 12:02pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872 "2024-05-07T12:02:25Z")\ **Posts on this page:** 13\ **Page:** 1
**Author:** ![AdiMoh08](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/adimoh08/32/4155_2.png) [@AdiMoh08](https://forum.generic-mapping-tools.org/u/AdiMoh08)\ **Post date:** [May 7, 2024, 12:02pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872/1 "2024-05-07T12:02:25Z")
Hi all, Recently I tried to plot fault geometry along with slip rate using GMT but it is not plotted. I am unable to understand where is the error. Please help. This is the script I am using: ```auto !/usr/bin/env sh gmt begin fault_grid.png gmt basemap -R70.2/80.3/28.2/36 -JM10i -Ba2f2 gmt makecpt -Cjet.cpt -G0.4/1 -Ic -T0/22 > fault.cpt gmt coast -R70.2/80.3/28.2/36 -JM10i -B5 -I1/1p,blue -N1/0.25p,black -W0.25p,white -Gyellow -Sblue gmt plot fault.txt -R70.2/80.3/28.2/36 -JM10i -W0.1p -L -Cfault.cpt gmt coast -R70.2/80.2/28.2/36 -JM10i -Bx4 -By2 -Da -W229/195/148 awk '{ if ($5==1 && $6==1) print $8, $9, $12, $17, $15, $20, $27, $1 }' jks0.vsum | gmt velo -R70.2/80.3/28.2/36 -JM10i -Se0.2/0.32/4 -W1p,black -L0.8p,black -A5p+e -Gblack #gmt psconvert $ps -A -Tg #rm $ps gmt.* $cpt gmt end show ``` File format of fault.txt ```auto > -Z 14.398 (Fault slip) 74.96000000 32.66000000 (Lon1,Lat1) 74.89000000 32.69600000 (Lon2,Lat2) 74.91838328 32.73829841 ....... 74.99014893 32.70155361 > -Z 14.290 74.99014893 32.70155361 74.91838328 32.73829841 74.94676657 32.78059681 75.02029785 32.74310723 > -Z 14.185 75.02029785 32.74310723 74.94676657 32.78059681 74.97514985 32.82289522 75.05044678 32.78466084 > -Z 14.084 75.05044678 32.78466084 74.97514985 32.82289522 75.00353313 32.86519362 75.08059570 32.82621445 > -Z 13.985 75.08059570 32.82621445 75.00353313 32.86519362 75.03191642 32.90749203 75.11074463 32.86776807 > -Z 13.891 75.11074463 32.86776807 75.03191642 32.90749203 75.06029970 32.94979044 75.14089356 32.90932168 > -Z 13.799 75.14089356 32.90932168 75.06029970 32.94979044 75.08868298 32.99208884 75.17104248 32.95087529 ``` The plot should somewhat should look like this which I am not getting. ![5c62bfb36da24](https://canada1.discourse-cdn.com/flex047/uploads/gmt/original/2X/7/7966c8d86770e35e38ed5a883ef12c33f5bd7475.png) ---
**Author:** ![Esteban82](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/esteban82/32/54_2.png) [@Esteban82](https://forum.generic-mapping-tools.org/u/Esteban82)\ **Post date:** [May 7, 2024, 12:14pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872/2 "2024-05-07T12:14:35Z")
> [@AdiMoh08](#): > > Recently I tried to plot fault geometry along with slip rate using GMT but it is not plotted. I am unable to understand where is the error. Please help. Hi Adi What is the desire result that would you like achieve? ---
**Author:** ![AdiMoh08](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/adimoh08/32/4155_2.png) [@AdiMoh08](https://forum.generic-mapping-tools.org/u/AdiMoh08)\ **Post date:** [May 7, 2024, 12:28pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872/3 "2024-05-07T12:28:49Z")
![7966c8d86770e35e38ed5a883ef12c33f5bd7475_2_504x500](https://canada1.discourse-cdn.com/flex047/uploads/gmt/original/2X/5/54dda0463bd462f3db9c2fe4c1319bb2082fe49f.png) This has to be the result. I have to plot fault slip magnitude along the fault geometry as shown in this image . ---
**Author:** ![AdiMoh08](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/adimoh08/32/4155_2.png) [@AdiMoh08](https://forum.generic-mapping-tools.org/u/AdiMoh08)\ **Post date:** [May 7, 2024, 12:30pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872/4 "2024-05-07T12:30:20Z")
Its not taking the color variation given in cpt ---
**Author:** ![Esteban82](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/esteban82/32/54_2.png) [@Esteban82](https://forum.generic-mapping-tools.org/u/Esteban82)\ **Post date:** [May 7, 2024, 12:45pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872/5 "2024-05-07T12:45:25Z")
Sorry, I didn’t read well. And what is map that you are gettiing? ---
**Author:** ![PlanetGus](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/planetgus/32/303_2.png) [@PlanetGus](https://forum.generic-mapping-tools.org/u/PlanetGus)\ **Post date:** [May 7, 2024, 1:32pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872/6 "2024-05-07T13:32:49Z")
Is it only the color that is not working ? Try `makecpt [...] -H > ...` ---
**Author:** ![AdiMoh08](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/adimoh08/32/4155_2.png) [@AdiMoh08](https://forum.generic-mapping-tools.org/u/AdiMoh08)\ **Post date:** [May 7, 2024, 3:46pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872/7 "2024-05-07T15:46:38Z")
![image_fault](https://canada1.discourse-cdn.com/flex047/uploads/gmt/original/2X/9/9b359ee7763d2fa819c0716efe051e0ed130e850.png) Whatever cpt I am giving its plotting the outline.No colour variation. its also not giving different colour to land and sea in “pscoast” ---
**Author:** ![Esteban82](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/esteban82/32/54_2.png) [@Esteban82](https://forum.generic-mapping-tools.org/u/Esteban82)\ **Post date:** [May 7, 2024, 4:25pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872/8 "2024-05-07T16:25:33Z")
I managed to get this figure with the following script ```auto #! /usr/bin/env sh gmt begin fault_grid png gmt basemap -R70.2/80.3/28.2/36 -JM10i -Ba2f2 gmt makecpt -Cjet -G0.4/1 -Ic -T0/22 -H > fault.cpt gmt coast -I1/1p,blue -N1/0.25p,black -W0.25p,white -Gyellow -Sblue gmt plot fault.txt -W0.1p -L -Cfault.cpt gmt coast -Da -W229/195/148 #awk '{ if ($5==1 && $6==1) print $8, $9, $12, $17, $15, $20, $27, $1 }' jks0.vsum | gmt velo -R70.2/80.3/28.2/36 - Se0.2/0.32/4 -W1p,black -L0.8p,black -A5p+e -Gblack gmt end ``` ![fault_grid](https://canada1.discourse-cdn.com/flex047/uploads/gmt/original/2X/1/16282881751e54046852ba90ef7b39e43ee66d49.png) ---
**Author:** ![AdiMoh08](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/adimoh08/32/4155_2.png) [@AdiMoh08](https://forum.generic-mapping-tools.org/u/AdiMoh08)\ **Post date:** [May 7, 2024, 4:48pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872/9 "2024-05-07T16:48:53Z")
Yes its plotting the pscoast properly now but its not taking the cpt . ---
**Author:** ![AdiMoh08](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/adimoh08/32/4155_2.png) [@AdiMoh08](https://forum.generic-mapping-tools.org/u/AdiMoh08)\ **Post date:** [May 7, 2024, 4:58pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872/10 "2024-05-07T16:58:45Z")
[fault.txt](https://forum.generic-mapping-tools.org/uploads/short-url/wj9LWsuBa5LYkCzL4IHPqUhwRBb.txt) (50.4 KB) I have attached the fault.txt file which I am using to plot the fault slip magnitude using cpt (not working) . please have a look. ---
**Author:** ![Esteban82](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/esteban82/32/54_2.png) [@Esteban82](https://forum.generic-mapping-tools.org/u/Esteban82)\ **Post date:** [May 7, 2024, 5:07pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872/11 "2024-05-07T17:07:06Z")
> [@AdiMoh08](#): > > ```auto > > -Z 14.398 (Fault slip) > 74.96000000 32.66000000 (Lon1,Lat1) > 74.89000000 32.69600000 (Lon2,Lat2) > 74.91838328 32.73829841 ....... > 74.99014893 32.70155361 > > ``` Ok, now I think I understand. You want the polygons to be painted with the value of the segment header. ---
**Author:** ![Esteban82](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/esteban82/32/54_2.png) [@Esteban82](https://forum.generic-mapping-tools.org/u/Esteban82)\ **Post date:** [May 7, 2024, 5:12pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872/12 "2024-05-07T17:12:37Z")
Just remove the white space in the segment header between the -Z and the value (e.g.: `-Z14.398`). [fault.txt](https://forum.generic-mapping-tools.org/uploads/short-url/yRiiFtC6hgmElUXeILqorrs0NOz.txt) (50,0 KB) ![fault_grid](https://canada1.discourse-cdn.com/flex047/uploads/gmt/original/2X/0/00545f9eb69edb22a804ea01a5f958ffa934c7b7.png) ---
**Author:** ![AdiMoh08](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/adimoh08/32/4155_2.png) [@AdiMoh08](https://forum.generic-mapping-tools.org/u/AdiMoh08)\ **Post date:** [May 7, 2024, 7:37pm UTC](https://forum.generic-mapping-tools.org/t/how-to-plot-fault-geometry-with-fault-slip-rate-in-gmt/4872/13 "2024-05-07T19:37:23Z")
Yes its working now.Thanks alot @Esteban82 and @PlanetGus for the help !! Can you please suggest how can i use gmt script to make such polygons with different cpts ?