From BEElkins at mtech.edu Fri Apr 1 12:19:11 2011 From: BEElkins at mtech.edu (Elkins, Bailey E) Date: Fri, 1 Apr 2011 13:19:11 -0600 Subject: [SAC-HELP] Sac-Matlab Interface Message-ID: <088DBEC7FEB96246815271B837D93FF401EAA268B2@mtmail8c.butte.campus> I'm trying to utilize SAC's Matlab interface on an Intel Mac (64 bit) machine. I am running SAC v 101.4 and MATLAB_R2010bSP1. After I read in files I receive the following error while trying to execute the mat command: SAC> r /Users/BMcCommon/Thesis/Data/SAC/CCW/2009.142.00.00.13.0350.Y7.CCW..ELN.D.SAC SAC> mat Error: dlsym(0x100710ec0, engPutArray): symbol not found ERROR 8003: Cannot link to a MATLAB function: engPutArray ERROR 8003: Cannot link to a MATLAB function: Has anyone else encountered a similar problem? The Mathworks website seems to indicate that the "engPutArray" had been replaced with " engGetVariable" in the newer versions of Matlab. (http://www.mathworks.com/matlabcentral/newsreader/view_thread/292187) Is there anyway to fix this or do I need to downgrade my version of Matlab? I would greatly appreciate any help. Best, Bailey Elkins Geophysical Engineering Department Montana Tech Butte, MT From civyfei at gmail.com Wed Apr 6 13:13:37 2011 From: civyfei at gmail.com (Ching-Wen Chen) Date: Wed, 6 Apr 2011 15:13:37 -0500 Subject: [SAC-HELP] problems with evalresp - sac 101.2 OSX In-Reply-To: <20101223163241.18525ufvb60xhcow@webmail.rm.ingv.it> References: <20101223163241.18525ufvb60xhcow@webmail.rm.ingv.it> Message-ID: Hi , I wonder if an older version of sac is still available for request. I have a 10.5.8 Intel core 2 duo Macbook and it seems not working for current version (101.4). Thanks. Best regards, ChingWen Chen -------------- next part -------------- An HTML attachment was scrubbed... URL: From snoke at vt.edu Wed Apr 6 14:38:20 2011 From: snoke at vt.edu (Arthur Snoke) Date: Wed, 6 Apr 2011 17:38:20 -0400 (EDT) Subject: [SAC-HELP] problems with evalresp - sac 101.2 OSX In-Reply-To: References: <20101223163241.18525ufvb60xhcow@webmail.rm.ingv.it> Message-ID: There are two binary tarballs for SAC on the Mac: sac-101.4-mac_osx_64.tar.gz works only on OS 10.6 sac-101.4-mac_osx.tar.gz was built on a 10.4 Intel Mac and is a fat binary that should work on both PPC and Intel up to 10.5. Or, if you have the compilers, you could get the source distribution. For using SAC with Evalresp, look at the transfer help file: . As it says in the HISTORY file, all 101.4 versions of SAC work with version 3.3.3 of Evalresp. On Wed, 6 Apr 2011, Ching-Wen Chen wrote: > Hi , > > I wonder if an older version of sac is still available for request. I > have a 10.5.8 Intel core 2 duo Macbook and it seems not working for > current version (101.4). Thanks. > > Best regards, > > ChingWen Chen From shutian33 at yahoo.ca Sat Apr 9 05:21:03 2011 From: shutian33 at yahoo.ca (Shutian Ma) Date: Sat, 9 Apr 2011 05:21:03 -0700 (PDT) Subject: [SAC-HELP] data conversion from WGSN to SAC In-Reply-To: <088DBEC7FEB96246815271B837D93FF401EAA268B2@mtmail8c.butte.campus> Message-ID: <964793.36961.qm@web88301.mail.re4.yahoo.com> Dear colleagues, We obtained some waveform records in WGSN format (in asc), which is not a popular format. We do not have a program to process those records and we did not find a program to convert them into SAC. We assembled a program and successfully converted those records into SAC (binary). We would like to know the binary format and assemble a program to convert WGSN files in binary into SAC binary. If some of you know WGSN binary format, let me know. Thanks in advance, Shutian From warrenk at cwu.EDU Tue Apr 19 13:44:59 2011 From: warrenk at cwu.EDU (Kathryn Warren) Date: Tue, 19 Apr 2011 13:44:59 -0700 Subject: [SAC-HELP] Spectrogram command help Message-ID: <4DAD91DB02000081000714A4@hermes.cwu.edu> I am trying to create spectrograms using seismic data associated with snow avalanches. For my most recent data the Spectrogram command keeps giving me errors Instability detected in spectral estimator, increasing ridge regression factor to 1.00000e-04. Instability detected in spectral estimator, increasing ridge regression factor to 1.00000e-03. Instability detected in spectral estimator, increasing ridge regression factor to 1.00000e-02. Instability detected in spectral estimator, increasing ridge regression factor to 1.00000e-01. Error: Could not stabilize high res. spectrogram or segmentation fault or crashing the SAC program. Any help on this issue would be greatly appreciated. Thanks, Katie -------------- next part -------------- An HTML attachment was scrubbed... URL: From darko.dudjak at rgn.hr Wed Apr 20 05:21:41 2011 From: darko.dudjak at rgn.hr (darko.dudjak at rgn.hr) Date: Wed, 20 Apr 2011 14:21:41 +0200 Subject: [SAC-HELP] macro error in ch origin time Message-ID: <851446A2F16B413AB97A0087E154668C@ostr10418> Hi all, I made a macro in which I want to change a header value for origin time: READ $FILES READ $ORIGIN_TIME CH O GMT 2006 136 $ORIGIN_TIME WRITE $FILES However, if there is a 0 (zero) in a sec input, I get next error SAC> macro tscf FILES? 13_16052006.1.sac ORIGIN_TIME? 15 41 01 79 ERROR interpreting command: READ 15 41 01 79 ILLEGAL OPTION: ^ ERROR 1016: Terminating execution of macro If I leave a zero out, error still occurs. If values for secs are i.e. 40 or 50, macro goes through just fine. What should I do for sac to apply entered time? Thank you, Darko Dudjak -------------- next part -------------- An HTML attachment was scrubbed... URL: From warrenk at cwu.EDU Thu Apr 21 09:45:25 2011 From: warrenk at cwu.EDU (Kathryn Warren) Date: Thu, 21 Apr 2011 09:45:25 -0700 Subject: [SAC-HELP] spectrogram help: Arthur Snoke Message-ID: <4DAFFCB50200008100071745@hermes.cwu.edu> Arthur, It would be great if you could take a look at my data. It in on my ftp site ftp://ftp.geology.cwu.edu/warrenk/Outgoing/ Data file: march9 Commands run on the data are... rtrend hp co 10 cut 74020 74270 hp co 10 .......These all run fine spectrogram window 2 slice 1 method mem order 100 noscaling ymin 0 ymax 100 color .......This command fails Thanks, Katie -------------- next part -------------- An HTML attachment was scrubbed... URL: From kkrol001 at ucr.edu Sun Apr 24 12:10:37 2011 From: kkrol001 at ucr.edu (Kayla Kroll) Date: Sun, 24 Apr 2011 12:10:37 -0700 Subject: [SAC-HELP] getfhv - t0 marker Message-ID: Hello, I would like to get the 't0' header information from a large number of events. The fortran code (f77) below works perfectly when I replace the 't0' in the getfhv call by 'a', however when written as shown, I get the following error: * ERROR: Undefined header field value. t0* call rsac1(fn2,array,nlen,beg,del,max,nerr) if (nerr.gt.0) then print*, nerr go to 8888 endif call getfhv('t0',t0,nerr) if(t0.lt.0) then go to 9999 endif I have also tried replacing 't0' with 't1' as it is in the manual, but that produces the same error. The data does exist in the headers: [aten:EventWaveforms/0406_0614/10145282] kakroll% sac SEISMIC ANALYSIS CODE [02/05/2008 (Version 101.1)] Copyright 1995 Regents of the University of California SAC> r *ZY*EHN*sac 10145282.ZY.COON.EHN.sac 10145282.ZY.MLDN.EHN.sac 10145282.ZY.OCCM.EHN.sac 10145282.ZY.OYSB.EHN.sac 10145282.ZY.PNTO.EHN.sac 10145282.ZY.SHCM.EHN.sac 10145282.ZY.USG.EHN.sac 10145282.ZY.WAGN.EHN.sac SAC> lh t0 FILE: 10145282.ZY.COON.EHN.sac - 1 ------------------------------ t0 = 5.096052e+00 FILE: 10145282.ZY.MLDN.EHN.sac - 2 ------------------------------ t0 = 3.196018e+00 FILE: 10145282.ZY.OCCM.EHN.sac - 3 ------------------------------ t0 = 3.546471e+00 Any help would be appreciated. Thank you, Kayla Kroll UC Riverside -------------- next part -------------- An HTML attachment was scrubbed... URL: From snoke at vt.edu Sun Apr 24 12:29:34 2011 From: snoke at vt.edu (Arthur Snoke) Date: Sun, 24 Apr 2011 15:29:34 -0400 (EDT) Subject: [SAC-HELP] getfhv - t0 marker In-Reply-To: References: Message-ID: I have not looked aat the code, but I have looked at the discussion about header variables in the manual, and I think your problem is that the origin time is not an explicit header variable. In fact, it says that you can introduce ten different times t0-t9. One way that might work for you is to use ch to define t0 to be the origin time and then work with the t0 you have defined. On Sun, 24 Apr 2011, Kayla Kroll wrote: > Hello, > > ?? I would like to get the 't0' header information from a large number of > events. The fortran code (f77) below works perfectly when I replace the 't0' > in the getfhv call by 'a', however when written as shown, I get the > following error: > > ?ERROR: Undefined header field value. t0 > > ?? ? ?call rsac1(fn2,array,nlen,beg,del,max,nerr) > ?? ? ? if (nerr.gt.0) then > ?? ? ? print*, nerr > ?? ? ? go to 8888 > ?? ? ? endif > ?? ? ?call getfhv('t0',t0,nerr) > ?? ? ? if(t0.lt.0) then > ?? ? ? go to 9999 > ?? ? ? endif > > > I have also tried replacing 't0' with 't1' as it is in the manual, but that > produces the same error. > > The data does exist in the headers: > > [aten:EventWaveforms/0406_0614/10145282] kakroll% sac > ?SEISMIC ANALYSIS CODE [02/05/2008 (Version 101.1)] > ?Copyright 1995 Regents of the University of California > > SAC> r *ZY*EHN*sac > 10145282.ZY.COON.EHN.sac 10145282.ZY.MLDN.EHN.sac 10145282.ZY.OCCM.EHN.sac > 10145282.ZY.OYSB.EHN.sac 10145282.ZY.PNTO.EHN.sac 10145282.ZY.SHCM.EHN.sac > 10145282.ZY.USG.EHN.sac 10145282.ZY.WAGN.EHN.sac > SAC> lh t0 > > > ??FILE: 10145282.ZY.COON.EHN.sac - 1 > ?------------------------------ > > ?? ? t0 = 5.096052e+00 > > ??FILE: 10145282.ZY.MLDN.EHN.sac - 2 > ?------------------------------ > > ?? ? t0 = 3.196018e+00 > > ??FILE: 10145282.ZY.OCCM.EHN.sac - 3 > ?------------------------------ > > ?? ? t0 = 3.546471e+00 > > > Any help would be appreciated. > > Thank you, > Kayla Kroll > UC Riverside > >