[Shake-dev] mapping program hangs then GMT Fatal Error

Mon Dec 15 03:12:29 UTC 2014

Hi all:

I am running shakemap revision 1362 on Ubuntu 14.04 for the example event
*9583161*.

mapping hangs then spits a GMT *error*:

osop at seismo1:/opt/shakemap/config$ /opt/shakemap/bin/../bin/mapping -event
*9583161* -timestamp -itopo -gsm -verbose
Running: /opt/shakemap/bin/../bin/echo 0.000000 19.046508 |
/usr/lib/gmt/bin/mapproject -I
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-JM16.51
Running: /opt/shakemap/bin/../bin/echo -113.900000 38.683000 |
/usr/lib/gmt/bin/mapproject
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-JM16.51
Running: /opt/shakemap/bin/../bin/echo -113.900000 38.683000 |
/usr/lib/gmt/bin/mapproject
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-JM25.4
Running: /opt/shakemap/bin/../bin/echo -113.900000 38.683000 |
/usr/lib/gmt/bin/mapproject
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-JQ0/16.51
Running: /usr/lib/gmt/bin/psbasemap -Bw -P
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-JM16.51 -X2.54 -Y2.54 -K
Running: /usr/lib/gmt/bin/psbasemap -Bw -P
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-JM16.51 -X2.54 -Y6.35 -K
Running: /usr/lib/gmt/bin/psxy
/opt/shakemap/bin/../lib/mapping/global_faults.txt -P
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-JM16.51 -W4/204/0/0 -O -K -m
Running: /usr/lib/gmt/bin/psxy
/opt/shakemap/bin/../lib/mapping/global_faults.txt -P
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-JM25.4 -W4/204/0/0 -O -K -m
Running: /usr/lib/gmt/bin/pscoast -P
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-JM16.51 -N1/5 -N2/5 -Df -W5 -S120/160/220 -Ia/2/120/160/220 -O -K
Running: /opt/shakemap/bin/../bin/echo -118.900000 34.563000 |
/usr/lib/gmt/bin/psxy -P
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-JM16.51 -Sa0.762 -W6/255/0/0 -L -O -K
Running: /opt/shakemap/bin/../bin/echo "0 -1.143 9 0 1 5 Map Version 1
Processed 2014-12-15 02:56:10 UTC " | /usr/lib/gmt/bin/pstext
-JX16.4727084188/16.51 -R0/16.4727084188/0/16.51 -N -O -K
Running: /usr/lib/gmt/bin/psbasemap -P
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-JM16.51 -Ba120mf60mWSen -Lfx3.81/1.78/34.5630/150+l -O
2014-12-14 21:56:10 : mapping: Creating topography files.

/opt/shakemap/bin/../bin/topo2grd 9583161
-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
2m regime=active 1
Using /opt/shakemap/dem. Final resolution fixed at .
/opt/shakemap/dem/dem2grd.pl
-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
tmp.topo2grd.14454
HELLO WORLD. I am /opt/shakemap/dem/dem2grd.pl
EXTRACTING w140n40.grd W:-123.90000000000001 E:-113.90000000000001
S:30.4333333333333 N:38.6833333333333
/usr/lib/gmt/bin//grdcut
-R-123.90000000000001/-113.90000000000001/30.4333333333333/38.6833333333333
w140n40.grd -Gtmp.14455.w140n40
/usr/lib/gmt/bin//grdcut
-R-123.90000000000001/-113.90000000000001/30.4333333333333/38.6833333333333
w140n40.grd -Gtmp.14455.w140n40
mv tmp.14455.w140n40 tmp.14455.w140
mv tmp.14455.w140 tmp.topo2grd.14454.grd
Done with pasting.
FINAL BOUNDS:
-R-123.90000000000001/-113.90000000000001/30.4333333333333/38.6833333333333
2014-12-14 21:56:10 : mapping: Topography file created.

Making topo intensity file :
/usr/lib/gmt/bin/grdgradient
/opt/shakemap/bin/../data/9583161/mapping/topo_data.grd -A45
-G/opt/shakemap/bin/../data/9583161/mapping/intens_data.grd -Ne0.6...done.
/usr/lib/gmt/bin/grdcut -V
/opt/shakemap/bin/../data/9583161/mapping/topo_data.grd
-G/opt/shakemap/bin/../data/9583161/mapping/Q7xr6eDZaP.grd
-R-123.9/-113.90000000000001/30.4416666666666/38.68333333333
GMT_grd_is_global: no!
GMT_grd_is_global: no!
grdcut: File spec: W E S N dx dy nx ny:
grdcut: Old: -123.9 -113.9 30.4333333333 38.6833333333 0.00833333333333
0.00833333333333 1200 990
grdcut: New: -123.9 -113.9 30.4416666667 38.6833333333 0.00833333333333
0.00833333333333 1200 989
/usr/lib/gmt/bin/grdsample -V
/opt/shakemap/bin/../data/9583161/mapping/Q7xr6eDZaP.grd
-G/opt/shakemap/bin/../data/9583161/mapping/topo.grd
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-I30c/30c -Ql
grdsample: Given domain implies y_inc = 0.00833165
GMT_grd_is_global: no!
GMT_boundcond_param_prep determined edgeinfo: gn = 0, gs = 0, nxp = 0, nyp
= 0
grdsample: New grid (-123.9/-113.9/30.443/38.683) nx = 1200 ny = 989 dx =
0.00833333333333 dy = 0.00833164812942 node_offset = 1
GMT_grd_is_global: no!
/usr/lib/gmt/bin/grdcut -V
/opt/shakemap/bin/../data/9583161/mapping/intens_data.grd
-G/opt/shakemap/bin/../data/9583161/mapping/lFn3LcANAe.grd
-R-123.9/-113.90000000000001/30.4416666666666/38.68333333333
GMT_grd_is_global: no!
GMT_grd_is_global: no!
grdcut: File spec: W E S N dx dy nx ny:
grdcut: Old: -123.9 -113.9 30.4333333333 38.6833333333 0.00833333333333
0.00833333333333 1200 990
grdcut: New: -123.9 -113.9 30.4416666667 38.6833333333 0.00833333333333
0.00833333333333 1200 989
/usr/lib/gmt/bin/grdsample -V
/opt/shakemap/bin/../data/9583161/mapping/lFn3LcANAe.grd
-G/opt/shakemap/bin/../data/9583161/mapping/inten.grd
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-I30c/30c -Ql
grdsample: Given domain implies y_inc = 0.00833165
GMT_grd_is_global: no!
GMT_boundcond_param_prep determined edgeinfo: gn = 0, gs = 0, nxp = 0, nyp
= 0
grdsample: New grid (-123.9/-113.9/30.443/38.683) nx = 1200 ny = 989 dx =
0.00833333333333 dy = 0.00833164812942 node_offset = 1
GMT_grd_is_global: no!
Sampling to a finer grid; old: 0.025000/0.024970 new: 30c/30c
Running: /usr/lib/gmt/bin/grdsample
/opt/shakemap/bin/../data/9583161/output/mi.grd
-G/opt/shakemap/bin/../data/9583161/mapping/tmp_mi.grd -I30c/30c
Running: /usr/lib/gmt/bin/grdview
/opt/shakemap/bin/../data/9583161/mapping/topo.grd -P
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000/10/3000
-JM16.51 -JZ1 -G/opt/shakemap/bin/../data/9583161/mapping/tmp_mi.grd
-C/opt/shakemap/bin/../lib/mapping/Ii.cpt
-I/opt/shakemap/bin/../data/9583161/mapping/inten.grd -Qi72 -O -K 2>
/dev/null
Running: /usr/lib/gmt/bin/grdview
/opt/shakemap/bin/../data/9583161/mapping/topo.grd -P
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000/10/3000
-JM25.4 -JZ1 -G/opt/shakemap/bin/../data/9583161/mapping/tmp_mi.grd
-C/opt/shakemap/bin/../lib/mapping/Ii_ntsc.cpt
-I/opt/shakemap/bin/../data/9583161/mapping/inten.grd -Qi72 -O -K 2>
/dev/null
Running: /usr/lib/gmt/bin/grdcontour
/opt/shakemap/bin/../data/9583161/mapping/_WrpwdfoQR.grd -P
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-JM16.51 -C/opt/shakemap/bin/../lib/mapping/mi_cont.txt
-Wafaint,125,0_100:0 -Wc5/255/255/255 -Af12/255/255/255
-D/opt/shakemap/bin/../data/9583161/mapping/cont_mi.xyz -Q20 -m -L1/10 -O -K

*mapping hangs for a long time here then spits this message:*

*GMT Fatal Error: grdcontour could not reallocate memory [8.00 Gb, n_items
= 1073741824]*
mapping: Error running command:
/usr/lib/gmt/bin/grdcontour
/opt/shakemap/bin/../data/9583161/mapping/_WrpwdfoQR.grd -P
-R-123.90000000000001/-113.90000000000001/30.44300000000000/38.68300000000000
-JM16.51 -C/opt/shakemap/bin/../lib/mapping/mi_cont.txt
-Wafaint,125,0_100:0 -Wc5/255/255/255 -Af12/255/255/255
-D/opt/shakemap/bin/../data/9583161/mapping/cont_mi.xyz -Q20 -m -L1/10 -O
-K|

Thank you for any insights.

Saludos,

*Branden Christensen*
Director, OSOP
Sign up for the OSOP mailing list:
http://www.osop.com.pa/about/osop-mailing-list/
Visit us at *AGU* next week: Booth 2136
LA SALSA VIVE...Gózala. Cali 2014.

OSOP
--
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20141215/0a82a65f/attachment.html>