From huang.zhch at gmail.com Sun Sep 5 23:52:19 2010 From: huang.zhch at gmail.com (Zhouchuan Huang) Date: Mon, 6 Sep 2010 15:52:19 +0900 Subject: [SAC-HELP] how to display distance along each seismogram Message-ID: *Dear everyone, * ** *When plotting using sac, some information such as the network (KNETWK), station (KSTNM), channels (KCMPNM) and time (KZDATE, KZTIME) will be plotted together with each seismogram. * *I wonder if I can add more information, such as distance (DIST) on each seismogram simultaneously.* ** *Thank you very much.* ** *Best wishes,* ** *Zhouchuan* -- Zhouchuan Huang School of Earth Sciences and Engineering Nanjing University 22 Hankou Road Nanjing 210093 P.R.CHINA. Tel: +8613675137130 E-mail: huang.zhch at gmail.com -------------- next part -------------- An HTML attachment was scrubbed... URL: From mplasencia at ogs.trieste.it Mon Sep 6 00:52:34 2010 From: mplasencia at ogs.trieste.it (Milton P. Plasencia Linares) Date: Mon, 06 Sep 2010 09:52:34 +0200 Subject: [SAC-HELP] how to display distance along each seismogram In-Reply-To: References: Message-ID: <20100906095234.168659l2rqq7etvm@webmail.ogs.trieste.it> Hi, you can use the "fileid" command. sac> fileid type list "here all header list that you want" I wait this help you, Cheers, Milton ********************************** Milton P. PLASENCIA LINARES Dipartimento Centro di Ricerche Sismologiche Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS Borgo Grotta Gigante 42/C (34010) Sgonico - TRIESTE - ITALIA Tel: +39-040-2140136 Fax: +39-040-327307 E-mail: mplasencia at ogs.trieste.it ASAIN (Antarctic Seismographic Argentinean Italian Network) ********************************* Quoting "Zhouchuan Huang" : > *Dear everyone, * > ** > *When plotting using sac, some information such as the network (KNETWK), > station (KSTNM), channels (KCMPNM) and time (KZDATE, KZTIME) will be plotted > together with each seismogram. * > *I wonder if I can add more information, such as distance (DIST) on each > seismogram simultaneously.* > ** > *Thank you very much.* > ** > *Best wishes,* > ** > *Zhouchuan* > > -- > Zhouchuan Huang > School of Earth Sciences and Engineering > Nanjing University > 22 Hankou Road > Nanjing 210093 > P.R.CHINA. > Tel: +8613675137130 > E-mail: huang.zhch at gmail.com > > -- > This email was Anti Virus checked by INOGS Antivirus filter. > ---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program. From ndikumeric at yahoo.com Fri Sep 10 05:28:01 2010 From: ndikumeric at yahoo.com (ndoh ndikum) Date: Fri, 10 Sep 2010 05:28:01 -0700 (PDT) Subject: [SAC-HELP] converting sac data!! Message-ID: <181040.10033.qm@web51706.mail.re2.yahoo.com> please, I am new in this forum. But i am a user of SAC for some time now. I seek to know if i could have some software to convert data in sac format (output from SAC) into SEED and AH format. Also, I wish to know if there are softwares which can accept data in sac format as input to produce group velocity maps and inversed structures. hope you help me out. Ndikum. --------------------------------------------------- NDIKUM ERIC NDOH, C/O C.M.F.I. YAOUNDE, P.O. Box 6090, CAMEROON, WEST AFRICA . TEL:(237) 75 17 91 66 -------------- next part -------------- An HTML attachment was scrubbed... URL: From chad at iris.washington.edu Fri Sep 10 08:01:33 2010 From: chad at iris.washington.edu (Chad Trabant) Date: Fri, 10 Sep 2010 08:01:33 -0700 Subject: [SAC-HELP] converting sac data!! In-Reply-To: <181040.10033.qm@web51706.mail.re2.yahoo.com> References: <181040.10033.qm@web51706.mail.re2.yahoo.com> Message-ID: Hello Ndikum, The sac2mseed program will convert SAC data into miniSEED, the source code for the current version can be found here: http://www.iris.edu/pub/programs/converters/ You will need to compile the source code, instructions are included. This converter is free to use but not an officially supported IRIS software package. The miniSEED data will contain the SAC time series, but will not contain any of the metadata in the SAC headers because there is no place to store that information in miniSEED. I do not know of any converters for SAC to full SEED or AH. regards, Chad Trabant On Sep 10, 2010, at 5:28 AM, ndoh ndikum wrote: > please, I am new in this forum. But i am a user of SAC for some time now. > I seek to know if i could have some software to convert data in sac format > (output from SAC) into SEED and AH format. > Also, I wish to know if there are softwares which can accept data in sac > format as input to produce group velocity maps and inversed structures. > hope you help me out. > Ndikum. > > --------------------------------------------------- > NDIKUM ERIC NDOH, > C/O C.M.F.I. YAOUNDE, > P.O. Box 6090, > CAMEROON, > WEST AFRICA . > > TEL:(237) 75 17 91 66 > > > > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help -------------- next part -------------- An HTML attachment was scrubbed... URL: From jattanayake at gmail.com Mon Sep 13 09:07:04 2010 From: jattanayake at gmail.com (Januka Attanayake) Date: Mon, 13 Sep 2010 12:07:04 -0400 Subject: [SAC-HELP] TRANSFER freqlimits Message-ID: Greetings, Could someone please clarify the physical meaning of freqlimits when using TRANSFER function in SAC? I have read the documentation but for me it still is a little fuzzy. Does it mean that when we use certain freqlimits, for example freqlimits 0.01 0.015 100 200, the result we see is a narrow band passed record with corner frequencies of, in this example 0.15Hz and 100Hz? If that were the case and if I wanted to have the broadest bandpass possible, what would be the limits to use? When I use lower freqlimits like 0.001 0.005 the resulting record seems to have a trend and sometimes amplified long period signals. Therefore, I am forced to narrow the freqlimits to 0.01 0.15 (lower end) 100 200 (higher end). I believe if freqlimits is acting as a bandpass I am loosing some lower and higher frequency information. Any feed back is much appreciated! Thank you! Januka - using SAC 101.3b -------------- next part -------------- An HTML attachment was scrubbed... URL: From snoke at vt.edu Mon Sep 13 11:28:43 2010 From: snoke at vt.edu (Arthur Snoke) Date: Mon, 13 Sep 2010 14:28:43 -0400 (EDT) Subject: [SAC-HELP] TRANSFER freqlimits In-Reply-To: References: Message-ID: Januka, I suggest looking at the help file for transfer that was included with 101.4 and is posted on the IRIS SAC Web site at . I think/hope the description is clearer than it was earlier, although the functionality is probably the same. Dealing with your example: freqlimits 0.01 0.015 100 200. Assuming your Nyquist is above 200.0, what you have here is a filter that has unit weight between f2 and f3: 0.015 and 100.0. That is, f2 and f3 are not corners as in a "normal" bandpass. There is a cosine taper from unity to 0 between 100.0 and 200.0 and between 0.015 and 0.01. The response is "zero" for frequencies below 0.01 and above 200.0. I put quotes around the zero, because unless your npts is a power of 2, because (apparently) of the nature of the FFT, the response will not be zero. I hope this helps. If not, let me know. Arthur On Mon, 13 Sep 2010, Januka Attanayake wrote: > Greetings, > Could someone please clarify the physical meaning of freqlimits when > using TRANSFER function in SAC? I have read the documentation but for me > it still is a little fuzzy. Does it mean that when we use certain > freqlimits, for example freqlimits 0.01 0.015 100 200, the result we see > is a narrow band passed record with corner frequencies of, in this > example 0.15Hz and 100Hz? If that were the case and if I wanted to have > the broadest bandpass possible, what would be the limits to use? When I > use lower freqlimits like 0.001 0.005 the resulting record seems to have > a trend and sometimes amplified long period signals. Therefore, I am > forced to narrow the freqlimits to 0.01 0.15 (lower end) 100 200 (higher > end). I believe if freqlimits is acting as a bandpass?I am loosing some > lower and higher frequency information. Any feed back is much > appreciated! > > Thank you! > > Januka - using SAC 101.3b > > From jattanayake at gmail.com Tue Sep 14 07:06:40 2010 From: jattanayake at gmail.com (Januka Attanayake) Date: Tue, 14 Sep 2010 10:06:40 -0400 Subject: [SAC-HELP] TRANSFER freqlimits In-Reply-To: References: Message-ID: Thanks Arthur, the process is clear now. However, one of the questions I raised still remains. That is how would you determine what low/high pass frequencies (f1 f2 f3 f4) would you use? Documentation says "for a Nyquist frequency of 0.5, set f3=10. and f4=20". When processing seismograms in bulk with different instruments (1) you wouldn't know what the Nyquist freq. is (2) If the answer to (1) is " yes you don't" how would you determine f1,f2,f3,f4 for a given Nyquist freq on general terms. Thanks again ! Januka. On Mon, Sep 13, 2010 at 2:28 PM, Arthur Snoke wrote: > Januka, > > I suggest looking at the help file for transfer that was included with > 101.4 and is posted on the IRIS SAC Web site at < > http://www.iris.edu/software/sac/commands/transfer.html>. I think/hope > the description is clearer than it was earlier, although the functionality > is probably the same. > > Dealing with your example: freqlimits 0.01 0.015 100 200. Assuming your > Nyquist is above 200.0, what you have here is a filter that has unit weight > between f2 and f3: 0.015 and 100.0. That is, f2 and f3 are not corners as > in a "normal" bandpass. There is a cosine taper from unity to 0 between > 100.0 and 200.0 and between 0.015 and 0.01. The response is "zero" for > frequencies below 0.01 and above 200.0. I put quotes around the zero, > because unless your npts is a power of 2, because (apparently) of the nature > of the FFT, the response will not be zero. > > I hope this helps. If not, let me know. > > Arthur > > > On Mon, 13 Sep 2010, Januka Attanayake wrote: > > Greetings, >> Could someone please clarify the physical meaning of freqlimits when using >> TRANSFER function in SAC? I have read the documentation but for me it still >> is a little fuzzy. Does it mean that when we use certain freqlimits, for >> example freqlimits 0.01 0.015 100 200, the result we see is a narrow band >> passed record with corner frequencies of, in this example 0.15Hz and 100Hz? >> If that were the case and if I wanted to have the broadest bandpass >> possible, what would be the limits to use? When I use lower freqlimits like >> 0.001 0.005 the resulting record seems to have a trend and sometimes >> amplified long period signals. Therefore, I am forced to narrow the >> freqlimits to 0.01 0.15 (lower end) 100 200 (higher end). I believe if >> freqlimits is acting as a bandpass I am loosing some lower and higher >> frequency information. Any feed back is much appreciated! >> >> Thank you! >> >> Januka - using SAC 101.3b >> >> -- "Nothing can start to exist nor can cease to exist, only transformation is possible" Januka Attanayake Solid Earth Geophysics & Economics University of Connecticut Beach Hall U-2045 354,Mansfield Rd; Storrs.CT 06269 Tel: +1 860 486 0475 (Office) +1 860 486 3914 (Lab - P405) URL: http://sites.google.com/site/janukaattanayake/ -------------- next part -------------- An HTML attachment was scrubbed... URL: From jattanayake at gmail.com Tue Sep 14 07:06:40 2010 From: jattanayake at gmail.com (Januka Attanayake) Date: Tue, 14 Sep 2010 10:06:40 -0400 Subject: [SAC-HELP] TRANSFER freqlimits In-Reply-To: References: Message-ID: Thanks Arthur, the process is clear now. However, one of the questions I raised still remains. That is how would you determine what low/high pass frequencies (f1 f2 f3 f4) would you use? Documentation says "for a Nyquist frequency of 0.5, set f3=10. and f4=20". When processing seismograms in bulk with different instruments (1) you wouldn't know what the Nyquist freq. is (2) If the answer to (1) is " yes you don't" how would you determine f1,f2,f3,f4 for a given Nyquist freq on general terms. Thanks again ! Januka. On Mon, Sep 13, 2010 at 2:28 PM, Arthur Snoke wrote: > Januka, > > I suggest looking at the help file for transfer that was included with > 101.4 and is posted on the IRIS SAC Web site at < > http://www.iris.edu/software/sac/commands/transfer.html>. I think/hope > the description is clearer than it was earlier, although the functionality > is probably the same. > > Dealing with your example: freqlimits 0.01 0.015 100 200. Assuming your > Nyquist is above 200.0, what you have here is a filter that has unit weight > between f2 and f3: 0.015 and 100.0. That is, f2 and f3 are not corners as > in a "normal" bandpass. There is a cosine taper from unity to 0 between > 100.0 and 200.0 and between 0.015 and 0.01. The response is "zero" for > frequencies below 0.01 and above 200.0. I put quotes around the zero, > because unless your npts is a power of 2, because (apparently) of the nature > of the FFT, the response will not be zero. > > I hope this helps. If not, let me know. > > Arthur > > > On Mon, 13 Sep 2010, Januka Attanayake wrote: > > Greetings, >> Could someone please clarify the physical meaning of freqlimits when using >> TRANSFER function in SAC? I have read the documentation but for me it still >> is a little fuzzy. Does it mean that when we use certain freqlimits, for >> example freqlimits 0.01 0.015 100 200, the result we see is a narrow band >> passed record with corner frequencies of, in this example 0.15Hz and 100Hz? >> If that were the case and if I wanted to have the broadest bandpass >> possible, what would be the limits to use? When I use lower freqlimits like >> 0.001 0.005 the resulting record seems to have a trend and sometimes >> amplified long period signals. Therefore, I am forced to narrow the >> freqlimits to 0.01 0.15 (lower end) 100 200 (higher end). I believe if >> freqlimits is acting as a bandpass I am loosing some lower and higher >> frequency information. Any feed back is much appreciated! >> >> Thank you! >> >> Januka - using SAC 101.3b >> >> -- "Nothing can start to exist nor can cease to exist, only transformation is possible" Januka Attanayake Solid Earth Geophysics & Economics University of Connecticut Beach Hall U-2045 354,Mansfield Rd; Storrs.CT 06269 Tel: +1 860 486 0475 (Office) +1 860 486 3914 (Lab - P405) URL: http://sites.google.com/site/janukaattanayake/ -------------- next part -------------- An HTML attachment was scrubbed... URL: From snoke at vt.edu Tue Sep 14 07:54:46 2010 From: snoke at vt.edu (Arthur Snoke) Date: Tue, 14 Sep 2010 10:54:46 -0400 (EDT) Subject: [SAC-HELP] TRANSFER freqlimits In-Reply-To: References: Message-ID: Comments interspersed. On Tue, 14 Sep 2010, Januka Attanayake wrote: > Thanks Arthur, the process is clear now. However, one of the questions I > raised still remains. That is how would you determine what low/high pass > frequencies (f1 f2 f3 f4) would you use? Documentation says "for a > Nyquist frequency of 0.5, set f3=10. and f4=20". When processing > seismograms in bulk with different instruments (1) you wouldn't know > what the Nyquist freq. is (2) If the answer to (1) is " yes you don't" > how would you determine f1,f2,f3,f4 for a given Nyquist freq on general > terms. In the help file that comes with 101.4 and which can bde accessed from the IRIS SAC manual pages, it says: If you want to do a low-pass filter but have no filtering at low frequencies, one way is to set f1=-2 and f2=-1. If you want to do a high-pass fiter but have no filtering at the high frequencies, for a Nyquist frequency of 0.5, set f3=10. and f4=20. The point of this passage is how do you use freqlimits if you want to have no filtering at one end of the spectrum but you do at the other. The program as written says that if f2 > f1, but f2 is less than df, there is no filtering at low frequencies (a negative frequency is always less than any df.). Similarly, for a Nyquist of 0.5 Hz, if f4 > f3 and f3 is greater than 0.5, there is no filtering at high frequencies (10 is much greater than 0.5). Note that it is important that f2 > f1, or the program wil think this is dealing with the high-frequency end. You ask how one knows the Nyquist if one has many instruments, etc. DF and the Nyquist are related to the time series itself. Here I read in a time series and do an lh: SAC> r N11A.lhz SAC> lh FILE: N11A.lhz - 1 -------------- NPTS = 3101 B = -9.980005e+01 E = 3.000200e+03 IFTYPE = TIME SERIES FILE LEVEN = TRUE DELTA = 1.000000e+00 DEPMIN = -1.928006e+03 DEPMAX = 1.540401e+03 Now I do an fft and a lh SAC> fft DC level after DFT is -724.33 SAC> lh FILE: N11A.lhz - 1 -------------- NPTS = 4096 B = 0.000000e+00 E = 5.000000e-01 IFTYPE = SPECTRAL FILE-AMPL/PHASE LEVEN = TRUE DELTA = 2.441406e-04 DEPMIN = -1.928006e+03 DEPMAX = 1.540401e+03 In the first lh, we see that there are 3101 points with a sampling rate of 1 sps for a total time, then, of (NPTS-1)*DELTA =3100 seconds. After the fft, there are 4096 points because the fft adds zeroes to the time serie to the next power of 2. The lowest frequency is 0, and E is now the Nyquist of 0.5 Hz (0.5/DELTA). Now, DELTA = 2.441406e-04 is DF, the frequency interval. I hope this helps. If you still have questions ... From jattanayake at yahoo.com Tue Sep 14 08:31:44 2010 From: jattanayake at yahoo.com (Januka Attanayake) Date: Tue, 14 Sep 2010 08:31:44 -0700 (PDT) Subject: [SAC-HELP] TRANSFER freqlimits In-Reply-To: Message-ID: <252222.84834.qm@web55508.mail.re4.yahoo.com> GREAT ! Thanks again Arthur, I haven't got a better explanation than the one that you just gave me! May be you should include this in SAC documentation. Serenity isn't freedom from the storm, but peace within the storm Januka Attanayake? Solid Earth Geophysics & Economics The University of Connecticut Beach Hall - U2045 354, Mansfield Rd; Storrs, CT 06269 Tel :? 860 486 0475 (Office) ????????? 860 486 3914 (Lab - P405) URL:?http://sites.google.com/site/janukaattanayake/ --- On Tue, 9/14/10, Arthur Snoke wrote: From: Arthur Snoke Subject: Re: [SAC-HELP] TRANSFER freqlimits To: "Januka Attanayake" Cc: sac-help at iris.washington.edu Date: Tuesday, September 14, 2010, 10:54 AM Comments interspersed. On Tue, 14 Sep 2010, Januka Attanayake wrote: > Thanks Arthur, the process is clear now. However, one of the questions I raised still remains. That is how would you determine what low/high pass frequencies (f1 f2 f3 f4) would you use? Documentation says "for a Nyquist frequency of 0.5, set f3=10. and f4=20". When processing seismograms in bulk with different instruments (1) you wouldn't know what the Nyquist freq. is (2) If the answer to (1) is " yes you don't" how would you determine f1,f2,f3,f4 for a given Nyquist freq on general terms. In the help file that comes with 101.4 and which can bde accessed from the IRIS SAC manual pages, it says: If you want to do a low-pass filter but have no filtering at low frequencies, one way is to set f1=-2 and f2=-1. If you want to do a high-pass fiter but have no filtering at the high frequencies, for a Nyquist frequency of 0.5, set f3=10. and f4=20. The point of this passage is how do you use freqlimits if you want to have no filtering at one end of the spectrum but you do at the other.? The program as written says that if f2 > f1, but f2 is less than df, there is no filtering at low frequencies (a negative frequency is always less than any df.).? Similarly, for a Nyquist of 0.5 Hz, if f4 > f3 and f3 is greater than 0.5, there is no filtering at high frequencies (10 is much greater than 0.5).? Note that it is important that f2 > f1, or the program wil think this is dealing with the high-frequency end. You ask how one knows the Nyquist if one has many instruments, etc.? DF and the Nyquist are related to the time series itself.? Here I read in a time series and do an lh: SAC> r N11A.lhz SAC> lh ? FILE: N11A.lhz - 1 -------------- ? ? ???NPTS = 3101? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? B = -9.980005e+01 ? ? ? ? ? E = 3.000200e+03? ? ? ? ? ? ? ? ???IFTYPE = TIME SERIES FILE ? ? ? LEVEN = TRUE? ? ? ? ? ? ? ? ? ? ? ? ? ? DELTA = 1.000000e+00 ? ???DEPMIN = -1.928006e+03? ? ? ? ? ? ? ? ? DEPMAX = 1.540401e+03 Now I do an fft and a lh SAC> fft DC level after DFT is -724.33 SAC> lh ? FILE: N11A.lhz - 1 -------------- ? ? ???NPTS = 4096? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? B = 0.000000e+00 ? ? ? ? ? E = 5.000000e-01? ? ? ? ? ? ? ? ???IFTYPE = SPECTRAL FILE-AMPL/PHASE ? ? ? LEVEN = TRUE? ? ? ? ? ? ? ? ? ? ? ? ? ? DELTA = 2.441406e-04 ? ???DEPMIN = -1.928006e+03? ? ? ? ? ? ? ? ? DEPMAX = 1.540401e+03 In the first lh, we see that there are 3101 points with a sampling rate of 1 sps for a total time, then, of (NPTS-1)*DELTA =3100 seconds.? After the fft, there are 4096 points because the fft adds zeroes to the time serie to the next power of 2.? The lowest frequency is 0, and E is now the Nyquist of 0.5 Hz (0.5/DELTA).? Now, DELTA = 2.441406e-04 is DF, the frequency interval. I hope this helps.? If you still have questions ... _______________________________________________ sac-help mailing list sac-help at iris.washington.edu http://www.iris.washington.edu/mailman/listinfo/sac-help -------------- next part -------------- An HTML attachment was scrubbed... URL: From srichandprajapati at yahoo.co.in Thu Sep 23 11:15:31 2010 From: srichandprajapati at yahoo.co.in (srichand prajapati) Date: Thu, 23 Sep 2010 23:45:31 +0530 (IST) Subject: [SAC-HELP] Welcome to the "sac-help" mailing list In-Reply-To: Message-ID: <560967.95111.qm@web94008.mail.in2.yahoo.com> Sir, ? I try to installed the SAC on Cygwin (Cygwin based SAC)?but i am not able to installed it. I already registered on IRIS site to get SAC software. when i run ./configure, its work well but i try to compile by using ?make it show some error in src directory. please help to installed SAC on Cygwin plateform. regards srichand --- On Thu, 23/9/10, sac-help-request at iris.washington.edu wrote: From: sac-help-request at iris.washington.edu Subject: Welcome to the "sac-help" mailing list To: srichandprajapati at yahoo.co.in Date: Thursday, 23 September, 2010, 11:43 PM Welcome to the sac-help at iris.washington.edu mailing list! To post to this list, send your email to: ? sac-help at iris.washington.edu General information about the mailing list is at: ? http://www.iris.washington.edu/mailman/listinfo/sac-help If you ever want to unsubscribe or change your options (eg, switch to or from digest mode, change your password, etc.), visit your subscription page at: ? http://www.iris.washington.edu/mailman/options/sac-help/srichandprajapati%40yahoo.co.in You can also make such adjustments via email by sending a message to: ? sac-help-request at iris.washington.edu with the word `help' in the subject or body (don't include the quotes), and you will get back a message with instructions. You must know your password to change your options (including changing the password, itself) or to unsubscribe.? It is: ? shubh Normally, Mailman will remind you of your iris.washington.edu mailing list passwords once every month, although you can disable this if you prefer.? This reminder will also include instructions on how to unsubscribe or change your account options.? There is also a button on your options page that will email your current password to you. -------------- next part -------------- An HTML attachment was scrubbed... URL: From chad at iris.washington.edu Thu Sep 23 12:28:39 2010 From: chad at iris.washington.edu (Chad Trabant) Date: Thu, 23 Sep 2010 12:28:39 -0700 Subject: [SAC-HELP] Welcome to the "sac-help" mailing list In-Reply-To: <560967.95111.qm@web94008.mail.in2.yahoo.com> References: <560967.95111.qm@web94008.mail.in2.yahoo.com> Message-ID: Dear srichand, It would be helpful if you shared the error messages with the list members, can you include them. I personally have no experience with Cygwin but someone on the list might recognize the problem. regards, Chad On Sep 23, 2010, at 11:15 AM, srichand prajapati wrote: > Sir, > I try to installed the SAC on Cygwin (Cygwin based SAC) but i am not able to installed it. I already registered on IRIS site to get SAC software. > when i run ./configure, its work well > but i try to compile by using make > it show some error in src directory. > please help to installed SAC on Cygwin plateform. > regards > srichand > > > --- On Thu, 23/9/10, sac-help-request at iris.washington.edu wrote: > > From: sac-help-request at iris.washington.edu > Subject: Welcome to the "sac-help" mailing list > To: srichandprajapati at yahoo.co.in > Date: Thursday, 23 September, 2010, 11:43 PM > > Welcome to the sac-help at iris.washington.edu mailing list! > > To post to this list, send your email to: > > sac-help at iris.washington.edu > > General information about the mailing list is at: > > http://www.iris.washington.edu/mailman/listinfo/sac-help > > If you ever want to unsubscribe or change your options (eg, switch to > or from digest mode, change your password, etc.), visit your > subscription page at: > > http://www.iris.washington.edu/mailman/options/sac-help/srichandprajapati%40yahoo.co.in > > > You can also make such adjustments via email by sending a message to: > > sac-help-request at iris.washington.edu > > with the word `help' in the subject or body (don't include the > quotes), and you will get back a message with instructions. > > You must know your password to change your options (including changing > the password, itself) or to unsubscribe. It is: > > shubh > > Normally, Mailman will remind you of your iris.washington.edu mailing > list passwords once every month, although you can disable this if you > prefer. This reminder will also include instructions on how to > unsubscribe or change your account options. There is also a button on > your options page that will email your current password to you. > > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help -------------- next part -------------- An HTML attachment was scrubbed... URL: From wikango at yahoo.fr Fri Sep 24 04:02:20 2010 From: wikango at yahoo.fr (Samir L.) Date: Fri, 24 Sep 2010 11:02:20 +0000 (GMT) Subject: [SAC-HELP] Polezero file Message-ID: <503449.45221.qm@web27907.mail.ukl.yahoo.com> Dear all, I work currently with a python script which invokes SAC to process some data. Since my database usually involves mseed files, I used "ms2sac" to convert them to sac files. However, I found out that during the processing steps Polezero files are also needed. Could someone tell me what are these files and how I can get them? Thank you in advance, Cheers, Wik. -------------- next part -------------- An HTML attachment was scrubbed... URL: From sheila at blacknest.gov.uk Fri Sep 24 04:51:03 2010 From: sheila at blacknest.gov.uk (Sheila Peacock) Date: Fri, 24 Sep 2010 12:51:03 +0100 Subject: [SAC-HELP] Polezero file In-Reply-To: <503449.45221.qm@web27907.mail.ukl.yahoo.com> References: <503449.45221.qm@web27907.mail.ukl.yahoo.com> Message-ID: <4C9C90A7.4080905@blacknest.gov.uk> Dear Wik, Have you a "dataless SEED" file accompanying your miniSEED file? If so, then you can extract the poles and zeros from it in SAC format with "rdseed", which you can download from the IRIS website. Before running "rdseed", declare the environment variable ALT_RESPONSE_FILE to be the name of your dataless SEED file (so in "bash" you say "export ALT_RESPONSE_FILE=datalessSEEDfilename"). Then run "rdseed" without any parameters, i.e. in interactive mode. Give it the name of your miniSEED file as the data file to use, and select option "d" (dump). When it prompts you "output poles & zeros", answer y (for "yes"). It will create poles-and-zeros files with names beginning "SAC_PZs_" (it will also create SAC files of your data and "response" files with names beginning "RESP" containing filter information plus the poles and zeros again). The "SAC_PZs" files are in the right format for "transfer from polezero..." in SAC. Good luck, Sheila Peacock. From panxu.zhang at usc.edu Fri Sep 24 20:36:46 2010 From: panxu.zhang at usc.edu (Panxu Zhang) Date: Fri, 24 Sep 2010 20:36:46 -0700 Subject: [SAC-HELP] begin time in SAC header is incorrect Message-ID: Hi all, I download some earthquakes from IRIS. I set time range as 100 before first P and 600 after first P. But the begin time(B in SAC header) in SAC file is zero. How could I get the correct time? (assume the time when earthquake happened is zero, then B should be a number like 600 or something) I think of a method that, first calculate the first P time as T1MARKER with TAUP, then change header B as (T1-100). But SAC cannot recognize the command---"ch b (t1-100)". How could I get it run? Thank you in advance! Panxu From snoke at vt.edu Sat Sep 25 05:06:26 2010 From: snoke at vt.edu (Arthur Snoke) Date: Sat, 25 Sep 2010 08:06:26 -0400 (EDT) Subject: [SAC-HELP] begin time in SAC header is incorrect In-Reply-To: References: Message-ID: The default setting for SAC data files is to have IZTYPE set to IB, the begin time (B). Alternatives of possible interest to you are IO, origin time, or IA, first arrival time. Check out the SAC Header Table at . To change the value for IZTYPE, you use the command CH within SAC after reading in the file of interest. On Fri, 24 Sep 2010, Panxu Zhang wrote: > Hi all, > > I download some earthquakes from IRIS. I set time range as 100 before > first P and 600 after first P. But the begin time(B in SAC header) in > SAC file is zero. How could I get the correct time? (assume the time > when earthquake happened is zero, then B should be a number like 600 > or something) > > I think of a method that, first calculate the first P time as T1MARKER > with TAUP, then change header B as (T1-100). But SAC cannot recognize > the command---"ch b (t1-100)". How could I get it run? > > Thank you in advance! > > Panxu > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help > From m.tchelo at gmail.com Tue Sep 28 02:12:49 2010 From: m.tchelo at gmail.com (Marcelo Bianchi) Date: Tue, 28 Sep 2010 11:12:49 +0200 Subject: [SAC-HELP] bug on transfer evalresp Message-ID: Dear All, Today I was using the transfer command to remove the instrument from resp files and found something strange. The transfer evalresp couldn?t find the RESP file because of a time calculation error. My file started at 01:45:05.280 and the resp file, here is the output from SAC and the header from the RESP file: SAC> lh kzdate kztime FILE: 2008.330.01.45.05.2800.PD.EW01..HHZ.D.SAC - 1 ----------------------------------------------- kzdate = NOV 25 (330), 2008 kztime = 01:45:05.280 SAC> B052F22 Start date: 2008,330,01:45:05.2800 B052F23 End date: 2008,330,02:00:11.9200 when I try to run the transfer I got the folowing message: SAC> transfer from evalresp to none freqlimits 0.01 0.5 2. 5. Extracting evresp response for EW01, HHZ... [PD.EW01.*.HHZ] WARNING: no response found for NET=PD,STA=EW01,LOCID=*,CHAN=HHZ,DATE=2008,330,1:45:5 No transfer function applied! ERROR 2118: No transfer function applied. SAC> Pay special attention to the error output I noticed it was not considering the miliseconds, it was trying to find the date equal to 'DATE=2008,330,1:45:5'. As a test case, If I change the RESP file, removing the miliseconds part, like show below: B052F22 Start date: 2008,330,01:45:05 B052F23 End date: 2008,330,02:00:11.9200 It works: SAC> transfer from evalresp to none freqlimits 0.01 0.5 2. 5. Extracting evresp response for EW01, HHZ... Station (EW01 ), Channel (HHZ ) Waveform multiplied by 1.000000 after deconvolution. SAC> As I could read from the Seed manual it should consider the miliseconds for searching for the resp file don't ? Does anyone noticed this problem also ? This bug was noticied on the folowing SAC version: SEISMIC ANALYSIS CODE [06/07/2010 (Version 101.4)] Copyright 1995 Regents of the University of California SAC> Regards, Marcelo -- Mobile ~ +47 9080 6225 http://sites.google.com/site/foo4funreborn/ -------------- next part -------------- An HTML attachment was scrubbed... URL: From crotwell at seis.sc.edu Tue Sep 28 13:27:59 2010 From: crotwell at seis.sc.edu (Philip Crotwell) Date: Tue, 28 Sep 2010 16:27:59 -0400 Subject: [SAC-HELP] sac header enumerated values Message-ID: Hi I have a java based seismic file reading package called seisFile that handles SAC files amongst others (http://www.seis.sc.edu/seisFile.html). As part of that I have all the "enumerated header field values" as constants (ie static final fields). At the time I created this, they were all unique. I was looking at the sac docs today and noticed that there have been new values added (51-86). I was going to add them, but they are no longer unique. For example there are 2 IO "enumerated headers", 11 and 82. Is there a reason that these names were not kept unique? http://www.iris.edu/software/sac/manual/file_format.html thanks Philip From chad at iris.washington.edu Tue Sep 28 13:46:11 2010 From: chad at iris.washington.edu (Chad Trabant) Date: Tue, 28 Sep 2010 13:46:11 -0700 Subject: [SAC-HELP] sac header enumerated values In-Reply-To: References: Message-ID: <64CE4048-1ADC-4086-AC8E-C741EA424D72@iris.washington.edu> Hi Philip, I suspect there is a documentation error, maybe just the transcription to the HTML page. The values on that page don't quite match the enumerated values in the SacHeader.h that ships with 101.4 source code, the deviation starts . The defined value for 82 is "INU" (Nuclear cavity collapse), there is another IO-like value for "IO_" at number 84, there is also an "IOOTHER" at 44, at least that one matches the header. Somewhere around number 60 they deviate, looks like IPDEQ (59) and IPDEW (60) were included at some point in one list and not the other. Chad On Sep 28, 2010, at 1:27 PM, Philip Crotwell wrote: > Hi > > I have a java based seismic file reading package called seisFile that > handles SAC files amongst others > (http://www.seis.sc.edu/seisFile.html). > > As part of that I have all the "enumerated header field values" as > constants (ie static final fields). At the time I created this, they > were all unique. I was looking at the sac docs today and noticed that > there have been new values added (51-86). I was going to add them, but > they are no longer unique. For example there are 2 IO "enumerated > headers", 11 and 82. > > Is there a reason that these names were not kept unique? > > http://www.iris.edu/software/sac/manual/file_format.html > > thanks > Philip > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help From crotwell at seis.sc.edu Tue Sep 28 14:57:03 2010 From: crotwell at seis.sc.edu (Philip Crotwell) Date: Tue, 28 Sep 2010 17:57:03 -0400 Subject: [SAC-HELP] sac header enumerated values In-Reply-To: <64CE4048-1ADC-4086-AC8E-C741EA424D72@iris.washington.edu> References: <64CE4048-1ADC-4086-AC8E-C741EA424D72@iris.washington.edu> Message-ID: OK, so I will assume utils/sac.h from the 101.4 tarball should be considered definitive? I have to say "IO_" is a pretty dumb name for a constant, but I am going to make a real effort to use IODOR = 97 /* Odors */ as soon as I get a chance! :) Philip On Tue, Sep 28, 2010 at 4:46 PM, Chad Trabant wrote: > > Hi Philip, > > I suspect there is a documentation error, maybe just the transcription to the HTML page. ?The values on that page don't quite match the enumerated values in the SacHeader.h that ships with 101.4 source code, the deviation starts . ?The defined value for 82 is "INU" (Nuclear cavity collapse), there is another IO-like value for "IO_" at number 84, there is also an "IOOTHER" at 44, at least that one matches the header. ?Somewhere around number 60 they deviate, looks like IPDEQ (59) and IPDEW (60) were included at some point in one list and not the other. > > Chad > > > On Sep 28, 2010, at 1:27 PM, Philip Crotwell wrote: > >> Hi >> >> I have a java based seismic file reading package called seisFile that >> handles SAC files amongst others >> (http://www.seis.sc.edu/seisFile.html). >> >> As part of that I have all the "enumerated header field values" as >> constants (ie static final fields). At the time I created this, they >> were all unique. I was looking at the sac docs today and noticed that >> there have been new values added (51-86). I was going to add them, but >> they are no longer unique. For example there are 2 IO "enumerated >> headers", 11 and 82. >> >> Is there a reason that these names were not kept unique? >> >> http://www.iris.edu/software/sac/manual/file_format.html >> >> thanks >> Philip >> _______________________________________________ >> sac-help mailing list >> sac-help at iris.washington.edu >> http://www.iris.washington.edu/mailman/listinfo/sac-help > > From chad at iris.washington.edu Tue Sep 28 16:11:30 2010 From: chad at iris.washington.edu (Chad Trabant) Date: Tue, 28 Sep 2010 16:11:30 -0700 Subject: [SAC-HELP] bug on transfer evalresp In-Reply-To: References: Message-ID: Hi Marcelo, I don't think this is a very common error, most response definitions are not so precise to match the data to the millisecond. In fact most response information we have does not including any sub-second time component. I did some quick testing and there does seem to be a bug in evalresp where responses are not matched when the sub-second times are specified. We'll look into that. The practical solution is to not include such high precision times in your RESP files. cheers, Chad PS. Where did you get that RESP? On Sep 28, 2010, at 2:12 AM, Marcelo Bianchi wrote: > Dear All, > > Today I was using the transfer command to remove the instrument from resp files and found something strange. The transfer evalresp couldn?t find the RESP file because of a time calculation error. > > My file started at 01:45:05.280 and the resp file, here is the output from SAC and the header from the RESP file: > > SAC> lh kzdate kztime > > > FILE: 2008.330.01.45.05.2800.PD.EW01..HHZ.D.SAC - 1 > ----------------------------------------------- > > kzdate = NOV 25 (330), 2008 > kztime = 01:45:05.280 > SAC> > > > B052F22 Start date: 2008,330,01:45:05.2800 > B052F23 End date: 2008,330,02:00:11.9200 > > when I try to run the transfer I got the folowing message: > > SAC> transfer from evalresp to none freqlimits 0.01 0.5 2. 5. > Extracting evresp response for EW01, HHZ... > [PD.EW01.*.HHZ] WARNING: no response found for NET=PD,STA=EW01,LOCID=*,CHAN=HHZ,DATE=2008,330,1:45:5 > No transfer function applied! > ERROR 2118: No transfer function applied. > SAC> > > Pay special attention to the error output I noticed it was not considering the miliseconds, it was trying to find the date equal to 'DATE=2008,330,1:45:5'. As a test case, If I change the RESP file, removing the miliseconds part, like show below: > > B052F22 Start date: 2008,330,01:45:05 > B052F23 End date: 2008,330,02:00:11.9200 > > It works: > > SAC> transfer from evalresp to none freqlimits 0.01 0.5 2. 5. > Extracting evresp response for EW01, HHZ... > > Station (EW01 ), Channel (HHZ ) > Waveform multiplied by 1.000000 after deconvolution. > SAC> > > As I could read from the Seed manual it should consider the miliseconds for searching for the resp file don't ? Does anyone noticed this problem also ? This bug was noticied on the folowing SAC version: > > SEISMIC ANALYSIS CODE [06/07/2010 (Version 101.4)] > Copyright 1995 Regents of the University of California > > SAC> > > Regards, > > Marcelo > -- > Mobile ~ +47 9080 6225 > http://sites.google.com/site/foo4funreborn/ > > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help -------------- next part -------------- An HTML attachment was scrubbed... URL: From m.tchelo at gmail.com Wed Sep 29 05:55:36 2010 From: m.tchelo at gmail.com (Marcelo Bianchi) Date: Wed, 29 Sep 2010 14:55:36 +0200 Subject: [SAC-HELP] bug on transfer evalresp In-Reply-To: References: Message-ID: Hi Chad, Actually I got this error from a dataset that I get from GEOFON. I requested data for P wave windows. And the responses I got (I don't know why) were matched to the windows I got. I will see if I could send you a package of seed that I results on the error for you to test since it's a restricted dataset. Lets see If I can get the permission to do it. The quick solution that I am doing is parsing my RESP files to sed filter like: sed -n '/^B052F22/s/\..*$/' InputResp > CorrectedResp regards, Marcelo Bianchi -- Mobile ~ +47 9080 6225 http://sites.google.com/site/foo4funreborn/ 2010/9/29 Chad Trabant > > Hi Marcelo, > > I don't think this is a very common error, most response definitions are > not so precise to match the data to the millisecond. In fact most response > information we have does not including any sub-second time component. > > I did some quick testing and there does seem to be a bug in evalresp where > responses are not matched when the sub-second times are specified. We'll > look into that. The practical solution is to not include such high > precision times in your RESP files. > > cheers, > Chad > > PS. Where did you get that RESP? > > On Sep 28, 2010, at 2:12 AM, Marcelo Bianchi wrote: > > Dear All, > > Today I was using the transfer command to remove the instrument from resp > files and found something strange. The transfer evalresp couldn?t find the > RESP file because of a time calculation error. > > My file started at 01:45:05.280 and the resp file, here is the output from > SAC and the header from the RESP file: > > SAC> lh kzdate kztime > > > FILE: 2008.330.01.45.05.2800.PD.EW01..HHZ.D.SAC - 1 > ----------------------------------------------- > > kzdate = NOV 25 (330), 2008 > kztime = 01:45:05.280 > SAC> > > > B052F22 Start date: 2008,330,01:45:05.2800 > B052F23 End date: 2008,330,02:00:11.9200 > > when I try to run the transfer I got the folowing message: > > SAC> transfer from evalresp to none freqlimits 0.01 0.5 2. 5. > Extracting evresp response for EW01, HHZ... > [PD.EW01.*.HHZ] WARNING: no response found for > NET=PD,STA=EW01,LOCID=*,CHAN=HHZ,DATE=2008,330,1:45:5 > No transfer function applied! > ERROR 2118: No transfer function applied. > SAC> > > Pay special attention to the error output I noticed it was not considering > the miliseconds, it was trying to find the date equal to > 'DATE=2008,330,1:45:5'. As a test case, If I change the RESP file, removing > the miliseconds part, like show below: > > B052F22 Start date: 2008,330,01:45:05 > B052F23 End date: 2008,330,02:00:11.9200 > > It works: > > SAC> transfer from evalresp to none freqlimits 0.01 0.5 2. 5. > Extracting evresp response for EW01, HHZ... > > Station (EW01 ), Channel (HHZ ) > Waveform multiplied by 1.000000 after deconvolution. > SAC> > > As I could read from the Seed manual it should consider the miliseconds for > searching for the resp file don't ? Does anyone noticed this problem also ? > This bug was noticied on the folowing SAC version: > > SEISMIC ANALYSIS CODE [06/07/2010 (Version 101.4)] > Copyright 1995 Regents of the University of California > > SAC> > > Regards, > > Marcelo > -- > Mobile ~ +47 9080 6225 > http://sites.google.com/site/foo4funreborn/ > > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help > > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From montugg08 at gmail.com Thu Sep 30 21:44:11 2010 From: montugg08 at gmail.com (MONTU GOGOI) Date: Fri, 1 Oct 2010 10:14:11 +0530 Subject: [SAC-HELP] Sac tutorial Message-ID: Respected sir, Most humbly and respectfully I would like to inform you that I,Mr Montu Gogoi a studend of computational seismology(M.Tech)and I am a new user of the software SAC(seismic analysis code).I try to learn this software but in many steps I haved faced many problems and confusions.As a seismology student it is compulsory for me and I must learn this.So I kindly request you that if you have some materials and tutorials then please send to me. Yours faithfully Mr Montu Gogoi Tezpur University Assam(India) -------------- next part -------------- An HTML attachment was scrubbed... URL: From Onur.Tan at mam.gov.tr Thu Sep 30 23:13:34 2010 From: Onur.Tan at mam.gov.tr (Onur Tan) Date: Fri, 1 Oct 2010 09:13:34 +0300 Subject: [SAC-HELP] Sac tutorial In-Reply-To: Message-ID: Dear Montu The best way is the manual web pages : http://www.iris.edu/software/sac/manual.html On the other hand, the SAC team should give the manuals as zipped file for off-line usage. best regards onur Dr. Onur TAN ---------------------------------------------- 40.7866N 29.4500E --------- T?B?TAK Marmara Ara?t?rma Merkezi, Yer ve Deniz Bilimleri Enstit?s? TUBITAK Marmara Research Center, Earth and Marine Sciences Institute Gebze - Kocaeli - TURKEY On 01.10.2010 07:44, "MONTU GOGOI" wrote: Respected sir, Most humbly and respectfully I would like to inform you that I,Mr Montu Gogoi a studend of computational seismology(M.Tech)and I am a new user of the software SAC(seismic analysis code).I try to learn this software but in many steps I haved faced many problems and confusions.As a seismology student it is compulsory for me and I must learn this.So I kindly request you that if you have some materials and tutorials then please send to me. Yours faithfully Mr Montu Gogoi Tezpur University Assam(India)