<HTML><BODY style="word-wrap: break-word; -khtml-nbsp-mode: space; -khtml-line-break: after-white-space; ">Hi Zehra,<DIV><BR class="khtml-block-placeholder"></DIV><DIV>At its simplest, a fault file is just pairs of points separated by spaces or commas:</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>latitude1, longitude1</DIV><DIV>latitude2, longitude2</DIV><DIV>...</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>orÊ</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>latitude1 longitude1</DIV><DIV>latitude2 longitude2</DIV><DIV>...</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>A third value may be added to each line for depth.</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>For the dist_rjb() function (i.e., Joyner-Boore distance) each point is considered to be connected to the next one in the file.Ê Depth is ignored as JB projects the points to the surface anyway.Ê If the first and last point are the same, the fault is treated as a polygon (i.e., the fault plane projected to the surface), with all points inside the polygon at zero distance. Ê</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>dist_hypo() just computes the hypocentral distance (i.e., the distance to a (subsurface) point only).</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>dist_rrupt() computes the distance to a multi-planar 3-D fault surface.Ê It is complicated to explain, though; see the Distance.pm module for this function's documentation and fault definition syntax.</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>If you need some other distance measure, the functions in Distance.pm should provide you with helpful examples.</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>There is currently no function that allows for separate, disconnected fault segments, though there really should be.Ê It wouldn't be that hard to modify Fault.pm and dist_rjb() (in Distance.pm) to handle that situation.Ê I guess nobody has felt the need for it yet.</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>I home this helps.</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>Bruce</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV><BR><DIV><DIV>On Mar 6, 2007, at 12:25 AM, zehra cagnan wrote:</DIV><BR class="Apple-interchange-newline"><BLOCKQUOTE type="cite">Dear Dr Worden,<BR>I started to prepare Shakemap input files for past major earthquakes in Turkey.Especially I am interested with obtaining the instrumental intensity distribution for the Kocaeli earthquake. For cases like Kocaeli earthquake I will have to include a _fault.txt file. I am not clear about this file's format and whether I can include more than one segment for strike slip type fault mechanisms or not. I was wondering if you could send me an example _fault.txt fõle for this purpose. <BR>Thank you very much for all your help.<BR>Best Regards,<BR>Zehra Cagnan<BR><BR><BR><DIV><SPAN class="gmail_quote">On 2/16/07, <B class="gmail_sendername">Bruce Worden</B> <<A href="mailto:cbworden@usgs.gov">cbworden@usgs.gov </A>> wrote:</SPAN><BLOCKQUOTE class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><DIV style="">Hi Zehra,<DIV><BR></DIV><DIV>I'm glad to hear you got it working -- most of the issues with GMT are a result of the differences in the way the headers behave between Intel and Sun hardware.Ê You may have hit on a combination of GMT and ShakeMap that just happens to work.Ê I would recommend upgrading to ShakeMap V3.2, but that can wait until the official version is released (hopefully next week -- the code is done, I'm just waiting for a new distribution site).</DIV><SPAN class="sg"><DIV><BR></DIV><DIV>Bruce</DIV></SPAN><DIV> <SPAN class="e" id="q_110cd133ecb2c0d8_2"><DIV><BR></DIV><DIV><DIV><DIV>On Feb 16, 2007, at 2:00 PM, zehra cagnan wrote:</DIV><BR><BLOCKQUOTE type="cite"><DIV>Hi Bruce</DIV> <DIV>Thanks a lot for all your help. I solved the problems by installing GMT 4.1 instead of GMT 3.4. The grind, mapping, genex modules are working fine under Linux OS. </DIV> <DIV>Best regards</DIV> <DIV>zehra</DIV> <DIV><BR>Ê</DIV> <DIV><SPAN class="gmail_quote">On 2/16/07, <B class="gmail_sendername"> Bruce Worden</B> <<A href="mailto:cbworden@usgs.gov" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">cbworden@usgs.gov</A>> wrote:</SPAN> <BLOCKQUOTE class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0px 0px 0px 0.8ex; padding-left: 1ex;"> <DIV style="">Hmmm, that makes no sense to me.Ê grdinfo can read the header, but another program can't seek past it.Ê Can you run grind and mapping again, this time using the -verbose flag with each program (you'll have to run them individually from the command line since "shake" won't pass the verbose flag to its sub-programs -- use the -dryrun flag to shake to get a listing of how it would run the sub-programs, then just add -verbose to each in turn). <DIV><BR>Ê</DIV> <DIV>Also send me "mmi.grd" so I can see what is going on inside it. <DIV><BR>Ê</DIV> <DIV>Bruce</DIV> <DIV><SPAN> <DIV><BR>Ê</DIV> <DIV><BR> <DIV> <DIV>On Feb 16, 2007, at 12:40 AM, zehra cagnan wrote: </DIV><BR> <BLOCKQUOTE type="cite"> <DIV><FONT size="2"><P>Hi Bruce</P><P>this is the output of grdinfo.</P><P>thanks</P><P>zehra</P><P>/home/ozden/Desktop/SHAKEMAP/shake_home/data/9583161/output/mmi.grd: Title:</P><P>/home/ozden/Desktop/SHAKEMAP/shake_home/data/9583161/output/mmi.grd: Command: /h ome/ozden/Desktop/SHAKEMAP/Kurulum/GMT/bin/nearneighbor - R-120.150000/-117.65000 0/33.729666/35.396333 -G/home/ozden/Desktop/SHAKEMAP/shake_home/bin/../data/9583 161/output/mmi.grd=1 - I0.0166666666666667/0.0166666666666667 -N1 -S0.01666666666 66667</P><P>/home/ozden/Desktop/SHAKEMAP/shake_home/data/9583161/output/mmi.grd: Remark: </P><P>/home/ozden/Desktop/SHAKEMAP/shake_home/data/9583161/output/mmi.grd: Normal node registration used</P><P>/home/ozden/Desktop/SHAKEMAP/shake_home/data/9583161/output/mmi.grd: grdfile for mat # 1</P><P>/home/ozden/Desktop/SHAKEMAP/shake_home/data/9583161/output/mmi.grd: x_min: -120 .15 x_max: - 117.65 x_inc: 0.0166667 units: user_x_unit nx: 151</P><P>/home/ozden/Desktop/SHAKEMAP/shake_home/data/9583161/output/mmi.grd: y_min: 33.7 297 y_max: 35.3963 y_inc: 0.0166667 units: user_y_unit ny: 101</P><P>/home/ozden/Desktop/SHAKEMAP/shake_home/data/9583161/output/mmi.grd: z_min: 1 z_ max: 4.46 units: user_z_unit</P><P>/home/ozden/Desktop/SHAKEMAP/shake_home/data/9583161/output/mmi.grd: scale_facto r: 1 add_offset: 0</P></FONT></DIV> <DIV>Ê</DIV> <DIV>Ê</DIV> <DIV>Ê</DIV> <DIV>Ê</DIV> <DIV><BR><BR>Ê</DIV> <DIV> <SPAN class="gmail_quote">On 2/16/07, <B class="gmail_sendername">Bruce Worden</B> <<A href="mailto:cbworden@usgs.gov" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">cbworden@usgs.gov</A>> wrote: </SPAN> <BLOCKQUOTE class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0px 0px 0px 0.8ex; padding-left: 1ex;"> <DIV style=""> <DIV><BR>Ê</DIV>What did the output of grdinfo look like? <DIV><BR>Ê </DIV> <DIV>Bruce</DIV> <DIV><BR> <DIV> <DIV><SPAN> <DIV>On Feb 15, 2007, at 2:13 PM, zehra cagnan wrote:</DIV><BR></SPAN></DIV> <BLOCKQUOTE type="cite"> <DIV><SPAN> <DIV>Hi Bruce,</DIV> <DIV>I checked the macros and ShakeConfig.pm files. Architecture is set to CISC.</DIV> <DIV>Ê</DIV> <DIV>mmi.grd,pga.grd,pgv.grd they are each around 60 KB. On the other hand I have approximately 10000 rows in the grid.xml file .</DIV> <DIV>zehra</DIV> <DIV><BR><BR> Ê</DIV> <DIV><SPAN class="gmail_quote">On 2/15/07, <B class="gmail_sendername">Bruce Worden</B> <<A href="mailto:cbworden@usgs.gov" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">cbworden@usgs.gov </A>> wrote: </SPAN> <BLOCKQUOTE class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0px 0px 0px 0.8ex; padding-left: 1ex;"> <DIV style="">Hi Zehra, <DIV><BR>Ê</DIV> <DIV>On your second question -- there are some issues in getting to the lists from outside the gps domain.Ê We are going to move shake-dev to a USGS site fairly soon, so you should have easier access (as well as searchable archives).Ê Have you tried this URL: </DIV> <DIV><BR>Ê</DIV> <DIV><A href="http://scsnmail.gps.caltech.edu/mailman/listinfo/shake-dev" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">http://scsnmail.gps.caltech.edu/mailman/listinfo/shake-dev </A></DIV> <DIV><BR>Ê</DIV> <DIV>I don't know if that will work better for you.Ê I am inside the domain, so I see it all just fine.</DIV> <DIV><BR>Ê</DIV> <DIV>On your primary problem, your first error is:</DIV><SPAN> <DIV><BR>Ê</DIV> <DIV style="margin: 0px;">GMT Fatal Error: Error seeking past header file /home/ozden/Desktop/SHAKEMAP/shake_home/bin/../data/9583161/output/mmi.grd!</DIV> <DIV><BR>Ê</DIV></SPAN> <DIV>Which tells me that the .grd files are probably not being created properly.Ê Have a look at that file and see how big it is.Ê You can also use 'grdinfo' to examine the contents -- but remember to append '=1' to the filename so grdinfo can read it ( i.e., 'grdinfo mmi.grd=1').Ê If it gags, 'grind' isn't writing the file correctly.Ê If nx * ny * 4 isn't within a couple of dozen bytes of the file size, there is some other problem.</DIV> <DIV><BR> Ê</DIV> <DIV>Check your 'include/macros' file -- your 'ARCHITECTURE' macro should be set to 'CISC'.Ê If it isn't, you'll need to change it and run "make all" in the shakemap home directory again.Ê Then make sure that the '$architecture' variable in lib/ShakeConfig.pm is set to the proper thing.Ê If it isn't, delete the file and run 'make all' again. </DIV><SPAN> <DIV><BR>Ê</DIV> <DIV>Bruce</DIV></SPAN> <DIV><SPAN> <DIV><BR>Ê</DIV> <DIV><BR>Ê</DIV> <DIV> <DIV> <DIV>On Feb 15, 2007, at 5:11 AM, zehra cagnan wrote:</DIV><BR> <BLOCKQUOTE type="cite">Hi Bruce<BR>I followed the shakemap manual and installed GMT version 3.4 thinking that installing GMT version 4.1 would create problems.<BR>Thank you very much for the help you provided. I solved the problem below. I uncommented certain lines within mapping as they were applicable to verion 4.1 of GMT.<BR><BR>The grind module seems to work with no problems but I am still experiencing problems with the mapping module. I am getting the error message below. Somehow mapping does not create the tmpmmi.grd file. I am seeing the mmi.grd file under the data/mapping directory but tmpmmi.grd file is not being created.<BR><BR>Some of the files I need to complete regional configuration is not ready yet so I am using california data and trying to do a test run just to understand the program better. <BR>I commented out certain parts within shake.conf so that only grind,tag,mapping,genex modules will be called in the order I wrote above.<BR><BR>My second question is regarding the shakemap list. I am trying to access the website of the list thinking that I might be able access the list archieves via there (may be some of the questions I am asking now were answered there before). I am using my email address as my user id and the password you sent me as my password but could not access the web site. is my user id something different than <A href="mailto:zcagnan@gmail.com" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">zcagnan@gmail.com</A><BR><BR>thanks<BR>zehra<BR><BR><BR><BR>grind: Type: pga; bias 0.815; absolute deviation was 0.019171 is 0.017089<BR>grind: Type: pgv; bias 0.775; absolute deviation was 0.006004 is 0.005444<BR>grind: 1 stations flagged out this iteration<BR>grind: Type: pga; bias 0.815; absolute deviation was 0.019171 is 0.017089<BR>grind: Type: pgv; bias 0.775; absolute deviation was 0.006004 is 0.005444 <BR>grind: 0 stations flagged out this iteration<BR>GMT Fatal Error: Error seeking past header file /home/ozden/Desktop/SHAKEMAP/shake_home/bin/../data/9583161/output/mmi.grd! <BR>grdview: No such file or directory [/home/ozden/Desktop/SHAKEMAP/shake_home/bin/../data/9583161/mapping/tmp_mmi.grd] <BR>mapping: Error running command:<BR>/home/ozden/Desktop/SHAKEMAP/Kurulum/GMT/bin/grdview /home/ozden/Desktop/SHAKEMAP/shake_home/bin/../data/9583161/mapping/topo.grd -P - R-120.15/-117.65/33.729666/35.396333/10/3000 -JM6.5 -JZ1 -G/home/ozden/Desktop/SHAKEMAP/shake_home/bin/../data/9583161/mapping/tmp_mmi.grd -C/home/ozden/Desktop/SHAKEMAP/shake_home/bin/../lib/mapping/Ii.cpt -I/home/ozden/Desktop/SHAKEMAP/shake_home/bin/../data/9583161/mapping/inten.grd -Qi72 -O -K| <BR><BR>shake: Error in mapping: 256<BR><BR><BR><BR><BR><BR> <DIV><SPAN class="gmail_quote">On 2/14/07, <B class="gmail_sendername">Bruce Worden</B> <<A href="mailto:cbworden@usgs.gov" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)"> cbworden@usgs.gov</A>> wrote: </SPAN> <BLOCKQUOTE class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi Zehra,<BR><BR>What version of GMT are you using?ÊÊThere are problems with this <BR>version of ShakeMap and certain versions of GMT.ÊÊThe next version of <BR>ShakeMap (V3.2, which you can get a pre-release version of at the FTP<BR>site) should be more forgiving of GMT's idiosyncrasies.<BR><BR>But in the version you have: <BR>Edit 'mapping' and search for the string "PLOT_DEGREE_FORMAT".ÊÊThere <BR>is some commentary there about fixing the first error, below.<BR><BR>The topography error is odd -- can you verify that you have the <BR>topography files "cal_neva_9c[i].grd" (or whatever topo files you<BR>have configured in mapping.conf) and not just links?ÊÊIf they are<BR>links, you can get the actual files at the FTP site.<BR><BR>Also, I am adding you to the shake-dev mailing list (you should get a <BR>confirmation email shortly -- make sure your spam filter doesn't get <BR>it) so you can ask your questions there, and maybe get more help.<BR><BR>I hope this helps.<BR><BR>Bruce<BR><BR>On Feb 14, 2007, at 2:59 AM, zehra cagnan wrote: <BR><BR>> Okey<BR>> I solved the previous problem. The binary form of the vsgrid.txt<BR>> that we were using was not generated on a linux machine. That was<BR>> the problem.<BR>> Now I am having problems with the GMT. <BR>><BR>> The error message that i am getting is below. Any help you can <BR>> provide?<BR>><BR>> grind: Type: pga; bias 0.815; absolute deviation was 0.019171 is<BR>> 0.017089<BR>> grind: Type: pgv; bias 0.775; absolute deviation was 0.006004 is<BR>> 0.005444<BR>> grind: 1 stations flagged out this iteration <BR>> grind: Type: pga; bias 0.815; absolute deviation was 0.019171 is<BR>> 0.017089<BR>> grind: Type: pgv; bias 0.775; absolute deviation was 0.006004 is<BR>> 0.005444<BR>> grind: 0 stations flagged out this iteration <BR>> gmtset: GMT SYNTAX ERROR in GMT_setparameter:ÊÊUnrecognized keyword<BR>> PLOT_DEGREE_FORMAT<BR>> gmtset: GMT SYNTAX ERROR:ÊÊPLOT_DEGREE_FORMAT given illegal value (D)! <BR>> gmtset:ÊÊ1 conversion errors<BR> > mapping: WARNING: no topography data found; will continue without<BR>> topo.<BR>> GMT Fatal Error: Error seeking past header file /home/ozden/Desktop/<BR>> SHAKEMAP/shake_home/bin/../data/9583161/output/mmi.grd! <BR>> mapping: WARNING: can't find topo data: making mmi without topography<BR>> grdimage: No such file or directory [/home/ozden/Desktop/SHAKEMAP/<BR>> shake_home/bin/../data/9583161/mapping/tmp_mmi.grd]<BR> > mapping: Error running command:<BR>> /home/ozden/Desktop/SHAKEMAP/Kurulum/GMT/bin/grdimage /home/ozden/<BR>> Desktop/SHAKEMAP/shake_home/bin/../data/9583161/mapping/tmp_mmi.grd<BR>> -P -R-120.15/-117.65/33.729666/35.396333/10/3000 - JM6.5 -C/home/<BR>> ozden/Desktop/SHAKEMAP/shake_home/bin/../lib/mapping/Ii.cpt -O -K|<BR>><BR>><BR><BR></BLOCKQUOTE></DIV><BR></BLOCKQUOTE></DIV><BR>Ê</DIV></SPAN></DIV></DIV></BLOCKQUOTE></DIV><BR></SPAN></DIV> <DIV style="margin: 0px;">_______________________________________________</DIV> <DIV style="margin: 0px;">Shake-dev mailing list</DIV> <DIV style="margin: 0px;"><A href="mailto:Shake-dev@scsnmail.gps.caltech.edu" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)"> Shake-dev@scsnmail.gps.caltech.edu</A></DIV> <DIV style="margin: 0px;"><A href="http://unix1.gps.caltech.edu/mailman/listinfo/shake-dev" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">http://unix1.gps.caltech.edu/mailman/listinfo/shake-dev </A> </DIV></BLOCKQUOTE></DIV><BR>Ê</DIV></DIV></BLOCKQUOTE></DIV><BR> <DIV style="margin: 0px;">_______________________________________________</DIV> <DIV style="margin: 0px;">Shake-dev mailing list</DIV> <DIV style="margin: 0px;"> <A href="mailto:Shake-dev@scsnmail.gps.caltech.edu" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">Shake-dev@scsnmail.gps.caltech.edu</A></DIV> <DIV style="margin: 0px;"><A href="http://unix1.gps.caltech.edu/mailman/listinfo/shake-dev" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)"> http://unix1.gps.caltech.edu/mailman/listinfo/shake-dev</A> </DIV></BLOCKQUOTE></DIV><BR>Ê</DIV></SPAN></DIV></DIV></DIV></BLOCKQUOTE></DIV><BR></BLOCKQUOTE></DIV><BR></DIV></SPAN></DIV></DIV></BLOCKQUOTE></DIV><BR></BLOCKQUOTE></DIV><BR></DIV></BODY></HTML>