[SAC-HELP] printed comment in transfer command

Arthur Snoke snoke at vt.edu
Tue Jul 26 13:57:45 PDT 2011


Neil (and Sheila),

Brian and I have run several experiments and find that indeed, there is an 
error in runs with v101.4 and earlier when a single transfer call is 
applied to more than one file in memory if the header variable SCALE is 
defined, the NWFID is not defined, and the file is the first in the list 
of files.  The bug can occur when only one file has been read.  Reading in 
one of your files twice and applying your transfer call, one gets 
different depmin and depmax for the two (identical) files:

  SEISMIC ANALYSIS CODE [06/07/2010 (Version 101.4)]
  Copyright 1995 Regents of the University of California

SAC> r ./20020511_000000.EKB1.SHZ.sac ./20020511_000000.EKB1.SHZ.sac
SAC> rmean
SAC> taper
SAC> transfer from polezero s ./EKA-SZ.pz to none freq 1 2 3 4
  Extracting polezero response for EKB1    , SHZ     ...

Station (EKB1    ), Channel (SHZ     )
         Waveform multiplied by 1.000000 after deconvolution.
  Extracting polezero response for EKB1    , SHZ     ...

Station (EKB1    ), Channel (SHZ     )
         Nominal Calper not set. Setting Nominal Calper to -999.000000 to 
match Wfdisc Calper
         Sensor Calper not set. Setting Sensor Calper to 1.000000 to match 
Wfdisc Calper
         Response appears to be Normalized.
         Waveform appears to be unscaled.
         Wfdisc calib value is (0.072250).
         Waveform multiplied by 0.072250 after deconvolution.
SAC> lh scale depmin depmax


   FILE: ./20020511_000000.EKB1.SHZ.sac - 1
  ------------------------------------

      scale = 1.000000e+00                   depmin = -1.396522e+01
     depmax = 1.375477e+01

   FILE: ./20020511_000000.EKB1.SHZ.sac - 2
  ------------------------------------

      scale = 1.000000e+00                   depmin = -1.008987e+00
     depmax = 9.937820e-01
SAC>

The SCALE was indeed used for processing CSS files for which CALIB was 
used and for which NWFID (Css 3.0 wavform id) is defined.  For v101.5 
SCALE will truly not be used (at least in TRANSFER).

For now, if one undefines SCALE, one should get consistent results:

SAC> r ./20020511_000000.EKB1.SHZ.sac ./20020511_000000.EKB1.SHZ.sac
SAC> ch SCALE undef

    dot dot dot

SAC> lh scale depmin depmax


   FILE: ./20020511_000000.EKB1.SHZ.sac - 1
  ------------------------------------

      scale = 1.000000e+00                   depmin = -1.396522e+01
     depmax = 1.375477e+01

   FILE: ./20020511_000000.EKB1.SHZ.sac - 2
  ------------------------------------

      scale = 1.000000e+00                   depmin = -1.396522e+01
     depmax = 1.375477e+01
SAC>

We will put some warnings in the help file for READCSS to alert users that 
use CALIB that it is not recognized in SAC v101.5.

Arthur

p.s. Thanks, Sheila for helping us find an error we hopefully have solved.





More information about the sac-help mailing list