[Shake-dev] mapping program hangs then GMT Fatal Error
Bruce Worden
cbworden at gmail.com
Wed Jan 7 20:36:26 UTC 2015
Hi Branden,
That suggests an issue with grdcontour on your system. For some reason GMT and Ubuntu seem to have issues (this is something that Paul Wessel has commented on in the past). My only suggestion is to take one of the data files that's causing the problem and try to isolate the error. If you can't find a workaround, submit it as a GMT bug.
Bruce
On Jan 7, 2015, at 11:50 AM, Branden Christensen <branden.christensen at osop.com.pa> wrote:
> 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/4bd1c466/attachment.html>
More information about the Shake-dev
mailing list