[Shake-dev] Help

[Bruce Worden]

Hi Lily,

Yes, if you want topography on your maps, you need to install and configure a GMT grid of topography. A Google search will lead you to various sources (use the terms "GMT" and "topography" in your search). If you can't find something that works for you, many people use the SRTM30 Digital Elevation Model (DEM). Attached is a program that will convert the DEM files to GMT .grd files. You can then use grdcut and grdpaste to stitch together a grid file that covers your region of interest. 

Also, how are you providing your Vs30 grid? If you have a grid, great, but if not, you can use the USGS Vs30 server: 

	http://earthquake.usgs.gov/hazards/apps/vs30/

Note that the GMT grid you get may be have pixel-node registration, and ShakeMap requires gridline-node registration. You can fix the file using grdsample:

	grdsample your_vs30_grid.grd -Gnew_file_name.grd -T

Make sure your topography grid covers at least the area of your Vs30 grid.

Finally, you can combine both of these operations into one step (sort of) by getting the SRTM30 data, converting it to .grd with the attached program, configuring include/macros by setting DEMDIR to point to your SRTM30 directory, and configuring grind.conf so that "qtm_file" is unspecified. This will make a Vs30 grid on the fly. You will also then want to run the program "mapping" with the option -gsm, which will use the DEM data for topography.

Attached is dem2grd.pl. You will need to edit it and set a couple of paths at the top of the file to point to the right places.

Hope this helps,
Bruce





On Nov 7, 2012, at 11:51 AM, Lily Seidman wrote:

> Dear ShakeMap Users,
> 
> I am currently a new user of ShakeMap and am trying to customize ShakeMap (after having downloaded it from Fedora and Virtual Machine) for Chile.  I have created an event.xml file in which to test my customizations on.  When I run "shake" for the event, I get lines that say "no topo data found" or "can't find topo data."  Is this because the current topo data is for California and not Chile?  If so, what do I do to customize the topo data?  Also, do I need to do the same (i.e. customize) with the basindepths.xy file so that it is for Chile and not the LA basin?
> 
> Also, I have a few questions about turning the Chilean services data files into xml files.  Do you suggest I use a XML parser and, if so, which do you suggest?
> 
> Thank you so much in advance for your help and suggestions.
> 
> Saludos,
> Lily
> 
> -- 
> Lily Seidman (Smith '11)
> Fulbright Fellow
> LESeidman at gmail.com
> 
> _______________________________________________
> 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/20121108/7ece48c0/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dem2grd.pl
Type: text/x-perl-script
Size: 4283 bytes
Desc: not available
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20121108/7ece48c0/attachment.bin>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20121108/7ece48c0/attachment-0001.html>