[Shake-dev] error runing grind

[Bruce Worden]

Hi Baishan,

If you use -qtm, you must specify a "qtm_file" in grind.conf. For the  
test event this will be 'usgs_cgs_geology_60s.grd'. You should also  
probably change "ampfactor_file" to 'site_corr_cdmg.dat' as well.

For your general use, you'll want to develop a Vs30 grid for your  
region, as well as the corresponding site_corr_*.dat. If you don't  
have a Vs30 grid, you can approximate one from topography using the  
Allen & Wald technique, which is built into the program topo2grd.  
That's what grind was attempting to do (on the fly) in the messages  
you included, but the digital elevation model is not generally  
distributed with ShakeMap (it's 3GB). If you obtain the DEM, you'll  
need to set the $DEMdir macro in the include/macros file to point to  
your DEM directory. The topo2grd program requires the DEM files to be  
in a fairly specific format. You can get DEM tiles from:

http://www.dgadv.com/srtm30/

Or other places by searching on SRTM30.

We use a command like:

my $command = "$gmtdir/xyz2grd $file -G$section.grd -R$lon/$lat -I30c - 
Ddeg/deg/m/1/0 -L -F -ZTLhw -N-9999 -V";

To convert the DEM file to a GMT grd. The program topo2grd can then  
cut a chunk out of this file and convert it to Vs30.

None of this is really intended to be part of the normal ShakeMap  
distribution. topo2grd exists for the use of the global ShakeMap group  
at NEIC, so if you are going to try to make it work, be prepared to do  
some hacking on your own. You're much better off making a static Vs30  
GMT grid file for your region and configuring its use in grind.conf.

Bruce


On Nov 3, 2010, at 5:04 PM, Baishan Peng wrote:

> Hi,
>
> I am stuck with my newly installed shakmap-3.5,
>
> 1. the first error message when I tried to run shake is:
>
> [shake at hutl13447 bin]$ ./shake -event 3366146
> transfer: ----- Starting Transfer at 11/04/2010 09:40:14 -----
> transfer: ----- Transfer finished at 11/04/2010 09:40:14 -----
> grind: Predictive map: no data in '/home/shake/shakemap-3.5/bin/../ 
> data/3366146/input'
> cp: cannot stat `/tmp.topo2grd.5251.grd': No such file or directory
> grdsample: Could not find file [/tmp.topo2grd.5251.grd.1]
> /home/shake/shakemap-3.5/bin/../bin/topo2grd 3366146  
> 169.174350/175.174350/-45.727310/-41.377310 0.025/0.025  
> regime=activeUsing . Final resolution fixed at 0.025.
> Done with pasting.
> FINAL BOUNDS: -R///
> cp /tmp.topo2grd.5251.grd /tmp.topo2grd.5251.grd.1
> /usr/bin//grdsample /tmp.topo2grd.5251.grd.1 -G/tmp.topo2grd. 
> 5251.grd -I0.025
> FINAL BOUNDS: -R///
> Could not find /tmp.topo2grd.5251.grd, aborting.
> grind: Error: Error in topo/sitecorr map creation, exiting.
> shake: Error in grind: 256
>
> It looks like it has something to do with  grind.conf, so I changed  
> the topobin value in the config file to:
> topobin : /home/shake/shakemap-3.5/bin/topo2grd 3366146  
> 165.0/180.0/-47.0/-34.0 regime=stable
>
>
> 2. and then I got the following error:
> [shake at hutl13447 bin]$ ./shake -event 3366146
> transfer: ----- Starting Transfer at 11/04/2010 11:07:53 -----
> transfer: ----- Transfer finished at 11/04/2010 11:07:53 -----
> grind: Predictive map: no data in '/home/shake/shakemap-3.5/bin/../ 
> data/3366146/input'
> Invalid res regime=stable at /home/shake/shakemap-3.5/bin/topo2grd  
> line 349.
> grind: Error: Error in topo/sitecorr map creation, exiting.
> /home/shake/shakemap-3.5/bin/topo2grd 3366146  
> 165.0/180.0/-47.0/-34.0 regime=stableshake: Error in grind: 256
>
> It complains about the last statment 'regime=stable' which is just  
> the same as that in the documentation(seems it expects a value for  
> resolution).
>
> Am I doing anything wrong with the configuration, or the program  
> itself?
>
> Please help,
>
> Thanks,
>
> Baishan
>
>
>
> Notice: This email and any attachments are confidential. If received  
> in error please destroy and immediately notify us. Do not copy or  
> disclose the contents.
>
>
> _______________________________________________
> 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/20101104/cd9123ec/attachment.html>