Thread: printed comment in transfer command

Started: 2011-07-25 21:58:44
Last activity: 2011-07-26 23:57:45
Topics: SAC Help
Neil Selby
2011-07-25 21:58:44
Hello, hope this gets through!


this is an example of the behaviour I described to Sheila:


[neil@braeriach scratch]$ /sharedprograms/sac/v101.4/64/sac/bin/sac
SEISMIC ANALYSIS CODE [06/07/2010 (Version 101.4)]
Copyright 1995 Regents of the University of California

SAC> r *SHZ*sac
20020511_000000.EKB10.SHZ.sac 20020511_000000.EKB1.SHZ.sac
20020511_000000.EKB2.SHZ.sac 20020511_000000.EKB3.SHZ.sac
20020511_000000.EKB4.SHZ.sac 20020511_000000.EKB5.SHZ.sac
20020511_000000.EKB6.SHZ.sac 20020511_000000.EKB7.SHZ.sac
20020511_000000.EKB8.SHZ.sac 20020511_000000.EKB9.SHZ.sac
20020511_000000.EKR10.SHZ.sac 20020511_000000.EKR1.SHZ.sac
20020511_000000.EKR2.SHZ.sac 20020511_000000.EKR3.SHZ.sac
20020511_000000.EKR4.SHZ.sac 20020511_000000.EKR5.SHZ.sac
20020511_000000.EKR6.SHZ.sac 20020511_000000.EKR7.SHZ.sac
20020511_000000.EKR8.SHZ.sac 20020511_000000.EKR9.SHZ.sac
SAC> lh scale


FILE: 20020511_000000.EKB10.SHZ.sac - 1
-----------------------------------

scale = 7.225000e-02

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

scale = 7.225000e-02

FILE: 20020511_000000.EKB2.SHZ.sac - 3
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKB3.SHZ.sac - 4
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKB4.SHZ.sac - 5
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKB5.SHZ.sac - 6
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKB6.SHZ.sac - 7
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKB7.SHZ.sac - 8
----------------------------------

scale = 7.225000e-02
Waiting


FILE: 20020511_000000.EKB8.SHZ.sac - 9
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKB9.SHZ.sac - 10
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKR10.SHZ.sac - 11
-----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKR1.SHZ.sac - 12
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKR2.SHZ.sac - 13
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKR3.SHZ.sac - 14
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKR4.SHZ.sac - 15
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKR5.SHZ.sac - 16
----------------------------------

scale = 7.225000e-02
Waiting


FILE: 20020511_000000.EKR6.SHZ.sac - 17
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKR7.SHZ.sac - 18
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKR8.SHZ.sac - 19
----------------------------------

scale = 7.225000e-02

FILE: 20020511_000000.EKR9.SHZ.sac - 20
----------------------------------

scale = 7.225000e-02
SAC>
SAC> rmean
SAC> taper
SAC> trans from pol s EKA-SZ.pz to none freq 1 2 3 4
Extracting polezero response for EKB10 , SHZ ...

Station (EKB10 ), 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.
Extracting polezero response for EKB2 , SHZ ...

Station (EKB2 ), 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.
Extracting polezero response for EKB3 , SHZ ...

Station (EKB3 ), 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.
Extracting polezero response for EKB4 , SHZ ...

Station (EKB4 ), 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.
Extracting polezero response for EKB5 , SHZ ...

Station (EKB5 ), 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.
Extracting polezero response for EKB6 , SHZ ...

Station (EKB6 ), 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.
Extracting polezero response for EKB7 , SHZ ...

Station (EKB7 ), 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.
Extracting polezero response for EKB8 , SHZ ...

Station (EKB8 ), 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.
Extracting polezero response for EKB9 , SHZ ...

Station (EKB9 ), 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.
Extracting polezero response for EKR10 , SHZ ...

Station (EKR10 ), 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.
Extracting polezero response for EKR1 , SHZ ...

Station (EKR1 ), 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.
Extracting polezero response for EKR2 , SHZ ...

Station (EKR2 ), 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.
Extracting polezero response for EKR3 , SHZ ...

Station (EKR3 ), 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.
Extracting polezero response for EKR4 , SHZ ...

Station (EKR4 ), 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.
Extracting polezero response for EKR5 , SHZ ...

Station (EKR5 ), 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.
Extracting polezero response for EKR6 , SHZ ...

Station (EKR6 ), 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.
Extracting polezero response for EKR7 , SHZ ...

Station (EKR7 ), 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.
Extracting polezero response for EKR8 , SHZ ...

Station (EKR8 ), 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.
Extracting polezero response for EKR9 , SHZ ...

Station (EKR9 ), 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>



plotting the data shows that all but the first channel has been
multiplied by the "scale" factor 0.072250. If you don't want to multiply
by scale (which you often don't) then care needs to be taken. I assume
that the above behaviour is not intended.

Also I don't think it's fair to say that the "scale" value is unused -
it's unused when rdseed generated sac files are being used, certainly.

Neil Selby
AWE Blacknest

  • Arthur Snoke
    2011-07-26 23:57:45
    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.



02:59:58 v.01697673