├── hw1 ├── output │ ├── count.gif │ ├── count.mp4 │ ├── test1.pdf │ └── test1.png ├── test_1.sh ├── test_2.sh └── README.md ├── beachballs ├── README.md ├── meca.dat ├── beachball-same-size.sh ├── beachball-custom-colors.sh ├── beachball-color-by-depth.sh ├── beachball.sh ├── beachball-cross-section.sh ├── japan-beachballs.sh └── japan-focal.dat ├── grids ├── README.md ├── contours.sh ├── images.sh ├── contour-antarctica.sh ├── images-bonus-shading.sh └── images-antarctica.sh ├── basic ├── world.sh ├── twoworlds.sh └── italy.sh ├── .github └── ISSUE_TEMPLATE │ ├── homework.md │ └── install.md ├── animation ├── earth.sh ├── quakes.sh └── get_quakes.sh ├── README.md └── INSTALL.md /hw1/output/count.gif: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/GenericMappingTools/2019-unavco-course/HEAD/hw1/output/count.gif -------------------------------------------------------------------------------- /hw1/output/count.mp4: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/GenericMappingTools/2019-unavco-course/HEAD/hw1/output/count.mp4 -------------------------------------------------------------------------------- /hw1/output/test1.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/GenericMappingTools/2019-unavco-course/HEAD/hw1/output/test1.pdf -------------------------------------------------------------------------------- /hw1/output/test1.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/GenericMappingTools/2019-unavco-course/HEAD/hw1/output/test1.png -------------------------------------------------------------------------------- /beachballs/README.md: -------------------------------------------------------------------------------- 1 | # Plotting beachballs 2 | 3 | These scripts are examples of plotting beachballs on maps and cross-sections. 4 | -------------------------------------------------------------------------------- /grids/README.md: -------------------------------------------------------------------------------- 1 | # Plotting grids 2 | 3 | These scripts are examples of plotting grids with contours and pseudo-color images. 4 | -------------------------------------------------------------------------------- /hw1/test_1.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # Test script 1: Make sure your GMT installation works 3 | gmt begin test1 4 | gmt grdimage @earth_relief_20m -R0/40/0/40 -JM15c -I+d -U 5 | gmt coast -W1p -B 6 | gmt end show 7 | -------------------------------------------------------------------------------- /basic/world.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # To make more than one plot per session you need to use gmt figure 3 | gmt begin 4 | gmt figure world1 5 | gmt coast -Rg -JH0/15c -Gpurple -Bg 6 | gmt figure world2 7 | gmt coast -Rg -JR0/15c -Ggreen -Bg 8 | gmt end show 9 | -------------------------------------------------------------------------------- /beachballs/meca.dat: -------------------------------------------------------------------------------- 1 | # lon lat depth mrr mtt mpp mrt mrp mtp iexp name 2 | 114 32.9 24 0.06 -2.48 2.42 -0.34 0.52 -5.10 25 X Y strike-slip fault 3 | 115 33.0 5 -1.55 0.13 1.42 0.24 -0.29 -0.38 24 X Y normal fault 4 | 116 33.1 13 2.41 -0.15 -2.26 1.03 -1.81 0.48 22 X Y reverse fault 5 | -------------------------------------------------------------------------------- /basic/twoworlds.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # Show how subplot works. Here a 2 by 1. Change to 3 and add one more coast -c command! 3 | gmt begin 4 | # See subplot man page for the mysterious -F option 5 | gmt subplot begin 2x1 -Ff16c/25c -M0 -A 6 | gmt coast -Rg -JH0 -Gpurple -Bg 7 | gmt coast -Rg -JR0 -Ggreen -Bg -c 8 | gmt subplot end 9 | gmt end show 10 | -------------------------------------------------------------------------------- /hw1/test_2.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # Test script 2: Test that ffmpeg is present and able to make a MP4 or WEBM movie 3 | cat << EOF > main.sh 4 | gmt begin 5 | echo \${MOVIE_FRAME} | gmt text -R0/1/0/1 -JX10c -F+f200p+cCM -B0 -X0 -Y0 6 | gmt end 7 | EOF 8 | gmt movie main.sh -Ncount -T25 -C10cx10cx30 -Fmp4 -Fwebm -Gpink -D4 -Z 9 | rm -f main.sh 10 | echo "Open either count.mp4 or count.webm in an application that can view movies" 11 | -------------------------------------------------------------------------------- /beachballs/beachball-same-size.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # 3 | # Plot beachballs of three basic faults (strike-slip, normal and reverse faults) 4 | # for the same size 5 | # 6 | 7 | gmt begin beachball-same-size 8 | # By default, the beachball size is proportional to the magnitude, 9 | # with the size for magnitude=5 set by -S. 10 | # Use -M if you want to use the same size for any magnitude. 11 | gmt meca meca.dat -JM15c -R113/117/32/34 -Baf -Sd2c -M 12 | gmt end show 13 | -------------------------------------------------------------------------------- /beachballs/beachball-custom-colors.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # 3 | # Plot beachballs of three basic faults (strike-slip, normal and reverse faults), 4 | # with the compressional and extensive quadrants filled by different colors. 5 | # 6 | 7 | gmt begin beachball-custom-colors 8 | # -G select the filling of the compressional quadrants 9 | # -E select the filling of the extensive quadrants 10 | gmt meca meca.dat -JM15c -R113/117/32/34 -Baf -Sd2c -Gred -Elightblue 11 | gmt end show 12 | -------------------------------------------------------------------------------- /.github/ISSUE_TEMPLATE/homework.md: -------------------------------------------------------------------------------- 1 | --- 2 | name: Homework help 3 | about: Questions about homework assignments 4 | title: '' 5 | labels: homework 6 | assignees: '' 7 | 8 | --- 9 | 10 | 16 | **Operating system**: 17 | 18 | **Homework number**: 19 | 20 | **Question**: 21 | -------------------------------------------------------------------------------- /.github/ISSUE_TEMPLATE/install.md: -------------------------------------------------------------------------------- 1 | --- 2 | name: Install help 3 | about: Ask for help installing the required software 4 | title: '' 5 | labels: install 6 | assignees: '' 7 | 8 | --- 9 | 10 | 11 | 12 | **Operating system**: 13 | 14 | 15 | **The commands you ran or things you installed**: 16 | 17 | 18 | **Which step is giving you trouble**: 19 | 20 | 21 | **What is the exact error message (copy and paste it)**: 22 | -------------------------------------------------------------------------------- /animation/earth.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # Make a very simple spinning globe movie using coast 3 | # Run one of the gmt movie commands commented out below 4 | gmt begin 5 | gmt coast -Rg -JG${MOVIE_FRAME}/30/12c -Gburlywood -Bg -Stomato -Xc -Yc 6 | gmt end 7 | # 1. Test that script works by just making a single frame (55) as PDF, no movie generated 8 | # gmt movie earth.sh -C360p -T91 -M55,pdf -Fnone -Nearth 9 | # 2. Run the movie and build a MP4 animation 10 | # gmt movie earth.sh -C360p -T91 -Fmp4 -Nearth -Le -------------------------------------------------------------------------------- /basic/italy.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # Make our first GMT map of Italy, of all places 3 | gmt begin italy pdf,png # Starting our new gmt modern mode session, calling plot italy and ask for pdf and png 4 | #Lay down painted continent with national borders on a Mercator map 5 | gmt coast -R5/20/35/50 -Wthin -Gbeige -Df -Sazure -B -N1/thick,red -JM15c 6 | # Show where Italy is in the world via a map inset 7 | gmt inset begin -DjTR+w4c+o0.2c -M0 -F+gwhite+pthick 8 | gmt coast -Rg -JG70W/15S/4c -Gbeige -Bg -EIT+gred 9 | gmt inset end 10 | gmt end show # end will finish the plots and open them in a viewer 11 | -------------------------------------------------------------------------------- /grids/contours.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # 3 | # Create a contour plot of Earth relief at 60 arc-minute resolution 4 | 5 | gmt begin contours pdf 6 | # Use basemap to plot the map frame and set the plot region and projection. 7 | # -JN = Robinson projection 8 | gmt basemap -JN15c -B -Rg 9 | # Without any options will use the defaults. 10 | # Set the regular (C) and annotated (A) contours and the pen for regular (Wc) and 11 | # annotated (Wa). +f sets the font property of the annotations (in this case, the 12 | # size to 6 points). 13 | gmt grdcontour @earth_relief_60m \ 14 | -C1000 -A4000+f6p \ 15 | -Wathin,black,-- -Wcthinnest,black 16 | gmt end show 17 | -------------------------------------------------------------------------------- /beachballs/beachball-color-by-depth.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # 3 | # Plot beachballs of three basic faults, with the compressive quadrants color 4 | # determined by earthquake depth 5 | # 6 | 7 | gmt begin beachball-color-by-depth 8 | # Customize the jet CPT 9 | gmt makecpt -Cjet -T0/30/10 10 | # Plot beachballs using the "current" CPT, by using -Z with no arguments. 11 | gmt meca meca.dat -JM15c -R113/117/32/34 -Baf -BWSen -Sd2c -Z 12 | # Plot the colorbar in the Middle Right (MR), shift it 1.5 cm to the right. 13 | # Set the label for the colorbar with -Bx+l, 14 | # and move the label to the left of the colorbar (+ml) 15 | gmt colorbar -DJMR+ml+o1.5c/0c -Bx+l"Depth (km)" -Np 16 | gmt end show 17 | -------------------------------------------------------------------------------- /animation/quakes.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # Make sure you run get_quakes.sh first to obtain some data and create the CPT 3 | # and times.txt data file needed by movie. 4 | # This script (below) creates a single frame, laying down a coast map and then 5 | # plotting only the events that should be visible for the given time (-T). 6 | gmt begin 7 | gmt coast -Rg -JG${MOVIE_COL1}/30/12c -Gburlywood -Bg -Sgray -Xc -Yc 8 | gmt events quakes.txt -Cquakes.cpt -T${MOVIE_COL0} -SE- -Es+r2+d6 -Ms5+c0.5 -Mi1+c-0.6 -Mt+c0 --TIME_UNIT=d 9 | gmt end 10 | # 1. Test that script works by just making a single frame (50) as PDF, no movie generated 11 | # gmt movie quakes.sh -C360p -Ttimes.txt -M50,pdf -Fnone -Nquakes 12 | # 2. Run the movie and build a MP4 animation at 360p resolution 13 | # gmt movie quakes.sh -C360p -Ttimes.txt -Fmp4 -Nquakes 14 | -------------------------------------------------------------------------------- /grids/images.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # 3 | # Create a pseudo-color image plot of global Earth relief 4 | 5 | gmt begin images pdf 6 | # Use basemap to create a frame and set the region and projection (in this case, 7 | # Mollweide) 8 | gmt basemap -Rg -JW20c -B 9 | # Plot Earth relief by mapping values to colors using the "etopo1" CPT (colormap). 10 | # Set automatic shading effects 11 | gmt grdimage @earth_relief_30m -Cetopo1 -I+d 12 | # Add a colorbar anchored to the bottom center (BC) in the outside of the plot (J) 13 | # -C tells it to use the "current" colormap (this will not be necessary in the near 14 | # future). Use -B to configure the frame of the colorbar (just like the map frame). 15 | # The number is the interval between annotations. 16 | gmt colorbar -DJBC -C -B2000 17 | gmt end show 18 | -------------------------------------------------------------------------------- /beachballs/beachball.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # 3 | # Plot beachballs of three basic faults (strike-slip, normal and reverse faults) 4 | # 5 | 6 | gmt begin beachball 7 | # meca supports focal mechanisms in different conventions. 8 | # Here -Sd means the focal mechanisms in input file meca.dat are given 9 | # in Global CMT moment tensor format, and only the double couple part 10 | # is plotted. 11 | # The beachball size is proportional to the magnitude, with the size 12 | # for magnitude=5 is set to be 2 cm. 13 | # 14 | gmt meca meca.dat -JM15c -R113/117/32/34 -Baf -Sd2c 15 | 16 | # Adjust the attributions of beachball event titles: 17 | # - set the font size (12p) 18 | # - offset from beachball (6p) 19 | # - place the event title below the beachball (+u) 20 | # 21 | #gmt meca meca.dat -JM15c -R113/117/32/34 -Baf -Sd2c/12p/6p+u 22 | gmt end show 23 | -------------------------------------------------------------------------------- /beachballs/beachball-cross-section.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # 3 | # Plot cross-section of beachballs 4 | # 5 | 6 | gmt begin beachball-cross-section 7 | # The module name "coupe" comes from the French verb "to cut". The best translation is a (vertical) cross section. 8 | # coupe uses the same input file format and same -S options as meca. 9 | # The cross-section is selected by specifly the locations (longitude and latitude) of a starting point (113/33), 10 | # and a ending point (117/33). The cross-section plane is vertical (dip angle=90), with the width set to be 100 km, 11 | # and depth to be 0-30 km. 12 | # The plot region (usually set by -R) is automatically determined and set if +f is used. 13 | # 14 | # To reverse the Y axis, set the figure heigth to a negative value (-10 cm) 15 | gmt coupe meca.dat -JX15c/-10c -Bxaf+l"Distance (km)" -Byaf+l"Depth (km)" -BWSen \ 16 | -Sd2c -Aa113/33/117/33/90/100/0/30+f -Q 17 | gmt end show 18 | -------------------------------------------------------------------------------- /grids/contour-antarctica.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # 3 | # Create a contour map of Antarctica with different color and style for continents and 4 | # oceans 5 | 6 | gmt begin contour_antarctica pdf 7 | # Create the map frame and set the region and projection. 8 | # There was a bug with the Antactica ISO country code so we have to specify the 9 | # reigon manually. 10 | # The projection is Azimuthal orthographic centered at the South pole (0/-90) 11 | gmt basemap -R0/360/-90/-60 -JG0/-90/15c -B 12 | # Restrict to the positive contours by passing -Lp (P would include 0) 13 | gmt grdcontour @earth_relief_30m -Lp \ 14 | -C500 -Wcthinnest,black \ 15 | -A2000+f6p -Wathin,black 16 | # Restrict to the negative contours by passing -Ln (N would include 0) 17 | gmt grdcontour @earth_relief_30m -Ln \ 18 | -C500 -Wcthinnest,lightblue \ 19 | -A2000+f6p -Wathin,lightblue 20 | # Plot the coastlines using a thick black line 21 | gmt coast -Wthick,black 22 | gmt end show 23 | -------------------------------------------------------------------------------- /animation/get_quakes.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # This script is a helper script to obtain NEIC seismicity from the USGS 3 | # You can edit the dates, times, magnitudes, etc. 4 | SITE="https://earthquake.usgs.gov/fdsnws/event/1/query.csv" 5 | TIME="starttime=2018-01-01%2000:00:00&endtime=2018-12-31%2000:00:00" 6 | MAG="minmagnitude=5" 7 | ORDER="orderby=time-asc" 8 | URL="${SITE}?${TIME}&${MAG}&${ORDER}" 9 | # We get the data via the URL and reformat using gmt convert. Here we use the 10 | # -i option to shuffle the columns, plus scaling the magintue (col 4) by 50 to 11 | # yield a fake symbol size in km 12 | gmt convert ${URL} -i2,1,3,4+s50,0 -hi1 > quakes.txt 13 | # Make standard r,g,b seismicity depth color table 14 | gmt makecpt -Cred,green,blue -T0,70,300,10000 > quakes.cpt 15 | # Create a file with one record per day in 2018 in first column, with 16 | # longitudes spanning the Pacific (from 160 to 240 evenly) in the second column. 17 | # We will use this file to make the quakes movie in quakes.sh 18 | gmt math -T2018-01-01T/2018-12-31T/1 --TIME_UNIT=d TNORM 40 MUL 200 ADD = times.txt 19 | -------------------------------------------------------------------------------- /grids/images-bonus-shading.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # 3 | # BONUS: Create an image plot of the geoid but shade it using topography to give a 4 | # slight impression of it on the smooth geoid map. 5 | 6 | gmt begin geoid 7 | # Make a global Mollweide map 8 | gmt basemap -Rg -JW25c -B 9 | 10 | # We can't use the automatic -I+d since we want to use a different data grid to 11 | # create the shadow. For this, we have to go old school. 12 | # Shading is done by calculating a directional derivative of the data. That is what 13 | # grdgradient does. It will create a new grid file (-G) with the derivative along 14 | # the 45 degree azimuth (-A) and normalize it to [-1,1] (-N), which is needed for 15 | # shading. 16 | gmt grdgradient @earth_relief_10m -A45 -N -Grelief_gradient.nc 17 | 18 | # The geoid data is in the remote file @osu91a1f_16.nc. This is a 60 arc-minute 19 | # resolution grid. To use shading, the data and gradient grids have to match. We 20 | # will resample the geoid grid at 10 arc-minutes (-I) to match the gradient grid. 21 | # Resampling to smaller grid spacing is OK but doing the opposite may result in 22 | # aliasing! 23 | gmt grdsample @osu91a1f_16.nc -I10m -Ggeoid.nc 24 | 25 | # Finally, make a image plot the same way we did before but passing the gradient 26 | # grid to -I instead of the defaults (+d). The polar CPT is good for showing 27 | # diverging data (which has negative and positive). 28 | gmt grdimage geoid.nc -Cpolar -Irelief_gradient.nc 29 | 30 | # Make a colorbar with arguments we've seen before 31 | gmt colorbar -DJBC+w10c+o0/0.5c -B25 -Bx+l"geoid height [m]" 32 | gmt end show 33 | -------------------------------------------------------------------------------- /grids/images-antarctica.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # 3 | # Create a pseudo-color image plot of Antactica using different colormaps (CPTs) for 4 | # land and ocean. 5 | 6 | gmt begin images_antarctica pdf 7 | # Same as the contour example 8 | gmt basemap -R0/360/-90/-60 -JG0/-90/20c -B 9 | 10 | # Plot the oceans: 11 | 12 | # Customize the abyss CPT to use only negative data in the given range. We'll use 13 | # this for the oceans. 14 | gmt makecpt -Cabyss -T-5000/0 15 | # Another coast superpower: clip the next plotting commands to include only "wet 16 | # areas". Whatever is plot between this and "coast -Q" will have the continental 17 | # part cut out. Remember that -S is used to set the fill for oceans. 18 | gmt coast -Sc 19 | # Plot using the "current" CPT by using -C with no arguments. Will take the one 20 | # created by makecpt. Don't use shading for this part. 21 | gmt grdimage @earth_relief_05m -C 22 | gmt coast -Q 23 | # Plot the colorbar in the bottom center but shift it 6 cm to the left and 0.5 cm 24 | # down (+o) and set the width to 8 cm (+w). Set the label for the x-axis of the 25 | # colorbar using -Bx+l (l for label). This uses the current colorbar which is now 26 | # the modified abyss 27 | gmt colorbar -DJBC+o-6c/0.5c+w8c -B1500 -Bx+l"bathymetry [m]" -C 28 | 29 | # Plot the land: 30 | 31 | # Customize the grays CPT (which goes from black to white) so that the land 32 | # topography goes from gray to white by shifting the black part down to -6000 33 | gmt makecpt -Cgray -T-6000/6000 34 | # Clip the land portion of the plot. Remember that -G can set the fill for land. 35 | gmt coast -Gc 36 | # Plot using the "current" CPT, which is not the modified gray 37 | gmt grdimage @earth_relief_05m -C -I+d 38 | gmt coast -Q 39 | # Now the current colorbar is our modified gray, so this call will draw that instead 40 | # of the abyss. Shift it to the right (+o6c instead of +o-6c) instead. But since we 41 | # cheated when making the CPT by stretching it out to -6000, we can tell colorbar to 42 | # crop between 0 and 6000 (-G) 43 | gmt colorbar -DJBC+o6c/0.5c+w8c -B1500 -Bx+l"topography [m]" -C -G0/6000 44 | 45 | gmt coast -Wthick,black 46 | 47 | # Add a light contour to the land portion to make it look nicer. Use a lower 48 | # resolution grid here to make it smoother. 49 | gmt grdcontour @earth_relief_30m -Lp -C500 -Wcthinnest -A2000+f6p -Wathin 50 | gmt end show 51 | 52 | -------------------------------------------------------------------------------- /hw1/README.md: -------------------------------------------------------------------------------- 1 | # Homework 1 2 | 3 | **Goal**: Install GMT and check that everything is correctly setup in your computer. 4 | 5 | Follow these steps to setup your computer with GMT and other software needed for the 6 | workshop and check that everything is working. You'll run some GMT scripts and check 7 | that you get the correct output. 8 | 9 | **If you have any questions/problems**: 10 | 11 | * Open an *Issue* on this repository by 12 | [clicking here](https://github.com/GenericMappingTools/2019-unavco-course/issues/new/choose). 13 | You'll need a GitHub account to do this. 14 | * Alternatively, post to our [Gitter chatroom](https://gitter.im/GenericMappingTools/2019-unavco-course) 15 | 16 | ## Tasks 17 | 18 | 0. Fill out the following survey so that we can better prepare for the workshop: 19 | https://forms.gle/1e8kBsN22KSjhioN9 20 | 1. Follow the [install instructions](../INSTALL.md) to setup your computer. 21 | 2. Open a terminal (Mac: open the "Terminal" app; Windows: open "Git Bash"). 22 | The following steps should be done in the terminal. 23 | To run a command, type it out and press *Enter*. 24 | 3. Run this command to download the workshop material using [git](https://en.wikipedia.org/wiki/Git): 25 | 26 | ``` 27 | git clone https://github.com/GenericMappingTools/2019-unavco-course.git 28 | ``` 29 | 30 | This will create a folder called `2019-unavco-course` in your computer. 31 | 4. Run this command to navigate to the `hw1` folder in `2019-unavco-course` and figure 32 | out where it is placed on your computer (each line is a command): 33 | 34 | ``` 35 | cd 2019-unavco-course/hw1 36 | pwd 37 | ``` 38 | 39 | 5. Run the [`test_1.sh`](test_1.sh) script: 40 | 41 | ``` 42 | ./test_1.sh 43 | ``` 44 | 45 | A window should pop up with a colored relief map of northern Africa. It should look 46 | like this: 47 | 48 | ![`2019-unavco-course/hw1/output/test1.pdf`](output/test1.png) 49 | 6. Run the [`test_2.sh`](test_2.sh) script: 50 | 51 | ``` 52 | ./test_2.sh 53 | ``` 54 | 55 | This should produce 2 files called `count.mp4` and `count.webm` in the 56 | `2019-unavco-course/hw1` folder. To open the mp4 movie (replace with `.webm` 57 | if the mp4 doesn't work): 58 | 59 | * Mac: run `open count.mp4` 60 | * Linux: run `xdg-open count.mp4` 61 | * Windows: run `explorer count.mp4` 62 | 63 | The result should be an animation with numbers counting from 0 to 24 that looks like 64 | this: 65 | 66 | ![`2019-unavco-course/hw1/output/count.mp4`](output/count.gif) 67 | -------------------------------------------------------------------------------- /beachballs/japan-beachballs.sh: -------------------------------------------------------------------------------- 1 | #!/usr/bin/env bash 2 | # 3 | # Plot real-world beachballs on maps and cross-sections 4 | # 5 | # The focal mechanism information is obtained from Global CMT project 6 | # - region: 125E-150E, 30N-47N 7 | # - time range: 2010-01-01 to 2019-07-01 8 | # - magnitude: >=5.5 9 | 10 | gmt begin 11 | # Figure 1: Beachballs on maps 12 | gmt figure japan-beachballs-mapview 13 | 14 | # plot basemap, setting map width, region and frames 15 | gmt basemap -R122/147/30/48 -JM15c -BWSen -Baf 16 | # Plot Earth relief by mapping values to colors using the "globe" CPT 17 | # Set automatic shading effects 18 | gmt grdimage @earth_relief_02m -Cglobe -I+d 19 | # Add the colorbar to the Top Right (TR) and shift to the right by 2 cm. 20 | # The colorbar width is set to be 6 cm. 21 | gmt colorbar -DjTR+o-2c/0c+w6c+ml -Ba2000+l"Earth relief (m)" 22 | 23 | # Customize the jet CPT for earthquake depth from 0 to 700 24 | gmt makecpt -Cjet -T0/700/100 25 | # Plot beachballs using the "current" CPT, by using -Z with no arguments. 26 | # Use -i0-11 to skip the trailing "event title" in the 13rd column (column number starts at 0). 27 | gmt meca japan-focal.dat -Sd0.5c -Z -i0-11 28 | # Add the colorbar to the Bottom Right (BR) and shift to the right by 2 cm. 29 | # The colorbar is reversed by giving a negative length (+w-6c). 30 | gmt colorbar -DjBR+o-2c/0c+w-6c+ml -Ba100+l"Focal depth (km)" 31 | 32 | # Prepare data file for plotting and label the cross-section line 33 | echo 126 42 A > xsection.dat 34 | echo 146 40 B >> xsection.dat 35 | # plot the cross-section line, with a 2 point, red pen 36 | # the plot column only needs the longitude and latitude of the input file, and ignores the third column 37 | gmt plot xsection.dat -W2p,red 38 | # label the cross-section by "A" and "B", with a font size of 15 points 39 | gmt text xsection.dat -F+f15p 40 | 41 | 42 | # Figure 2: Beachball cross-sections 43 | gmt figure japan-beachballs-sideview 44 | # Customize the jet CPT for earthquake depth from 0 -700 45 | gmt makecpt -Cjet -T0/700/100 46 | # The cross-section is selected by specifly the locations (longitude and latitude) of a starting point (126/42), 47 | # and a ending point (146/40). The cross-section plane is vertical (dip angle=90), with the width set to be 500 km, 48 | # and depth to be 0-70 km. 49 | # 50 | # To reverse the Y axis, set the figure heigth to a negative value (-10 cm) 51 | # Use -i0-11 to skip the trailing "event title" in the 13rd column (column number starts at 0). 52 | gmt coupe japan-focal.dat -Sd0.5c -Aa126/42/146/40/90/500/0/700+f -Q -Z -i0-11 \ 53 | -JX15c/-10c -R0/2000/0/700 -BWSen -Bxaf+l"Distance (km)" -Byaf+l"Depth (km)" 54 | # Add the colorbar to the Bottom Right (BR) and shift to the right by 2 cm. 55 | # The colorbar is reversed by giving a negative length (+w-10c). 56 | gmt colorbar -DjBR+o-2c/0c+w-10c+ml -Np -Ba100+l"Focal depth (km)" 57 | 58 | # cleanup temporary files 59 | rm xsection.dat 60 | gmt end show 61 | -------------------------------------------------------------------------------- /README.md: -------------------------------------------------------------------------------- 1 | # The Generic Mapping Tools for Geodesy 2 | 3 | Material for the 2019 UNAVCO Short Course. 4 | 5 | **Instructors:** 6 | [Paul Wessel](http://www.soest.hawaii.edu/wessel/), 7 | [Joaquim Luis](http://w3.ualg.pt/~jluis/), 8 | [Leonardo Uieda](http://www.leouieda.com), 9 | [Dongdong Tian](https://msu.edu/~tiandong/), 10 | [Xiaohua Xu](https://www.researchgate.net/profile/Xiaohua_Xu7) 11 | 12 | **When**: 13 | July 22-23 (Monday-Tuesday), 2019. 9:00 am - 5:00 pm (beach time 5-6 pm) 14 | 15 | **Where**: 16 | Scripps Institution of Oceanography, La Jolla, CA 17 | 18 | **Bring**: 19 | Participants must **bring their own computers** (Mac, Windows, or Linux) with the 20 | required **software already installed** (see [Installing](#installing) below). 21 | 22 | This 2-day course will cover the use of Generic Mapping Tools (GMT) in geodesy 23 | using UNIX shell scripting. Lectures and exercises will be given to teach the 24 | basic conventions of using GMT, such as plotting grids, images, and vector data 25 | (points, lines, polygons). Labs will include both processing and mapping of 26 | various data sets relevant to geodesy. 27 | 28 | ## Installing 29 | 30 | Participants will need to install some software for the workshop. 31 | Please follow 32 | [these instructions](INSTALL.md) 33 | **before the workshop**. 34 | 35 | ## Need help? 36 | 37 | If you have any trouble installing or questions about the course/homeworks, 38 | there are a few ways to ask for help: 39 | 40 | * Open an *Issue* on this repository by 41 | [clicking here](https://github.com/GenericMappingTools/2019-unavco-course/issues/new/choose). 42 | You'll need a GitHub account to do this. 43 | * Post to our [Gitter chatroom](https://gitter.im/GenericMappingTools/2019-unavco-course) 44 | 45 | ## Video, Slides and Notes 46 | 47 | The workshop video recording can be found on Youtube: 48 | https://www.youtube.com/playlist?list=PLzmugeDoplFMwr1W7dgmwf26p83ie6ZEN 49 | 50 | The slides with notes for the course are available in this 51 | Google Drive folder: 52 | https://drive.google.com/drive/folders/1zElyVNOYKK_Uea9FWCTlvYig75KZlzjZ?usp=sharing 53 | 54 | We also have shared notes on this Google Docs: 55 | https://docs.google.com/document/d/1TzWVXUwHRVdQki1kKNr07CU9-0bUPRsLCQTdVyutlNg/edit?usp=sharing 56 | 57 | 58 | ## Schedule 59 | 60 | | **Day 1** | | | 61 | |:------------|:----------|:-------------| 62 | | 08:30-09:00 | *Breakfast* | 30 min | 63 | | 09:00-09:15 | Introductions | 15 min | 64 | | 09:15-09:30 | Overview | 15 min | 65 | | 09:30-10:00 | Install help | 30 min | 66 | | 10:00-10:15 | *Break* | 15 min | 67 | | 10:15-12:00 | GMT Basics | 105 min | 68 | | 12:00-01:30 | *Lunch* | 90 min | 69 | | 01:30-03:15 | Symbols and lines | 105 min | 70 | | 03:15-03:30 | *Break* | 15 min | 71 | | 03:30-04:45 | Grids and images | 105 min | 72 | | 04:45-05:00 | Wrap up | 15 min | 73 | | 06:00-07:30 | Pizza (Munch room) | 90 min | 74 | | **Day 2** | | | 75 | | 08:30-09:00 | *Breakfast* | 30 min | 76 | | 09:00-10:30 | Seismology | 90 min | 77 | | 10:30-10:45 | *Break* | 15 min | 78 | | 10:45-12:00 | Geodesy | 75 min | 79 | | 12:00-01:30 | *Lunch* | 90 min | 80 | | 01:30-03:00 | Animations | 90 min | 81 | | 03:00-03:15 | *Break* | 15 min | 82 | | 03:15-04:45 | Final projects | 90 min | 83 | | 04:45-05:00 | Wrap up | 15 min | 84 | 85 | ## License 86 | 87 | Creative Commons License
This content is licensed under a 88 | Creative Commons Attribution 4.0 International License. 89 | -------------------------------------------------------------------------------- /INSTALL.md: -------------------------------------------------------------------------------- 1 | # Installing the required software 2 | 3 | 4 | > 🚨 **IMPORTANT NOTE** 🚨 5 | > Installing GMT and its dependencies can be a difficult task even for experienced 6 | > users (including the developers themselves). So **don't be discouraged** if you 7 | > have any trouble or questions. [We're here to help you work through it!](#need-help) 8 | 9 | 10 | Please follow these instructions to install **GMT 6.0.0rc3** (rc = release candidate) 11 | and other software that will be used. Here is a summary of what you'll need: 12 | 13 | * GMT 6.0.0rc3 and its dependencies 14 | * [Ghostscript](https://ghostscript.com/): Used by GMT to convert postscript (GMT's 15 | native graphics format) into PNG, PDF, etc. 16 | * [GDAL](https://gdal.org/): Library used by GMT for reading and writing raster and 17 | vector geospatial data formats (pre-installed on Windows). 18 | * [ffmpeg](https://ffmpeg.org/): Video creating and conversion software used to make GMT 19 | animations. 20 | * A text editor of your choice (see [instructions below](#install-a-text-editor)). 21 | 22 | **We will use features of GMT that are being introduced in GMT 6.** Showing up to the 23 | workshop with a GMT 5 installation could result in delays for you or not being able to 24 | fully participate in the tasks we have planned for the short course! 25 | 26 | See specific instructions for your operating system below: 27 | 28 | * [Linux](#linux) 29 | * [Mac](#mac) 30 | * [Windows](#windows) 31 | 32 | 33 | ## Need help? 34 | 35 | If have any trouble installing, there are a few you ways to ask for help: 36 | 37 | * Open an *Issue* on this repository by 38 | [clicking here](https://github.com/GenericMappingTools/2019-unavco-course/issues/new/choose). 39 | You'll need a GitHub account to do this. 40 | * Post to our [Gitter chatroom](https://gitter.im/GenericMappingTools/2019-unavco-course) 41 | 42 | Key things you need to tell us so that we can help you: 43 | 44 | 1. Your operating system 45 | 2. The commands you ran or things you installed 46 | 3. Which step is giving you trouble 47 | 4. What is the exact error message (copy and paste it) 48 | 49 | 50 | ## Linux 51 | 52 | There are 2 options for installing GMT 6.0.0rc3: 53 | 54 | 1. Use the [conda package manager](https://conda.io/en/latest/) to install 55 | binaries for GMT and all dependencies (**recommended**) 56 | 2. Build GMT from source and install dependencies from your system package manager (see 57 | [instructions below](#building-gmt-from-source)) 58 | 59 | Unfortunately, you cannot use your system package manager to install GMT because this 60 | course requires GMT 6, which is not available on most distributions. 61 | 62 | Most Linux distributions come with a text editor installed. You may also 63 | [install Visual Studio Code](#install-a-text-editor) if you want but it is not required. 64 | 65 | ### Installing with conda 66 | 67 | 1. Download and install the [Python **3.7** **64-bit** version of Miniconda](https://conda.io/en/latest/miniconda.html). 68 | This will give you access to the conda package manager. **Make sure you select to 69 | have conda added to your `PATH` when asked by the installer**. If you have 70 | the Anaconda Python distribution installed, you won't need to do this step. 71 | 2. Install GMT 6.0.0rc3 and its dependencies (including ghostscript and gdal) by running 72 | the following in a terminal: 73 | 74 | ``` 75 | conda install gmt=6.0.0rc3 -c conda-forge/label/dev -c conda-forge 76 | ``` 77 | 78 | 3. Install ffmpeg using conda by running the following in a terminal: 79 | 80 | ``` 81 | conda install ffmpeg -c conda-forge 82 | ``` 83 | 84 | 85 | ## Mac 86 | 87 | There are 3 options for installing GMT 6.0.0rc3 on Mac: 88 | 89 | 1. Use the GMT bundle installer and install dependencies separately using Homebrew 90 | (**recommended**) 91 | 2. Use the [conda package manager](https://conda.io/en/latest/) to install 92 | binaries for GMT and all dependencies 93 | (follow the same [instructions for conda on Linux](#installing-with-conda)) 94 | 3. Build GMT from source and install dependencies separately (see 95 | [instructions below](#building-gmt-from-source)) 96 | 97 | Make sure you also [have a text editor installed](#install-a-text-editor). 98 | 99 | ### Using the bundle installer and Homebrew 100 | 101 | 1. Download and install [Homebrew](https://brew.sh/) (a package manager for Mac that we 102 | can use to install our dependencies) 103 | 2. Download and install the bundle `gmt-6.0.0rc3-darwin-x86_64.dmg` from the 104 | [GMT 6.0.0rc3 release page](https://github.com/GenericMappingTools/gmt/releases/tag/6.0.0rc3) 105 | 3. Install Ghostscript, GDAL, and ffmpeg by running the following in a terminal (open 106 | the "Terminal" app, type this in, and hit *Enter*): 107 | 108 | ``` 109 | brew update 110 | brew install gdal ghostscript ffmpeg 111 | ``` 112 | 113 | 114 | ## Windows 115 | 116 | Before installing GMT, you will need access to a UNIX style terminal with a 117 | [bash shell](https://en.wikipedia.org/wiki/Bash_%28Unix_shell%29) for GMT scripting. 118 | **Don't worry if you don't know what this means. We'll go through it in the first day of 119 | the workshop.** 120 | Linux and Mac usually come with one installed but Windows users aren't so lucky. 121 | To install a terminal with a bash shell on Windows: 122 | 123 | 1. Download and install Git-for-Windows `Git-2.22.0-64-bit.exe` from 124 | [their official website](https://gitforwindows.org/). 125 | There are 126 | [video instructions for installing](https://www.youtube.com/watch?v=339AEqk9c-8) 127 | provided by [Software Carpentry](https://software-carpentry.org/) 128 | 2. In this course, when you see/hear "Terminal" you should open the "Git Bash" program. 129 | 130 | To install GMT, we'll use the Windows installer: 131 | 132 | 1. Download and install Ghostscript `gs926aw64.exe` 133 | [from the official website](https://github.com/ArtifexSoftware/ghostpdl-downloads/releases/tag/gs926) 134 | 2. Download and install the GMT Windows installer `gmt-6.0.0_2019.07.11-win64.exe` 135 | from [this link](http://w3.ualg.pt/~jluis/mirone/gmt-6.0.0_2019.07.11-win64.exe). 136 | It comes with GDAL and ffmpeg pre installed. 137 | (if downloading from the above is too slow you can try this 138 | [DropBox link](https://www.dropbox.com/s/ncgs0d9h5uqtep5/gmt-6.0.0_2019.07.11-win64.exe?dl=0) ) 139 | 140 | Make sure you also [have a text editor installed](#install-a-text-editor). 141 | 142 | 143 | ## Building GMT from source 144 | 145 | Follow the ["Building from source" instructions](https://github.com/GenericMappingTools/gmt/blob/master/BUILDING.md) 146 | from the GMT repository. 147 | 148 | 149 | ## Install a text editor 150 | 151 | You'll need a [text editor](https://en.wikipedia.org/wiki/Text_editor) to write GMT 152 | scripts. 153 | Any editor is fine and may use your preferred one. 154 | If you don't have one, we recommend installing the "Stable" version of [Visual Studio 155 | Code](https://code.visualstudio.com/) on Linux, Mac, or Windows. 156 | Installers and instructions are available on their website. 157 | -------------------------------------------------------------------------------- /beachballs/japan-focal.dat: -------------------------------------------------------------------------------- 1 | 143.60 38.00 19 -0.76 -1.34 2.10 0.28 2.17 1.45 24 X Y 201001170604A 2 | 131.08 30.91 39 0.20 -1.18 0.98 0.50 0.93 1.28 24 X Y 201001250715A 3 | 130.66 42.48 579 0.27 0.19 -0.46 0.23 -2.50 -0.59 26 X Y 201002180113A 4 | 147.78 44.08 52 1.19 -3.40 2.20 1.90 -1.97 -0.02 24 X Y 201003061331A 5 | 141.45 37.62 75 -2.11 0.25 1.86 0.52 0.67 0.79 24 X Y 201003131246A 6 | 141.98 37.70 46 5.55 -0.28 -5.27 1.83 4.97 -2.07 25 X Y 201003140808A 7 | 138.45 43.28 12 5.07 -0.51 -4.57 -0.23 -1.15 -1.67 24 X Y 201003300102A 8 | 141.76 37.42 33 0.88 -0.98 0.09 0.02 -0.30 -0.16 25 X Y 201006130332A 9 | 148.94 44.34 31 1.39 -0.37 -1.02 0.46 1.45 -1.03 25 X Y 201006180223A 10 | 141.84 30.61 18 3.86 0.40 -4.26 0.19 3.76 0.58 24 X Y 201006281207A 11 | 142.77 39.68 35 2.15 0.23 -2.37 0.47 2.35 -0.23 25 X Y 201007042155A 12 | 142.32 41.47 61 5.50 -0.10 -5.40 1.76 3.23 -2.03 24 X Y 201009130547A 13 | 140.01 37.31 12 1.87 0.03 -1.89 -0.32 0.01 0.37 24 X Y 201009290759A 14 | 143.15 42.14 62 2.06 -0.72 -1.34 0.30 1.35 -1.44 24 X Y 201010141358A 15 | 146.87 46.82 362 0.71 0.18 -0.89 0.35 -1.57 -0.81 24 X Y 201011021810A 16 | 143.26 40.89 24 2.79 0.03 -2.82 1.84 3.46 -0.63 24 X Y 201012060730A 17 | 142.29 31.49 12 -2.88 -0.21 3.09 0.25 0.81 -0.56 24 X Y 201101050057A 18 | 142.78 38.56 14 0.49 -0.02 -0.47 0.42 1.05 -0.13 27 X Y 201103090245A 19 | 142.80 38.33 29 0.95 -0.05 -0.90 0.51 1.15 -0.29 25 X Y 201103091816A 20 | 143.50 38.47 18 4.61 -0.24 -4.36 3.09 6.49 -1.45 24 X Y 201103091844A 21 | 142.91 38.29 23 0.70 -0.02 -0.68 0.27 0.79 -0.23 25 X Y 201103092122A 22 | 142.82 38.27 23 3.79 -0.28 -3.51 1.90 4.42 -1.04 25 X Y 201103092124B 23 | 143.61 38.53 17 2.08 -0.12 -1.96 1.50 2.83 -0.64 24 X Y 201103100808A 24 | 143.05 37.52 20 1.73 -0.28 -1.45 2.12 4.55 -0.66 29 X Y 201103110546A 25 | 141.38 35.92 29 4.39 -0.31 -4.07 3.17 6.34 -1.97 27 X Y 201103110615A 26 | 144.63 38.27 21 -3.03 0.03 3.00 -0.41 -0.34 0.47 27 X Y 201103110625A 27 | 138.60 37.06 12 2.35 -1.43 -0.92 -2.03 -0.51 -1.96 25 X Y 201103111859A 28 | 143.01 39.43 30 -2.84 2.25 0.58 0.70 1.07 0.45 25 X Y 201103111902A 29 | 139.09 40.37 14 -0.00 1.92 -1.92 0.42 -0.47 1.48 25 X Y 201103111946A 30 | 142.68 38.98 12 -1.93 1.08 0.85 0.41 0.03 -0.73 25 X Y 201103112011A 31 | 142.85 38.81 19 -6.09 2.84 3.25 0.43 1.23 -1.45 24 X Y 201103120134A 32 | 142.84 37.57 14 -4.95 1.21 3.74 -3.62 -2.70 2.37 25 X Y 201103120147A 33 | 143.71 37.63 15 -0.16 -0.57 0.73 0.10 0.30 0.33 25 X Y 201103120247A 34 | 141.54 36.03 35 -1.09 0.66 0.42 -0.07 -0.03 0.00 25 X Y 201103120311A 35 | 143.52 40.07 29 3.95 -0.88 -3.06 1.74 3.37 -0.31 24 X Y 201103120452A 36 | 142.35 39.30 37 -1.82 1.10 0.72 0.50 1.01 0.55 24 X Y 201103120618A 37 | 142.28 39.04 36 -3.10 1.27 1.83 -0.42 3.23 -1.15 24 X Y 201103121053A 38 | 143.60 37.71 21 -0.60 -0.43 1.02 0.19 0.41 0.29 25 X Y 201103121253A 39 | 141.53 37.23 46 1.67 -0.34 -1.34 0.93 1.22 -0.65 25 X Y 201103121315A 40 | 142.72 38.76 26 -2.30 0.18 2.12 -0.80 -1.31 -0.69 24 X Y 201103121403A 41 | 142.52 39.52 15 0.88 0.14 -1.02 0.15 5.80 0.14 24 X Y 201103121443A 42 | 142.94 36.47 12 2.24 -0.46 -1.77 0.88 1.67 -0.95 24 X Y 201103121719A 43 | 141.99 37.56 15 -0.86 0.20 0.67 -0.03 0.20 -0.26 25 X Y 201103122212A 44 | 141.94 38.07 23 3.05 -0.71 -2.34 0.69 4.59 -1.27 24 X Y 201103122324A 45 | 141.81 35.74 12 -1.59 0.29 1.30 0.89 -1.37 -0.52 25 X Y 201103130126A 46 | 142.24 36.35 12 5.87 -3.75 -2.12 -1.93 0.61 -1.06 24 X Y 201103130223A 47 | 143.55 39.60 22 2.20 -0.09 -2.11 1.65 3.01 -0.55 24 X Y 201103130756A 48 | 142.20 38.90 50 2.36 0.07 -2.43 0.87 1.68 -0.51 24 X Y 201103130952A 49 | 142.34 37.42 15 -4.88 1.40 3.48 4.55 3.72 -1.85 24 X Y 201103131137A 50 | 141.30 35.28 15 0.66 -0.11 -0.56 2.78 2.14 -0.80 24 X Y 201103131755A 51 | 141.08 36.43 12 -1.55 -2.98 4.53 1.22 0.30 1.24 24 X Y 201103140102A 52 | 142.62 37.77 20 -0.54 -0.34 0.88 -0.69 -0.81 0.13 25 X Y 201103140612A 53 | 142.42 37.23 16 2.19 -0.64 -1.54 2.38 2.31 -1.15 24 X Y 201103141759A 54 | 142.29 37.44 12 -5.90 0.98 4.92 1.95 2.84 1.93 24 X Y 201103150949A 55 | 142.13 37.67 28 0.20 0.35 -0.55 -1.26 0.96 0.19 25 X Y 201103151327A 56 | 138.58 35.27 19 0.38 -1.03 0.65 0.04 0.20 -0.66 25 X Y 201103151331A 57 | 143.50 40.36 23 0.77 -0.03 -0.73 0.44 0.95 -0.16 25 X Y 201103151523A 58 | 141.16 35.17 22 1.87 -2.86 0.99 1.79 1.16 0.13 24 X Y 201103152030A 59 | 140.83 35.78 12 -5.92 -0.87 6.79 0.10 2.54 -0.42 24 X Y 201103160352A 60 | 142.33 39.87 48 4.40 -0.18 -4.22 1.64 3.29 -1.16 24 X Y 201103160629A 61 | 142.55 40.17 41 1.29 0.07 -1.36 0.56 1.21 -0.29 25 X Y 201103170413A 62 | 143.46 37.79 12 -1.96 0.24 1.72 0.81 1.95 0.58 24 X Y 201103170612A 63 | 140.89 35.58 29 2.39 -1.78 -0.61 0.59 1.03 -0.94 24 X Y 201103171232A 64 | 141.42 36.80 39 -6.67 3.07 3.61 -0.82 -0.07 2.22 24 X Y 201103171254A 65 | 143.47 37.86 12 -1.49 0.06 1.43 0.29 1.22 0.34 24 X Y 201103180323A 66 | 142.36 39.24 25 0.09 0.09 -0.18 1.00 2.86 -0.15 24 X Y 201103182333A 67 | 143.29 39.67 22 4.65 -0.06 -4.59 2.44 5.90 -0.92 24 X Y 201103190122A 68 | 140.48 36.86 12 -6.88 3.01 3.87 1.33 2.14 -2.71 24 X Y 201103190956A 69 | 141.67 37.73 49 2.16 -0.38 -1.78 0.63 1.62 -0.57 24 X Y 201103200555A 70 | 142.14 39.39 51 5.07 -0.25 -4.82 2.12 3.67 -1.02 24 X Y 201103201203A 71 | 140.99 35.18 18 1.29 -2.53 1.24 3.08 1.39 0.59 24 X Y 201103220338A 72 | 144.00 37.11 12 -4.82 1.86 2.96 0.14 1.61 2.57 25 X Y 201103220718A 73 | 141.78 37.38 30 -0.20 0.44 -0.24 -2.45 0.52 0.65 25 X Y 201103220919A 74 | 143.81 39.84 13 2.14 0.05 -2.20 1.12 2.86 -0.57 25 X Y 201103220944A 75 | 143.34 39.90 26 1.91 -0.25 -1.67 1.01 2.43 -0.31 24 X Y 201103221121A 76 | 143.29 36.83 12 -2.80 1.12 1.68 -0.04 -0.10 1.30 24 X Y 201103221201A 77 | 141.79 36.12 34 3.42 -1.27 -2.15 1.17 1.96 -1.12 24 X Y 201103221204A 78 | 141.66 35.77 19 4.57 0.18 -4.75 4.13 1.88 -1.63 24 X Y 201103221350A 79 | 141.56 35.88 15 3.37 -3.18 -0.19 -0.64 -1.13 0.40 24 X Y 201103221503A 80 | 140.76 37.09 12 -3.65 0.14 3.51 0.71 1.72 1.14 24 X Y 201103222212A 81 | 142.34 39.08 37 0.58 0.07 -0.64 0.31 0.69 -0.06 25 X Y 201103240821A 82 | 142.17 38.78 46 1.88 -0.04 -1.83 0.82 1.60 -0.49 25 X Y 201103251136A 83 | 142.09 38.40 22 -2.28 1.88 0.39 1.85 0.32 -0.64 25 X Y 201103272223A 84 | 142.37 37.47 18 -1.64 0.13 1.51 1.55 0.90 0.77 25 X Y 201103291054A 85 | 142.63 36.09 12 5.51 -2.11 -3.40 -0.07 0.86 -2.08 24 X Y 201103300529A 86 | 142.05 38.97 50 1.07 -0.00 -1.07 0.42 0.84 -0.25 25 X Y 201103310715A 87 | 142.17 39.38 48 0.76 -0.03 -0.73 0.33 0.57 -0.16 25 X Y 201104011157A 88 | 143.54 40.27 22 1.95 -0.10 -1.85 0.88 2.09 -0.41 24 X Y 201104011703A 89 | 141.85 38.32 53 5.99 -1.43 -4.55 -1.11 -1.25 -1.65 26 X Y 201104071432A 90 | 131.84 30.04 12 -1.04 -0.05 1.09 0.25 0.35 0.18 25 X Y 201104091257A 91 | 140.54 37.00 12 -1.27 0.76 0.51 -0.01 0.44 -0.33 26 X Y 201104110816A 92 | 140.53 37.01 12 -1.86 1.95 -0.09 0.45 0.96 0.01 24 X Y 201104111142A 93 | 144.50 37.70 15 -2.03 1.37 0.66 0.11 0.47 0.81 24 X Y 201104111753A 94 | 140.85 35.42 16 -0.27 -1.88 2.15 -1.71 0.16 -0.48 25 X Y 201104112308A 95 | 140.63 37.08 18 0.24 0.22 -0.46 0.55 0.11 -0.73 25 X Y 201104120507A 96 | 142.13 39.39 52 2.45 -0.18 -2.26 1.04 1.88 -0.64 24 X Y 201104121937A 97 | 143.59 39.63 15 0.81 0.04 -0.85 0.42 1.12 -0.18 25 X Y 201104131957A 98 | 143.48 39.64 24 1.40 -0.17 -1.23 0.72 1.45 -0.31 24 X Y 201104132032A 99 | 142.25 35.50 12 2.36 -0.34 -2.02 0.91 2.38 -0.79 24 X Y 201104140608A 100 | 140.00 36.40 81 3.34 0.44 -3.77 3.17 2.95 -1.13 24 X Y 201104160219A 101 | 143.92 36.95 12 -1.75 1.41 0.34 0.42 0.91 0.81 24 X Y 201104162036A 102 | 143.85 40.30 15 1.60 -0.09 -1.51 1.24 3.23 -0.53 24 X Y 201104210039A 103 | 143.87 40.33 15 1.47 0.17 -1.65 1.90 6.47 -0.08 24 X Y 201104210154A 104 | 140.78 35.56 42 1.81 0.09 -1.90 -0.17 1.42 -0.28 25 X Y 201104211337A 105 | 141.49 37.55 53 2.74 -0.30 -2.43 0.78 1.88 -0.85 24 X Y 201104211611A 106 | 142.91 39.21 30 1.05 -1.00 -0.05 -0.30 0.21 0.19 25 X Y 201104231012A 107 | 141.85 37.50 43 2.98 -0.38 -2.60 -0.05 -0.38 -0.73 24 X Y 201104280927A 108 | 144.17 38.24 14 -1.12 -0.39 1.51 0.11 0.37 0.84 25 X Y 201105051458A 109 | 142.53 40.29 42 3.85 0.09 -3.94 1.32 3.16 -0.84 24 X Y 201105072052A 110 | 143.67 37.85 13 -3.13 1.17 1.96 0.30 0.61 1.84 24 X Y 201105092015A 111 | 130.98 43.61 564 0.74 2.16 -2.91 1.53 -3.24 2.33 24 X Y 201105101526A 112 | 141.67 37.31 41 1.23 -0.14 -1.09 0.40 1.06 -0.40 25 X Y 201105132335A 113 | 141.22 35.74 37 4.38 -0.50 -3.88 1.02 3.63 -1.34 24 X Y 201105200046A 114 | 140.79 35.64 49 2.10 -0.06 -2.04 -0.15 1.32 -0.39 24 X Y 201105212206A 115 | 143.59 37.59 17 -0.86 -2.17 3.02 0.05 1.05 1.36 24 X Y 201105221446A 116 | 143.42 39.75 25 3.45 -0.04 -3.41 3.08 6.37 -0.90 24 X Y 201105240340A 117 | 143.96 37.35 21 -1.41 0.85 0.56 -0.37 0.41 0.99 25 X Y 201106030005A 118 | 143.58 37.73 20 -1.51 -4.76 6.27 2.16 2.81 1.33 24 X Y 201106141306A 119 | 141.76 37.68 16 -1.89 3.02 -1.13 -1.90 5.86 -1.17 24 X Y 201106181131A 120 | 142.89 40.13 32 2.28 -0.12 -2.16 1.35 2.75 -0.52 24 X Y 201106221429A 121 | 142.51 39.99 41 1.05 0.05 -1.09 0.12 1.12 -0.03 26 X Y 201106222150A 122 | 141.91 36.33 23 3.51 -0.49 -3.02 1.21 3.07 -1.12 24 X Y 201107061515A 123 | 141.11 37.18 55 2.31 -0.59 -1.72 -0.04 0.05 -0.68 24 X Y 201107071835A 124 | 143.33 37.98 22 0.49 -3.25 2.76 0.41 1.36 2.12 26 X Y 201107100057A 125 | 140.16 36.13 67 2.07 0.27 -2.33 -0.12 1.29 -0.29 24 X Y 201107151201A 126 | 142.10 38.96 47 2.82 -0.12 -2.70 1.20 2.50 -0.71 25 X Y 201107230434A 127 | 141.66 37.70 48 2.67 -0.41 -2.26 0.95 2.38 -0.88 25 X Y 201107241851A 128 | 140.97 35.18 15 0.84 -1.12 0.28 2.73 1.06 0.11 24 X Y 201107251154A 129 | 143.60 40.31 19 1.85 -0.03 -1.81 1.08 2.83 -0.39 24 X Y 201107280901A 130 | 141.16 37.00 53 3.86 0.09 -3.95 0.84 -1.41 -1.22 25 X Y 201107301853A 131 | 142.31 39.87 48 3.36 -0.15 -3.22 1.37 2.65 -0.83 24 X Y 201108011344A 132 | 138.39 34.65 19 0.94 -0.97 0.04 0.42 -0.15 0.13 25 X Y 201108011458A 133 | 141.25 37.02 50 0.66 -0.10 -0.56 0.37 0.43 -0.27 25 X Y 201108111822A 134 | 143.84 36.80 12 -2.18 1.24 0.94 0.22 0.58 1.06 25 X Y 201108171144A 135 | 141.71 37.69 48 3.19 -0.07 -3.12 -0.17 -1.19 -0.61 25 X Y 201108190536A 136 | 141.99 36.01 18 0.55 -0.10 -0.46 0.32 0.69 -0.22 25 X Y 201108221123A 137 | 141.54 36.22 28 1.75 -0.19 -1.55 0.73 2.08 -0.74 25 X Y 201109150800A 138 | 143.22 40.20 22 0.71 -0.00 -0.71 0.44 0.94 -0.13 26 X Y 201109161926A 139 | 143.82 40.25 32 3.11 -0.44 -2.67 1.20 2.61 -0.74 24 X Y 201109162011A 140 | 143.28 40.28 28 0.92 -0.05 -0.87 0.33 1.07 -0.16 25 X Y 201109162108A 141 | 143.55 40.11 28 3.78 -0.05 -3.73 1.84 4.54 -0.88 24 X Y 201109162136A 142 | 143.51 40.30 24 3.18 -0.39 -2.79 1.25 3.00 -0.64 24 X Y 201109162240A 143 | 142.93 40.33 33 2.89 -0.06 -2.83 1.89 2.93 -0.78 24 X Y 201109170734A 144 | 143.13 39.93 24 2.38 -0.16 -2.22 1.07 2.41 -0.54 24 X Y 201109180704A 145 | 141.48 37.56 53 3.25 -0.40 -2.85 0.84 2.38 -0.93 24 X Y 201110100245A 146 | 142.52 43.83 191 -0.76 0.93 -0.17 1.41 -0.25 1.40 25 X Y 201110210802A 147 | 141.65 37.32 41 1.35 -0.16 -1.19 0.48 1.27 -0.46 25 X Y 201111231924A 148 | 142.94 41.79 43 1.54 -0.36 -1.19 0.58 1.42 -0.68 25 X Y 201111241025A 149 | 144.60 46.93 390 0.69 2.34 -3.03 -2.90 -3.32 -2.19 24 X Y 201112091942A 150 | 141.89 31.73 12 3.20 0.20 -3.40 -0.33 1.26 -0.25 24 X Y 201112151512A 151 | 138.17 31.61 352 -0.37 -0.11 0.48 -0.26 -1.77 -0.52 26 X Y 201201010527A 152 | 141.25 37.06 16 -2.81 1.30 1.51 1.18 1.81 1.38 24 X Y 201201120320A 153 | 142.51 40.23 44 2.12 0.03 -2.15 0.91 1.74 -0.42 24 X Y 201201280022A 154 | 141.68 36.19 27 1.42 -0.32 -1.10 0.71 1.69 -0.75 24 X Y 201202140327A 155 | 141.73 36.17 26 4.56 -0.81 -3.75 2.48 5.22 -1.82 24 X Y 201202140621A 156 | 141.12 35.47 23 -1.95 0.56 1.39 2.77 0.62 1.15 24 X Y 201202291432A 157 | 147.87 45.09 125 0.22 0.97 -1.19 -1.16 -2.63 0.09 24 X Y 201203121232A 158 | 144.93 40.88 13 -3.50 3.11 0.38 -0.04 0.57 -0.13 26 X Y 201203140908A 159 | 144.90 40.91 19 -1.90 0.50 1.40 0.56 -0.16 0.52 25 X Y 201203141049A 160 | 140.89 35.83 23 -1.17 0.20 0.97 0.46 0.04 0.23 25 X Y 201203141205A 161 | 139.52 35.86 102 1.57 0.20 -1.77 -1.15 0.76 0.66 24 X Y 201203151920A 162 | 142.22 39.76 14 0.84 -0.66 -0.19 -0.66 -0.19 0.86 25 X Y 201203271100A 163 | 141.23 37.10 52 4.21 -0.69 -3.52 1.07 3.07 -1.56 24 X Y 201204011404A 164 | 141.66 37.48 14 -0.21 0.61 -0.40 0.69 0.51 1.72 24 X Y 201204121450A 165 | 141.35 36.98 12 -4.32 0.08 4.24 0.25 -0.15 1.62 24 X Y 201204131010A 166 | 140.80 35.68 45 1.92 0.08 -2.00 0.33 1.19 -0.37 24 X Y 201204242022A 167 | 140.71 35.62 48 5.09 0.26 -5.35 -0.51 4.20 -0.41 24 X Y 201204291028A 168 | 143.60 39.67 17 3.47 -0.01 -3.46 2.26 5.53 -0.61 24 X Y 201205191905A 169 | 143.58 39.61 17 2.06 0.03 -2.09 1.45 3.44 -0.34 25 X Y 201205200719A 170 | 143.41 39.72 15 0.68 -0.04 -0.64 0.84 2.06 -0.18 24 X Y 201205200820A 171 | 142.22 41.39 54 0.90 -0.03 -0.88 0.37 0.61 -0.23 25 X Y 201205231502A 172 | 141.21 34.96 15 0.34 -1.99 1.65 0.37 -0.48 -0.05 25 X Y 201206051931A 173 | 142.15 38.92 47 2.65 -0.06 -2.59 1.01 2.32 -0.58 25 X Y 201206172032A 174 | 142.86 42.39 50 4.85 -3.85 -1.00 -6.30 2.79 -0.46 24 X Y 201208251416A 175 | 143.30 39.86 21 0.73 0.02 -0.74 0.60 1.01 -0.20 25 X Y 201210012221A 176 | 141.91 38.35 51 2.39 -0.11 -2.28 1.18 1.82 -0.68 24 X Y 201210251032A 177 | 143.69 37.86 17 -1.58 -1.09 2.67 0.30 1.62 1.32 24 X Y 201211050430A 178 | 144.09 38.01 58 0.51 0.44 -0.95 0.08 -0.14 -0.16 27 X Y 201212070818A 179 | 143.83 37.77 20 -0.65 -0.12 0.76 0.04 0.12 0.31 27 X Y 201212070818B 180 | 143.90 37.94 25 -2.36 0.16 2.20 0.30 -0.08 1.04 25 X Y 201212070831A 181 | 142.24 38.79 44 1.63 -0.11 -1.53 0.71 1.33 -0.44 24 X Y 201212291459A 182 | 142.51 40.18 44 1.57 0.05 -1.62 0.73 1.26 -0.35 24 X Y 201301080751A 183 | 143.24 42.85 105 -0.33 0.67 -0.35 3.06 -0.16 0.21 26 X Y 201302021417A 184 | 141.00 35.22 25 0.60 -1.06 0.47 1.45 0.88 0.05 24 X Y 201302191227A 185 | 139.39 36.89 13 -0.55 3.26 -2.71 -0.50 -0.12 -5.55 24 X Y 201302250723A 186 | 143.23 39.56 23 0.60 -0.01 -0.59 0.39 0.78 -0.15 25 X Y 201304011853A 187 | 143.43 39.50 21 1.18 -0.04 -1.14 0.83 1.78 -0.28 24 X Y 201304020208A 188 | 131.02 42.77 572 0.66 -0.25 -0.41 0.10 -3.09 0.45 25 X Y 201304051300A 189 | 131.05 42.82 575 2.38 -0.86 -1.52 0.29 -5.18 1.01 24 X Y 201304060029A 190 | 134.69 34.43 18 0.57 0.02 -0.59 -0.04 -0.58 0.12 25 X Y 201304122033A 191 | 139.37 34.00 12 -0.34 -1.59 1.94 -0.44 -0.44 -4.66 24 X Y 201304170857A 192 | 141.70 38.51 56 0.77 -0.07 -0.70 0.20 0.59 -0.18 25 X Y 201304171203A 193 | 138.88 30.02 430 -1.34 1.58 -0.25 0.15 -1.02 0.56 25 X Y 201304210322A 194 | 141.19 35.73 34 2.58 -0.20 -2.38 0.83 2.37 -0.85 24 X Y 201304291301A 195 | 141.72 37.71 50 1.00 -0.14 -0.85 0.35 0.82 -0.36 25 X Y 201305180548A 196 | 141.83 38.23 54 0.67 0.09 -0.76 -0.02 -0.14 -0.10 25 X Y 201308040328A 197 | 145.22 46.97 378 1.17 3.11 -4.28 -4.68 -3.47 -1.78 24 X Y 201308041556A 198 | 133.65 42.19 455 -0.05 -0.21 0.26 0.06 -3.99 0.49 24 X Y 201309020251A 199 | 138.79 30.02 412 -2.80 2.56 0.23 2.41 -6.17 -1.64 25 X Y 201309040018A 200 | 144.66 37.17 25 -6.29 0.05 6.24 -1.38 -0.64 0.23 26 X Y 201310251710A 201 | 144.69 37.21 34 -5.89 1.63 4.26 -0.10 -0.24 0.54 24 X Y 201310252127A 202 | 144.53 37.21 32 -3.17 -0.17 3.34 0.73 0.84 0.69 24 X Y 201310271813A 203 | 140.21 35.95 65 2.88 0.24 -3.12 -0.27 1.63 -0.56 24 X Y 201311092237A 204 | 142.26 41.51 56 1.96 -0.07 -1.89 0.64 1.19 -0.77 24 X Y 201311150200A 205 | 140.26 35.57 58 2.13 0.22 -2.35 -0.24 0.48 -0.05 24 X Y 201311161144A 206 | 136.86 34.35 340 3.07 -0.44 -2.62 1.33 -2.82 0.74 24 X Y 201311181910A 207 | 149.24 44.53 28 0.33 -0.06 -0.28 1.55 0.45 -0.82 25 X Y 201312081724A 208 | 140.79 35.60 41 2.46 -0.22 -2.24 0.35 1.46 -1.03 24 X Y 201312140406A 209 | 140.79 35.61 45 2.26 -0.13 -2.13 -0.14 1.54 -0.36 24 X Y 201312210134A 210 | 142.23 35.52 12 1.50 -0.26 -1.25 0.44 1.15 -0.60 24 X Y 201312230657A 211 | 131.82 33.62 83 1.07 -0.37 -0.71 -1.06 3.00 -1.66 25 X Y 201403131706A 212 | 141.92 36.35 23 1.44 -0.15 -1.29 0.44 1.17 -0.50 24 X Y 201403290153A 213 | 144.22 37.91 12 -3.13 2.63 0.50 0.02 0.81 0.93 24 X Y 201405020915A 214 | 139.48 34.89 165 -0.35 -0.37 0.73 0.63 0.99 0.22 25 X Y 201405042018A 215 | 141.07 39.45 89 -1.22 -0.06 1.28 -0.90 -2.38 0.28 24 X Y 201406141731A 216 | 141.71 36.64 20 -3.24 0.31 2.93 0.20 0.22 0.61 24 X Y 201406151819A 217 | 141.21 37.12 53 3.22 -0.51 -2.71 0.82 2.22 -1.28 24 X Y 201406152014A 218 | 142.18 39.73 52 3.23 -0.01 -3.21 0.82 2.67 -0.65 24 X Y 201407042242A 219 | 142.39 36.97 12 -5.42 -0.83 6.25 -2.13 -4.31 2.13 25 X Y 201407111922A 220 | 148.94 44.63 67 2.11 0.43 -2.54 1.72 1.36 -0.69 25 X Y 201407201832A 221 | 142.37 41.25 49 1.32 0.03 -1.36 0.58 1.03 -0.35 25 X Y 201408100343A 222 | 132.11 32.07 28 5.19 -2.02 -3.17 2.94 5.41 -3.19 24 X Y 201408281914A 223 | 139.85 36.07 51 1.76 -1.20 -0.55 0.85 1.06 -0.75 24 X Y 201409160328A 224 | 142.68 40.23 37 2.76 -0.25 -2.51 1.15 1.78 -0.81 24 X Y 201410030057A 225 | 143.25 41.06 24 0.86 -0.05 -0.81 0.48 1.02 -0.24 25 X Y 201410110235A 226 | 143.42 41.01 31 1.33 -0.06 -1.28 0.59 1.43 -0.33 24 X Y 201410110520A 227 | 140.60 31.86 71 -2.21 1.43 0.78 -1.09 0.03 1.22 24 X Y 201410152353A 228 | 141.69 37.37 42 2.49 -0.27 -2.23 0.94 1.98 -0.80 24 X Y 201411200151A 229 | 137.88 36.76 12 2.03 -0.72 -1.31 -1.22 -0.25 -2.07 25 X Y 201411221308A 230 | 141.69 37.43 44 0.80 -0.14 -0.66 0.23 0.67 -0.30 25 X Y 201412200929A 231 | 143.22 39.78 22 0.86 -0.02 -0.84 0.62 1.16 -0.22 26 X Y 201502162306B 232 | 143.83 39.85 12 1.36 0.07 -1.43 0.66 1.77 -0.41 25 X Y 201502200425A 233 | 143.74 39.83 12 0.58 0.02 -0.60 0.36 0.78 -0.17 25 X Y 201502211013A 234 | 143.32 39.65 22 3.57 -0.20 -3.37 1.65 4.06 -0.82 24 X Y 201502240228A 235 | 141.93 30.99 13 0.65 0.04 -0.69 -0.09 0.35 0.13 25 X Y 201502250701A 236 | 139.94 31.47 12 4.44 -2.73 -1.72 -1.06 -2.81 -0.83 24 X Y 201505021650A 237 | 142.05 31.09 12 0.72 -0.01 -0.72 -0.04 0.25 0.00 25 X Y 201505102125A 238 | 142.29 38.89 44 1.37 -0.03 -1.34 0.68 1.30 -0.38 26 X Y 201505122112A 239 | 143.00 30.80 17 -2.62 0.67 1.95 -0.25 -0.42 -1.39 25 X Y 201505301849A 240 | 142.24 41.49 54 1.19 -0.06 -1.12 0.39 0.86 -0.37 25 X Y 201506080601A 241 | 143.64 39.66 22 3.16 -0.18 -2.98 1.32 2.78 -0.63 24 X Y 201506100833A 242 | 143.54 39.65 19 2.23 0.00 -2.23 0.64 3.23 -0.20 24 X Y 201506110445A 243 | 143.52 39.61 25 1.80 -0.06 -1.74 1.49 2.91 -0.41 24 X Y 201506110451A 244 | 148.12 43.93 54 2.18 0.07 -2.25 2.58 1.30 -0.73 25 X Y 201507070510A 245 | 141.44 40.43 87 2.13 -0.55 -1.58 0.53 -1.43 -1.59 24 X Y 201507091832A 246 | 131.77 33.04 58 0.14 -0.82 0.68 0.15 1.10 1.68 24 X Y 201507121752A 247 | 141.80 31.16 12 0.60 0.14 -0.74 -0.03 0.28 -0.12 25 X Y 201509011524A 248 | 142.72 41.89 51 1.66 -0.45 -1.21 0.54 1.04 -0.69 24 X Y 201509121338A 249 | 128.67 31.13 12 -0.37 0.80 -0.44 0.08 0.38 1.53 26 X Y 201511132051B 250 | 128.69 31.39 12 -0.90 2.04 -1.13 2.00 0.60 4.56 24 X Y 201511141920A 251 | 146.56 43.33 76 -0.43 1.27 -0.84 2.00 1.34 -0.27 24 X Y 201511280251A 252 | 129.75 44.83 587 2.43 0.80 -3.23 -3.50 -3.66 -4.05 24 X Y 201601020422A 253 | 132.83 30.57 13 -0.38 -0.44 0.83 -0.05 0.36 0.25 25 X Y 201601050221A 254 | 141.01 44.46 254 -2.07 0.93 1.14 1.69 0.24 -1.80 25 X Y 201601111708A 255 | 142.78 42.00 54 0.93 -0.19 -0.74 0.29 0.73 -0.38 26 X Y 201601140325B 256 | 138.67 30.43 422 -4.00 4.11 -0.11 2.93 -1.29 1.38 24 X Y 201602141809A 257 | 136.45 33.22 19 0.54 -0.31 -0.23 0.56 0.49 -0.26 25 X Y 201604010239A 258 | 130.74 32.77 15 -0.33 2.06 -1.74 0.10 -0.47 1.08 25 X Y 201604141226A 259 | 130.74 32.81 19 -0.98 1.73 -0.75 0.29 -0.45 1.23 24 X Y 201604141307A 260 | 130.69 32.71 12 0.01 0.97 -0.98 0.22 -0.12 0.58 25 X Y 201604141503A 261 | 130.77 32.84 13 -1.66 5.06 -3.40 0.38 0.94 0.82 26 X Y 201604151625A 262 | 131.12 33.00 13 0.10 2.21 -2.31 0.16 -0.51 0.09 24 X Y 201604181142A 263 | 141.77 37.81 48 0.70 -0.06 -0.64 0.28 0.62 -0.22 25 X Y 201604201219A 264 | 147.65 43.66 74 -0.90 2.99 -2.09 2.66 1.87 -2.63 24 X Y 201605030000A 265 | 140.07 30.98 144 0.84 -0.12 -0.72 -0.44 -2.24 -1.05 24 X Y 201605210606A 266 | 142.24 37.08 22 2.58 -0.22 -2.35 1.03 2.59 -0.92 24 X Y 201606262257A 267 | 141.73 37.37 39 0.54 0.31 -0.85 -0.23 1.76 -0.65 24 X Y 201608150704A 268 | 143.91 40.33 12 0.68 -0.01 -0.67 0.41 1.03 -0.17 25 X Y 201608200901A 269 | 143.79 40.30 15 0.41 -0.04 -0.37 0.43 0.84 -0.14 25 X Y 201608201558A 270 | 137.78 30.69 473 -0.28 0.66 -0.38 0.46 -0.38 0.41 25 X Y 201608251704A 271 | 137.96 31.78 343 -0.24 -0.30 0.54 -1.47 -1.37 -0.04 24 X Y 201609031226A 272 | 141.13 36.35 40 1.79 -0.09 -1.70 0.51 1.29 -0.73 24 X Y 201609091153A 273 | 129.12 35.79 20 0.23 1.55 -1.78 0.72 -0.26 1.39 24 X Y 201609121132A 274 | 142.09 30.53 12 1.02 -0.02 -0.99 -0.47 0.40 0.20 25 X Y 201609201621A 275 | 141.74 34.47 12 1.97 -0.10 -1.87 -0.10 0.95 0.30 25 X Y 201609230014A 276 | 133.72 35.40 16 -0.19 1.94 -1.74 -0.15 -0.86 -2.07 25 X Y 201610210507B 277 | 148.26 44.07 35 0.64 -0.25 -0.39 0.17 0.46 -0.42 25 X Y 201610232025A 278 | 141.76 38.46 54 1.65 -0.06 -1.59 0.44 1.37 -0.41 25 X Y 201611112142A 279 | 141.46 37.31 12 -2.97 1.67 1.30 0.80 0.75 1.50 26 X Y 201611212059A 280 | 141.28 37.19 12 -6.41 0.84 5.58 1.50 2.18 1.97 24 X Y 201611232123A 281 | 140.52 36.82 12 -0.77 0.26 0.52 0.17 0.21 -0.32 25 X Y 201612281238A 282 | 140.74 35.64 45 1.84 -0.16 -1.68 0.03 1.13 -0.40 24 X Y 201702190919A 283 | 141.45 37.56 53 4.74 -0.40 -4.34 0.96 3.16 -1.30 24 X Y 201702280749A 284 | 142.16 41.49 55 3.56 -0.43 -3.13 1.23 2.32 -1.27 24 X Y 201703071608A 285 | 131.43 30.82 26 -4.89 -0.02 4.91 1.84 2.89 0.71 24 X Y 201704291232A 286 | 147.85 43.88 55 1.56 1.09 -2.65 1.55 0.79 -0.85 24 X Y 201705221100A 287 | 147.98 43.77 65 -0.50 1.78 -1.27 1.52 1.27 -0.44 24 X Y 201706281207A 288 | 131.56 40.89 556 -0.08 0.67 -0.59 -0.15 -0.47 -0.54 25 X Y 201707121948A 289 | 141.70 37.32 42 6.22 -0.83 -5.39 2.09 4.80 -1.88 24 X Y 201707200011A 290 | 143.47 40.09 24 2.09 -0.20 -1.88 0.83 1.84 -0.35 24 X Y 201707231535A 291 | 142.96 41.80 45 2.93 -0.73 -2.21 0.82 1.87 -1.04 24 X Y 201709100844A 292 | 144.69 38.02 17 -1.83 0.17 1.66 0.43 1.75 0.95 25 X Y 201709201637A 293 | 142.53 40.31 44 0.72 -0.01 -0.71 0.28 0.56 -0.17 25 X Y 201709262022A 294 | 144.00 37.48 14 -1.84 -0.18 2.01 0.25 0.88 1.43 25 X Y 201710060759A 295 | 141.25 37.10 53 4.48 -0.77 -3.71 1.14 2.92 -1.58 24 X Y 201710061456A 296 | 141.52 32.45 19 1.07 0.08 -1.15 -0.07 0.46 -0.03 25 X Y 201711090742A 297 | 145.01 38.07 14 -1.10 0.42 0.69 0.09 0.08 0.48 25 X Y 201711122224A 298 | 129.21 36.13 12 1.12 1.08 -2.21 -0.62 -0.82 -0.62 24 X Y 201711150529A 299 | 140.55 32.29 84 2.24 -0.04 -2.20 0.12 -6.24 -4.15 24 X Y 201711160943A 300 | 146.34 45.31 17 1.91 -0.81 -1.11 -0.87 -1.08 -0.95 24 X Y 201711292232A 301 | 132.14 30.77 12 -2.14 -0.02 2.17 0.92 0.31 0.40 24 X Y 201712201340A 302 | 142.55 41.07 44 3.97 0.16 -4.13 1.79 3.32 -1.12 25 X Y 201801241051A 303 | 141.71 37.55 38 1.75 -0.97 -0.78 3.54 1.00 -0.98 24 X Y 201802251628A 304 | 147.96 44.98 109 -0.51 0.53 -0.02 -0.89 -1.94 0.27 24 X Y 201803041952A 305 | 140.75 32.58 66 2.60 0.51 -3.11 0.66 -2.19 -2.91 24 X Y 201803251437A 306 | 132.48 35.17 18 -0.15 4.28 -4.13 -0.23 -1.18 -2.05 24 X Y 201804081632A 307 | 145.73 43.34 92 -1.39 1.36 0.04 1.94 0.66 0.76 24 X Y 201804240853A 308 | 145.49 42.72 31 3.97 -2.21 -1.75 2.08 3.75 -2.66 24 X Y 201805171842A 309 | 148.05 43.94 15 0.72 0.63 -1.35 -1.43 -1.45 0.87 24 X Y 201806020948A 310 | 131.69 31.23 37 1.97 -0.25 -1.73 0.58 1.05 -0.55 24 X Y 201806111954A 311 | 135.47 34.83 17 1.65 1.84 -3.50 0.22 -1.13 -0.87 24 X Y 201806172258A 312 | 140.56 35.11 45 0.22 0.64 -0.86 -0.11 0.30 -0.48 25 X Y 201807071123A 313 | 141.20 37.21 12 -1.74 1.22 0.51 0.26 0.62 0.85 24 X Y 201807310842A 314 | 144.01 37.87 30 -2.68 -0.20 2.88 -1.18 0.21 1.52 24 X Y 201808071512A 315 | 142.04 42.75 31 0.77 0.01 -0.78 -0.11 0.83 0.42 26 X Y 201809051807A 316 | 142.01 42.78 36 2.30 -0.28 -2.02 -0.03 -0.45 0.45 24 X Y 201809052111A 317 | 141.94 38.31 50 4.54 -0.12 -4.43 1.86 3.17 -1.20 24 X Y 201810251836A 318 | 145.68 44.59 12 3.83 -1.85 -1.98 -1.18 -1.39 -2.12 24 X Y 201810260304A 319 | 140.37 31.67 105 -0.81 1.20 -0.39 -0.03 -1.95 -0.32 24 X Y 201810290127A 320 | 145.68 44.60 12 1.06 -0.13 -0.93 -0.36 -0.41 -0.41 25 X Y 201811041926A 321 | 130.11 30.43 129 0.83 -0.23 -0.60 0.40 2.47 -0.30 24 X Y 201811201909A 322 | 143.18 41.47 12 1.53 -2.12 0.59 -0.35 -1.00 -0.31 24 X Y 201811280223A 323 | 141.83 33.98 20 3.90 -0.43 -3.46 0.51 1.33 -0.10 24 X Y 201812240018A 324 | 141.96 33.94 22 3.47 -0.57 -2.90 0.15 1.37 0.01 24 X Y 201812240043A 325 | --------------------------------------------------------------------------------