[Shake-dev] Help

Sat Nov 10 19:21:57 UTC 2012

Thank you so much, Mr. Worden and Mr. Burgoa, for the great information!  I
will be sure to try your suggestions out next week and will be in touch on
how it goes!

Thank you again and hope you have a great weekend!

Saludos,
Lily

On Thu, Nov 8, 2012 at 8:18 PM, Billy edson Burgoa Rosso <
billyburgoa at yahoo.com> wrote:

> 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
> ****************************************
>
>
>
> _______________________________________________
> Shake-dev mailing list
> Shake-dev at geohazards.usgs.gov
> https://geohazards.usgs.gov/mailman/listinfo/shake-dev
>
>

-- 
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/20121110/8aa19f23/attachment.html>