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

Wed Jan 7 20:41:51 UTC 2015

Bruce,

More GMT fun :)

I will try to install GMT 4.5.11 from sources and get back to you.

*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 3:36 PM, Bruce Worden <cbworden at gmail.com> wrote:

> 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/c6c2e18d/attachment-0001.html>