Seriestime GMT_part2

Q&A
1

hello this is my code

gmt begin series_times pdf,png

gmt set GMT_THEME cookbook

gmt set  FONT_ANNOT_PRIMARY +9p

gmt set  FORMAT_DATE_MAP=-o    

gmt set GMT_LANGUAGE IT

gmt set IO_COL_SEPARATOR=tab

gmt info C:\programs\Instance_Sicilia_A.csv

gmt plot C:\programs\Instance_Sicilia_A.csv -i3,4 -Sc0.15c -Gred -R2005-01-01T/2020-01-30T/0.0/6.5  -JX10i   -Bxafg+l"Time" -Byafg+l"Magnitude" -fi3T,2d

gmt end show

i want understand why my points no its align for time label.

this is what i received with this code

Screenshot (747)
Screenshot (747)1920Ă—1080 140 KB

and than this is want i want…

series_times
series_times2248Ă—2223 62.4 KB

in .csv i have
35911 13.3917 38.5727 2012-07-11T10:07:40.96Z 2.9 11.0
long lat time mag depth

2

The -f option looks suspicious, why is there a 2d there? You do not have any columns in cm or inches as far as I can tell.

1 Like
3

I increment my code with this
gmt set FONT_ANNOT_PRIMARY +9p

GMT set FORMAT_DATA_IN 2005-01-01T

gmt set FORMAT_DATA_OUT 2005-01-01T

gmt set FORMAT_DATE_MAP=-yyyy

gmt set TIME_EPOCH 2005-01-01T

and in gmt plot i remove +2d…but i dont receive my result…

4

What’s the problem exactly ?

The only obvious difference between your plot and “this is what I want” is the y-axis limits.

-R2005-12-01T/2021-01-31T/1/4.1

no?
(don’t hesitate to use tools such as deepl to better explain and translate what you really want)

Best,
G.

5

“What I really want” is an example to make you understand.

The code is the reality of my plot.

I want the points to be perfectly aligned as in the plot example.

I honestly don’t understand why they don’t line up…non-aligned is meaningless.

Surely there is something wrong with the code, but I don’t understand what.
I have quoted the first line of the .csv file so you can see how the columns are combined.

6

Do you mean vertically aligned ?

Your data have time records down to 0.01 seconds (eg. 2012-07-11T10:07:40.96)

What GMT is doing is placing this time coordinate as precisely as resolution allows.

So a point in July 7th and a point the very next day (July 8th) won’t be aligned.

If I understood correctly and you need the points to be aligned, I suggest

  • using a different module (such as histogram)
    -or-
  • specify the entry to be rounded (everything at July 1st for anything in July for example)

Is it sensible or am I off track ?

7

aligned as in the second screenshot you see.

In this way I will have a time series that allows me to see in that year what earthquakes, and therefore what magnitude, I had.

8

Can you try these lines ?

gmt begin series_times png

gmt math C:\programs\Instance_Sicilia_A.csv --FORMAT_CLOCK_OUT=- --FORMAT_DATE_OUT=yyyy-mm = test.txt
gmt plot  test.txt -i3,4 -Sc0.15c -Gred -R2005-01-01T/2020-01-30T/0.0/6.5  -JX10i   -Bxafg+l"Time" -Byafg+l"Magnitude"

gmt end show

First one converts your csv file into text file with rounded timestamps,
Second line is your plot

9

Screenshot (749)
Screenshot (749)1920Ă—1080 131 KB
Uploading: Screenshot (749).png…

I obtain this.
It is the same

I dont change anythings

10

Well, I have “aligned” them per month. If you want even coarser resolution (aligned per year), just remove mm in --FORMAT_DATE_OUT=yyyy-mm

11

Screenshot (750)
Screenshot (750)1920Ă—1080 113 KB

YEAH!!!
It’s fantastic.

do you know how I can add the years 2005 or 2007, which don’t appear yet? there seems to be an interval of 2…

12

-Bxa1Yf1Yg1Y+l"Time" instead of automatic -Bxafg maybe ?

Like I said last time, a controls annotations … without further instruction GMT defines the interval automatically, else it uses your steps (here 1y = annotation every year)
f is for ticks (frame) and g for grid-line

You should play with the tutorial : https://docs.generic-mapping-tools.org/dev/tutorial/session-1.html#logarithmic-projection

13

Thank you very much.
Yes I know, I should play it more…only then I fail that I try several times…and I ask you to help me out of desperation. I’m new to this world… I’m getting the hang of it little by little.
Thank you for your prompt and sincere help.

1 Like
14

Same at the beginning of my PhD too :slight_smile:

1 Like
15

Hi, I wanted to ask you for confirmation…about this line of code which is similar to everything you have seen above.

-Bx15ua1Yfg+l “Time”

in this case, am I telling it to select my data every 15 days making the axis stay with a year to year interval?

Because if I change and put everything with 15u it’s a mess.

I need to have a bigger window, 1 or 2 weeks…for the time series.

As in this example…

4-Figure2-1
4-Figure2-1916Ă—1216 81.2 KB

16

You want to “hide” your data ? Or do you want take an average value per 2 weeks window ?

The order for -B parameters is [letter]value. In the code you shared the “15u” is doing nothing … I think.

17

No I want to represent all the data…the time window must be enlarged.

Instead of having plotted are years… I need a window of 2 weeks and not every year.

The same in the screenshot.

18

You mean the grid lines ?

-Ba1Yfg2u (or U?)

This will annotate the calendar year (2020, 2021 …)
Use automatic ticks
Put gridlines every 2 weeks.

If you just want to “zoom in”
It’s -R that you need to change (-R2000-01-01T/2000-02-01T would plot only January 2000)

19

Screenshot (751)
Screenshot (751)1920Ă—1080 268 KB
Uploading: Screenshot (751).png…

this is the result.
it very terrible…I had already done this option, but I was not enthusiastic about it because it shows almost everything grey.

isn’t there a way to do what I wrote but not make it like this? But like in the screenshot?

20

The screenshot uses gridlines every years. There’s nothing apparent about a 2 weeks period.

It’d be simply -Ba1Yf1Yg1Y, and you can change the line style with --MAP_GRID_PEN_PRIMARY=thin,lightgray,dotted for example