Thread: irisFetch.Traces returning duplicate values for unique SNCL-time

Started: 2013-02-08 00:01:55
Last activity: 2013-02-08 00:33:15
Topics: Web Services
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'
net =CI
sta='PASC'
sta =PASC
loc='00'
loc =00
chan='BHZ'
chan =BHZ
starttime
starttime =2013-02-06 01:12:35
endtime
endtime =2013-02-06 03:12:35

primary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
primary =
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'
net =IC
sta='ENH'
sta =ENH

secondary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
secondary =
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
ans = 20
secondary.network
ans =IC
ans =IC

etc.

The only field where the values aren't duplicated is sampleCount:
secondary.sampleCount
ans = 75712
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'
    net =CI
    sta='PASC'
    sta =PASC
    loc='00'
    loc =00
    chan='BHZ'
    chan =BHZ
    starttime
    starttime =2013-02-06 01:12:35
    endtime
    endtime =2013-02-06 03:12:35

    primary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
    primary =
    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'
    net =IC
    sta='ENH'
    sta =ENH

    secondary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
    secondary =
    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
    ans = 20
    secondary.network
    ans =IC
    ans =IC

    etc.

    The only field where the values aren't duplicated is sampleCount:
    secondary.sampleCount
    ans = 75712
    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




    • Hi 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'
      net =CI
      sta='PASC'
      sta =PASC
      loc='00'
      loc =00
      chan='BHZ'
      chan =BHZ
      starttime
      starttime =2013-02-06 01:12:35
      endtime
      endtime =2013-02-06 03:12:35

      primary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
      primary =
      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'
      net =IC
      sta='ENH'
      sta =ENH

      secondary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)
      secondary =
      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
      ans = 20
      secondary.network
      ans =IC
      ans =IC

      etc.

      The only field where the values aren't duplicated is sampleCount:
      secondary.sampleCount
      ans = 75712
      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



  • Hi Mary,
    You've found a gap in the data that is to big to be automatically bridged.

    Here...
    datestr([primary.startTime])

    ans =

    06-Feb-2013 01:12:35
    06-Feb-2013 02:22:21

    datestr([primary.endTime])

    ans =

    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'
    net =CI
    sta='PASC'
    sta =PASC
    loc='00'
    loc =00
    chan='BHZ'
    chan =BHZ
    starttime
    starttime =2013-02-06 01:12:35
    endtime
    endtime =2013-02-06 03:12:35

    primary = irisFetch.Traces(net,sta,loc,chan,starttime,endtime)


    • Thanks 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])

      ans =

      06-Feb-2013 01:12:35
      06-Feb-2013 02:22:21

      datestr([primary.endTime])

      ans =

      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'
      net =CI
      sta='PASC'
      sta =PASC
      loc='00'
      loc =00
      chan='BHZ'
      chan =BHZ
      starttime
      starttime =2013-02-06 01:12:35
      endtime
      endtime =2013-02-06 03:12:35

      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



13:10:16 v.22510d55