From re.rezaei at gmail.com Fri Oct 2 01:46:21 2009 From: re.rezaei at gmail.com (reza rezaei) Date: Fri, 2 Oct 2009 12:16:21 +0330 Subject: [SAC-HELP] Depmax or Depmin period? Message-ID: <5feb258a0910020146k16ab7012u8c62c61a2f9bb162@mail.gmail.com> Hi every one Is it possible to read related period of Depmax or Depmin in Sac? Markptp can read period of peak to peak not peak or trough points. How can I read their periods? I would be pleased if you help me. Sincerely, Reza, From jattanayake at gmail.com Mon Oct 12 13:42:58 2009 From: jattanayake at gmail.com (Januka Attanayake) Date: Mon, 12 Oct 2009 16:42:58 -0400 Subject: [SAC-HELP] Synchronize issue Message-ID: <984f6f900910121342h4ee94523mb7cc704ec3fcfa0d@mail.gmail.com> Greetings All, I am trying to use the "synchorize" command in SAC 101.2 and have been getting the following error. Any advise is much appreciated. Thank you. Januka. SAC> r *.VEL.SAC AQU_BHZ.VEL.SAC BFO_BHZ.VEL.SAC BNI_BHZ.VEL.SAC DIVS_BHZ.VEL.SAC EIL_BHZ.VEL.SAC GRA1_BHZ.VEL.SAC GRFO_BHZ.VEL.SAC IDI_BHZ.VEL.SAC KHC_BHZ.VEL.SAC TIP_BHZ.VEL.SAC TRI_BHZ.VEL.SAC TUE_BHZ.VEL.SAC VTS_BHZ.VEL.SAC SAC> synch ERROR 0: Resulting error from shift exceeded tolerance (1e-6): 0.00040234 Please report bug to: sac-help at iris.washington.edu ERROR 0: Resulting error from shift exceeded tolerance (1e-6): 0.00036975 Please report bug to: sac-help at iris.washington.edu ERROR 0: Resulting error from shift exceeded tolerance (1e-6): 0.00024768 Please report bug to: sac-help at iris.washington.edu ERROR 0: Resulting error from shift exceeded tolerance (1e-6): 0.00048779 Please report bug to: sac-help at iris.washington.edu ERROR 0: Resulting error from shift exceeded tolerance (1e-6): 0.00014087 Please report bug to: sac-help at iris.washington.edu ERROR 0: Resulting error from shift exceeded tolerance (1e-6): 0.0005 Please report bug to: sac-help at iris.washington.edu ERROR 0: Resulting error from shift exceeded tolerance (1e-6): 0.00036975 Please report bug to: sac-help at iris.washington.edu ERROR 0: Resulting error from shift exceeded tolerance (1e-6): 0.0003728 Please report bug to: sac-help at iris.washington.edu ERROR 0: Resulting error from shift exceeded tolerance (1e-6): 0.00035754 Please report bug to: sac-help at iris.washington.edu ERROR 0: Resulting error from shift exceeded tolerance (1e-6): 0.00025586 Please report bug to: sac-help at iris.washington.edu SAC> p1 -------------- next part -------------- An HTML attachment was scrubbed... URL: From zm1984 at mail.igcas.ac.cn Mon Oct 12 13:53:36 2009 From: zm1984 at mail.igcas.ac.cn (ÕÔÃ÷) Date: 12 Oct 2009 20:53:36 -0000 Subject: [SAC-HELP] sac-help Digest, Vol 50, Issue 1 Message-ID: <20091012205336.17181.eqmail@mail.igcas.ac.cn> From: "??" Subject: ?????? = = = = = = = = = = = = = = = = = = = = > > ????????? > ?? > > > ?????????? > ????????zm1984 at mail.igcas.ac.cn ************************************************ Mr. Ming Zhao Institute of Geology and Geophysics Chinese Academy of Sciences Beijing 100029 China Tel: 0086-010-88258647(Office) ************************************************ From savage at uri.edu Tue Oct 13 05:37:31 2009 From: savage at uri.edu (Brian Savage) Date: Tue, 13 Oct 2009 08:37:31 -0400 Subject: [SAC-HELP] Synchronize issue In-Reply-To: <984f6f900910121342h4ee94523mb7cc704ec3fcfa0d@mail.gmail.com> References: <984f6f900910121342h4ee94523mb7cc704ec3fcfa0d@mail.gmail.com> Message-ID: Dear Januka, You can safely ignore these messages and they will be removed in upcoming version of SAC. Brian Savage savage at uri.edu On Oct 12, 2009, at 4:42 PM , Januka Attanayake wrote: > Greetings All, > I am trying to use the "synchorize" command in SAC 101.2 and have > been getting the following error. Any advise is much appreciated. > > Thank you. > > Januka. > > SAC> r *.VEL.SAC > AQU_BHZ.VEL.SAC BFO_BHZ.VEL.SAC BNI_BHZ.VEL.SAC DIVS_BHZ.VEL.SAC > EIL_BHZ.VEL.SAC GRA1_BHZ.VEL.SAC GRFO_BHZ.VEL.SAC IDI_BHZ.VEL.SAC > KHC_BHZ.VEL.SAC TIP_BHZ.VEL.SAC TRI_BHZ.VEL.SAC TUE_BHZ.VEL.SAC > VTS_BHZ.VEL.SAC > SAC> synch > ERROR 0: Resulting error from shift exceeded tolerance (1e-6): > 0.00040234 > Please report bug to: sac-help at iris.washington.edu > ERROR 0: Resulting error from shift exceeded tolerance (1e-6): > 0.00036975 > Please report bug to: sac-help at iris.washington.edu > ERROR 0: Resulting error from shift exceeded tolerance (1e-6): > 0.00024768 > Please report bug to: sac-help at iris.washington.edu > ERROR 0: Resulting error from shift exceeded tolerance (1e-6): > 0.00048779 > Please report bug to: sac-help at iris.washington.edu > ERROR 0: Resulting error from shift exceeded tolerance (1e-6): > 0.00014087 > Please report bug to: sac-help at iris.washington.edu > ERROR 0: Resulting error from shift exceeded tolerance (1e-6): > 0.0005 > Please report bug to: sac-help at iris.washington.edu > ERROR 0: Resulting error from shift exceeded tolerance (1e-6): > 0.00036975 > Please report bug to: sac-help at iris.washington.edu > ERROR 0: Resulting error from shift exceeded tolerance (1e-6): > 0.0003728 > Please report bug to: sac-help at iris.washington.edu > ERROR 0: Resulting error from shift exceeded tolerance (1e-6): > 0.00035754 > Please report bug to: sac-help at iris.washington.edu > ERROR 0: Resulting error from shift exceeded tolerance (1e-6): > 0.00025586 > Please report bug to: sac-help at iris.washington.edu > SAC> p1 > > _______________________________________________ > 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 =?ISO-8859-15?Q?=22langlaude_philippe_=2D_CETE_M=E9diterr=2E=2FLABO_NI?= Tue Oct 20 07:01:47 2009 From: =?ISO-8859-15?Q?=22langlaude_philippe_=2D_CETE_M=E9diterr=2E=2FLABO_NI?= (=?ISO-8859-15?Q?=22langlaude_philippe_=2D_CETE_M=E9diterr=2E=2FLABO_NI?=) Date: Tue, 20 Oct 2009 16:01:47 +0200 Subject: [SAC-HELP] Problems under cygwin Message-ID: <4ADDC2CB.4050103@developpement-durable.gouv.fr> Hello, I have just installed a version of SAC for Cygwin. It runs without any problems but when I use the plot function, I cannot close the graphic window which is created. I am obliged to kill SAC. On the other hand, I have some scripts which uses SAC and when I run this scripts I have this error : "readline: readline_callback_read_char() called with no handler!" Would anybody have a solution/explaination ? Thanks in advance. -- Philippe Langlaude Equipe de Recherche "Risque Sismique" CETE M?diterran?e Laboratoire de Nice 56 bd Stalingrad 06359 Nice cedex 4 France Tel: 06 70 83 57 85 / (33) 04 92 00 81 82 http://www.cete-mediterranee.fr philippe.langlaude at developpement-durable.gouv.fr From savage at uri.edu Tue Oct 20 08:51:40 2009 From: savage at uri.edu (Brian Savage) Date: Tue, 20 Oct 2009 11:51:40 -0400 Subject: [SAC-HELP] Problems under cygwin In-Reply-To: <4ADDC2CB.4050103@developpement-durable.gouv.fr> References: <4ADDC2CB.4050103@developpement-durable.gouv.fr> Message-ID: <70219D10-B4DE-449A-8D24-1656B7485AAB@uri.edu> Philippe A couple of things are happening here. First, the closing of the graphics window by using the "close window" button on the window is disabled as it was crashing sac. We need to find a solution to this issue. Second, sac normally uses the editline library (a compatible, but less restrictive library), but has the ability to use the readline library. We do not test normally test with the readline library. If you could send a simplified version of your script to me, I can test it with the readline library on windows to attempt to resolve this issue. If you expect your scripts to quit sac, you might want to make sure your scripts have the quit command at the end. This has been an issue previously. Brian savage at uri.edu On Oct 20, 2009, at 10:01 AM , langlaude philippe - CETE M?diterr./ LABO NICE/RS wrote: > Hello, > > I have just installed a version of SAC for Cygwin. It runs without > any problems but when I use the plot function, I cannot close the > graphic window which is created. I am obliged to kill SAC. On the > other hand, I have some scripts which uses SAC and when I run this > scripts I have this error : > > "readline: readline_callback_read_char() called with no handler!" > > Would anybody have a solution/explaination ? > > Thanks in advance. > -- > Philippe Langlaude > Equipe de Recherche "Risque Sismique" > CETE M?diterran?e Laboratoire de Nice 56 bd Stalingrad > 06359 Nice cedex 4 France > Tel: 06 70 83 57 85 / (33) 04 92 00 81 82 > http://www.cete-mediterranee.fr > philippe.langlaude at developpement-durable.gouv.fr > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help From pku132 at 163.com Sat Oct 24 05:21:09 2009 From: pku132 at 163.com (Lee) Date: Sat, 24 Oct 2009 20:21:9 +0800 Subject: [SAC-HELP] missing header and asynchronous 3C SAC Message-ID: <4AE2F137.0F2AF7.06288@m50-132.163.com> Hi, there, I received *.mseed.tar package as well as the related dataless file from IRIS DMC's ftp and extracted many mseed files from the package. Then I used rdseed_v4.8 to convert mseed to SAC with the command: rdseed -d -f PATH-TO-MSEEDFILE -g PATH-TO-DATALESSFILE That worked but there is no event information in the SAC header and many 3-component SAC files don't start from the same time. So I used the summary file generated by JWEED to make them synchronous as I did for SEED: rdseed -d -f PATH-TO-MSEEDFILE -g PATH-TO-DATALESSFILE -x PATH-TO-SUMMARYFILE Unfortunately, this doesn't work for mseed file. How can I deal with this problem? Good Luck, Lee -------------- next part -------------- An HTML attachment was scrubbed... URL: From dmikesell at cgiss.boisestate.edu Wed Oct 21 09:19:27 2009 From: dmikesell at cgiss.boisestate.edu (Dylan Mikesell) Date: Wed, 21 Oct 2009 10:19:27 -0600 Subject: [SAC-HELP] problem with readsac? In-Reply-To: References: <4ADF2713.7060505@cgiss.boisestate.edu> <4ADF2AE5.4060600@cgiss.boisestate.edu> Message-ID: <4ADF348F.5020900@cgiss.boisestate.edu> Hi SAC users, I am having some issues with timing...please see my steps below and I have attached the miniseed file in case anyone wants to try this themselves. I am losing millisecond information when I use the command ms2sac to convert data. I wonder if anyone has seen this before and how to fix it? Thanks, Dylan Mikesell 1) I read miniseed time information dmikesell at stonefly:$ mseedhdr 07.222.15.32.56.9541.1.m Net: FST: 2007:222:15:32:56.8800 LST: 2007:222:16:02:56.876000 n: 450000 time: 1800.000000 delta: 4.000000e-03 sps: 250.000000 This command reads the time stamp for every data point. You can see the first is 2007:222:15:32:56.8800 and the last is 2007:222:16:02:56.876000...everything looks correct...freq = 250 Hz, npts in 30 minutes, etc. 2) I convert to SAC dmikesell at stonefly:$ ms2sac 07.222.15.32.56.9541.1.m 07.222.15.32.56.9541.1.sac Warning: unable to open leap second file: /usr/local/lib/leapseconds Warning - unable to open coord file: /usr/contrib/data/bdsn/stat.db.coord The warnings don't matter...they are just options for adding coordinates to the headers and smoothing timing errors. I don't even care about time gaps right now...I just want the starting time to be correct. 3) read sac headers in two ways A) dmikesell at stonefly:$ sachdr 07.222.15.32.56.9541.1.sac nzyear = -687407104 nzjday = -570425344 nzhour = 251658240 nzmin = 536870912 nzsec = 956301312 nzmsec = 0 ...........this is now zero! B) read sac headers using SAC program dmikesell at stonefly:$ sac SEISMIC ANALYSIS CODE [04/07/2009 (Version 101.3)] Copyright 1995 Regents of the University of California SAC> r 07.222.15.32.56.9541.1.sac SAC> lh FILE: 07.222.15.32.56.9541.1.sac - 1 -------------------------------- NPTS = 450000 B = -1.200000e-01 E = 1.799876e+03 IFTYPE = TIME SERIES FILE LEVEN = TRUE DELTA = 4.000000e-03 IDEP = UNKNOWN DEPMIN = -3.049800e+04 DEPMAX = 2.401200e+04 DEPMEN = -1.041911e+03 KZDATE = AUG 10 (222), 2007 KZTIME = 15:32:57.000 IZTYPE = EVENT ORIGIN TIME KSTNM = 09541 CMPINC = -1.225500e+04 KEVNM = 072221531C1 IEVTYP = EARTHQUAKE It is obvious the KZTIME is not the same as the FST time in the miniseed header. However, the KZTIME was rounded up correctly. I need the milliseconds to be accurate though because my stations are very close. Does anyone know how to fix this issue? Many thanks! Dylan Mikesell -------------- next part -------------- A non-text attachment was scrubbed... Name: 07.222.15.32.56.9541.1.m Type: text/x-objcsrc Size: 517120 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: dmikesell.vcf Type: text/x-vcard Size: 360 bytes Desc: not available URL: