From seismo_003 at rediffmail.com Sun Jan 1 04:44:29 2012 From: seismo_003 at rediffmail.com (Sunil Roy) Date: 1 Jan 2012 12:44:29 -0000 Subject: [SAC-HELP] =?utf-8?q?Now_able_to_write_sac_header_in_C?= Message-ID: <20120101124429.17583.qmail@f4mail-235-247.rediffmail.com> Hi SAC user, I am trying to write the complete sac header for evenly space file using wsac0() in c program. But i am always getting error or segmentation fault. The problem is that I am not able to understand what input it is taking to write the complete header information using wsac0() function. what is the 2nd and 3rd input in wsac0(). The code to write the sac file is given here. What is the problem in the code or I would be thankful for your advise to correct the code. int main(int argc, char *argv[]) {   FILE *fp1;   char **fi;   char fi1[400];   static int nrf=0;   int i;   nrf=count(argv[1]);   fp1=fopen(argv[1],"r");     fscanf(fp1,"%s",fi1);     fi=malloc ( strlen(fi1) + 1 );     strcpy(fi,fi1);   fclose(fp1);   /* Define variables to be used in the call to rsac1() and getfhv() */     static  int max = MAX, nlen, nerr, n1, n2, npts=0;     static float xarray[ MAX ] , yarray [ MAX ], beg , del , delta , b , o, e, stla=0, stlo=0, stel=0, evla,evlo,evdp, mag, cmpaz;   int leven=1;   float minimum=1.0, maximum=6.0;     static char kname[ 31 ], out[31], *kstnm, *kcmpnm;     strcpy(kname,fi[i]);     rsac1(kname, yarray, & nlen, & beg, & del, & max, & nerr, strlen( kname ) ) ;     /* Check the Error status */     if ( nerr != 0 ) {       fprintf(stderr, "Error reading SAC file: %s\n", kname);       exit(-1);     }     newhdr();     setnhv("NPTS", & npts, & nerr, strlen("NPTS"));     setfhv("DELTA", & del, &nerr, strlen("DELTA"));     setfhv("B", & b, & nerr, strlen("B"));     setfhv("E", & e, & nerr, strlen("E"));     setfhv("O", & o, & nerr, strlen("O"));     setfhv("STLA", & stla, & nerr, strlen("STLA"));     setfhv("STLO", & stlo, & nerr, strlen("STLO"));     setfhv("STEL", & stel, & nerr, strlen("STEL"));     setfhv("EVLA", & evla, & nerr, strlen("EVLA"));     setfhv("EVLO", & evlo, & nerr, strlen("EVLO"));     setfhv("EVDP", & evdp, & nerr, strlen("EVDP"));     setkhv("KSTNM", kstnm, & nerr, strlen("KSTNM"),strlen(kstnm));     setfhv("CMPAZ", &cmpaz, &nerr, strlen("CMPAZ"));     setkhv("KCMPNM", kcmpnm, & nerr, strlen("KCMPNM"),strlen(kcmpnm));     setfhv("MAG", & mag, & nerr, strlen("MAG"));       if(nerr != 0) {       fprintf(stderr, "Error in reading SAC File: %s\n", kname);       exit(-1);       }     sprintf(out,"%s.out",fi);     strcpy(kname,out);    wsac0 (kname, b, yarray, &nerr, strlen( out ));     if ( nerr != 0 ) {       fprintf(stderr, "Error getting header variable: delta\n");       exit(-1); }   free(fi);   return 0; } Error Message: xx3.c: In function ?main?: xx3.c:105:4: error: incompatible type for argument 2 of ?wsac0? ../include/sacio.h:104:6: note: expected ?float *? but argument is of type ?float? What will be the input in 2nd and 3rd place of wsac0() for evenly space file Thanking you, Sunil Roy Research Scholar National Geophysical Research Institute Hyderbad, India -------------- next part -------------- An HTML attachment was scrubbed... URL: From snoke at vt.edu Sun Jan 1 15:05:51 2012 From: snoke at vt.edu (Arthur Snoke) Date: Sun, 1 Jan 2012 18:05:51 -0500 (EST) Subject: [SAC-HELP] Now able to write sac header in C In-Reply-To: <20120101124429.17583.qmail@f4mail-235-247.rediffmail.com> References: <20120101124429.17583.qmail@f4mail-235-247.rediffmail.com> Message-ID: I had not answered because my calls to wsac0 have been using Fortran so I hoped that someone would respond. This being a vacation time for much of the world, that has not happened. So I will try. I looked at the help pages and was surprised to find that between v101.3 and 101.4 discussionas about C and Fortran subroutines were left out of the user_manual sections on input_output and appendix. Here is a short section from the old input_output that addresses some of your questions? WSAC0 C: void wsac0(kname, xarray, yarray, nerr, kname_s) char *kname; float *xarray; float *yarray; long int *nerr; long int kname_s; FORTRAN: call wsac0(kname, xarray, yarray, nerr) character*8 kname real*4 xarray(*) real*4 yarray(*) integer*4 nerr purpose: To write a SAC file to disk using current header values. Input variables: kname: Name of disk file to write (padded at the end with blanks or '\0's). xarray: contains the independent variable. yarray: contains the dependent variable. kname_s: (C only) string length of kname Output variables: nerr: error status; always 0 if there was no error. Subroutine wsac0 can write SAC files with even spacing (LEVEN true) in which case xarray is just a dummy file. Or it can write unevenly0spaced or spectral files (LEVEN false). Perhaps you need to include the line setlhv ( "LEVEN" , &true , &nerr , 5 ) ; Let me know if this helps -- or if you need more information. We will look into adding more documentation about subroutines to the HELP file. We cannot change files in the distribution until the next release, but it is easy to alter the online HELP files at IRIS. Arthur From branden.christensen at osop.com.pa Mon Jan 2 12:09:13 2012 From: branden.christensen at osop.com.pa (Branden Christensen) Date: Mon, 2 Jan 2012 15:09:13 -0500 Subject: [SAC-HELP] Employment Opportunity: Electronics Guru and Field System Engineer Message-ID: *Electronics Guru and Field System Engineer, OSOP, Panama* OSOP, located in the highlands of Western Panama, has an opening for a full-time electronics specialist and field system engineer. We are looking for someone who can design circuit boards for geophysical applications, is good with their hands and can install geophysical equipment. Knowledge of geophysics is a plus. *Company Description* OSOP provides vertically integrated solutions and turn-key systems for geophysical monitoring. OSOP develops, supports, installs, maintains and configures instrumentation and software for early warning civil defense systems to mitigate the effects of natural disasters. OSOP promotes a dynamic, fast-paced and entrepreneurial work environment. *Position Description* The successful applicant will be in charge of the OSOP eletronics and field teams, manage network installation and instrument development projects from cradle to grave, have free-range of the electronics laboratory and machine shop, and will be granted the freedom to innovate. The position will require world travel to exotic and often tropical environments. Resilience, independence, adaptability, efficiency and tenacity are desired qualities. Before applying applicants should consider if they can adapt to an international team in a tropical land at the base of a volcano within a county known for its endless beaches. *Essential Job Functions* * Assists in the design, fabrication, testing, documentation and repair of hardware for geophysical projects * Independently designs and installs geophysical networks: logistics, site reconnaissance, site selection, installation, construction, maintenance, and data transmission. * Independently identifies, troubleshoots and solves technical issues related to hardware, software, and data transmission. * Independently trains OSOP field personnel. * Travel domestically or internationally for extended periods, often on short notice. Location: Volcano Baru, Panama Deadline: January 31st, 2012 Perks: Paid flight to Panama, interim housing, company transportation (including use of car on weekends), international travel * * *Qualifications* * Formal education in electronics. * Experience in circuitry and instrument design. * Programming of embedded systems and micro-processors. * Some knowledge of geophysical methodologies and instrumentation, data communications, and power systems. * Proven ability to plan and coordinate projects. * Experience in construction. * Functional English skills. Spanish is helpful but not a requirement. Applicants not gung-ho about pushing the envelope need not apply. *Selection Procedure* Applications will be evaluated on their experience in network and hardware design, capacity for critical thought and lateral thinking skills. *Still interested?* Send a C.V. to Branden Christensen (branden . christensen @ osop . com . pa) and Angel Rodriguez (angel @ volcanbaru . com). Tailor the C.V. to your exposure to and proficiency in: - Electronics and instrument design - Programming of embedded systems and micro-processors. - Geophysical network configuration and management - Volcano, earthquake and/or geophysical institute experience -------------- next part -------------- An HTML attachment was scrubbed... URL: From branden.christensen at osop.com.pa Mon Jan 2 12:10:25 2012 From: branden.christensen at osop.com.pa (Branden Christensen) Date: Mon, 2 Jan 2012 15:10:25 -0500 Subject: [SAC-HELP] Employment Opportunity: Geophysical Software Developer and Programmer Message-ID: *Geophysical Software Developer and Programmer, OSOP, Panama* * * OSOP, located in the highlands of Western Panama, has an opening for a full-time geophysical software developer and programmer. We are looking for an individual who is versed in geophysical signal processing and at-home in the Linux environment. *Company Description* OSOP provides vertically integrated solutions and turn-key systems for geophysical monitoring. OSOP develops, supports, installs, maintains and configures instrumentation and software for early warning civil defense systems to mitigate the effects of natural disasters. OSOP promotes a dynamic, fast-paced and entrepreneurial work environment. * * *Position Description* The successful applicant will be in charge of the OSOP programming team, manage software development projects from cradle to grave, install software for institutions across the globe and be granted the freedom to innovate. Duties will depend on the successful candidate's skill set and are likely to include: software development, installation, configuration and employee training/ supervision. The position will require regular domestic and international travel. Resilience, independence, adaptability, efficiency and tenacity are desired qualities. Before applying applicants should consider if they can adapt to an international team in a tropical land at the base of a volcano within a county known for its endless beaches. Location: Volcano Baru, Panama Deadline: January 31st, 2012 Perks: Paid flight to Panama, interim housing, company transportation (including use of car on weekends), international travel *Qualifications* * High-level programming aptitude in C++ and Python * Experience programming in the Linux environment * Demonstrated understanding of signal-processing * Experience with RDBMS * Either an M.S./ PhD in Geophysics, Mathematics or Physics or advanced programming skills in C++ and Python and a demonstrated interest in solving problems related to signal processing applications. * Functional English skills. Spanish is helpful but not a requirement. * Capable of regular domestic and international travel Applicants not gung-ho about pushing the envelope need not apply. *Selection Procedure* Applications will be evaluated on their proficiency in programming, signal processing, industry-wide software packages, capacity for critical thought and lateral thinking skills. *Still interested?* Send a C.V. to Branden Christensen (branden . christensen @ osop . com . pa) and Angel Rodriguez (angel @ volcanbaru . com). Tailor the C.V. to your exposure to and proficiency in: - Programming languages - UNIX/ Linux - Signal Processing - Geophysics, Physics, Mathematics - Antelope, Earlybird, Earthworm, SEISAN, SeisComp3 - RDBMS - Communication protocols: TCP/IP, UDP - Basic Linux system administration - Geophysical network configuration and management - Volcano, earthquake and/or geophysical institute experience - Electronics and instrument design -------------- next part -------------- An HTML attachment was scrubbed... URL: From branden.christensen at osop.com.pa Mon Jan 2 13:54:29 2012 From: branden.christensen at osop.com.pa (Branden Christensen) Date: Mon, 2 Jan 2012 16:54:29 -0500 Subject: [SAC-HELP] Announcing the First SeisComp Workshop for the Americas Message-ID: *The Observatorio Sismico del Occidente de Panama (OSOP) is pleased to announce the First SeisComp Workshop for the Americas * * * Date: Sunday, June 3rd to Saturday, June 9th Location: OSOP, Volcan Baru, Panama Application Deadline: March 15th, 2012 Slots: 15 (The first 10 participants will receive free room) Cost: 1,900 USD (includes ice breaker, coffee, snacks, lunch, a gala dinner and RUMBA. Transportation not included.) 2,500 after deadline. *Company Description* OSOP develops, supports, installs, maintains and configures instrumentation and software for early warning civil defense systems to mitigate the effects of natural disasters. OSOP's vertically integrated solutions and turn-key SeisComp-Earthworm-Earlybird-Seisan systems are being installed world-wide. *Workshop Description* This short-course is comprehensive and ambitious in its scope: Day -1: Introduction to Linux for Seismologists Day 1: Introduction to SeisComp | Metadata Management (PDCC) | SeisComp Installation Part 1 Day 2: SeisComp Installation Part 2 Day 3: Integrating SeisComp, Earthworm and Seisan Day 4: SWOT analysis | The day-to-day use of SeisComp in observatories | Data Quality Control (PQLX) Day 5: Advanced Topics as decided by the participants | Q&A Session Day 6: Instrument Calibration and Sensor Orientation This course is targeted at people who need to install, integrate and maintain different seismological software, technicians who install field equipment and need to maintain response databases for seismic networks as well as those who are on the fence about installing integrated software packages which include SeisComp. * Lectures will be given in English and Spanish. * Certificates will be given in conjunction with the University of Panama. * Formal invitations from the University of Panama are available upon request. * Instructors: Wilfried Strauss, PhD, Branden Christensen, MSc, Angel Rodriguez, MSc and others. * Classroom: OSOP has moved to a new, larger facility with an educational wing. * Presentations: Please note that you do not have to bring a presentation to the workshop. You may, however, bring a poster of any nature so that you may hash out ideas with staff and participants during down-time. * We can help with transportation arrangements. *How to apply* Send an email to Branden Christensen (branden . christensen @ osop . com . pa) and Angel Rodriguez (a ngel @ volcanbaru . com) expressing your interest. Please be sure to include your full name, phone number(s), title, organization and country of residence. -------------- next part -------------- An HTML attachment was scrubbed... URL: From rob at iris.washington.edu Thu Jan 5 14:09:04 2012 From: rob at iris.washington.edu (Robert Casey) Date: Thu, 5 Jan 2012 14:09:04 -0800 Subject: [SAC-HELP] Fwd: SAC Software Request v101.5 -Linux 64 References: <3E07ED2C-BB00-4031-B608-B6F0A90A563F@iris.washington.edu> Message-ID: <74D7F4E3-0906-4BDC-AA29-14D1897C0B96@iris.washington.edu> Hi Margaret- This question gets to the heart of the SAC builds themselves, which we do not have control over. We only broker the distribution. I need to forward your question to the SAC HELP list and see if the list members can help you address this. Please consider subscribing to the SAC Help list through the IRIS website where you're sure to get the best expert answers to your SAC questions. -Rob > > Begin forwarded message: > > From: "Doll, Margaret Ann" > Date: January 5, 2012 1:22:07 PM PST > To: "Gale J. Cox" > Subject: Re: SAC Software Request v101.5 -Linux 64 > > Gale, > > Thanks for your quick response. > > The version of sac that you sent me needs GLIBC_2.7 or 2.11. My > system is using glibc-2.5; I have a ROCKS 5.4.3 cluster based on > CentOS release 5.6 (Final). Would you send me a version of sac that > would run using the glibc2.5 library. > > Thanks. > > By the way in the README for sac-101.5, you have > > For bash, edit ~/.bashrc adding the lines > export SACHOME=/usr/local/sac > source ${SACHOME}/bin/sacinit.csh > > > The source line should use sacinit.sh >> >> On Thu, Jan 5, 2012 at 1:08 PM, Gale J. Cox >> wrote: >> Attached is the binary compressed tar file of the SAC v101.5 >> software for >> Linux (64 bit): sac-101.5-linux_x86_64.tar.gz. The first step is >> to expand >> that file. If your version of tar has the "z" option, you can do >> this in >> one step: >> [....snip.....] > >> If you need help ... >> >> 1. If you have SAC working, use the help command from within SAC. >> You can >> print hard-copies of the help files using printhelp. HTML >> versions of the >> help files can be found at URL > >. >> 2. Some features of SAC are described in the CHANGES file in this >> directory, >> and changes in new versions are summarized in HISTORY. >> 3. URL has links to other >> resources and >> information about updates. >> 4. A place to share problems and find solutions is the sac-help >> mailing >> list: . You must be a subscriber >> to post to this listserv. For further information, go to >> URL . >> You do not have to be a subscriber to access the list archives. >> > -------------- next part -------------- An HTML attachment was scrubbed... URL: From wangwtustc at gmail.com Sat Jan 7 23:49:43 2012 From: wangwtustc at gmail.com (weitao wang) Date: Sun, 8 Jan 2012 15:49:43 +0800 Subject: [SAC-HELP] sac101.5 32bit problem: can not change(update) header using sac101.5 Message-ID: Hi everyone, I recently update sac into sac 101.5 on centos6.2 and centos 5.5, both 32bit and encountered problems while changing headers using "ch ", as described bellow. cat <test.sh funcgen seis w demo.sac rh demo.sac ch t1 12.0 wh r demo.sac lh t1 q EOF sac < test.sh As shown in my linux, the header t1 do not exist in demo.sac and moreover it seems the waveform has been changed. Is this a bug or it is due to the setup of my linux env? Any advices ? thanks! PS. I use the pre-compiled binary distribution of sac101.32bit. wangwt -------------- next part -------------- An HTML attachment was scrubbed... URL: From fabioludias at gmail.com Sun Jan 8 10:57:24 2012 From: fabioludias at gmail.com (=?ISO-8859-1?Q?F=E1bio_Dias?=) Date: Sun, 8 Jan 2012 16:57:24 -0200 Subject: [SAC-HELP] sac101.5 32bit problem: can not change(update) header using sac101.5 In-Reply-To: References: Message-ID: I have the same problem using 'ch' version 101.5. I came back to version 101.4. Fabio 2012/1/8 weitao wang > Hi everyone, > > I recently update sac into sac 101.5 on centos6.2 and centos 5.5, both > 32bit and encountered problems while changing headers using "ch ", as > described bellow. > > cat <test.sh > funcgen seis > w demo.sac > rh demo.sac > ch t1 12.0 > wh > r demo.sac > lh t1 > q > EOF > > sac < test.sh > > As shown in my linux, the header t1 do not exist in demo.sac and moreover > it seems the waveform has been changed. > > Is this a bug or it is due to the setup of my linux env? > > Any advices ? thanks! > > PS. I use the pre-compiled binary distribution of sac101.32bit. > > wangwt > > _______________________________________________ > 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 savage at uri.edu Sun Jan 8 12:20:29 2012 From: savage at uri.edu (Brian Savage) Date: Sun, 8 Jan 2012 15:20:29 -0500 Subject: [SAC-HELP] sac101.5 32bit problem: can not change(update) header using sac101.5 In-Reply-To: References: Message-ID: Dear wangwt and Fabio, There is a bug in sac version 101.5 command to write headers. It only affects the command writehdr (wh). We should be able to create and upload to IRIS a patched version of version 101.5 in the coming days (week). In the meantime, your options are - Version 101.4 - Using the read/write command sequence instead of readhdr/writehdr If you compiled the binary yourself by requesting the source file from IRIS, the attached patch file can be applied to remedy the writehdr bug. wangwt, Thank you for submitting the bug report and a simple command sequence It was quite helpful to identifying source of the the bug. Thanks again. B -------------- next part -------------- A non-text attachment was scrubbed... Name: write_header.patch Type: application/octet-stream Size: 305 bytes Desc: not available URL: -------------- next part -------------- On Jan 8, 2012, at 1:57 PM, F?bio Dias wrote: > I have the same problem using 'ch' version 101.5. I came back to version 101.4. > > Fabio > > 2012/1/8 weitao wang > Hi everyone, > > I recently update sac into sac 101.5 on centos6.2 and centos 5.5, both 32bit and encountered problems while changing headers using "ch ", as described bellow. > > cat <test.sh > funcgen seis > w demo.sac > rh demo.sac > ch t1 12.0 > wh > r demo.sac > lh t1 > q > EOF > > sac < test.sh > > As shown in my linux, the header t1 do not exist in demo.sac and moreover it seems the waveform has been changed. > > Is this a bug or it is due to the setup of my linux env? > > Any advices ? thanks! > > PS. I use the pre-compiled binary distribution of sac101.32bit. > > wangwt > > _______________________________________________ > 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 From jattanayake at gmail.com Mon Jan 16 19:50:44 2012 From: jattanayake at gmail.com (Januka Attanayake) Date: Mon, 16 Jan 2012 22:50:44 -0500 Subject: [SAC-HELP] GCARC header differences Message-ID: Hello, I have converted velocity seismograms to displacement and I see that header GCARC are not the same. I haven't checked other headers at this time. In one case, the difference was ~7 degrees. Is there a reason for this? How do I mitigate this problem. Thanks for your time! Januka. -------------- next part -------------- An HTML attachment was scrubbed... URL: From C.Holden at gns.cri.nz Thu Jan 19 20:44:02 2012 From: C.Holden at gns.cri.nz (Caroline Holden) Date: Fri, 20 Jan 2012 17:44:02 +1300 Subject: [SAC-HELP] polezero file for Milne-Shaw historical instrument Message-ID: Hello, could you please let me know what would be the sac polezero file for a Milne-Shaw instrument with the following characteristics: Seismometer period Ts=12 Damping ratio 20:1 Magnification 250 Thank you very much for your help! Caroline PS: my guess is: Zero 2 poles 2 -0.3613 0.3790 -0.3613 0.3790 Constant 250*2*pi Notice: This email and any attachments are confidential. If received in error please destroy and immediately notify us. Do not copy or disclose the contents. -------------- next part -------------- An HTML attachment was scrubbed... URL: From sgao at mst.edu Sat Jan 21 05:43:15 2012 From: sgao at mst.edu (Gao, Stephen S.) Date: Sat, 21 Jan 2012 07:43:15 -0600 Subject: [SAC-HELP] lost headers? Message-ID: <4012DAA935478B428921E410B28D9BA7584437@mst-vmail02.srv.mst.edu> Hi, My understanding was that wsac1 keeps all the header values from the most recent rsac1 call? However, in the following example, file tmp1.z lost all its event/station coordinates and related headers. Any hints? Thanks. Steve Gao ------------------------------------ nin1='tmp.z' call rsac1(nin1, x, npts, beg, delta, max, nerr) nin1='tmp1.z' call wsac1(nin1, x, npts, beg, delta, nerr) From snoke at vt.edu Sat Jan 21 17:51:47 2012 From: snoke at vt.edu (Arthur Snoke) Date: Sat, 21 Jan 2012 20:51:47 -0500 (EST) Subject: [SAC-HELP] lost headers? In-Reply-To: <4012DAA935478B428921E410B28D9BA7584437@mst-vmail02.srv.mst.edu> References: <4012DAA935478B428921E410B28D9BA7584437@mst-vmail02.srv.mst.edu> Message-ID: > My understanding was that wsac1 keeps all the header values from the > most recent rsac1 call? However, in the following example, file tmp1.z > lost all its event/station coordinates and related headers. No. wsac1 and wsac2 write "minimum headers." wsac0 writes full headers and includes updates to header values for a SAC file that had been read in and then processed. We found that the help command for the last couple of versions was incomplete so did not explain and give examples for the various read and write routines. That will be corrected in v101.5b -- coming soon. WSAC1 writes evenly spaced files WSAC2 writes unevenly spaced and spectral files WSAC0 writes either format but has more comprehensive header files than the other two WSAC1 and WSAC2 write SAC files with a minimum header contains only those variables needed to be able to read the file: B, E, DELTA, LEVEN, and NPTS. For calls to WSAC0, If it is a new file, it requires a call to subroutine NEWHDR supplemented by additional header variables to be set using the SETXXX routines (see examples below). If it is writing to a file that is based on one that had been read in previously in the program, one should not call NEWHDR. Before writing such a file using WSAC0, SAC updates the header variables such as DEPMAX and BAZ. As shown in the examples below, the type of SAC data file that gets written depends on header variables that must be set: IFTYPE and LEVEN. IFTYPE has the following values: From sheila at blacknest.gov.uk Mon Jan 23 03:18:12 2012 From: sheila at blacknest.gov.uk (Sheila Peacock) Date: Mon, 23 Jan 2012 11:18:12 +0000 Subject: [SAC-HELP] pole-zero for Milne-Shaw seismometer Message-ID: <4F1D41F4.7010703@blacknest.gov.uk> Caroline Holden wrote, on 20th January 2012: Hello, could you please let me know what would be the sac polezero file for a Milne-Shaw instrument with the following characteristics: Seismometer period Ts=12 Damping ratio 20:1 Magnification 250 Thank you very much for your help! Caroline PS: my guess is: Zero 2 poles 2 -0.3613 0.3790 -0.3613 0.3790 Constant 250*2*pi I looked up in Scherbaum ("Of Poles and Zeros", Kluwer, 2001) and found the following. First of all, the second pole must be the complex conjugate of the first, so the second 0.3790 should be negative. Secondly, I'm not sure how the "damping ratio" above fits into the equations on p. 57 of Scherbaum. Does it correspond to the amplitude ratio of two consecutive maxima, or of two consecutive extreme values i.e. maximum and the following minimum? The logarithmic decrement, Lambda, is given by ln (ratio of two consecutive maxima), or 2 x ln (ratio of two consecutive extremes). On p. 64 of Scherbaum the poles for an underdamped system are given by p1,2 = -(h +/- i sqrt (1-h^2)) x omega0 where h is the damping coefficient and omega0 is given by 2 x pi / T0, and T0 is your 12 s. On page 57 the logarithmic decrement Lambda is given as Lambda = 2 x pi x h / sqrt(1-h^2). I reproduced your two poles using these formulae and with the second definition of Lambda, i.e. that your "damping ratio" 20:1 is the ratio of two consecutive extremes of the response. I got a damping coefficient h = 0.69 which is close to the value of 0.7 said by my senior colleague to be favoured by operators of this type of seismometer. The two zeros are necessary to convert the response in acceleration to displacement. If the magnification of 250 is the magnification of displacement then the constant required by SAC is equal to the ratio of the zeros polynomial to the poles polynomial at 12 s (Scherbaum p. 40), multiplied by the magnification. Here it is 1.38x250 = 345.0 (I use a program for doing the polynomial ratio - it needs the ability to handle complex numbers, which is in fortran or Matlab). There should be a further factor introduced by the digitiser, with units of digital counts per nanometre (or similar), which will need to be incorporated into the constant. I hope that's some help - and that it's correct, which is why I'm sending it to the list, so that others can tell me if it isn't. I find it very easy to get mixed up about whether to multiply or divide the constant by the magnification. Regards, Sheila Peacock AWE Blacknest. From C.Holden at gns.cri.nz Mon Jan 23 19:17:47 2012 From: C.Holden at gns.cri.nz (Caroline Holden) Date: Tue, 24 Jan 2012 16:17:47 +1300 Subject: [SAC-HELP] pole-zero for Milne-Shaw seismometer In-Reply-To: <4F1D41F4.7010703@blacknest.gov.uk> Message-ID: Thank you very much Sheila, that's very clear now, only one last point, regarding the constant, I read somewhere that we need to multiply it by 2*pi to follow sac conventions ? A0 then becomes 1.38*250*2*pi Kind regards, Caroline ------------------------------------------ GNS Science - Te Pu Ao 1 Fairway Drive, PO Box 30368, Lower Hutt, New Zealand Web: www.gns.cri.nz Email: c.holden at gns.cri.nz Sheila Peacock Sent by: sac-help-bounces at iris.washington.edu 24/01/2012 00:34 To "sac-help at iris.washington.edu" cc Subject Re: [SAC-HELP] pole-zero for Milne-Shaw seismometer Caroline Holden wrote, on 20th January 2012: Hello, could you please let me know what would be the sac polezero file for a Milne-Shaw instrument with the following characteristics: Seismometer period Ts=12 Damping ratio 20:1 Magnification 250 Thank you very much for your help! Caroline PS: my guess is: Zero 2 poles 2 -0.3613 0.3790 -0.3613 0.3790 Constant 250*2*pi I looked up in Scherbaum ("Of Poles and Zeros", Kluwer, 2001) and found the following. First of all, the second pole must be the complex conjugate of the first, so the second 0.3790 should be negative. Secondly, I'm not sure how the "damping ratio" above fits into the equations on p. 57 of Scherbaum. Does it correspond to the amplitude ratio of two consecutive maxima, or of two consecutive extreme values i.e. maximum and the following minimum? The logarithmic decrement, Lambda, is given by ln (ratio of two consecutive maxima), or 2 x ln (ratio of two consecutive extremes). On p. 64 of Scherbaum the poles for an underdamped system are given by p1,2 = -(h +/- i sqrt (1-h^2)) x omega0 where h is the damping coefficient and omega0 is given by 2 x pi / T0, and T0 is your 12 s. On page 57 the logarithmic decrement Lambda is given as Lambda = 2 x pi x h / sqrt(1-h^2). I reproduced your two poles using these formulae and with the second definition of Lambda, i.e. that your "damping ratio" 20:1 is the ratio of two consecutive extremes of the response. I got a damping coefficient h = 0.69 which is close to the value of 0.7 said by my senior colleague to be favoured by operators of this type of seismometer. The two zeros are necessary to convert the response in acceleration to displacement. If the magnification of 250 is the magnification of displacement then the constant required by SAC is equal to the ratio of the zeros polynomial to the poles polynomial at 12 s (Scherbaum p. 40), multiplied by the magnification. Here it is 1.38x250 = 345.0 (I use a program for doing the polynomial ratio - it needs the ability to handle complex numbers, which is in fortran or Matlab). There should be a further factor introduced by the digitiser, with units of digital counts per nanometre (or similar), which will need to be incorporated into the constant. I hope that's some help - and that it's correct, which is why I'm sending it to the list, so that others can tell me if it isn't. I find it very easy to get mixed up about whether to multiply or divide the constant by the magnification. Regards, Sheila Peacock AWE Blacknest. _______________________________________________ sac-help mailing list sac-help at iris.washington.edu http://www.iris.washington.edu/mailman/listinfo/sac-help Notice: This email and any attachments are confidential. If received in error please destroy and immediately notify us. Do not copy or disclose the contents. -------------- next part -------------- An HTML attachment was scrubbed... URL: From jar at email.unc.edu Tue Jan 24 14:01:18 2012 From: jar at email.unc.edu (jose rial) Date: Tue, 24 Jan 2012 17:01:18 -0500 Subject: [SAC-HELP] sac problems with MacOSx 10.7 Message-ID: <4F1F2A2E.108@email.unc.edu> Hi; I'm running SAC on Mac OSX 10.7.2 and thought it installed correctly and seems to work OK, it crashes a lot. It crashes sometimes as in the example below after seemingly routine commands, like read and plot. Sometimes it crashes for a segmentation fault, without apparent motive, just when plot is invoked. SAC> lp co .001 n 4 p 2 SAC> p SAC> xlim 8000 20000 SAC> p SAC> ls 159.18E 72.82S.txt TALDICE_18O_TALDICE1_8-25ka.txt HP [Converted].jpg cucu.sac HP.ai taldice.sac Stenni et al_Nature Geoscience_2011.pdf untitled [Converted].jpg TALDICE-d18Orecord_8-25kyr.xls untitled.ai SAC> r taldice.sac SAC> p Assertion failed: (dpy->xcb->pending_requests_tail->next == NULL), function append_pending_request, file xcb_io.c, line 129. Other times it crashes due to an apparent segmentation fault. If sac is launched again it fails right off the bat, so I have to restart my computer in order for it to work for another, say half hour, until it crashes again... Does anyone have similar problems? Thanks, JRial -------------- next part -------------- A non-text attachment was scrubbed... Name: jar.vcf Type: text/x-vcard Size: 273 bytes Desc: not available URL: From snoke at vt.edu Tue Jan 24 14:32:11 2012 From: snoke at vt.edu (Arthur Snoke) Date: Tue, 24 Jan 2012 17:32:11 -0500 (EST) Subject: [SAC-HELP] sac problems with MacOSx 10.7 In-Reply-To: <4F1F2A2E.108@email.unc.edu> References: <4F1F2A2E.108@email.unc.edu> Message-ID: Alas ... the developers do not have access to a Mac OS 10.7 computer so were not able to test it. We assumed -- perhaps erroneously -- that the changes from 10.6 would not affect SAC. I am assuming you using the binary. True? If so, could you try building from sources and let me know. Within a day or so v101.5b will be out, so I suggest you wait -- but request the source distribution. On Tue, 24 Jan 2012, jose rial wrote: > Hi; > I'm running SAC on Mac OSX 10.7.2 and thought it installed correctly and > seems to work OK, it crashes a lot. It > crashes sometimes as in the example below after seemingly routine commands, > like read and plot. Sometimes it crashes for > a segmentation fault, without apparent motive, just when plot is invoked. > > > > SAC> lp co .001 n 4 p 2 > SAC> p > SAC> xlim 8000 20000 > SAC> p > SAC> ls > 159.18E 72.82S.txt TALDICE_18O_TALDICE1_8-25ka.txt > HP [Converted].jpg cucu.sac > HP.ai taldice.sac > Stenni et al_Nature Geoscience_2011.pdf untitled [Converted].jpg > TALDICE-d18Orecord_8-25kyr.xls untitled.ai > SAC> r taldice.sac > SAC> p > Assertion failed: (dpy->xcb->pending_requests_tail->next == NULL), function > append_pending_request, file xcb_io.c, line 129. > > Other times it crashes due to an apparent segmentation fault. > > If sac is launched again it fails right off the bat, so I have to restart my > computer in order for it to work for another, say half hour, until it crashes > again... > > Does anyone have similar problems? > > Thanks, > > JRial > From snoke at vt.edu Wed Jan 25 12:13:49 2012 From: snoke at vt.edu (Arthur Snoke) Date: Wed, 25 Jan 2012 15:13:49 -0500 (EST) Subject: [SAC-HELP] sac problems with MacOSx 10.7 Message-ID: Jose, I got temporary access to a OS10.7.2 machine but had no problems. I had sent out a "beta" version of v101.5, and I think someone wrote back that it had worked okay for them on OS 10.7.2. Still, it is not working for you, so we need to figure out the reason. Please repeat the steps I do below and make a file of the screen interactions and send it to me. Important questions are the version of gcc you are using. (The one here is the defalt that comes with xcode -- which was put on this morning.) Also, send me (off list) your SAC file so I can try truly to reproduce your tests. My tar file: sac-101.5b-mac_64.tar.gz so my results are for 101.5b, but I would be amazed if the results are any different from those for 101.5 as only one change was made in the source code. Most of the changes are in the help files. ==================================================== uname -a Darwin kingmac.local 11.2.0 Darwin Kernel Version 11.2.0: Tue Aug 9 20:56:15 PDT 2011; root:xnu-1699.24.8~1/RELEASE_I386 i386 [gs4039:~/Desktop/jan25/tests] arthur% gcc --version i686-apple-darwin11-llvm-gcc-4.2 (GCC) 4.2.1 (Based on Apple Inc. build 5658) (LLVM build 2336.1.00) Copyright (C) 2007 Free Software Foundation, Inc. This is free software; see the source for copying conditions. There is NO warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. sac SEISMIC ANALYSIS CODE [01/13/2012 (Version 101.5b)] Copyright 1995 Regents of the University of California SAC> fg seismo SAC> p SAC> quit [gs4039:~/Desktop/jan25/tests] arthur% sac SEISMIC ANALYSIS CODE [01/13/2012 (Version 101.5b)] Copyright 1995 Regents of the University of California SAC> fg seismo SAC> p1 SAC> save seismo.png save file seismo.png [PNG] SAC> save seismo.ps save file seismo.ps [PS] SAC> save seismo.pdf save file seismo.pdf [PDF]which Xcode SAC> quit [gs4039:~/Desktop/jan25] arthur% otool -L sac/bin/sac sac/bin/sac: /usr/X11/lib/libSM.6.dylib (compatibility version 7.0.0, current version 7.0.0) /usr/X11/lib/libICE.6.dylib (compatibility version 10.0.0, current version 10.0.0) /usr/X11/lib/libXpm.4.dylib (compatibility version 16.0.0, current version 16.0.0) /usr/lib/libz.1.dylib (compatibility version 1.0.0, current version 1.2.3) /usr/X11/lib/libX11.6.dylib (compatibility version 9.0.0, current version 9.0.0) /usr/lib/libSystem.B.dylib (compatibility version 1.0.0, current version 125.2.11) /usr/lib/libncurses.5.4.dylib (compatibility version 5.4.0, current version 5.4.0) [gs4039:~/Desktop/jan25] arthur% strings sac/bin/sac | grep version Header version number is incorrect. sgf_to_ps: No SGF file for conversion from sgf_to_ps: No PS file for conversion to Matlab routines in the default version of SAC are not available VARS version number is incorrect. Unable to open file (%s) for version check. Failed to determine version number! This version of SAC does not support database queries. libpng version 1.2.44 - June 26, 2010 libpng version 1.2.44 - June 26, 2010 zlib failed to initialize compressor -- version error Incompatible libpng version in application and library [gs4039:~/Desktop/jan25] arthur% ============================================= I will lose access to this machine tomorrow afternoon (the system was built for a student). Arthur On Tue, 24 Jan 2012, jose rial wrote: > Hi; > > I'm running SAC on Mac OSX 10.7.2 and thought it installed correctly and > seems to work OK, it crashes a lot. It crashes sometimes as in the > example below after seemingly routine commands, like read and plot. > Sometimes it crashes for a segmentation fault, without apparent motive, > just when plot is invoked. > > > > SAC> lp co .001 n 4 p 2 > SAC> p > SAC> xlim 8000 20000 > SAC> p > SAC> ls > 159.18E 72.82S.txt TALDICE_18O_TALDICE1_8-25ka.txt > HP [Converted].jpg cucu.sac > HP.ai taldice.sac > Stenni et al_Nature Geoscience_2011.pdf untitled [Converted].jpg > TALDICE-d18Orecord_8-25kyr.xls untitled.ai > SAC> r taldice.sac > SAC> p > Assertion failed: (dpy->xcb->pending_requests_tail->next == NULL), > function append_pending_request, file xcb_io.c, line 129. > > Other times it crashes due to an apparent segmentation fault. > > If sac is launched again it fails right off the bat, so I have to > restart my computer in order for it to work for another, say half hour, > until it crashes again... > > Does anyone have similar problems? > > Thanks, > > JRial From jar at email.unc.edu Thu Jan 26 08:27:42 2012 From: jar at email.unc.edu (Jose Rial) Date: Thu, 26 Jan 2012 11:27:42 -0500 Subject: [SAC-HELP] problems with SAC in Lion Message-ID: <4F217EFE.8000209@email.unc.edu> Has anyone had problems running SAC on Lion? I have, plenty... Lately, it has been running but I know it'll crash either when I want to plot anything using p. I get tons of unrequested comments as I am processing data, such as setcolor3:0 0 [default] This appears after almost every command. Then, it suddenly, and for no apparent reason, crashes, like this SAC> r curt curtnoint SAC> p2 SAC> r curt SAC> p sac(1081) malloc: *** error for object 0x1008107a0: incorrect checksum for freed object - object was probably modified after being freed. *** set a breakpoint in malloc_error_break to debug Abort Also, why am getting this annoying message about aliasing no matter whether there is potential or not? SAC> interpolate delta 10 WARNING potential for aliasing. new delta: 10.000000 data delta: -12345.000000 I have installed Xcode, and am running sac 101.5b on a MacBook using MacOS 10.7.2 Thanks for the help! Very frustrating exercise.. JRial -------------- next part -------------- A non-text attachment was scrubbed... Name: jar.vcf Type: text/x-vcard Size: 222 bytes Desc: not available URL: From maytalsadeh at post.tau.ac.il Sun Jan 29 07:17:23 2012 From: maytalsadeh at post.tau.ac.il (maytalsadeh at post.tau.ac.il) Date: Sun, 29 Jan 2012 17:17:23 +0200 Subject: [SAC-HELP] (no subject) Message-ID: <20120129171723.9851947qqs9uc7ib@webmail.tau.ac.il> Hi, I am very new to sac (and a bit to seismology in general) and have started experimenting with sac. one of my header lines is: IDEP = JSOP I didn't find any reference for JSOP. Any ideas what that could be? Thanks From kadmiel at bgu.ac.il Sun Jan 29 08:08:05 2012 From: kadmiel at bgu.ac.il (Shahar Shani-Kadmiel) Date: Sun, 29 Jan 2012 18:08:05 +0200 Subject: [SAC-HELP] (no subject) In-Reply-To: <20120129171723.9851947qqs9uc7ib@webmail.tau.ac.il> References: <20120129171723.9851947qqs9uc7ib@webmail.tau.ac.il> Message-ID: <155AAEF2-1FE4-4E34-B14D-E45666D0BB94@bgu.ac.il> Hi Tooly, The only thing that pops to mind is "Joint Seismic Observation Program" but I cannot imagine the connection. I suggest you take a look here- http://www.iris.edu/software/sac/manual/file_format.html for a description of header fields. Let us know if you find anything. On Jan 29, 2012, at 5:17 PM, maytalsadeh at post.tau.ac.il wrote: > Hi, > > I am very new to sac (and a bit to seismology in general) and have started experimenting with sac. > > one of my header lines is: > > IDEP = JSOP > > I didn't find any reference for JSOP. > > Any ideas what that could be? > > > Thanks > > _______________________________________________ > 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 snoke at vt.edu Sun Jan 29 11:42:59 2012 From: snoke at vt.edu (Arthur Snoke) Date: Sun, 29 Jan 2012 14:42:59 -0500 (EST) Subject: [SAC-HELP] (no subject) In-Reply-To: <20120129171723.9851947qqs9uc7ib@webmail.tau.ac.il> References: <20120129171723.9851947qqs9uc7ib@webmail.tau.ac.il> Message-ID: > one of my header lines is: > > IDEP = JSOP > > I didn't find any reference for JSOP. > > Any ideas what that could be? Could it be a typo? >From DEP I Type of dependent variable: IUNKN (Unknown) IDISP (Displacement in nm) IVEL (Velocity in nm/sec) IVOLTS (Velocity in volts) IACC (Acceleration in nm/sec/sec) From Onur.Tan at mam.gov.tr Sun Jan 29 12:19:48 2012 From: Onur.Tan at mam.gov.tr (Onur Tan) Date: Sun, 29 Jan 2012 22:19:48 +0200 Subject: [SAC-HELP] (no subject) In-Reply-To: Message-ID: I've search it in google, and found Joint Seismic Observation Period (JSOP) ??? this is may be a earthquake data set about Africa and/or Mediterranean 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 29.01.2012 21:42, "Arthur Snoke" wrote: >> one of my header lines is: >> >> IDEP = JSOP >> >> I didn't find any reference for JSOP. >> >> Any ideas what that could be? > >Could it be a typo? > >>From > >DEP I >Type of dependent variable: > IUNKN (Unknown) > IDISP (Displacement in nm) > IVEL (Velocity in nm/sec) > IVOLTS (Velocity in volts) > IACC (Acceleration in nm/sec/sec) >_______________________________________________ >sac-help mailing list >sac-help at iris.washington.edu >http://www.iris.washington.edu/mailman/listinfo/sac-help Bu e-posta mesaj? ve onunla iletilen t?m ekler g?nderildi?i ki?i ya da kuruma ?zel olup, gizli ve ?zel bilgiler i?erebilece?i gibi gizlilik y?k?ml?l??? de ta??yor olabilir. Bu mesajda ve ekindeki dosyalarda bulunan t?m fikir ve g?r??ler sadece adres sahip(ler)ine ait olup, T?B?TAK MAM bu e-posta i?eri?indeki bilgilerin kullan?lmas? nedeniyle hi? kimseye kar?? sorumlu tutulamaz Mesaj?n yetkili al?c?s? veya al?c?s?na iletmekten sorumlu ki?i de?ilseniz, mesaj i?eri?ini ya da eklerini kullanmay?n?z, kopyalamay?n?z, yaymay?n?z, ba?ka ki?ilere y?nlendirmeyiniz ve mesaj? g?nderen ki?iyi derhal uyararak bu mesaj? siliniz. Kurumumuz size, mesaj?n ve bilgilerinin de?i?ikli?e u?ramamas?, b?t?nl???n?n ve gizlili?in korunmas? konusunda garanti vermemekte olup, e-posta i?eri?ine yetkisiz olarak yap?lan m?dahale, vir?s i?ermesi ve/veya bilgisayar sisteminize verebilece?i herhangi bir zarardan da sorumlu de?ildir. This e-mail message, including any attachments, is intended only for the use of the individual or entity to whom it is addressed and may contain confidential, private information as well as the exemption from disclosure. All ideas and opinions in this message and attached files can only belong to the addressee(s) and TUBITAK MRC shall have no liability to any person with regard to the use of the information contained in this message If you are not the intended addressee(s) or responsible person to inform the addressee(s), you are hereby notified that; any use, dissemination, distribution, or copying of this message and attached files is strictly prohibited. Delete the message and notify the sender immediately. TUBITAK MRC do not warrant for the accuracy, completeness of the contents of this email and/or the preservation of confidentiality and shall not be liable for the unauthorized changes made to this message, viruses and/or any damages caused in anyway to your computer system