Background
psuedo Spectral Accelerations (pSA) values are calculated in a cython script called rspectra.
It was found that the RotD values for pSA intensity measures were smaller than expected. Upon investigation it was determined that the spectral displacements were being returned instead of the pSA values.
This did not effect components such as the basic components and geom as they were calculated and returned separately.
Changes made
In order to get back the correct pSA values, the variable u was multiplied buy the square of the circular frequency (ω) before being returned. As the maximum pSA value for each period can be determined from this data, this return value has been removed from the script.
The maximum pSA values are now calculated outside of the rspectra script. Any future issues with pSA values should be reflected here.
Further changes in the same pull request include name changes to clarify that the RotD component can be taken of any intensity measure that takes an acceleration waveform as an input.
Verification
This change was tested using a copy of a Hossack broadband binary that was available to the tester.
Because the method of calculating the geom component was changed this is also included in the below table
While only one station is presented, this held true for all stations in the resulting intensity measure files.
Source | station | component | pSA_0.02 | pSA_0.05 | pSA_0.1 | pSA_0.2 | pSA_0.3 | pSA_0.4 | pSA_0.5 | pSA_0.75 | pSA_1.0 | pSA_2.0 | pSA_3.0 | pSA_4.0 | pSA_5.0 | pSA_7.5 | pSA_10.0 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Original rspectra | WRAV | geom | 0.06013336675273 | 0.089115745322454 | 0.125133181113855 | 0.121005303566589 | 0.144512778030697 | 0.103801783325256 | 0.080523663769336 | 0.067167579184118 | 0.028598548143978 | 0.01028351531141 | 0.004587755557469 | 0.001993342513191 | 0.001102070004063 | 0.000443169802414 | 0.000250627974779 |
Original rspectra | WRAV | rotd50 | 6.2335278555554E-07 | 5.72149302246633E-06 | 2.96559252638543E-05 | 0.000149028614283 | 0.000344619064028 | 0.00042019840063 | 0.000535445191359 | 0.000950297921658 | 0.000755494339119 | 0.001070922696095 | 0.001041079503484 | 0.000849442741323 | 0.000697894240638 | 0.000575025983932 | 0.000597973586903 |
Original rspectra with ω2 | WRAV | geom | 0.06013336675273 | 0.089115745322454 | 0.125133181113855 | 0.121005303566589 | 0.144512778030697 | 0.103801783325256 | 0.080523663769336 | 0.067167579184118 | 0.028598548143978 | 0.01028351531141 | 0.004587755557469 | 0.001993342513191 | 0.001102070004063 | 0.000443169802414 | 0.000250627974779 |
Original rspectra with ω2 | WRAV | rotd50 | 0.061522453957503 | 0.090350196344537 | 0.117076900201006 | 0.147085346741202 | 0.151166836934552 | 0.103679799604779 | 0.084554115474825 | 0.066695570132952 | 0.02982572101747 | 0.010569583354605 | 0.004566685710875 | 0.002095915954559 | 0.001102070411023 | 0.000403575394258 | 0.000236070509801 |
New rspectra | WRAV | geom | 0.06013336675273 | 0.089115745322454 | 0.125133181113855 | 0.121005303566589 | 0.144512778030697 | 0.103801783325256 | 0.080523663769336 | 0.067167579184118 | 0.028598548143978 | 0.01028351531141 | 0.004587755557469 | 0.001993342513191 | 0.001102070004063 | 0.000443169802414 | 0.000250627974779 |
New rspectra | WRAV | rotd50 | 0.061522453957503 | 0.090350196344537 | 0.117076900201006 | 0.147085346741202 | 0.151166836934552 | 0.103679799604779 | 0.084554115474825 | 0.066695570132952 | 0.02982572101747 | 0.010569583354605 | 0.004566685710875 | 0.002095915954559 | 0.001102070411023 | 0.000403575394258 | 0.000236070509801 |