Hello,
I'm running a matlab amplitude check on GSN stations for the recent Santa
Cruz Islands event. I've done this for past event of M>=7 with no problem.
Today my routine is failing because the structure returned by irisFetch.Traces
for IC.ENH.00.BHZ (for example) is a structure of 1x2 arrays having duplicate
values rather than a structure of single values. I'm not getting this behavior with
CI.PASC.00.BHZ.
Here's what I'm used to:
network: 'CI'
station: 'PASC'
location: '00'
channel: 'BHZ'
quality: 'M'
latitude: 34.1714
longitude: -118.1852
elevation: 341
depth: 0
azimuth: 0
dip: -90
sensitivity: 5.9199e+09
sensitivityFrequency: 0.0300
instrument: 'STS-1'
sensitivityUnits: 'M/S'
data: [288000x1 single]
sampleCount: 288000
sampleRate: 40
startTime: 7.3527e+05
endTime: 7.3527e+05
sacpz: [1x1 struct]
Here's how IC.ENH differs:
1x2 struct array with fields:
network
station
location
channel
quality
latitude
longitude
elevation
depth
azimuth
dip
sensitivity
sensitivityFrequency
instrument
sensitivityUnits
data
sampleCount
sampleRate
startTime
endTime
sacpz
ans = 20
ans =IC
etc.
The only field where the values aren't duplicated is sampleCount:
ans = 60270
Is this what you would expect? If so, can you explain when/why it occurs? I'm wondering
if this is an error or something I need to catch and handle.
Thanks a bunch,
Mary Templeton
I'm running a matlab amplitude check on GSN stations for the recent Santa
Cruz Islands event. I've done this for past event of M>=7 with no problem.
Today my routine is failing because the structure returned by irisFetch.Traces
for IC.ENH.00.BHZ (for example) is a structure of 1x2 arrays having duplicate
values rather than a structure of single values. I'm not getting this behavior with
CI.PASC.00.BHZ.
Here's what I'm used to:
net =CInet='CI'
sta =PASCsta='PASC'
loc =00loc='00'
chan =BHZchan='BHZ'
starttime =2013-02-06 01:12:35starttime
endtime =2013-02-06 03:12:35endtime
primary =primary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
network: 'CI'
station: 'PASC'
location: '00'
channel: 'BHZ'
quality: 'M'
latitude: 34.1714
longitude: -118.1852
elevation: 341
depth: 0
azimuth: 0
dip: -90
sensitivity: 5.9199e+09
sensitivityFrequency: 0.0300
instrument: 'STS-1'
sensitivityUnits: 'M/S'
data: [288000x1 single]
sampleCount: 288000
sampleRate: 40
startTime: 7.3527e+05
endTime: 7.3527e+05
sacpz: [1x1 struct]
Here's how IC.ENH differs:
net =ICnet='IC'
sta =ENHsta='ENH'
secondary =secondary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
1x2 struct array with fields:
network
station
location
channel
quality
latitude
longitude
elevation
depth
azimuth
dip
sensitivity
sensitivityFrequency
instrument
sensitivityUnits
data
sampleCount
sampleRate
startTime
endTime
sacpz
ans = 20secondary.sampleRate
ans = 20
ans =ICsecondary.network
ans =IC
etc.
The only field where the values aren't duplicated is sampleCount:
ans = 75712secondary.sampleCount
ans = 60270
Is this what you would expect? If so, can you explain when/why it occurs? I'm wondering
if this is an error or something I need to catch and handle.
Thanks a bunch,
Mary Templeton
-
Hi Mary,
In the examples you have given, in the case of CI_PASC you are specifying a channel and location code. With IC_ENH you are not specifying either, you wouldn't you expect more than one record (or in Matlab language, structure of 1x2 arrays)? Maybe try specifying a channel and/or loc code and see what you get back.
Best regards,
- Rob N
On Feb 7, 2013, at 3:01 PM, Mary Templeton wrote:
Hello,
_____________________________________________
I'm running a matlab amplitude check on GSN stations for the recent Santa
Cruz Islands event. I've done this for past event of M>=7 with no problem.
Today my routine is failing because the structure returned by irisFetch.Traces
for IC.ENH.00.BHZ (for example) is a structure of 1x2 arrays having duplicate
values rather than a structure of single values. I'm not getting this behavior with
CI.PASC.00.BHZ.
Here's what I'm used to:
net='CI'
sta='PASC'
loc='00'
chan='BHZ'
starttime
endtime
primary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
network: 'CI'
station: 'PASC'
location: '00'
channel: 'BHZ'
quality: 'M'
latitude: 34.1714
longitude: -118.1852
elevation: 341
depth: 0
azimuth: 0
dip: -90
sensitivity: 5.9199e+09
sensitivityFrequency: 0.0300
instrument: 'STS-1'
sensitivityUnits: 'M/S'
data: [288000x1 single]
sampleCount: 288000
sampleRate: 40
startTime: 7.3527e+05
endTime: 7.3527e+05
sacpz: [1x1 struct]
Here's how IC.ENH differs:
net='IC'
sta='ENH'
secondary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
1x2 struct array with fields:
network
station
location
channel
quality
latitude
longitude
elevation
depth
azimuth
dip
sensitivity
sensitivityFrequency
instrument
sensitivityUnits
data
sampleCount
sampleRate
startTime
endTime
sacpz
secondary.sampleRate
ans = 20
secondary.network
ans =IC
etc.
The only field where the values aren't duplicated is sampleCount:
secondary.sampleCount
ans = 60270
Is this what you would expect? If so, can you explain when/why it occurs? I'm wondering
if this is an error or something I need to catch and handle.
Thanks a bunch,
Mary Templeton
_______________________________________________
webservices mailing list
webservices<at>iris.washington.edu
http://www.iris.washington.edu/mailman/listinfo/webservices
Rob Newman
Technology Architect, IRIS DMC
[e] rnewman<at>iris.washington.edu
-
Mary Templeton2013-02-08 00:28:27Hi Rob,
My apologies - it's still chan='BHZ' and loc='00'. I didn't reset them,
so I didn't relist them. They're appropriate for IC.ENH.
At any rate, the two chan and loc values returned are duplicates -
{00, 00} and {BHZ,BHZ}
--mt
On Feb 7, 2013, at 4:14 PM, Rob Newman wrote:
Hi Mary,
In the examples you have given, in the case of CI_PASC you are specifying a channel and location code. With IC_ENH you are not specifying either, you wouldn't you expect more than one record (or in Matlab language, structure of 1x2 arrays)? Maybe try specifying a channel and/or loc code and see what you get back.
Best regards,
- Rob N
On Feb 7, 2013, at 3:01 PM, Mary Templeton wrote:
Hello,
_____________________________________________
I'm running a matlab amplitude check on GSN stations for the recent Santa
Cruz Islands event. I've done this for past event of M>=7 with no problem.
Today my routine is failing because the structure returned by irisFetch.Traces
for IC.ENH.00.BHZ (for example) is a structure of 1x2 arrays having duplicate
values rather than a structure of single values. I'm not getting this behavior with
CI.PASC.00.BHZ.
Here's what I'm used to:
net='CI'
sta='PASC'
loc='00'
chan='BHZ'
starttime
endtime
primary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
network: 'CI'
station: 'PASC'
location: '00'
channel: 'BHZ'
quality: 'M'
latitude: 34.1714
longitude: -118.1852
elevation: 341
depth: 0
azimuth: 0
dip: -90
sensitivity: 5.9199e+09
sensitivityFrequency: 0.0300
instrument: 'STS-1'
sensitivityUnits: 'M/S'
data: [288000x1 single]
sampleCount: 288000
sampleRate: 40
startTime: 7.3527e+05
endTime: 7.3527e+05
sacpz: [1x1 struct]
Here's how IC.ENH differs:
net='IC'
sta='ENH'
secondary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
1x2 struct array with fields:
network
station
location
channel
quality
latitude
longitude
elevation
depth
azimuth
dip
sensitivity
sensitivityFrequency
instrument
sensitivityUnits
data
sampleCount
sampleRate
startTime
endTime
sacpz
secondary.sampleRate
ans = 20
secondary.network
ans =IC
etc.
The only field where the values aren't duplicated is sampleCount:
secondary.sampleCount
ans = 60270
Is this what you would expect? If so, can you explain when/why it occurs? I'm wondering
if this is an error or something I need to catch and handle.
Thanks a bunch,
Mary Templeton
_______________________________________________
webservices mailing list
webservices<at>iris.washington.edu
http://www.iris.washington.edu/mailman/listinfo/webservices
Rob Newman
Technology Architect, IRIS DMC
[e] rnewman<at>iris.washington.edu
_______________________________________________
webservices mailing list
webservices<at>iris.washington.edu
http://www.iris.washington.edu/mailman/listinfo/webservices
-
-
Celso Reyes2013-02-07 23:24:59Hi Mary,
You've found a gap in the data that is to big to be automatically bridged.
Here...
datestr([primary.startTime])
06-Feb-2013 01:12:35
06-Feb-2013 02:22:21
datestr([primary.endTime])
06-Feb-2013 02:15:40
06-Feb-2013 03:12:34
So, it looks like there's a huge gap between 02:15:40 and 02:22:21
Because there was no data there, irisFetch doesn't know how to connect the two, so it gives you two traces instead.
This is likely a condition you'll have to keep an eye out for in your processing codes; as I bet it has other implications for processing.
Cheers,
Celso
On Feb 7, 2013, at 3:01 PM, Mary Templeton wrote:
net='CI'
sta='PASC'
loc='00'
chan='BHZ'
starttime
endtime
primary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
-
Mary Templeton2013-02-08 00:33:15Thanks Celso,
That makes sense. I'll see what I can do to catch it.
cheers,
Mary
On Feb 7, 2013, at 4:24 PM, Celso Reyes wrote:
Hi Mary,
You've found a gap in the data that is to big to be automatically bridged.
Here...
datestr([primary.startTime])
06-Feb-2013 01:12:35
06-Feb-2013 02:22:21
datestr([primary.endTime])
06-Feb-2013 02:15:40
06-Feb-2013 03:12:34
So, it looks like there's a huge gap between 02:15:40 and 02:22:21
Because there was no data there, irisFetch doesn't know how to connect the two, so it gives you two traces instead.
This is likely a condition you'll have to keep an eye out for in your processing codes; as I bet it has other implications for processing.
Cheers,
Celso
On Feb 7, 2013, at 3:01 PM, Mary Templeton wrote:
net='CI'
sta='PASC'
loc='00'
chan='BHZ'
starttime
endtime
primary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
webservices mailing list
webservices<at>iris.washington.edu
http://www.iris.washington.edu/mailman/listinfo/webservices
-