[Shake-dev] developing new GMPE module for ShakeMap v3.5

[Vince Quitoriano]

Hi Georgia,

I would look at MA2005.pm, at least the way it implements 
siteCorrectList and _makeAmps(). I bet your problem between those two 
functions.

--
Vince Quitoriano (vinceq at usgs.gov)
USGS Did You Feel It?	http://earthquake.usgs.gov/earthquakes/dyfi
USGS ShakeMaps:		http://earthquake.usgs.gov/earthquakes/shakemap

On 12/16/10 4:25 AM, Georgia Cua wrote:
> Hi Vince,
>
> Thanks for the hint.
>
> I replaced the subroutine siteCorrectList in my new module with what you
> attached (I'm not sure whether that was what you were suggesting), and
> have a different error now
>
> shake at phantom:~/opt/shakemap/bin$ ./shake -event newtest
> grind: Predictive map: no data in
> '/rdsk/projects/shake/opt/shakemap-3.5-vanilla/bin/../data/newtest/input'
> grind: Using fixed boundary 5.6 45.6 10.8 47.8
> grind: 47.8:45.6:10.8:5.6
> grind:
>
> arg <lats> must be an array ref at
> /rdsk/projects/shake/opt/shakemap-3.5-vanilla/bin/../perl/lib/Shake/GMPE/SwissNGA.pm
> line 309
> shake: Error in grind: 65280
>
> (originally, the error had to do with <lons>, now it is <lats>)
>
> Do you have a GMPE module that uses Rjb and HAZUS for site corrections?
> Maybe I could use that as my starting template instead? (Of the modules
> included in the distribution, only MA2005.pm seemed to use HAZUS, and it
> does not use Rjb.)
>
> Thanks,
> Georgia
>
>
>
> On Dec 15, 2010, at 4:43 PM, Vince Quitoriano wrote:
>
>>
>> I've run into this problem before when using the site correction function
>> in HazusPGV.pm instead of rolling my own. In HazusPGV.pm, the
>> functions siteCorrectList and siteCorrectGrid will throw an exception when
>> they get no points, so you need to catch that possibility before you call
>> them from within your own module. If you look at the other GMPEs that
>> use HAZUS, they all say something like:
>>
>> sub siteCorrectList {
>> my ($self, $type, $direction, $amps, $lons, $lats, $vels, $outdir,
>> $gmtbin, $verbose) = @_;
>> return [] if (!defined $amps or @$amps <= 0); #### Avoid exception
>>
>> #-------------------------------------------------------------------------
>> # Call the HazusPGV module to do the site corrections
>> #-------------------------------------------------------------------------
>> return $self->{hazus}->siteCorrectList(@_);
>> }
>>
>> Hope that helps,
>> --
>> Vince Quitoriano (vinceq at usgs.gov <mailto:vinceq at usgs.gov>)
>> USGS Did You Feel It?http://earthquake.usgs.gov/earthquakes/dyfi
>> USGS ShakeMaps:http://earthquake.usgs.gov/earthquakes/shakemap
>>
>> On 12/15/10 5:57 AM, Georgia Cua wrote:
>>> Dear Bruce,
>>>
>>> I'm trying to code a new GMPE module for our ShakeMap v3.5 installation
>>> in Switzerland. It has the functional form of Cua and Heaton
>>> (unpublished) constrained with Swiss weak motion data and NGA strong
>>> motion data, and should be valid in the magnitude range 2 <= M <= 8.
>>>
>>> I have written the new module (SwissNGA.pm) based on Small.pm (because
>>> it uses Rjb as a distance measure) and MA2005.pm (because it does not
>>> have a GMPE-spefic site amplification approach, and should thus use the
>>> Hazus module).
>>>
>>> Essentially the new module is Small.pm with a few lines/subroutines
>>> replaced with portions of MA2005.pm (from Constructor, and subroutines
>>> siteCorrectGrd, siteCorrectList), and of course, changing the
>>> calculation of the regression equation.
>>>
>>> In theory, should this result in a working module, or is the basic
>>> approach flawed?
>>>
>>> I have installed the new module and tried running it on a scenario event
>>> (just event information, no PGM observations), resulting in the
>>> following error
>>>
>>> shake at phantom:~/opt/shakemap/bin$ ./shake -event newtest -verbose
>>> shake: Running:
>>> '/rdsk/projects/shake/opt/shakemap-3.5-vanilla/bin/../bin/grind -event
>>> newtest -qtm -lonspan 5.5 -boundcheck -xml -lonspan 5.5 -psa'
>>> grind: Predictive map: no data in
>>> '/rdsk/projects/shake/opt/shakemap-3.5-vanilla/bin/../data/newtest/input'
>>> grind: Using fixed boundary 5.6 45.6 10.8 47.8
>>> grind: 47.8:45.6:10.8:5.6
>>> grind:
>>>
>>> arg <lons> must be an array ref at
>>> /rdsk/projects/shake/opt/shakemap-3.5-vanilla/bin/../perl/lib/Shake/GMPE/SwissNGA.pm
>>> line 306
>>> shake: Error in grind: 65280
>>>
>>> Is the approach of basing the new module primarily on Small.pm (with
>>> portions from MA2005.pm) a valid approach, or are there some subtleties
>>> that I am completely bulldozing over? Do you have any more detailed
>>> documentation/instructions for writing new modules?
>>>
>>> Sorry if this is not the type of question to be sending to the developer
>>> list!
>>>
>>> Thanks,
>>> Georgia
>>>
>>> PS. the new module SwissNGA.pm is attached
>>>
>>>
>>>
>>> _______________________________________________
>>> Shake-dev mailing list
>>> Shake-dev at geohazards.usgs.gov
>>> https://geohazards.usgs.gov/mailman/listinfo/shake-dev
>> _______________________________________________
>> Shake-dev mailing list
>> Shake-dev at geohazards.usgs.gov
>> https://geohazards.usgs.gov/mailman/listinfo/shake-dev
>
>
>
> -----------------------------
> Dr. Georgia Cua
> Swiss Seismological Service
> Institute of Geophysics
> ETH Zurich, NO H61
> Sonneggstr 5
> CH-8092 Zurich
> Switzerland
>
> Tel: +41-44-633-7574
> Fax: +41-44-633-1065
>
> georgia.cua at sed.ethz.ch <mailto:georgia.cua at sed.ethz.ch>
> -----------------------------
>
>
>