[Shake-dev] Help

Thu Nov 8 23:18:58 UTC 2012

Hello Lily,

I have been working with shakemap and I would like to help you.

The stdout message means there are no topo data for your region. The topo data is grd file and ShakeMap uses this topographic information to obtain the Vs30 (shear waves velocity beneath 30 m earth surface) and also for drawing maps taking into account the topography.

Now, where do you find a topographic data for your region??...well, you can 
download with a resolution of 30 seconds of arc (a grid with a 
resolution of ~1x1 km) in this site:

http://www1.gsi.go.jp/geowww/globalmap-gsi/gtopo30/gtopo30.html

select
 the area where is Chile located. The compress file, that You'll 
download, contains a DEM file, in your case is called W100S10.DEM. You 
have to use the GMT program grdraster in order to convert this DEM file 
to grd one. Please follow the next steps to convert the DEM file to grd:

============================================================================================================
1. create a temporal directory and put into this the DEM file.

2. Inside your temp directory create a file called grdraster.info with this line and save it:

    1 "SRTM topography data - Chile region"     "m"    -R-100/-60/-60/-10  -I0.5m   P i 1   0  -9999  W100S10.DEM   B

3.
 With the grdraster.info file and the DEM file for you region inside the
 temporal file that you created, you can obtain the grd file with this 
command:

[billy at myFedora17:~ ]$ grdraster 1   -R-100/-60/-60/-10  -I0.5m -Gw100s10.grd -V

4.
 After executing the last command please verify the elevation values: 
this information is in last line that you'll read in the std output. If 
you see unrealistic values in the elevation, please delete the B value 
in the line that contains the grdraster.info file.

5. In order to
 verify if your grd file was created properly please execute the next 
command that will draw a map taking into account the elevation:

[billy at myFedora17:~ ]$ grdimage w100s10.grd -R-100/-60/-60/-10 -C/usr/local/share/cpt/GMT_relief.cpt -JM6.5i -P -V > test.ps
[billy at myFedora17:~ ]$ gs test.ps
==============================================================================================================

If all is ok, you will have to set up some information in a program and a config file.

In the first case where Vs30 is calculated based on the topographic data, the program called topo2grd (which is located in $SHAKE_HOME/bin) needs to know where the grd files are to do this job. I suggest you to edit the program topo2grd in the line:

my $topo_dir = 'PATH_of_your_grd_files';

It is important to tell you that the path the directory that contains the grd files must be named DATA. In my case I have the data in this PATH:

/opt/gtopo30/DATA/

and in my topo2grd program I set this path:

my $topo_dir = '/opt/gtopo';

Note that I didn't set the complete path name because this program uses the path and adds a directory called DATA.

Now, if you want the topography to be drawn in your output image files, you will have to set the path and file name of your grd file in the maping.conf file which is located in the $SHAKE_HOME/config directory. There you'll have to edit this file and add this line (take a look to this file to understand more about it is needed to draw maps):

map_topo_hires = PATH_AND_FILE_NAME;

In your case you can add this information (here is necessary the complete path and grd file name):

map_topo_hires = /opt/gtopo30/DATA/w100s10.grd

regardering your last question, I think it is not a hard task to create the XML files, following the information that is in 3.6 section of the ShakeMap SoftwareGuide, which is located in the $SHAKE_HOME/doc directory as well as the example XML file that comes with ShakeMap as example date. Some realtime softwares such as Earthworm or SeisComP3 can produce automatically the XML file, you will only need to copy directory that is create by those programs in the directory $SHAKE_HOME/data. Note that the name of the directory that contains the XML files must have the name of the ID event. I don't know a parser....if you need help with Earthworm or SeisComP3 please let me know it, I can give you a hand.

I hope this will help you.

Kind regards,

Billy Edson Burgoa Rosso

________________________________
 De: "shake-dev-request at geohazards.usgs.gov" <shake-dev-request at geohazards.usgs.gov>
Para: shake-dev at geohazards.usgs.gov 
Enviado: Jueves, 8 de noviembre, 2012 8:00 A.M.
Asunto: Shake-dev Digest, Vol 54, Issue 1

Send Shake-dev mailing list submissions to
    shake-dev at geohazards.usgs.gov

To subscribe or unsubscribe via the World Wide Web, visit
    https://geohazards.usgs.gov/mailman/listinfo/shake-dev
or, via email,
 send a message with subject or body 'help' to
    shake-dev-request at geohazards.usgs.gov

You can reach the person managing the list at
    shake-dev-owner at geohazards.usgs.gov

When
 replying, please edit your Subject line so it is more specific
than "Re: Contents of Shake-dev digest..."

Today's Topics:

   1. Help (Lily Seidman)

----------------------------------------------------------------------

Message: 1
Date: Wed, 7 Nov 2012 16:51:26 -0300
From: Lily Seidman <leseidman at gmail.com>
To: shake-dev at geohazards.usgs.gov
Subject: [Shake-dev] Help
Message-ID:
    <CAOTEYN2i5v5ZuFcnGcc78xmL_tbSW9s9MRbdKmEFSefLPU8Crg at mail.gmail.com>
Content-Type: text/plain; charset="iso-8859-1"

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20121107/d8c33d3c/attachment-0001.html>

------------------------------

_______________________________________________
Shake-dev mailing list
Shake-dev at geohazards.usgs.gov
https://geohazards.usgs.gov/mailman/listinfo/shake-dev

End of Shake-dev Digest, Vol 54, Issue 1
****************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20121108/7b99c1f2/attachment-0001.html>