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

Wed Jan 7 19:50:49 UTC 2015

Hi Bruce:

I get the same result for all events.

Best,

*Branden Christensen*
Director, OSOP
Sign up for the OSOP mailing list:
http://www.osop.com.pa/about/osop-mailing-list/

OSOP
--.

On Wed, Jan 7, 2015 at 12:52 PM, Bruce Worden <cbworden at gmail.com> wrote:

> 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
>
>
>
> _______________________________________________
> 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/6525dbef/attachment-0001.html>