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

Bruce Worden cbworden at gmail.com
Wed Jan 7 17:52:03 UTC 2015 

Hi Branden,

Have you tried running ShakeMap on another event? Do you get the same result?

Bruce

On Jan 5, 2015, at 4:29 PM, Branden Christensen <branden.christensen at osop.com.pa> wrote:

> Hi all:
> 
> 
> Happy new year!
> 
> 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 and have a great start to the new year. 
> 
> 
> Branden
> --
> _______________________________________________
> Shake-dev mailing list
> Shake-dev at geohazards.usgs.gov
> https://geohazards.usgs.gov/mailman/listinfo/shake-dev

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20150107/09373f75/attachment.html>

More information about the Shake-dev mailing list