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

[Georgia Cua]

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








-----------------------------
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
-----------------------------



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20101215/cea56efd/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: SwissNGA.pm
Type: text/x-perl-script
Size: 13996 bytes
Desc: SwissNGA.pm
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20101215/cea56efd/attachment.bin>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20101215/cea56efd/attachment-0001.html>