Description: Calculate spectral intensity measures for synthetic seismograms provided in Surface-seismogram format. This code is adapted from "rspectra.f" provided by Doug Dreger (UC Berkeley). The CME version preserves the original spectral acceleration code and provides new capabilities of reading more than one seismogram plus pre-processing steps (conversion of seismogram units from velocity to acceleration, bandpass filtering).

Purpose:

Platform: This Fortran-90 code is coded to work on Solaris UNIX, Linux (HPCC Portland F90 and Intel "ifort").

Dependencies: Uses the CME Pathway 2 commandline/metadata parameter input library ("setparm").

Download: http://epicenter.usc.edu/cmeprograms/spec_accel.tar

--------------

Compilation: Use makefile within tar file.

Input information:

* Command line arguments for standalone execution.

* Ascii text file for Pathway 2 metadata.

#surfseis_rspectra in= out=

#
# in=                                       input Surf-seis file      [none]
# out=                                      output spectral file      [none]
# simulation_out_pointsX=                   # grid points in X        [none]
# simulation_out_pointsY=                   # grid points in Y        [none]
# simulation_out_timesamples=               # time points             [none]
# simulation_out_timeskip=                  sample rate (sec)         [none]
# surfseis_rspectra_output_type=     output: rd,rv,prv,aa,paa,magrat    [aa]
# surfseis_rspectra_output_format=          binary or text          [binary]
# surfseis_rspectra_galstog=                gals to g (1/980)     [.0010204]
# surfseis_rspectra_periodtype=             SCEC or TI                [SCEC]
# surfseis_rspectra_period=                 output period (sec)or "all" [2.]
# surfseis_rspectra_seismogram_units=                          [meterpersec]
#                                           vel:   "meterpersec"  "cmpersec"
#                                           accel: "meterpersec2" "cmpersec2"
# surfseis_rspectra_output_units=                                 [unitsofg]
#                                           rd:    "meter"        "cm"
#                                           rv,prv:"meterpersec"  "cmpersec"
#                                           aa,paa:"meterpersec2" "cmpersec2"
#                                           aa,paa:"unitsofg"     "percentg"
# surfseis_rspectra_apply_vel2accel=        apply: "yes" or "no"       [yes]
# surfseis_rspectra_apply_demean=           apply: "yes" or "no"       [yes]
# surfseis_rspectra_apply_filter=           apply: "yes" or "no"       [yes]
# surfseis_rspectra_apply_filter_lowHZ=     low HZ cutoff               [0.]
# surfseis_rspectra_apply_filter_highHZ=    high HZ cutoff         [Nyquist]
# surfseis_rspectra_numdampcoeffs=          #damping coeffs              [1]
# surfseis_rspectra_dampcoeffs1=            damping value #1          [0.05]
# surfseis_rspectra_dampcoeffs2=            damping value #2          [none]
#                                                          :
# surfseis_rspectra_dampcoeffs10=           damping value #10         [none]
#      More optional: output decimation in surface nodes:
# surfseis_rspectra_out_pointsXstart=       first seismogram in X        [1]
# surfseis_rspectra_out_pointsXend=         last  seismogram in X    [all X]
# surfseis_rspectra_out_pointsYstart=       first seismogram in Y        [1]
# surfseis_rspectra_out_pointsYend=         last  seismogram in Y    [all Y]
# surfseis_rspectra_out_pointsXdel=         increment in X               [1]
# surfseis_rspectra_out_pointsYdel=         increment in Y               [1]

Output information:

(A) binary or ascii format: surfseis_rspectra_output_format=      "binary" or "text"

(B) spectral intensity measure: surfseis_rspectra_output_type=     output: rd,rv,prv,aa,paa,magrat

(C) one or all periods: surfseis_rspectra_period=            output period (sec)or "all"

(D) output units: surfseis_rspectra_output_units=   rd: "meter"        "cm"
                                           rv,prv: "meterpersec"  "cmpersec"
                                           aa,paa: "meterpersec2" "cmpersec2"
                                           aa,paa: "unitsofg"     "percentg"