Introduction

The purpose of this page is to notify users of the Computer Program in Seismology (CPS) package of repairs or extensions made to the code. Whenever this page is updated, the down-loadable version of CPS is updated and provided to users.

Enhancements to Version 3.30 codes



2005, 2006, 2007, 2008, 2009, 2010, 2011, 2012, 2013, 2014 2015 , 2016, 2017, 2018, 2019, 2020, 2021, 2022 For each year, the changes are (listed in order from newest to oldest)




Noevmber 7, 2022

October 12, 2022

July 29,  2022

June 9, 2022

June 6, 2022

March 20, 2022

From
float tr, ti;
...
temp = sqrt(tr*tr + ti*ti);

To
double tr, ti;
...
temp = (float)sqrt(tr*tr + ti*ti);

The problem was that some spectra that should have been non-zero, e.g., tr and ti = 1.0e-20, were set to zero because of the multiplication inside the square root, even thought the square root is a double.

March 9, 2022

March 6, 2022



February 7, 2022

December 11, 2021

September 20, 2021

September 19, 2021

WVFMT96    1.0  288.   74.  -91.   4.15     1.000 0.186E-12     1.000     1.000 0.154E-12  54.5 
where the columns are
 1 WVFMT96 if created by wvf96.f
 WVFMTD96 if created by wvfmtd96
 The next four columns describe the best (major) double couple derived from the moment tensor
 2  depth - depth of the Greens functions
 3 stk - strike of best double couple
 4 dip - dip of best double couple
 5  rake - rake of best double couple
 6  Mw  - moment magnitude of best double couple
 7 Reduction of Variance: 1.0-total_err/total_sum_sq - raw data 0 zero is best
 8 Std Error of fit: sqrt(total_err*tmpsig) 
 9 xy/sqrt(xx*yy) - cosine of angle between observed and predicted
10  Weighted reduction of variance: 1.0-total_wtd/total_sum_sq_wt 
11  Weighted Std error of fit: sqrt(total_wtd*tmpsig) 
12  Rercent CLVD - 0 means pure double couple



The new output is
WVFMT961    1.0  288.   74.  -91.   4.15     1.000 0.186E-12     1.000     1.000 0.154E-12  54.5 -0.1010000E+23 -0.1060000E+23  0.3000003E+21  0.8799998E+22  0.2799999E+22 -0.2190001E+23
The program name in the first column changes and there are 6 additional columns


        1 WVFMT961  if
        created by the updated wvf96.f

          WVFMTD961  if
        created by the updated wvfmtd96

        Columns 2 - 12 are the 
        13  Mxx dyne-cm)
        14  Myy (dyne-cm)
        15  Mxy (dyne-cm)
        16  Mxz (dyne-cm)
        17  Myz (dyne-cm)
        18  Mzz (dyne-cm)

      



September 15, 2021

September 4, 2021

The reason for this option is that users may wish to use another program to plot the model parameters as a function of depth.

Examples of the plots created are given in SHWMOD96/index.html

August 31, 2021

June 2, 2021

February 1, 2021

January 31, 2021



October 29, 2020

October 1, 2020

'eml_H.xy' 0 0.01 'CI' 0.05  '1 sig'
where eml_H.xy has entries   
   2.5000000       3.3192308       0.0000000      0.40308103   
   2.7000000       3.3750000       0.0000000      0.29140610   
   2.9000001       3.4560604       0.0000000      0.28769836   

The result is shown in the next figure:

genplt

This figure was created by concatenating the CALPLOT files of three invocations of genplt:
#####
#   plot the local earthquake Mw vs ML
#####cat > lcmdfil << EOF
'local.xy' 4 0.01 'CI' 0.05 'Local'
EOF
genplt -XLEN 5 -YLEN 5 -X0 2 -Y0 1 -XMIN 1.0 -XMAX 7.5 -YMIN 1.0 -YMAX 7.5 -TX "ML (H)" -TY "Mw" -L lcmdfi
mv GENPLT.PLT LSLU.PLT
#####
# plot the SLU Mw vs ML
#####
cat > ecmdfil << EOF
'eml_H.xy' 0 0.01 'CI' 0.05 '1 sig'
EOF
genplt -XLEN 5 -YLEN 5 -X0 2 -Y0 1 -XMIN 1.0 -XMAX 7.5 -YMIN 1.0 -YMAX 7.5 -TX "ML (H)" -TY "Mw" -E ecmdfil
mv GENPLT.PLT ESLU.PLT
#####
# plot the SMSIM bounds as a polynomial. The first line says use the file SMSIMML.xy, plot in red, and fill
#    The second line says use the file SMSIMML.xy, use the color black, and only plot the outline of the polygon
#    The result is a shaded red area outlined in black
#####
cat > pcmdfil << EOF
'SMSIMML.xy' 2 1
'SMSIMML.xy' 1 0
EOF
#####
#    note that since there was no -LPOS "TR" on the command line, no legend is plotted.
genplt -XLEN 5 -YLEN 5 -X0 2 -Y0 1 -XMIN 1.0 -XMAX 7.5 -YMIN 1.0 -YMAX 7.5 -TX "ML (H)" -TY "Mw" -P pcmdfil
mv GENPLT.PLT SMSIM.PLT
#####
#      concatenate placing red on bottom, then results with error bars and then the local
#####
cat SMSIM.PLT ESLU.PLT LSLU.PLT > ALL.PLT



September 16, 2020

May 20, 2020

February 24, 2020



December 30, 2019

hprep96 -M model -d dfile 
hspec96
hpulse96 -V -p -l 1
hprep96 -M model -d dfile
rspec96
hpulse96 -V -p -l 1

February 19, 2019