[Shake-dev] shakemap topography

Fri Mar 9 10:51:22 UTC 2012

Hi all,

We are trying to run shakemap but we have some problems with topography. The example runs fine but, when we change the topography of the example and the region of the shakemap by the input data, it doesn't work.

We think that the problem is the topography file called Sispyr_m.grd. In fact, when shakemap runs the grdsample, the program says that the "grdsample:  Selected region exceeds the X-boundaries of the grid file!". And then, if you follow the command line you'll see the error: "mapping: WARNING: can't find topo data: making mi without topography".

We've already checked the config files so many times and we think that we've put it correct. We also tried to do the shakemap of our region (Pyrenees) using the topography of the example and it doesn't work too because we wanted to try if the problem was the topography file, but now we don't know what causes the error.

Our version of the Shakemap is the 3.5 and the release is 416.

We attached the command line and in cursive you could find where the problem appears the first time.

Thank you very much.

Best Regards

Jordi Roviró i Pou
Unitat de Risc Sísmic
Institut Geològic de Catalunya
Balmes, 209-211 · E-08006 Barcelona<http://mercuri.icc.cat/website/mob_nf/mob1/mob2/inici2.htm?&zoom=3&XY=429136:4583504&ident=Institut%20Geologic%20de%20Catalunya%3Cbr%3E%20c/Balmes,%20209-211%3Cbr%3E%20E-08006%20Barcelona%3Cbr%3E%20Telefon%2093%20553%2084%2030&titol=IGC&eina=2&base=m&adreca=Balmes,%20209>
Tel. 34 93 553 84 30 (2418)
jroviro at igc.cat<mailto:mail at igc.cat> www.igc.cat<http://www.igc.cat/>
[cid:image001.jpg at 01CCFDEA.F2DCC130]

---------------------------------------------------------------
Shakemap command line:

xsismica at nobel ./shake -event 20101219 -V
shake: Running: '/home/xsismica/shakemap_treball/bin/../bin/pending -event 20101219'
transfer: ----- Starting Transfer at 03/09/2012 10:20:54 -----
transfer: ----- Transfer finished at 03/09/2012 10:20:54 -----
shake: Done with 'pending'
shake: Running: '/home/xsismica/shakemap_treball/bin/../bin/grind -event 20101219 -V'
Map bound: 1.550000/4.050000/40.928750/42.791250
Using GMICE::Wald99 for PGM2I function
Using GMICE::Wald99 for Inverse PGM2I function
Using GMPE::Small for GMPE
Using GMPE::Small for GMPE
grind: Bias -0.81 -0.31
Outlier standard deviations: 3.
Outlier standard deviations: 3.
grind: 0 stations flagged out this iteration
Using GMPE::Small for GMPE
Using DefaultIPE for IPE
Using GMPE::Small for GMPE
Using GMICE::Wald99 for PGM2I function
Using GMICE::Wald99 for Inverse PGM2I function
Using DefaultIPE for IPE
Using GMPE::Small for GMPE
Using GMICE::Wald99 for PGM2I function
Using GMICE::Wald99 for Inverse PGM2I function
Intensity Bias=0.59
0 stations flagged out this iteration
Using DefaultIPE for IPE
Using GMPE::Small for GMPE
Using GMICE::Wald99 for PGM2I function
Using GMICE::Wald99 for Inverse PGM2I function
Uncertainty computation skipped (no Intensity > 6.0)
shake: Done with 'grind'
shake: Running: '/home/xsismica/shakemap_treball/bin/../bin/tag -event 20101219'
shake: Done with 'tag'
shake: Running: '/home/xsismica/shakemap_treball/bin/../bin/mapping -event 20101219 -timestamp -itopo -gsm -V'
Running: /home/xsismica/shakemap_treball/bin/../bin/echo 4.050000 42.791250 | /usr/local/bin//mapproject -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51
Running: /home/xsismica/shakemap_treball/bin/../bin/echo 4.050000 42.791250 | /usr/local/bin//mapproject -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM25.4
Running: /usr/local/bin//psbasemap -Ba60mf30m/a30mf30mWSen -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -X2.54 -Y2.54 -K
Running: /usr/local/bin//psbasemap -Ba60mf30m/a30mf30mWSen -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -X2.54 -Y6.35 -K
Running: /usr/local/bin//psxy /home/xsismica/shakemap_treball/bin/../lib/mapping/ca_roads.xy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -m -W6/153/153/153 -O -K
Running: /usr/local/bin//pscoast -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -N1/5 -N2/5 -Df -W5 -S120/160/220 -Ia/2/120/160/220 -O -K
Using GMT timezone.
Running: /home/xsismica/shakemap_treball/bin/../bin/echo 2.800000 41.860000 | /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -Sa0.762 -W6/255/0/0 -L -O -K
Running: /home/xsismica/shakemap_treball/bin/../bin/echo "0 -1.143 9 0 1 5 Map Version 1 Processed Fri Mar  9, 2012 10:20:56 AM GST " | /usr/local/bin//pstext -JX16.4553191047/16.51 -R0/16.4553191047/0/16.51 -N -O -K
Running: /usr/local/bin//psbasemap -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -Ba60mf30m/a30mf30mWSen -Lfx2.54/1.78/41.8600/50 -O

Using topo file /home/xsismica/shakemap_treball/bin/../lib/mapping/sispyr_m.grd, disabling -gsm.
Running: /usr/local/bin//grdsample /home/xsismica/shakemap_treball/bin/../lib/mapping/sispyr_m.grd -G/home/xsismica/shakemap_treball/bin/../data/20101219/mapping/topo.grd -I30c/30c -R1.5500000000/4.0500000000/40.9287500000/42.7912500000
grdsample:  Selected region exceeds the X-boundaries of the grid file!
Running: /usr/local/bin//grdsample /home/xsismica/shakemap_treball/bin/../lib/mapping/intensity.grd -G/home/xsismica/shakemap_treball/bin/../data/20101219/mapping/inten.grd -I30c/30c -R1.5500000000/4.0500000000/40.9287500000/42.7912500000
grdsample:  Selected region exceeds the X-boundaries of the grid file!
Sampling to a finer grid; old: 0.025000/0.025169 new: 30c/30c
Running: /usr/local/bin//grdsample /home/xsismica/shakemap_treball/bin/../data/20101219/output/mi.grd -G/home/xsismica/shakemap_treball/bin/../data/20101219/mapping/tmp_mi.grd -I30c/30c
mapping: WARNING: can't find topo data: making mi without topography

Running: /usr/local/bin//grdimage /home/xsismica/shakemap_treball/bin/../data/20101219/mapping/tmp_mi.grd -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000/10/3000 -JM16.51 -C/home/xsismica/shakemap_treball/bin/../lib/mapping/Ii.cpt -O -K
mapping: WARNING: can't find topo data: making mi without topography
Running: /usr/local/bin//grdimage /home/xsismica/shakemap_treball/bin/../data/20101219/mapping/tmp_mi.grd -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000/10/3000 -JM25.4 -C/home/xsismica/shakemap_treball/bin/../lib/mapping/Ii_ntsc.cpt -O -K
Running: /usr/local/bin//grdcontour /home/xsismica/shakemap_treball/bin/../data/20101219/output/pga.grd -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -C0.2 -W5/255/255/255 -Af12/255/255/255 -D/home/xsismica/shakemap_treball/bin/../data/20101219/mapping/cont_pga.dat -Q8 -m -L0.1/9999 -O -K
Running: /usr/local/bin//grdcontour /home/xsismica/shakemap_treball/bin/../data/20101219/output/pgv.grd -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -C0.04 -W5/255/255/255 -Af12/255/255/255 -D/home/xsismica/shakemap_treball/bin/../data/20101219/mapping/cont_pgv.dat -Q8 -m -L0.02/9999 -O -K
Running: /usr/local/bin//grdimage /home/xsismica/shakemap_treball/bin/../data/20101219/output/urat_pga.grd -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -C/home/xsismica/shakemap_treball/bin/../lib/mapping/sd.cpt -O -K
Running: /usr/local/bin//psscale -C/home/xsismica/shakemap_treball/bin/../lib/mapping/sd.cpt.legend.no1 -D6.7i/2.75i/5.00i/0.2i -O -K
Running: /usr/local/bin//psscale -Bp -C/home/xsismica/shakemap_treball/bin/../lib/mapping/sd.cpt.legend -D6.7i/2.75i/5.00i/0.2i -O -K
Running: /usr/local/bin//psscale -C/home/xsismica/shakemap_treball/bin/../lib/mapping/sd.cpt.legend.0 -D6.7i/0.1i/0.2i/0.2i -O -K
Running: echo 17.145 12.446 12 0 1 5 F  | /usr/local/bin//pstext -JX16.4553191047/16.51 -R0/16.4553191047/0/16.51 -N -O -K
Running: echo 17.145 1.524 12 0 1 5 A  | /usr/local/bin//pstext -JX16.4553191047/16.51 -R0/16.4553191047/0/16.51 -N -O -K
Running: echo 17.145 8.636 12 0 1 5 D  | /usr/local/bin//pstext -JX16.4553191047/16.51 -R0/16.4553191047/0/16.51 -N -O -K
Running: echo 17.145 4.318 12 0 1 5 C  | /usr/local/bin//pstext -JX16.4553191047/16.51 -R0/16.4553191047/0/16.51 -N -O -K
Running: echo 17.145 2.845 12 0 1 5 B  | /usr/local/bin//pstext -JX16.4553191047/16.51 -R0/16.4553191047/0/16.51 -N -O -K
Running: /usr/local/bin//mapproject -JM16.51 -R1.5500000000/4.0500000000/40.9287500000/42.7912500000
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -W1 -L -O -K
Done plotting stations.Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G255/255/0 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G0/224/224 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G102/102/255 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G255/102/102 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G0/0/255 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G255/255/0 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G102/255/102 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G0/0/255 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G255/120/0 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G0/0/255 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G102/255/102 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G255/0/255 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G102/102/102 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G204/0/0 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G0/153/0 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G255/255/0 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G0/224/224 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G102/102/255 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G255/102/102 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G0/0/255 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G255/255/0 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G102/255/102 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G0/0/255 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G255/120/0 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G0/0/255 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G102/255/102 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G255/0/255 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G102/102/102 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G204/0/0 -O -K
Running: /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -S -G0/153/0 -O -K
Running: /home/xsismica/shakemap_treball/bin/../bin/echo 2.800000 41.860000 | /usr/local/bin//mapproject -JM16.51 -R1.5500000000/4.0500000000/40.9287500000/42.7912500000
mapping: WARNING: couldn't find topo and intensity data; will continue without topo.
Making urat_pga plot
/usr/local/bin//grdcontour /home/xsismica/shakemap_treball/bin/../data/20101219/output/mi.grd -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -C6 -W20/0 -A- -Q50 -m -L0.5/9999 -O -K
Running: /usr/local/bin//grdcontour /home/xsismica/shakemap_treball/bin/../data/20101219/output/mi.grd -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM16.51 -C6 -W20/0 -A- -Q50 -m -L0.5/9999 -O -K
Printing grade - (-0.000)
Running: echo  '0.508 0.635 12 0 1 5 GRADE: - (Mean sigma: -)' | /usr/local/bin//pstext -JX19.05/16.4553191047 -R0/19.05/0/16.4553191047 -C.05i -W255o4/0 -N -O -K
Running: echo  '13.53 0.508 12 0 1 5 MI = 6.0' | /usr/local/bin//pstext -JX19.05/16.4553191047 -R0/19.05/0/16.4553191047 -N -O -K
Making pgv plot
Making pga plot
Making mi plot
Done with mi plot
Making tvmap
Running: /usr/local/bin//psbasemap -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM25.4 -Bwesn -X2.54 -Y2.54 -K
Running: /usr/local/bin//psxy /home/xsismica/shakemap_treball/bin/../lib/mapping/ca_roads.xy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM25.4 -m -W8/153/153/153 -O -K
Running: /usr/local/bin//pscoast -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM25.4 -Na/10/102/102/102 -Df -W10/102/102/102 -S120/160/220 -O -K
Running: /home/xsismica/shakemap_treball/bin/../bin/echo 2.800000 41.860000 | /usr/local/bin//psxy -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM25.4 -Sa0.762 -W6/255/0/0 -L -O -K
Running: /usr/local/bin//psbasemap -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM25.4 -Bwesn -Lfx5.54/2.54/41.8600/50m -O
Running: /usr/local/bin//psbasemap -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM25.4 -Ba60mf30m/a30mf30mWSen -O
Done with tvmap
Running: /usr/local/bin//grdimage /home/xsismica/shakemap_treball/bin/../data/20101219/mapping/tmp_mi.grd -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM25.4 -C/home/xsismica/shakemap_treball/bin/../lib/mapping/Ii_ntsc.cpt -K
Running: /usr/local/bin//pscoast -Df -P -R1.5500000000/4.0500000000/40.9287500000/42.7912500000 -JM25.4 -S0 -O
shake: Done with 'mapping'
shake: Running: '/home/xsismica/shakemap_treball/bin/../bin/genex -event 20101219 -zip -metadata -shape shape'
mp 2.9.12 - Peter N. Schweitzer (U.S. Geological Survey)
Info: input file = metadata.txt
Info: process date = 20120309
Info: process time = 10:20:59
No errors
shake: Done with 'genex'
shake: Running: '/home/xsismica/shakemap_treball/bin/../bin/setversion -event 20101219'
shake: Done with 'setversion'

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20120309/13646d2d/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.jpg
Type: image/jpeg
Size: 15856 bytes
Desc: image001.jpg
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20120309/13646d2d/attachment-0001.jpg>