Removed unused parts of GMT-script doing map projections.

parent 13f0be0a
......@@ -22,29 +22,17 @@
# @section DESCRIPTION
# Plots and projects the data.
##
# points of interest in lon/lat
poi_lonlat=../poi_lonlat.csv
# receivers in lon/lat
receivers_lonlat=../receivers_lonlat.csv
# create output directory
out_dir=out_pe
mkdir ${out_dir}
# points of interest in meters
poi=${out_dir}/poi.csv
# receivers in meters
receivers=${out_dir}/receivers.csv
# region of interest cut from the mesh
region_cut=-50012.5/50012.5/-50012.5/50012.5
region_cut=-50025/50025/-50025/50025
# sampling of the topo mesh (25m)
region_sampling=25/25
# sampling of the topo mesh (50m)
region_sampling=50/50
# projection center is the epicenter of the 2004 Parfield earthquake
# projection center is the epicenter of the "M 5.1 - 2km NW of Brea, CA"-event https://earthquake.usgs.gov/earthquakes/eventpage/ci15481673/executive
projection="tmerc/-117.916/33.933"
grids=(grids/n34w119/grdn34w119_13/w001001.adf grids/n34w118/grdn34w118_13/w001001.adf grids/n35w119/grdn35w119_13/w001001.adf grids/n35w118/grdn35w118_13/w001001.adf)
......@@ -64,18 +52,6 @@ gmt grdblend ${grids[*]} -R${blend_region} -G${tmp_dir}/blend.nc -V2
# plot topography in lon/lat
gmt grdimage -J${projection}/8c ${tmp_dir}/blend.nc -Ba -E75 -C${tmp_dir}/gmt_colors.cpt > ${tmp_dir}/blend.ps
# convert receiver coordinates
#cat ${receivers_lonlat} | cs2cs +proj=latlong +datum=WGS84 +to +proj=tmerc +units=m +axis=end +no_defs +ellps=GRS80 +k=0.9996 +lon_0=-120.365 +lat_0=35.817 -f %.12f > ${receivers}
# add receivers (symbols only)
#gmt psxy ${receivers} -R${blend_region} -Z -J${projection}/8c -Sd.1 -Gblack -O -K -V2 >> ${tmp_dir}/blend.ps
# add points of interests
#gmt psxy ${poi_lonlat} -R${blend_region} -Z -J${projection}/8c -Sd.25 -Gblack -O -K -V2 >> ${tmp_dir}/blend.ps
# add captions
#gmt pstext ${poi_lonlat} -R${blend_region} -J${projection}/8c -F+f7+jL -D0.2/0c -O -V2 >> ${tmp_dir}/blend.ps
# convert to pdf
ps2pdf ${tmp_dir}/blend.ps ${out_dir}/blend.pdf
......@@ -91,22 +67,9 @@ gmt grd2xyz ${tmp_dir}/topo.nc > ${tmp_dir}/topo.xyz
# replace all nans (under water) with zeroes, grdfill tends to segfault: manual sed
sed 's/NaN/0/g' ${tmp_dir}/topo.xyz > ${out_dir}/topo.xyz
# convert point of interests
#cat ${poi_lonlat} | cs2cs +proj=latlong +datum=WGS84 +to +proj=tmerc +units=m +axis=end +no_defs +ellps=GRS80 +k=0.9996 +lon_0=-120.365 +lat_0=35.817 -f %.12f > ${poi}
# plot topography
gmt grdimage -JX6.5i -C${tmp_dir}/colors.cpt ${tmp_dir}/topo.nc -Ba -C${tmp_dir}/gmt_colors.cpt -V2 > ${tmp_dir}/map.ps
# add receivers
#gmt psxy ${receivers} -Z -JX6.5i -R${region_cut} -Sd.1 -Gblack -O -K -V2 >> ${tmp_dir}/map.ps
#gmt pstext ${receivers} -Z -R${region_cut} -JX6.5i -F+f3+jL -D0.05/0c -O -K -V2 >> ${tmp_dir}/map.ps
# add points of interests
#gmt psxy ${poi} -Z -JX6.5i -R${region_cut} -Sd.25 -Gblack -O -K -V2 >> ${tmp_dir}/map.ps
#gmt pstext ${poi} -Z -R${region_cut} -JX6.5i -F+f7+jL -D0.2/0c -O -V2 >> ${tmp_dir}/map.ps
# convert to pdf
ps2pdf ${tmp_dir}/map.ps ${out_dir}/map.pdf
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment