From pku132 at 163.com Tue Jan 4 10:45:56 2011 From: pku132 at 163.com (=?UTF-8?Q?Inor=C2=A0Lee?=) Date: Wed, 5 Jan 2011 02:45:56 +0800 (CST) Subject: [SAC-HELP] where can i find the instruction for response file extracted from seed files by rdseed? Message-ID: <3d12b4e1.220.12d52598aab.Coremail.pku132@163.com> Hi, all, I want to know the meanings of each lines in the RESP.* files generated by RDSEED. Where can I find an instruction? Regards, Lee -------------- next part -------------- An HTML attachment was scrubbed... URL: From chad at iris.washington.edu Tue Jan 4 11:43:00 2011 From: chad at iris.washington.edu (Chad Trabant) Date: Tue, 4 Jan 2011 11:43:00 -0800 Subject: [SAC-HELP] where can i find the instruction for response file extracted from seed files by rdseed? In-Reply-To: <3d12b4e1.220.12d52598aab.Coremail.pku132@163.com> References: <3d12b4e1.220.12d52598aab.Coremail.pku132@163.com> Message-ID: Hi Lee, This is not a SAC question per-say, but since SAC reads RESP files perhaps it will be useful for others. There is no formal documentation for the RESP files themselves, but since they refer directly to the SEED format you can read the SEED manual to learn more about the values. The SEED manual is available here: http://www.iris.edu/manuals/SEEDManual_V2.4.pdf Inside a RESP file each non-comment line begins with a "blockette" and a "field" reference in the form of "BxxxFyy" where xxx is the blockette number and yy is the field number. A line beginning with B054F03 indicates it contains the information in SEED blockette type 54, field number 3. You can then lookup the definition of that field in the SEED manual. regards, Chad On Jan 4, 2011, at 10:45 AM, Inor Lee wrote: > Hi, all, > > I want to know the meanings of each lines in the RESP.* files generated by RDSEED. Where can I find an instruction? > > Regards, > Lee > > > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help From snoke at vt.edu Tue Jan 4 12:00:50 2011 From: snoke at vt.edu (Arthur Snoke) Date: Tue, 4 Jan 2011 15:00:50 -0500 (EST) Subject: [SAC-HELP] problem with interpolation In-Reply-To: <002901cb907e$73094ee0$262acf9b@maria2gngjfgil> References: <002901cb907e$73094ee0$262acf9b@maria2gngjfgil> Message-ID: I do not see any response to this message, so I will make a couple of comments. Looking at your plots, it looks like you have high amplitude information at high frequencies. The explanation to your problem may simply be that your data are aliased when you downsample. I suggest you try repeating your runs using decimate with the antialiasing filter on. Command interpolate has no antialiasing. On Tue, 30 Nov 2010, maria manakou wrote: > Hello, > > I face some problem when?applying?the command "interpolate delta XX" in > some?noise data. > > I have?some?1-hour?noise data recorded in CMG-5 guralp?sensors and I want to > interpolate them for further processing in?SAC software.?The raw?*.cgf files (the > format of the CMG-5 sensors) have?already transformed to *.sac files through > Scream V4.4 software provided by Guralp. > > The *.sac?data (un-interpolated data) have delta 0.01sec and I want to > interpolate them to delta 0.05sec.?I use the following commands in my SAC > 8/8/2001 V. 00.59.44 (Copyright 1995) software which runs in Fedora 9.0: > > r *.sac > qdp off > rmean > rtrend > p1 > interpolate delta 0.05 > w change .sac .inter > > When I plot the un-interpolated (*.sac) and interpolated data (*.inter) to see if > the interpolation works ok (Fig. 1 in the attach file) the result seems ok, but > when I zoom?in?different time windows the result is not at all good!?Specific > spikes in my un-interpolated noise data?seem to "move" in time in > the?interpolated data as you can see in the attach doc file.?This "moving" seems > to happen after some time (e.g. in fig. 2 the coincidence of?specific?spikes > between?un-interpolated and interpolated data is very good, but this does not > happen for the rest time windows,fig. 3-5). > > I cannot?understand why this happen. > If someone has faced this problem again and have any idea, please let me know. > > Best Regards, > maria manakou > > > > > > From snoke at vt.edu Tue Jan 4 12:05:31 2011 From: snoke at vt.edu (Arthur Snoke) Date: Tue, 4 Jan 2011 15:05:31 -0500 (EST) Subject: [SAC-HELP] where can i find the instruction for response file extracted from seed files by rdseed? In-Reply-To: <3d12b4e1.220.12d52598aab.Coremail.pku132@163.com> References: <3d12b4e1.220.12d52598aab.Coremail.pku132@163.com> Message-ID: I suggest you look at the SEED manual or one of the other manuals that can be downloaded from IRIS at . On Wed, 5 Jan 2011, Inor?Lee wrote: > Hi, all, > I want to know the meanings of each lines in the RESP.* files generated by > RDSEED. Where can I find an instruction? > > Regards, > Lee > > > > From sheila at blacknest.gov.uk Wed Jan 5 01:31:17 2011 From: sheila at blacknest.gov.uk (Sheila Peacock) Date: Wed, 05 Jan 2011 09:31:17 +0000 Subject: [SAC-HELP] RESP files from SEED data Message-ID: <4D243A65.90304@blacknest.gov.uk> I had a go some time ago to understand the meaning of RESP files and wrote some notes, which follow. Please note that they might not be correct, or might be incorrect for data from other stations or instruments. I was trying to establish the relationship between the parameters in RESP files and the parameters in headers of IMS (GSE)-format data from the CTBTO International Data Centre. Sheila Peacock Blacknest. Meaning of numbers in rdseed RESP files Here is a sample rdseed RESP file, which is output if you ask rdseed to output responses: * RESP.IU.FURI.00.BHE: "RESP" file for FURI component BHE The topmost part is like this: # << IRIS SEED Reader, Release 4.6 >> # # ======== CHANNEL RESPONSE DATA ======== B050F03 Station: FURI B050F16 Network: IU B052F03 Location: 00 B052F04 Channel: BHE B052F22 Start date: 1999,111,00:00:00 B052F23 End date: No Ending Time # ======================================= # + +--------------------------------------------+ + # + | Response (Poles & Zeros), FURI ch BHE | + # + +--------------------------------------------+ + # B053F03 Transfer function type: A [Laplace Transform (Rad/sec)] B053F04 Stage sequence number: 1 B053F05 Response in units lookup: M/S - Velocity in Meters Per Second B053F06 Response out units lookup: V - Volts B053F07 A0 normalization factor: 3948.58 B053F08 Normalization frequency: 0.02 B053F09 Number of zeroes: 2 B053F14 Number of poles: 4 # Complex zeroes: # i real imag real_error imag_error B053F10-13 0 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 B053F10-13 1 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 # Complex poles: # i real imag real_error imag_error B053F15-18 0 -1.234000E-02 1.234000E-02 0.000000E+00 0.000000E+00 B053F15-18 1 -1.234000E-02 -1.234000E-02 0.000000E+00 0.000000E+00 B053F15-18 2 -3.918000E+01 4.912000E+01 0.000000E+00 0.000000E+00 B053F15-18 3 -3.918000E+01 -4.912000E+01 0.000000E+00 0.000000E+00 # # + +---------------------------------------+ + # + | Channel Gain, FURI ch BHE | + # + +---------------------------------------+ + # B058F03 Stage sequence number: 1 B058F04 Gain: 2.296000E+03 B058F05 Frequency of gain: 2.000000E-02 HZ B058F06 Number of calibrations: 0 # # + +-------------------------------------------+ + # + | Response (Coefficients), FURI ch BHE | + # + +-------------------------------------------+ + # B054F03 Transfer function type: D B054F04 Stage sequence number: 2 B054F05 Response in units lookup: V - Volts B054F06 Response out units lookup: COUNTS - Digital Counts B054F07 Number of numerators: 0 B054F10 Number of denominators: 0 # # + +------------------------------+ + # + | Decimation, FURI ch BHE | + # + +------------------------------+ + # B057F03 Stage sequence number: 2 B057F04 Input sample rate: 5.120000E+03 B057F05 Decimation factor: 1 B057F06 Decimation offset: 0 B057F07 Estimated delay (seconds): 0.000000E+00 B057F08 Correction applied (seconds): 0.000000E+00 # # + +---------------------------------------+ + # + | Channel Gain, FURI ch BHE | + # + +---------------------------------------+ + # B058F03 Stage sequence number: 2 B058F04 Gain: 4.194300E+05 B058F05 Frequency of gain: 0.000000E+00 HZ B058F06 Number of calibrations: 0 # And the very end is like this: # + +---------------------------------------+ + # + | Channel Sensitivity, FURI ch BHE | + # + +---------------------------------------+ + # B058F03 Stage sequence number: 0 B058F04 Sensitivity: 9.630000E+08 B058F05 Frequency of sensitivity: 2.000000E-02 HZ B058F06 Number of calibrations: 0 # The corresponding poles-zeros file, which is output if you ask rdseed to put it out, is: ZEROS 3 POLES 4 -0.0123 0.0123 -0.0123 -0.0123 -39.1800 49.1200 -39.1800 -49.1200 CONSTANT 3.802483e+12 The "RESP" file represents the velocity response, and hence the poles and zeros listed there include two zeros while the pole-zero file immediately above has three since it represents the displacement response. The "Sensitivity" at the end of the "RESP" file is the product of all the gains in the file (all the numbers in lines containing "Gain": in this case, 9.63e8 = 2.296e3 V/(m/s) x 4.1943e5 counts/V). It is in the velocity domain and has units of counts per (metre per second) so to convert it to displacement sensitivity (counts per metre) you would MULTIPLY by 2pi times the calibration frequency (0.02 Hz in this case). This gives you 10^9 / the value of "calib" that would be found in GSE/IMS format files, CSS .wfdisc files and "static.instrument" IDC database table (8.263497 in this case). "Calib" is in units of nanometres per count. The number "A0 normalization factor" (value 3948.58 in this case) is said in the SEED manual to be "set such that when you evaluate the polynomial at the reference frequency the result will be one". It is the "poles-and-zeros constant" for the velocity poles-zeros in the "RESP" file. If you calculate the poles-zeros constant (a hefty job best done in fortran for the complex arithmetic required - beware of floating-point overflow, though!) from the displacement poles-zeros file output by rdseed, rather than take the "RESP" file value, you should get the "displacement" poles-and-zeros constant, which is equal to the A0 normalization factor value divided by (2pi times the calibration frequency) (in this case 0.02 Hz, which gives 31421.7). If you multiply the A0 normalization factor by the "sensitivity" (9.63e8 counts/(m/s)), you get the value of the poles-zeros constant at the end of the poles-zeros file produced by rdseed (3.802483e12 in this case). Note that because this constant consists of the product of a velocity sensitivity (a displacement sensitivity divided by 2pi x calibration frequency) and a velocity poles-zeros constant (a displacement poles-zeros constant multiplied by 2pi x calibration frequency), the factors of 2pi x calibration frequency cancel out and it is correctly a displacement poles-zeros constant (albeit multiplied by the sensitivity) and accompanies a set of displacement poles-zeros, as expected by SAC. Since the poles-zeros file written by rdseed contains a constant that has already been multiplied by the sensitivity and the sensitivity is equivalent to 1/(calibration factor "calib" in IDC parlance x 2pi x frequency), it is not necessary to multiply the seismogram by the calibration factor. From alessandro.pino at ingv.it Wed Jan 5 01:56:02 2011 From: alessandro.pino at ingv.it (alessandro.pino at ingv.it) Date: Wed, 05 Jan 2011 10:56:02 +0100 Subject: [SAC-HELP] The results of your email commands In-Reply-To: References: Message-ID: <20110105105602.61571cg0ic340ytc@webmail.rm.ingv.it> Hi Tsitsi, I think the simplest way could be to synchronize your files at first. This will fix a common origin time for all the seismograms, setting the beginning time in accordance, with beg=0 for the latest file and beg<0 for the others. Then you should write the seismograms on temporary files a read them again after setting the cut range from zero to the desired time. Ex.: r file1.sac file2.sac synch write change .sac .cut cut 0 1000 r file1.cut file2.cut w over Regards Alessandro Def. Quota Tsitsi Madziwa-Nussinov : >> >> Hi Sac Users >> > > >> I have 2 sac files (examples attached). Their start times are 09:49:50 and >> 09:46:43 >> I would like to cut off the start times so that both files would start at >> the same time, >> say, the start time of the 9:49:50 file. >> >> How can I make use of the synchronize and cut commands to achieve this? The >> end >> goal is be able to perform some math on these files with matched times for >> corresponding >> points? >> >> Thanks. >> >> >> > ---------------------------------------------------------------- Il contenuto di questa e-mail e' rivolto unicamente alle persone cui e' indirizzato, e puo' contenere informazioni la cui riservatezza e' tutelata. E' proibita la copia, la divulgazione o l'uso di questo messaggio o dell'informazione ivi contenuta da chiunque altro che non sia il destinatario indicato. Se avete ricevuto questa e-mail per errore, vogliate cortesemente comunicarlo immediatamente per telefono, fax o e-mail. Grazie. This e-mail is intended only for person or entity to which is addressed and may contain information that is privileged, confidential or otherwise protected from disclosure. Copying, dissemination or use of this e-mail or the information herein by anyone other than the intended recipient is prohibited. If you have received this e-mail by mistake, please notify us immediately by telephone, fax or e-mail. Thank you. From tmacdonald at usgs.gov Thu Jan 6 10:25:39 2011 From: tmacdonald at usgs.gov (Tim MacDonald) Date: Thu, 6 Jan 2011 10:25:39 -0800 Subject: [SAC-HELP] cutim issue for 101.4 Message-ID: I am having a problem getting cutim to work on the 101.4 release for Linux. I got a segmentation fault error when using cutim as follows: SAC> read * 2002.110.10.49.30.0000.CN.BBB..BHE.D.SAC 2002.110.10.49.30.0000.CN.BBB..BHN.D.SAC 2002.110.10.49.30.0000.CN.BBB..BHZ.D.SAC 2002.110.10.49.30.0000.CN.BIB..EHZ.D.SAC SAC> cutim 100 200 Segmentation fault I tracked the segmentation fault error to the gcCollect routine in gcCollect.c and made the following changes to the code: /* struct CSStree *tree = (struct CSStree *) dblList; struct wftagList *wtL = tree->wtHead ; struct sensorList *seL = tree->seHead ; struct sitechanList *scL = tree->scHead ; struct siteList *siL = tree->slHead ; struct instrumentList *inL = tree->inHead ; struct originList *orL = tree->orHead ; struct eventList *evL = tree->evHead ; if(!tree)return; */ struct CSStree *tree = (struct CSStree *) dblList; struct wftagList *wtL; struct sensorList *seL; struct sitechanList *scL; struct siteList *siL; struct instrumentList *inL; struct originList *orL; struct eventList *evL; if(!tree)return; wtL = tree->wtHead ; seL = tree->seHead ; scL = tree->scHead ; siL = tree->slHead ; inL = tree->inHead ; orL = tree->orHead ; evL = tree->evHead ; The changes cured the segmentation fault error but now I get error 1301 when I try to plot as follows: SAC> read * 2002.110.10.49.30.0000.CN.BBB..BHE.D.SAC 2002.110.10.49.30.0000.CN.BBB..BHN.D.SAC 2002.110.10.49.30.0000.CN.BBB..BHZ.D.SAC 2002.110.10.49.30.0000.CN.BIB..EHZ.D.SAC SAC> cutim 100 200 SAC> plot1 ERROR 1301: No data files read in. If I don't issue the cutim command the plot1 works fine. Also if I use the OSX binaries and the same data files everything works fine! Any ideas? Thanks! From alireza.alinaghi at gmail.com Wed Jan 12 05:43:10 2011 From: alireza.alinaghi at gmail.com (Alireza Alinaghi) Date: Wed, 12 Jan 2011 14:43:10 +0100 Subject: [SAC-HELP] Problem with depmec Message-ID: Hi I have been trying to use depmec in sac but after i saw matlab logo flashing by, I got the following error message in sac command window: *Allocating and filling Matlab workspace arrays... * *Copying to Matlab workspace ... * *Done.* *>> ??? Undefined function or variable 'sacdepmec'.* * * Could you please tell me what is the source of the problem and how i can fix it? I am using sac version 101.0 and matlab (7)R14 on Fedora 10 linux machine. Thank you in advance, Alireza -------------- next part -------------- An HTML attachment was scrubbed... URL: From geophysia at gmail.com Wed Jan 12 21:24:52 2011 From: geophysia at gmail.com (Sierra Boyd) Date: Wed, 12 Jan 2011 21:24:52 -0800 Subject: [SAC-HELP] distance and azimuth Message-ID: Hello, I would like to know how to get the distance and azimuth of a station with respect to the hypocentral info, evla and evlo. I read in a sac file that came from IRIS. I then change the values of evla and evlo. I then listhdr inclusive to look at all the header values. The distance and azimuth are not updated. However I can get the distance and azimuth of a station with respect to an event when I download miniSEED data from another source and convert it to SAC format. I then change the values of evla and evlo as above and the distance and azimuth are listed. Any clarification as to why my first attempt does not reveal an updated distance and azimuth would be greatly appreciated. Best regards, --sierra From abhey.bansal at gmail.com Wed Jan 12 23:16:44 2011 From: abhey.bansal at gmail.com (Abhey Bansal) Date: Thu, 13 Jan 2011 12:46:44 +0530 Subject: [SAC-HELP] installing SAC on UBUNTU Message-ID: Recently I was installing SAC on my laptop running on UBUNTU but when I modified my ~/.bashrc file adding following two lines export SACHOME=/usr/local/sac source ${SACHOME}/bin/sacinit.csh and after saving the file. I run following command source ~/.bashrc the following error message come. No command 'setenv' found, did you mean: Command 'netenv' from package 'netenv' (universe) Could any one please help me sorting out the problem. Thanking you Abhey -------------- next part -------------- An HTML attachment was scrubbed... URL: From Onur.Tan at mam.gov.tr Thu Jan 13 00:00:04 2011 From: Onur.Tan at mam.gov.tr (Onur Tan) Date: Thu, 13 Jan 2011 10:00:04 +0200 Subject: [SAC-HELP] installing SAC on UBUNTU In-Reply-To: Message-ID: Hi, You source Cshell file ?sacinit.csh? in bash. This can be problem. bash can not understand Cshell code. o. 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 13.01.2011 09:16, "Abhey Bansal" wrote: Recently I was installing SAC on my laptop running on UBUNTU but when I modified my ~/.bashrc file adding following two lines export SACHOME=/usr/local/sac source ${SACHOME}/bin/sacinit.csh and after saving the file. I run following command source ~/.bashrc the following error message come. No command 'setenv' found, did you mean: Command 'netenv' from package 'netenv' (universe) Could any one please help me sorting out the problem. Thanking you Abhey From Onur.Tan at mam.gov.tr Thu Jan 13 00:02:30 2011 From: Onur.Tan at mam.gov.tr (Onur Tan) Date: Thu, 13 Jan 2011 10:02:30 +0200 Subject: [SAC-HELP] installing SAC on UBUNTU In-Reply-To: Message-ID: Hi again, use sacinit.sh for bash 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 13.01.2011 09:16, "Abhey Bansal" wrote: Recently I was installing SAC on my laptop running on UBUNTU but when I modified my ~/.bashrc file adding following two lines export SACHOME=/usr/local/sac source ${SACHOME}/bin/sacinit.csh and after saving the file. I run following command source ~/.bashrc the following error message come. No command 'setenv' found, did you mean: Command 'netenv' from package 'netenv' (universe) Could any one please help me sorting out the problem. Thanking you Abhey From geophysia at gmail.com Thu Jan 20 18:43:01 2011 From: geophysia at gmail.com (Sierra Boyd) Date: Thu, 20 Jan 2011 18:43:01 -0800 Subject: [SAC-HELP] bandpass filter Message-ID: Hello, I have a question about the number of poles for a Butterworth filter when there are 2 passes. When n, the number of poles, is not specified, is the default value of 2 used? Does this mean there are 4 poles with 2 passes when n is not specified? Suppose n = 4 and p = 2. Are there 8 poles? Best regards, --Sierra From Fatih.Sevim at mam.gov.tr Mon Jan 24 06:52:24 2011 From: Fatih.Sevim at mam.gov.tr (Fatih Sevim) Date: Mon, 24 Jan 2011 16:52:24 +0200 Subject: [SAC-HELP] transfer command Message-ID: <9E88105FE4AC8C47BD6D00EB9831482F1E863D39FE@posta.mam.gov.tr> Hi, I have a pole zero files to avoid instrument response from signal and I use the command like " transfer from polezero subtype pz_file ". Is this the true way to obtain signal which doesn't include instrument response. Or I must add " to vel " end of the command line ? Best regards, Fatih. -------------- next part -------------- An HTML attachment was scrubbed... URL: From sheila at blacknest.gov.uk Mon Jan 24 07:12:23 2011 From: sheila at blacknest.gov.uk (Sheila Peacock) Date: Mon, 24 Jan 2011 15:12:23 +0000 Subject: [SAC-HELP] Fwd: transfer command Message-ID: <4D3D96D7.3030902@blacknest.gov.uk> To remove the instrument response completely you need "to none" at the end of the "transfer" command. You might need to use "freqlimits" in your "transfer" command to limit the low-frequency part of the signal; and I have always found it necessary to run "rmean", "rtrend" and "taper" on a seismogram before attempting to remove the instrument response. Sheila Peacock Blacknest. -------- Original Message -------- Subject: [SAC-HELP] transfer command Date: Mon, 24 Jan 2011 16:52:24 +0200 From: Fatih Sevim To: sac-help at iris.washington.edu Hi, I have a pole zero files to avoid instrument response from signal and I use the command like ? transfer from polezero subtype pz_file ?. Is this the true way to obtain signal which doesn?t include instrument response. Or I must add ? to vel ? end of the command line ? Best regards, Fatih. From shaharb at stanford.edu Sat Jan 29 10:38:30 2011 From: shaharb at stanford.edu (Shahar Barak) Date: Sat, 29 Jan 2011 10:38:30 -0800 Subject: [SAC-HELP] SAC - SORT command not valid In-Reply-To: References: Message-ID: > Dear sir/madam, > > I am trying to sort traces by distance. > I found a "sort" command in the user guide which seems to allow to sort > files according to a header field: > http://www.iris.edu/software/sac/commands/sort.html > > However, when I use it I get: "ERROR 1106: Not a valid SAC command" > > Thanks, > > Shahar Barak > > -- > Shahar Barak > ======================================= > *Department of Geophysics, Stanford University* > Mitchell Building , 397 Panama Mall, Stanford CA 94305-2215 > Office: (650) 725 1625 > Cell: (650) 276 6169 > Home: (650) 497 3818 > shaharb at stanford.edu * > * > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From georgesand at unb.br Mon Jan 31 02:09:24 2011 From: georgesand at unb.br (George Sand) Date: Mon, 31 Jan 2011 08:09:24 -0200 Subject: [SAC-HELP] sac-help Digest, Vol 64, Issue 2 In-Reply-To: References: Message-ID: <4D468A54.1060106@unb.br> Dear sierra For update Azimuth and distance (DIST and GCARC) you need to change evdp, too! cheers Prof. George Sand L. A. de Fran?a Observat?rio Sismol?gico Universidade de Bras?lia SG 13 Campus Universit?rio Darcy Ribeiro 70910-900 55 61 33404055 55 61 93033829 (pessoal) www.obsis.unb.br www.tremordeterra.blogspot.com "Never be afraid to try something new. Remember, amateurs built the ark, professionals built the Titanic. " sac-help-request at iris.washington.edu escreveu: > Send sac-help mailing list submissions to > sac-help at iris.washington.edu > > To subscribe or unsubscribe via the World Wide Web, visit > http://www.iris.washington.edu/mailman/listinfo/sac-help > or, via email, send a message with subject or body 'help' to > sac-help-request at iris.washington.edu > > You can reach the person managing the list at > sac-help-owner at iris.washington.edu > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of sac-help digest..." > > > Today's Topics: > > 1. cutim issue for 101.4 (Tim MacDonald) > 2. Problem with depmec (Alireza Alinaghi) > 3. distance and azimuth (Sierra Boyd) > 4. installing SAC on UBUNTU (Abhey Bansal) > 5. Re: installing SAC on UBUNTU (Onur Tan) > 6. Re: installing SAC on UBUNTU (Onur Tan) > 7. bandpass filter (Sierra Boyd) > 8. transfer command (Fatih Sevim) > 9. Fwd: transfer command (Sheila Peacock) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Thu, 6 Jan 2011 10:25:39 -0800 > From: Tim MacDonald > Subject: [SAC-HELP] cutim issue for 101.4 > To: sac-help at iris.washington.edu > Message-ID: > Content-Type: text/plain; charset=us-ascii > > I am having a problem getting cutim to work on the 101.4 release for Linux. I got a segmentation fault error when using cutim as follows: > > SAC> read * > 2002.110.10.49.30.0000.CN.BBB..BHE.D.SAC 2002.110.10.49.30.0000.CN.BBB..BHN.D.SAC 2002.110.10.49.30.0000.CN.BBB..BHZ.D.SAC 2002.110.10.49.30.0000.CN.BIB..EHZ.D.SAC > SAC> cutim 100 200 > Segmentation fault > > > I tracked the segmentation fault error to the gcCollect routine in gcCollect.c and made the following changes to the code: > > /* > struct CSStree *tree = (struct CSStree *) dblList; > struct wftagList *wtL = tree->wtHead ; > struct sensorList *seL = tree->seHead ; > struct sitechanList *scL = tree->scHead ; > struct siteList *siL = tree->slHead ; > struct instrumentList *inL = tree->inHead ; > struct originList *orL = tree->orHead ; > struct eventList *evL = tree->evHead ; > > if(!tree)return; > */ > struct CSStree *tree = (struct CSStree *) dblList; > struct wftagList *wtL; > struct sensorList *seL; > struct sitechanList *scL; > struct siteList *siL; > struct instrumentList *inL; > struct originList *orL; > struct eventList *evL; > > if(!tree)return; > > wtL = tree->wtHead ; > seL = tree->seHead ; > scL = tree->scHead ; > siL = tree->slHead ; > inL = tree->inHead ; > orL = tree->orHead ; > evL = tree->evHead ; > > The changes cured the segmentation fault error but now I get error 1301 when I try to plot as follows: > > SAC> read * > 2002.110.10.49.30.0000.CN.BBB..BHE.D.SAC 2002.110.10.49.30.0000.CN.BBB..BHN.D.SAC 2002.110.10.49.30.0000.CN.BBB..BHZ.D.SAC 2002.110.10.49.30.0000.CN.BIB..EHZ.D.SAC > SAC> cutim 100 200 > SAC> plot1 > ERROR 1301: No data files read in. > > If I don't issue the cutim command the plot1 works fine. Also if I use the OSX binaries and the same data files everything works fine! > > Any ideas? > Thanks! > > > ------------------------------ > > Message: 2 > Date: Wed, 12 Jan 2011 14:43:10 +0100 > From: Alireza Alinaghi > Subject: [SAC-HELP] Problem with depmec > To: sac-help at iris.washington.edu > Message-ID: > > Content-Type: text/plain; charset="iso-8859-1" > > Hi > > I have been trying to use depmec in sac but after i saw matlab logo flashing > by, I got the following error message in sac command window: > > *Allocating and filling Matlab workspace arrays... * > *Copying to Matlab workspace ... * > *Done.* > *>> ??? Undefined function or variable 'sacdepmec'.* > * > * > Could you please tell me what is the source of the problem and how i can fix > it? > > I am using sac version 101.0 and matlab (7)R14 on Fedora 10 linux machine. > > Thank you in advance, > > Alireza > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > > ------------------------------ > > Message: 3 > Date: Wed, 12 Jan 2011 21:24:52 -0800 > From: Sierra Boyd > Subject: [SAC-HELP] distance and azimuth > To: sac-help at iris.washington.edu > Message-ID: > > Content-Type: text/plain; charset=ISO-8859-1 > > Hello, > > I would like to know how to get the distance and azimuth of a station > with respect to the hypocentral info, evla and evlo. > > I read in a sac file that came from IRIS. I then change the values of > evla and evlo. I then listhdr inclusive to look at all the header > values. The distance and azimuth are not updated. > > However I can get the distance and azimuth of a station with respect > to an event when I download miniSEED data from another source and > convert it to SAC format. I then change the values of evla and evlo as > above and the distance and azimuth are listed. > > Any clarification as to why my first attempt does not reveal an > updated distance and azimuth would be greatly appreciated. > > Best regards, > --sierra > > > ------------------------------ > > Message: 4 > Date: Thu, 13 Jan 2011 12:46:44 +0530 > From: Abhey Bansal > Subject: [SAC-HELP] installing SAC on UBUNTU > To: sac-help at iris.washington.edu > Message-ID: > > Content-Type: text/plain; charset="iso-8859-1" > > Recently I was installing SAC on my laptop running on UBUNTU but when I > modified my ~/.bashrc file adding following two lines > export SACHOME=/usr/local/sac > source ${SACHOME}/bin/sacinit.csh > and after saving the file. I run following command > source ~/.bashrc > the following error message come. > No command 'setenv' found, did you mean: > Command 'netenv' from package 'netenv' (universe) > > Could any one please help me sorting out the problem. > Thanking you > Abhey > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > > ------------------------------ > > Message: 5 > Date: Thu, 13 Jan 2011 10:00:04 +0200 > From: Onur Tan > Subject: Re: [SAC-HELP] installing SAC on UBUNTU > To: Sac-Help > Message-ID: > Content-Type: text/plain; charset="utf-8" > > Hi, > > > You source Cshell file ?sacinit.csh? in bash. This can be problem. > bash can not understand Cshell code. > > o. > > 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 13.01.2011 09:16, "Abhey Bansal" wrote: > > Recently I was installing SAC on my laptop running on UBUNTU but when I modified my ~/.bashrc file adding following two lines > export SACHOME=/usr/local/sac > source ${SACHOME}/bin/sacinit.csh > and after saving the file. I run following command > source ~/.bashrc > the following error message come. > No command 'setenv' found, did you mean: > Command 'netenv' from package 'netenv' (universe) > > Could any one please help me sorting out the problem. > Thanking you > Abhey > > > > > > > > > ------------------------------ > > Message: 6 > Date: Thu, 13 Jan 2011 10:02:30 +0200 > From: Onur Tan > Subject: Re: [SAC-HELP] installing SAC on UBUNTU > To: Sac-Help > Message-ID: > Content-Type: text/plain; charset="utf-8" > > Hi again, > > use sacinit.sh for bash > > 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 13.01.2011 09:16, "Abhey Bansal" wrote: > > Recently I was installing SAC on my laptop running on UBUNTU but when I modified my ~/.bashrc file adding following two lines > export SACHOME=/usr/local/sac > source ${SACHOME}/bin/sacinit.csh > and after saving the file. I run following command > source ~/.bashrc > the following error message come. > No command 'setenv' found, did you mean: > Command 'netenv' from package 'netenv' (universe) > > Could any one please help me sorting out the problem. > Thanking you > Abhey > > > > > > > > > > ------------------------------ > > Message: 7 > Date: Thu, 20 Jan 2011 18:43:01 -0800 > From: Sierra Boyd > Subject: [SAC-HELP] bandpass filter > To: sac-help at iris.washington.edu > Message-ID: > > Content-Type: text/plain; charset=ISO-8859-1 > > Hello, > > I have a question about the number of poles for a Butterworth filter > when there are 2 passes. > > When n, the number of poles, is not specified, is the default value of > 2 used? Does this mean there are 4 poles with 2 passes when n is not > specified? > > Suppose n = 4 and p = 2. Are there 8 poles? > > Best regards, > --Sierra > > > ------------------------------ > > Message: 8 > Date: Mon, 24 Jan 2011 16:52:24 +0200 > From: Fatih Sevim > Subject: [SAC-HELP] transfer command > To: "sac-help at iris.washington.edu" > Message-ID: > <9E88105FE4AC8C47BD6D00EB9831482F1E863D39FE at posta.mam.gov.tr> > Content-Type: text/plain; charset="iso-8859-1" > > > Hi, > > I have a pole zero files to avoid instrument response from signal and I use the command like " transfer from polezero subtype pz_file ". Is this the true way to obtain signal which doesn't include instrument response. > Or I must add " to vel " end of the command line ? > > Best regards, > > Fatih. > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > > ------------------------------ > > Message: 9 > Date: Mon, 24 Jan 2011 15:12:23 +0000 > From: Sheila Peacock > Subject: [SAC-HELP] Fwd: transfer command > To: "sac-help at iris.washington.edu" > Message-ID: <4D3D96D7.3030902 at blacknest.gov.uk> > Content-Type: text/plain; charset=UTF-8; format=flowed > > To remove the instrument response completely you need "to none" at the end of the > "transfer" command. You might need to use "freqlimits" in your "transfer" command > to limit the low-frequency part of the signal; and I have always found it necessary > to run "rmean", "rtrend" and "taper" on a seismogram before attempting to remove > the instrument response. > > Sheila Peacock > Blacknest. > > > -------- Original Message -------- > Subject: [SAC-HELP] transfer command > Date: Mon, 24 Jan 2011 16:52:24 +0200 > From: Fatih Sevim > To: sac-help at iris.washington.edu > > > > Hi, > > I have a pole zero files to avoid instrument response from signal and I use the command like ? transfer from polezero subtype pz_file ?. Is this the true way to obtain signal which doesn?t include instrument response. > > Or I must add ? to vel ? end of the command line ? > > Best regards, > > Fatih. > > > > ------------------------------ > > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help > > > End of sac-help Digest, Vol 64, Issue 2 > *************************************** > >