[Shake-dev] shakemap topography

Fri Mar 9 15:32:05 UTC 2012

Hi Jordi,

Can you run the GMT program 'grdinfo' on your file Sispyr_m.grd and report the results?

Thanks,
Bruce

On Mar 9, 2012, at 2:51 AM, Roviró Pou, Jordi wrote:

> 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
> Tel. 34 93 553 84 30 (2418)
> jroviro at igc.cat www.igc.cat
> 
> <image001.jpg>
> 
>  
>  
> ---------------------------------------------------------------
> 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'
>  
> _______________________________________________
> 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/20120309/0917fa98/attachment-0001.html>