From mayenew at gmail.com Mon Feb 9 09:55:45 2009 From: mayenew at gmail.com (Melaku Ayenew) Date: Mon, 9 Feb 2009 10:55:45 -0700 Subject: [SAC-HELP] Help on PLOTPK Message-ID: hi all, I am just starting to learn sac ,and have some difficulties using the PLOTPK command, I know it produces a plot for picking of arrival time, by problem is how i pick those arrival times after i got the plot with a x-y cursor. Do i need to creat hypo file to store the picks, please be more explicit. thank you melaku Bogale New Mexico State University -------------- next part -------------- An HTML attachment was scrubbed... URL: From M2Sr at pge.com Mon Feb 9 14:41:32 2009 From: M2Sr at pge.com (Stanton, Megan) Date: Mon, 9 Feb 2009 14:41:32 -0800 Subject: [SAC-HELP] FW: Help on PLOTPK Message-ID: <1C99032ABE04AB42B56583CC44FA011A034AA58C@exchange16.Utility.pge.com> If you are looking for the actual commands to make picks look in the command listing for plotpk for all the options. Just place your cursor on where you want to make the pick and then type in the desired attributes. Usually you have (e)mergent or (i)mpulsive, (u)p or (d)own, quality (0,1,2,3,4), (p)wave or (s)wave. The letters or numbers in parenthesis are what you would type. Once you are done it will show you your pick and label it. Re-pick by moving the cursor and typing it in again. http://www.iris.edu/manuals/sac/SAC_Commands/plotpk There are 2 ways to store your picks. 1) you can use another program to create a pick file from your picks or 2) you can create a pick file with sac using ohpf. I use option 1 combined with sac2arc to make my pick file. Here is the macro for option 1 with comment lines preceded by #. #header info. fileid type list kevnm kstcmp #read all files found. r *.[E,H]*[E,N,Z].PG.-- #plot pick 3 files at a time, bell off. ppk p 3 bell off #write over existing wave file. writehdr over #end macro. q Here is the macro for option 2 with comment lines. #begin macro. :b #open hypo pick file 'picks'. ohpf picks #header info. fileid type list kevnm kstnm #do loop, find all files ending w E,N or Z. do ENZ wild c E,N,Z #read all files found. r *[E,N,Z]* #plot pick 3 files at a time, bell off. ppk p 3 bell off #write over existing wave file.w over #end do loop. enddo #end macro. :e HTH, Megan Megan Stanton Seismologist Pacific Gas & Electric Company Geosciences Dept. ________________________________ From: sac-help-bounces at iris.washington.edu [mailto:sac-help-bounces at iris.washington.edu] On Behalf Of Melaku Ayenew Sent: Monday, February 09, 2009 9:56 AM To: sac-help at iris.washington.edu Subject: [SAC-HELP] Help on PLOTPK hi all, I am just starting to learn sac ,and have some difficulties using the PLOTPK command, I know it produces a plot for picking of arrival time, by problem is how i pick those arrival times after i got the plot with a x-y cursor. Do i need to creat hypo file to store the picks, please be more explicit. thank you melaku Bogale New Mexico State University -------------- next part -------------- An HTML attachment was scrubbed... URL: From snoke at vt.edu Mon Feb 9 16:01:56 2009 From: snoke at vt.edu (Arthur Snoke) Date: Mon, 9 Feb 2009 19:01:56 -0500 (EST) Subject: [SAC-HELP] FW: Help on PLOTPK In-Reply-To: <1C99032ABE04AB42B56583CC44FA011A034AA58C@exchange16.Utility.pge.com> References: <1C99032ABE04AB42B56583CC44FA011A034AA58C@exchange16.Utility.pge.com> Message-ID: Also of potential use is the online plotpk table: SAC> help plotpktable From talalmg1 at yahoo.com Tue Feb 10 02:25:12 2009 From: talalmg1 at yahoo.com (Talal Merghelani) Date: Tue, 10 Feb 2009 02:25:12 -0800 (PST) Subject: [SAC-HELP] Help on reading/writing spectral file Message-ID: <524238.46352.qm@web37307.mail.mud.yahoo.com> Dear Sac user's, ? ? My problem is I can't read my file after saving it to the disk. For example, the original file is the amplitude file in x,y file type. Then, performing the division operation by other file,?and then FFT, will produce the spectral?data file type. Which I have to save it on my disk as a file(in binary or ASCII). The problem is, when, I save it using the command "WRITE SAC" it successfully saves my file, as binary, spectral file. But, when,?I use the command "WRITE SP", with no error. But, when I try to read this file it says "ERROR 1301: No data files read in.".???.?Even, when I use the command "READSP"?it gave me this message "ERROR? 108: File does not exist:". ? I don't know what is going wrong with me..??? ? Any Idea on what going on with my files, highly appreciated. ? ? My regards, Talal ? ? =======================My terminal================================ SAC>? r file1.am SAC>? p SAC> divf file2.am ?WARNING: Time overlap: file1.am file2.am ?? BEG1: OCT 14 (287), 2006? 06:38:01.000 ?? BEG2: OCT 14 (287), 2006? 06:38:15.000 SAC> fft am ?(10.6e)FFT default change: not removing the mean ?DC level after DFT is 125.91 SAC> w sac spectralfile.sac SAC> r spectralfile.sac SAC> psp Waiting SAC> wsp spectralfile.sac SAC> w alpha spectralfile.sac.txt SAC> sc ls file1.am???file2.am???spectralfile.sac??spectralfile.sac.am??spectralfile.sac.ph?? spectralfile.sac.txt? fort.19 SAC> r spectralfile.sac.am ?ERROR 1301: No data files read in. SAC> rsp spectralfile.sac.am ?ERROR? 108: File does not exist: SAC> rsp amph CHB012CHB0090610140638-1.sac.am ?ERROR? 108: File does not exist: SAC> r alpha spectralfile.sac.txt ?ERROR 1301: No data files read in. SAC> -------------- next part -------------- An HTML attachment was scrubbed... URL: From mohamed.salah at ist.utl.pt Wed Feb 11 06:44:38 2009 From: mohamed.salah at ist.utl.pt (Mohamed Kamaleldin Abdel Galeil Salah) Date: Wed, 11 Feb 2009 14:44:38 +0000 Subject: [SAC-HELP] How to remove the instrument response? Message-ID: <20090211144438.vagucohh4wkwkkss@mail1.ist.utl.pt> Hi all, I have some BB SAC files for a sensor CMG-40TD (Guralp systems) and I want to remove the instrument response from the seismogram. I prepared a file ?pz?'as follows: ZEROS 4 0 0 0 159 POLES 3 -50.0 0.0 -0.02356 + j0.02356 -0.02356 - j0.02356 CONSTANT -0.314 I have two questions: 1. Is such data ok for the instrument type given above? 2. when I run the command: transfer from polezero subtype pz to none, from within SAC, I get the following message: Extracting polezero response for MTE , BHZ... Station (MTE ), Channel (BHZ) Calper is not available. Cannot tell if response is normalized, so transfer function will be used without scaling. Calib is not available. Cannot tell if scale factor has been applied on input, so transfer function will be used without scaling. Waveform multiplied by 1.000000 after deconvolution. How I can fix this error? Thanks From sanne at berkeley.edu Wed Feb 11 15:36:14 2009 From: sanne at berkeley.edu (Sanne Cottaar) Date: Wed, 11 Feb 2009 15:36:14 -0800 (PST) Subject: [SAC-HELP] deconvolution errors Message-ID: <2722.128.32.102.163.1234395374.squirrel@calmail.berkeley.edu> Hi everyone, When I remove the instrument response, by poles&zeros or response files, I come across the following messages/errors: * Calper is not available. Cannot tell if response is normalized, so transfer function will be used without scaling. * Calib is not available. Cannot tell if scale factor has been applied on input, so transfer function will be used without scaling. * WARNING: FIR normalized Would you know what one of more of these messages mean? Thanks a lot for helping me out, Sanne -- Sanne Cottaar PhD student UC Berkeley sanne at berkeley.edu 510-529-8661 McCone hall 381 From sheila at blacknest.gov.uk Mon Feb 16 02:12:40 2009 From: sheila at blacknest.gov.uk (Sheila Peacock) Date: Mon, 16 Feb 2009 10:12:40 +0000 Subject: [SAC-HELP] How to remove the instrument response? In-Reply-To: <20090211144438.vagucohh4wkwkkss@mail1.ist.utl.pt> References: <20090211144438.vagucohh4wkwkkss@mail1.ist.utl.pt> Message-ID: <49993C18.3010309@blacknest.gov.uk> Dear Mohammed, Guralp data sheets list the poles and zeros in units of Hz, not radians, and usually for velocity or acceleration as appropriate for the instrument (see the details on the actual data sheet). For use in SAC you have to convert them to displacement in radians i.e. x2pi, then add one zero at (0,0) if they are in velocity or two if they are in acceleration. Regards, Sheila From Onur.Tan at mam.gov.tr Mon Feb 16 06:32:46 2009 From: Onur.Tan at mam.gov.tr (Onur Tan) Date: Mon, 16 Feb 2009 16:32:46 +0200 Subject: [SAC-HELP] How to remove the instrument response? In-Reply-To: <49993C18.3010309@blacknest.gov.uk> Message-ID: <9E88105FE4AC8C47BD6D00EB9831482F0C44B2D2@posta.mam.gov.tr> Dear Mohammed, You must also calculate the CONSTANT from the Guralp's A0 value . Read this: www.le.ac.uk/seis-uk/downloads/seisuk_data_manage.pdf best regards Onur Tan Dr. Onur TAN |------------------------------------------------------------ ------------- | | T?B?TAK Marmara Ara?t?rma Merkezi | TUBITAK Marmara Research Center | | Yer ve Deniz Bilimleri Enstit?s? | Earth and Marine Sciences Institute | | Gebze - Kocaeli | Gebze - Kocaeli - TURKEY | -------------------------------------------------------------------------- -----Original Message----- From: sac-help-bounces at iris.washington.edu [mailto:sac-help-bounces at iris.washington.edu] On Behalf Of Sheila Peacock Sent: Monday, February 16, 2009 12:13 PM To: Mohamed Kamaleldin Abdel Galeil Salah Cc: sac-help at iris.washington.edu Subject: Re: [SAC-HELP] How to remove the instrument response? Dear Mohammed, Guralp data sheets list the poles and zeros in units of Hz, not radians, and usually for velocity or acceleration as appropriate for the instrument (see the details on the actual data sheet). For use in SAC you have to convert them to displacement in radians i.e. x2pi, then add one zero at (0,0) if they are in velocity or two if they are in acceleration. Regards, Sheila _______________________________________________ sac-help mailing list sac-help at iris.washington.edu http://www.iris.washington.edu/mailman/listinfo/sac-help From grant.mott at gmail.com Fri Feb 20 10:43:07 2009 From: grant.mott at gmail.com (Grant Mott) Date: Fri, 20 Feb 2009 11:43:07 -0700 Subject: [SAC-HELP] Units and magnitude for dependent variables in SAC files Message-ID: <8b5bc2250902201043y566c3794yf197525d7d332a5@mail.gmail.com> I am a VERY new user for the SAC data files and working on a window system. Currently, I'm using a MATLAB function (rsac) to convert sac files to a usable data format (3 columns: time, dependent variable, header information). I've been able to figure out the time units used, but am having trouble deciphering the units, magnitude, and form of the dependent variable. How can I determine if the dependent variable is acceleration, displacement, velocity, etc? Also, is there a way that I can determine the magnitude of the dependent variable (i.e. x 10^-3, etc). Thanks for your help! Grant -------------- next part -------------- An HTML attachment was scrubbed... URL: From mayenew at gmail.com Mon Feb 23 10:33:58 2009 From: mayenew at gmail.com (Melaku Ayenew) Date: Mon, 23 Feb 2009 11:33:58 -0700 Subject: [SAC-HELP] Commands in sac for period measurement Message-ID: hi all, I am new to sac i was wandering if somebody can tell me What comand(s) in sac allows to measure the period for ms and ml amplitude . And how do I get the amplitude for mL and ms estimates. melaku New Mexico State University -------------- next part -------------- An HTML attachment was scrubbed... URL: From invite+pmd1imw_ at facebookmail.com Tue Feb 24 07:21:49 2009 From: invite+pmd1imw_ at facebookmail.com (Omar Bin Afif) Date: Tue, 24 Feb 2009 07:21:49 -0800 Subject: [SAC-HELP] Check out my Facebook profile Message-ID: <827534d8cd6b59daf9a27b5b2e79f784@localhost.localdomain> Hi Sac, I set up a Facebook profile where I can post my pictures, videos and events and I want to add you as a friend so you can see it. First, you need to join Facebook! Once you join, you can also create your own profile. Thanks, Omar To sign up for Facebook, follow the link below: http://www.facebook.com/p.php?i=1128824526&k=541X536RUYVM5FCHXC3TSS&r This e-mail may contain promotional materials. If you do not wish to receive future commercial mailings from Facebook, please click on the link below. Facebook's offices are located at 156 University Ave., Palo Alto, CA 94301. http://www.facebook.com/o.php?k=22fdd6&u=1478913827 -------------- next part -------------- An HTML attachment was scrubbed... URL: From nolet at geoazur.unice.fr Thu Feb 26 03:32:28 2009 From: nolet at geoazur.unice.fr (Guust Nolet) Date: Thu, 26 Feb 2009 12:32:28 +0100 Subject: [SAC-HELP] sacio.a on 64 bit machines Message-ID: <58F6693E-4480-4675-8620-1F0FB59E056F@geoazur.unice.fr> Can Brian or anyone else tell me what the plans are to provide the sacio.a library for 64 bit machines? This seems an urgent matter to me because most of us have tons of software that we cannot recompile with a 32 bit flag without getting into messy problems. best wishes, Guust nolet at geoazur.unice.fr Geosciences Azur 250, Rue Albert Einstein 06560 Sophia Antipolis France +33.4.92.94.26.32 http://geoweb.princeton.edu/people/faculty/nolet/ From savage at uri.edu Thu Feb 26 07:16:46 2009 From: savage at uri.edu (Brian Savage) Date: Thu, 26 Feb 2009 10:16:46 -0500 Subject: [SAC-HELP] sacio.a on 64 bit machines In-Reply-To: <58F6693E-4480-4675-8620-1F0FB59E056F@geoazur.unice.fr> References: <58F6693E-4480-4675-8620-1F0FB59E056F@geoazur.unice.fr> Message-ID: <9AC6BE00-7560-4146-A420-832019F43E13@uri.edu> Guust, We are planning to create a 64 bit friendly version of sacio.a soon. It is a large undertaking due to the history and age of SAC. Cheers Brian On Feb 26, 2009, at 6:32 AM , Guust Nolet wrote: > Can Brian or anyone else tell me what the plans are to provide the > sacio.a library for 64 bit machines? This seems an urgent matter to > me because most of us have tons of software that we cannot > recompile with a 32 bit flag without getting into messy problems. > best wishes, > Guust > nolet at geoazur.unice.fr > > Geosciences Azur > 250, Rue Albert Einstein > 06560 Sophia Antipolis > France > +33.4.92.94.26.32 > > http://geoweb.princeton.edu/people/faculty/nolet/ > > > > _______________________________________________ > 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 Thu Feb 26 13:16:03 2009 From: m.tchelo at gmail.com (Marcelo Bianchi) Date: Thu, 26 Feb 2009 18:16:03 -0300 Subject: [SAC-HELP] sacio.a on 64 bit machines In-Reply-To: <9AC6BE00-7560-4146-A420-832019F43E13@uri.edu> References: <58F6693E-4480-4675-8620-1F0FB59E056F@geoazur.unice.fr> <9AC6BE00-7560-4146-A420-832019F43E13@uri.edu> Message-ID: <89b338c10902261316i405598e3ic8b5299e273fded9@mail.gmail.com> As a quick and dirty solution maybe you could try to use the -m32 flag in the compiling time ... don't forgette to install the lib32 version of libraries. []'s marcelob 2009/2/26 Brian Savage : > Guust, > > We are planning to create a 64 bit friendly version of sacio.a soon. > > It is a large undertaking due to the history and age of SAC. > > Cheers > Brian > > > On Feb 26, 2009, at 6:32 AM , Guust Nolet wrote: > >> Can Brian or anyone else tell me what the plans are to provide the sacio.a >> library for 64 bit machines? This seems an urgent matter to me because most >> of us have tons of software that we cannot recompile with a 32 bit flag >> without getting into messy problems. >> best wishes, >> Guust >> nolet at geoazur.unice.fr >> >> Geosciences Azur >> 250, Rue Albert Einstein >> 06560 Sophia Antipolis >> France >> +33.4.92.94.26.32 >> >> http://geoweb.princeton.edu/people/faculty/nolet/ >> >> >> >> _______________________________________________ >> sac-help mailing list >> sac-help at iris.washington.edu >> http://www.iris.washington.edu/mailman/listinfo/sac-help > > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help > -- :: "This is all about life !" :: :: Find me @ http://www.foo4fun.net/ :: From sebastian.heimann at zmaw.de Fri Feb 27 01:45:18 2009 From: sebastian.heimann at zmaw.de (Sebastian Heimann) Date: Fri, 27 Feb 2009 10:45:18 +0100 Subject: [SAC-HELP] sacio.a on 64 bit machines In-Reply-To: <89b338c10902261316i405598e3ic8b5299e273fded9@mail.gmail.com> References: <58F6693E-4480-4675-8620-1F0FB59E056F@geoazur.unice.fr> <9AC6BE00-7560-4146-A420-832019F43E13@uri.edu> <89b338c10902261316i405598e3ic8b5299e273fded9@mail.gmail.com> Message-ID: <9F28FEC3-5652-43ED-AAAD-A8FBBEACE25E@zmaw.de> Hi! Because I had the same Problem, I wrote a little sacio module in Python, which does not depend on sacio.a. Of course this won't solve the general problem of the 32 bit sacio.a, but maybe it is of interest to those who only want to be able to read and write some SAC files from Python. BTW, -m32 also does not solve the problem, because in a large application, which uses lots of libraries, all of these libs would have to be compiled with -m32, which is not always easy. Also -m32 limits the amount of memory the app can address, which prohibits the most interesting applications. Cheers, Sebastian -------------- next part -------------- A non-text attachment was scrubbed... Name: pysacio.py Type: text/x-python-script Size: 9899 bytes Desc: not available URL: -------------- next part -------------- On 26.02.2009, at 22:16, Marcelo Bianchi wrote: > As a quick and dirty solution maybe you could try to use the -m32 flag > in the compiling time ... don't forgette to install the lib32 version > of libraries.