170 likes | 218 Views
Learn how to download, manage, and process seismic data using Jrdseed, PQL, SAC, and other tools. Explore examples and tutorials on how to work with SEED files, convert formats, and extract responses in this informative workshop.
E N D
Using SEED Using: Jrdseed PQL SAC RESP JPlotResp Data Management Workshop Kuala Lumpur, Malaysia
Data Examples download • Create a download folder DATA_EXAMPLES • http://www.iris.edu/pub/userdata/Rick_Benson/WKSHOP • Get dataless SEED file MY.dataless • Get miniseed data file like KKM.MY..LHE.2007.276.R • Get file calledfullSEED Data Management Workshop Petaling Jaya, Malaysia
Portable SEED File Reading Jrdseed - Java port of the IRIS-DMC’s ‘rdseed’ • Reads full SEED, dataless SEED and Mini-SEED+dataless • Writes SEED, SAC, AH, CSS, SEGY, SEED RESP • Time-window and other output selection possible • Runs anywhere a JVM is available • Based on generic JavaSeed classes available separately • Not as complete as ‘rdseed’ yet Data Management Workshop Petaling Jaya, Malaysia
Jrdseed To install from Internet: • Connect with browser to: • http://www.iris.edu/manuals/rdseed_request.htm • Get Jrdseed or Rdseed depending on platform • Put in Applications folder. Data Management Workshop Kuala Lumpur, Malaysia
Run Jrdseed • Two modes: command line or interactive Command Line, for scripting, etc Example: • Jrdseed -d -o 4 -f “seed_volume” • Will output a file called “mini.seed” in current working directory. (because we used -o 4) • NOTE: No alternate output filename supported, so if you run it iteratively, make sure you move the mini.seed file to a new name between iterations. Data Management Workshop Petaling Jaya, Malaysia
Interactive Jrdseed Options for output: << IRIS SEED Reader, Release 4.7.5 >> Input File (/dev/nrst0) or 'Quit' to Exit: seed Output File (stdout) : Volume # [(1)-N] : Options [acCsSpRtde] : a: Abbreviation dictionary c: contents, timeseries C: Station/channel comments s: summary of entire volume S: station list into a file “rdseed.stations” p: Print SAC Pole/Zero files SAC_PZs_TA_U05C_BHN__2004.341.00.00.00.0000_99999.9999.24.60.60.99999 R: Print RESP format files t: Listing “timeseries” of waveform data in volume d: Output Format [(1=SAC), 2=AH, 3=CSS, 4=mini seed, 5=seed, 6=sac ascii, 7=SEGY] Data Management Workshop Petaling Jaya, Malaysia
Jrdseed: conjoin w/ miniSEED and Dataless • By having separate miniseed and dataless, you can still run rdseed: • csh: setenv ALT_RESPONSE_FILE `pwd`/dataless • bash: export ALT_RESPONSE_FILE=`pwd`/dataless • Win: use the “-g” option • http://www.iris.edu/manuals/jrdseed.htm Data Management Workshop Petaling Jaya, Malaysia
Jrdseed Exercise #1 • Run jrdseed and create RESP files with the “R” option. << Java SEED Reader, Release 0.07 >> Input File (/dev/nrst0) or 'Quit' to Exit: MY.dataless Output File (stdout) : Volume # [(1)-N] : Options [acCsSpRtde] : R Station List (ALL) : Channel List (ALL) : Network List (ALL) : Loc Ids (ALL ["--" for spaces]) : Input File (/dev/nrst0) or 'Quit' to Exit: quit Jrdseed completed. • We will be using these when running JPlotResp in just a bit Data Management Workshop Petaling Jaya, Malaysia
Jrdseed exercise #2 • Run jrdseed and convert miniseed to SAC << Java SEED Reader, Release 0.07 >> Input File (/dev/nrst0) or 'Quit' to Exit: fullSEED Output File (stdout) : Volume # [(1)-N] : Options [acCsSpRtde] : d Station List (ALL) : Channel List (ALL) : Network List (ALL) : Loc Ids (ALL ["--" for spaces]) : Output Format [(1=SAC), 2=AH, 3=CSS, 4=mini seed, 5=seed, 6=sac ascii, 7=SEGY] : Select Data Type [(E=Everything), D=Data of Undetermined State, R=Raw waveform Data, Q=QC'd data] :E Start Time(s) YYYY,DDD,HH:MM:SS.FFFF : End Time(s) YYYY,DDD,HH:MM:SS.FFFF : Extract Responses [Y/(N)] : Writing IU.KIEV.00.BHE, 9267 samples (binary), starting 2005,196 08:34:39.3850 UT Writing IU.KIEV.00.BHN, 6372 samples (binary), starting 2005,196 08:37:58.1350 UT Writing IU.KIEV.00.BHZ, 7743 samples (binary), starting 2005,196 08:35:53.8350 UT Input File (/dev/nrst0) or 'Quit' to Exit: quit Jrdseed completed. • We will be using these to plot using SAC in a bit. Data Management Workshop Petaling Jaya, Malaysia
SAC IRIS distributes and maintains support for SAC “Seismic Analysis Code” software developed by LLNL Stay informed, ask questions, and get help using user-group list: http://norman.iris.washington.edu/mailman/listinfo sac-help and sac-dev To access and install: http://www.iris.edu/manuals/sac Data Management Workshop Petaling Jaya, Malaysia
PQL To acquire and install: http://www.passcal.nmt.edu/announcements/pqlII.htm Data Management Workshop Petaling Jaya, Malaysia
PQLX- Next generation Data Management Workshop Petaling Jaya, Malaysia
RESP Format: Response files • # • # ======== CHANNEL RESPONSE DATA ======== • B050F03 Station: BTDF • B050F16 Network: MS • B052F03 Location: ?? • B052F04 Channel: BHE • B052F22 Start date: 1996,214,00:00:00 • B052F23 End date: No Ending Time • # ======================================= • # + +--------------------------------------------+ + • # + | Response (Poles & Zeros), BTDF 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 Data Management Workshop Petaling Jaya, Malaysia
JPlotResp: Evaluate Response Locate and install from http://www.isti.com/JPlotResp/v1.57/install.htm Data Management Workshop Petaling Jaya, Malaysia
Plot Response curve Data Management Workshop Petaling Jaya, Malaysia
JPlotResp: get DHI RESP files Data Management Workshop Petaling Jaya, Malaysia