WANtaroHP (GMT: Sample)

Outline of this page

Following items are shown in this page.

Map of Japan

Script and figure by GMT

Map of world

Script and figure by GMT

World map

Great circle course

Schmidt net (Lambert azimuthal equal-area projection)

Script and figure by GMT

Seismicity map of East Asia

      Link to USGS Home Page 


  NEIC: Earthquake Search Results


                    U.  S.  G E O L O G I C A L  S U R V E Y

                     E A R T H Q U A K E  D A T A  B A S E



 FILE CREATED:  Sun Oct 30 05:53:38 2011
 Global Search   Earthquakes=    383175
 Catalog Used: PDE
 Date Range:   1900/01/01   to    2011/10/31
 Magnitude Range:   3.0  -  10.0
 Depth Range:    0  -  100
 Data Selection: Historical & Preliminary Data


 CAT    YEAR  MO DA  ORIG TIME   LAT    LONG  DEP  MAGNITUDE  IEM DTSVNWG DIST
                                                              NFO          km
                                                              TF

 PDE    1973  01 01 034609.80  -9.21  150.63  41  5.3 mbGS    ... .......
 PDE    1973  01 01 052229.80 -15.01 -173.96  33  5.0 mbGS    ... .......
 PDE    1973  01 01 114237.50 -35.51  -16.21  33  6.0 MsGS    ... .......
.......... .......... ..........
 PDE-Q  2011  10 28 084809.40  38.49   43.77  10  4.6 mbGS    ... .......
 PDE-Q  2011  10 28 090837.95 -40.63  126.42  10  4.7 mbGS    ... .......
 PDE-Q  2011  10 28 091846.27  35.53  -97.37   5  3.5 LgGS    ... .......

------------------------------------------------------------------------

*USGS National Earthquake Information Center*
USGS Privacy Statement  | Disclaimer
 	FirstGov 

Samples of explanatory drawing

Script and figure by GMT (Embedded pipe model)

Script and figure by GMT (Shell model)

Comparison of climate (temperature and rainfall)

Outline

• The Temperatures and rainfalls of 4-cities in Japan are shown in the graphs.
• 4 graphs are included in 1 page of image. This image was created using -X and -Y options of GMT.
• In this case, two bachfiles were made. One is control part for drawing, and the other is main part for drawing. In control part, to call main part is repeated, and 4 graphs are created.
• Input datas for 4-cities are quite same as them which are used in gnuplot drawing. So, data files are divided and arranged in the bachfile for control using awk script because of easy data reading by GMT.

-P option of GMT and rotation of image

• Since -P option of GMt is not used in this sample, A4 landscape image file is outputted.
• In case of no-use of -P option of GMT, it is necessary for 'eps' image to be rotated by 90 degree for creating 'png' image.
• Therefore the 'png' image shown below was created by rotating 'eps' image using ImageMagick command.

Script, data and figure by GMT

Time-sequence diagram (monthly variation of temperature and rainfall)

Outline

• This sample shows the variation of monthly average temperature and cumulative rainfall in each month.
• Temperature variation is shown by red solid line, and rainfall variation is shown by cyan bar chart with magenta outline.
• Although color name 'cyan' is valid, color name 'aqua' is invalid on GMT.

Script, data and figure by GMT

(Note) Input data for legend drawing is created in the batch file for GMT execution.

Time-sequence diagram (daily variation of river discharge and duration curve)

Outline

• This sample shows the variation of daily average discharge of the river.
• In addition, the discharge-duration curve and simple table are added in this graph.
• Data of table is obtained from sorted daily discharge in descending order.

Script, data and figure by GMT

(Note) Input data for legend drawing is created in the batch file for GMT execution.

Normal axes graph (Time history of earthquake acceleration)

Outline

• This sample shows earthquake accelerations in 3-directions (EW, NS and UD).
• The range of the vertical axis of each graph is determined as the same value as the maximum or minimum value of acceleration. In order to realize it, awk script is used to search maximum and minimum value of acceleration and to make batch file for passing the value to GMT. As far as I know, it is impossible to do it for gnuplot.

Input data format

2011/03/11 14:46:00 TCGH16 EW2 Max.acc=1196.697  | Comment
dt,0.010                                         | Sampling period (sec)
ndata,30000                                      | Number of data
-1.890                                           | Acceleration (gal)
-1.888                                           | ...
-1.888                                           | ...
-1.891                                           |
.....                                            |
.....                                            |

Script, data and figure by GMT

Log-Log graph (Acceleration response spectrum)

Outline

• This sample shows the superimposed graph of acceleration response spectrum of earthquake motions in 3 directions (EW, NS and UD).

Input data format

2011/03/11 14:46:00 TCGH16 EW2 Max.acc=1196.697             | Comment
nt,200                                                      | Number of data
period(sec)  , acc(gal)     , vel(kine)    , dis(cm)        |
2.000000e-002, 1.196764e+003, 2.924900e-001, 1.213241e-002  | Period and
2.063122e-002, 1.198466e+003, 3.113090e-001, 1.292889e-002  |    Response
2.128236e-002, 1.200244e+003, 3.314467e-001, 1.377853e-002  | .....
.....                                                       |
.....                                                       |

Script, data and figure by GMT

(Note) Input data for legend drawing is created in the batch file for GMT execution.

Log-Log graph (Tripartite response spectrum)

Outline

• This sample shows the graph of tripartite response spectrum of earthquake motions in 3 directions (EW, NS and UD).
• In this graph, the velocity response spectrum is drawn on the Log-Log plane firstly, next axes for acceleration and displacement responses are appended.

Input data format

2011/03/11 14:46:00 TCGH16 EW2 Max.acc=1196.697             | Comment
nt,200                                                      | Number of data
period(sec)  , acc(gal)     , vel(kine)    , dis(cm)        |
2.000000e-002, 1.196764e+003, 2.924900e-001, 1.213241e-002  | Period and
2.063122e-002, 1.198466e+003, 3.113090e-001, 1.292889e-002  |    Response
2.128236e-002, 1.200244e+003, 3.314467e-001, 1.377853e-002  | .....
.....                                                       |
.....                                                       |

Script, data and figure by GMT

(Note) Input data for legend drawing is created in the batch file for GMT execution.

Semi-log graph (Fourier spectrum ratio)

Outline

• This sample shows the superimposed graph of 2 Fourier spectra and ratio of them.
• The peak values of spectrum curve are searched using awk and they are shown in the graph by GMT.
• In this case, 2 vertical axes are used. One is for spectrum ratio as W-asix and another is for Fourier spectrum as E-axis.
• In output image 1, three lines overlap and it is hard to see.
• In output image 2, the range of E-axis is changed. And expressions of vertical axis scale are changed into exponent notation appending p after the command line option -B.

Input data format

i, frq(Hz)                 , sp1               , sp2               , spr
1, 0.0000000000000000      , 68.644469660094956, 93.489593441899203, 1.3619391905106006
2, 3.05175781250000000E-003, 68.646266058682215, 93.493209037733777, 1.3619562199903368
3, 6.10351562500000000E-003, 68.651654422441766, 93.504054788801866, 1.3620073044916456
.....
.....

Script, data and figure by GMT

(Note) Input data for legend drawing is created in the batch file for GMT execution.

Histogram (Maximum daily rainfall in a year, cumulative rainfall in a year)

Outline

• This sample shows histogram using the function pshistogram of GMT.
• In this case, since the ranges of axes are not specified (no '-R'), GMT defines the ranges of axes automatically.
• In addition, number of data, minimum and maximum values are shown in the graph. These values are obtained using the function minmax of GMT.

Script, data and figure by GMT