<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi Baishan,<div><br></div><div><span class="Apple-style-span" style="font-family: monospace; ">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.</span></div><div><font class="Apple-style-span" face="monospace"><br></font></div><div><font class="Apple-style-span" face="monospace">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:</font></div><div><font class="Apple-style-span" face="monospace"><br></font></div><div><font class="Apple-style-span" face="monospace"><a href="http://www.dgadv.com/srtm30/">http://www.dgadv.com/srtm30/</a></font></div><div><font class="Apple-style-span" face="monospace"><br></font></div><div><font class="Apple-style-span" face="monospace">Or other places by searching on SRTM30.</font></div><div><font class="Apple-style-span" face="monospace"><br></font></div><div><font class="Apple-style-span" face="monospace">We use a command like:</font></div><div><font class="Apple-style-span" face="monospace"><br></font></div><div><font class="Apple-style-span" face="monospace">my $command = "$gmtdir/xyz2grd $file -G$section.grd -R$lon/$lat -I30c -Ddeg/deg/m/1/0 -L -F -ZTLhw -N-9999 -V";</font></div><div><font class="Apple-style-span" face="monospace"><br></font></div><div><font class="Apple-style-span" face="monospace">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.</font></div><div><font class="Apple-style-span" face="monospace"><br></font></div><div><font class="Apple-style-span" face="monospace">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.</font></div><div><font class="Apple-style-span" face="monospace"><br></font></div><div><font class="Apple-style-span" face="monospace">Bruce</font></div><div><font class="Apple-style-span" face="monospace"><br></font></div><div><font class="Apple-style-span" face="monospace"><br></font><div><div>On Nov 3, 2010, at 5:04 PM, Baishan Peng wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"> <div> Hi,<br> <br> I am stuck with my newly installed shakmap-3.5, <br> <br> 1. the first error message when I tried to run shake is:<br> <br> [shake@hutl13447 bin]$ ./shake -event 3366146<br> transfer: ----- Starting Transfer at 11/04/2010 09:40:14 -----<br> transfer: ----- Transfer finished at 11/04/2010 09:40:14 -----<br> grind: Predictive map: no data in '/home/shake/shakemap-3.5/bin/../data/3366146/input'<br> cp: cannot stat `/tmp.topo2grd.5251.grd': No such file or directory<br> grdsample: Could not find file [/tmp.topo2grd.5251.grd.1]<br> /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.<br> Done with pasting.<br> FINAL BOUNDS: -R///<br> cp /tmp.topo2grd.5251.grd /tmp.topo2grd.5251.grd.1<br> /usr/bin//grdsample /tmp.topo2grd.5251.grd.1 -G/tmp.topo2grd.5251.grd -I0.025<br> FINAL BOUNDS: -R///<br> Could not find /tmp.topo2grd.5251.grd, aborting.<br> grind: Error: Error in topo/sitecorr map creation, exiting.<br> shake: Error in grind: 256<br> <br> It looks like it has something to do with grind.conf, so I changed the topobin value in the config file to:<br> <font color="#0000ff">topobin : /home/shake/shakemap-3.5/bin/topo2grd 3366146 165.0/180.0/-47.0/-34.0 regime=stable</font><br> <br> <br> 2. and then I got the following error:<br> [shake@hutl13447 bin]$ ./shake -event 3366146<br> transfer: ----- Starting Transfer at 11/04/2010 11:07:53 -----<br> transfer: ----- Transfer finished at 11/04/2010 11:07:53 -----<br> grind: Predictive map: no data in '/home/shake/shakemap-3.5/bin/../data/3366146/input'<br> <font color="#800000">Invalid res regime=stable</font> at /home/shake/shakemap-3.5/bin/topo2grd line 349.<br> grind: Error: Error in topo/sitecorr map creation, exiting.<br> /home/shake/shakemap-3.5/bin/topo2grd 3366146 165.0/180.0/-47.0/-34.0 regime=stableshake: Error in grind: 256<br> <br> 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).<br> <br> Am I doing anything wrong with the configuration, or the program itself?<br> <br> Please help,<br> <br> Thanks,<br> <br> Baishan<br> <br> <br> <br><p><span style="font-family:'Arial';font-size:8pt; color:#000000;">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.</span></p><div><span style="font-family:'Arial';font-size:8pt; color:#000000;"> </span><br class="webkit-block-placeholder"></div></div> _______________________________________________<br>Shake-dev mailing list<br><a href="mailto:Shake-dev@geohazards.usgs.gov">Shake-dev@geohazards.usgs.gov</a><br>https://geohazards.usgs.gov/mailman/listinfo/shake-dev<br></blockquote></div><br></div></body></html>