#
#                                            mlMagMeth2.props
#
#In database Applications table
#PROGID NAME         
#------ --------
#  9100 CISNml2      
#
#In database AppChannels table, config=1 for (highgain), config=2 for (lowgain) channels mapped to application
#PROGID NET      STA    SEE LO CONFIG ONDATE              OFFDATE            
#------ -------- ------ --- -- ------ ------------------- -------------------
#  9100 AZ       BZN    HGE -- 2      1900-01-01 00:00:00 3000-01-01 00:00:00
#  9100 AZ       BZN    HGN -- 2      1900-01-01 00:00:00 3000-01-01 00:00:00
#  9100 AZ       BZN    HHE -- 1      1900-01-01 00:00:00 3000-01-01 00:00:00
#  9100 AZ       BZN    HHN -- 1      1900-01-01 00:00:00 3000-01-01 00:00:00
#
#Set true to disable use of AppChannels table mapping of magmethod's progid (e.g. 9100)
disableAppChannelsMap = false
#
#When disableAppChannelsMap=true, use acceptChan seedchan as declared below for channel mag calc
acceptChan=HHE HHN HLN HLE HNN HNE HGE HGN
#
#When disableAppChannelsMap=true, use summaryChan seedchan declared below for summary mag calc
#to be used the seedchan template listed  must be a subset of acceptChan list template
#both HH and HL are used if included (unlike map option, where HL is culled when site HH is present)
summaryChan=HHE HHN HLN HLE HNN HNE HGE HGN
#
#Type of Butterworth secondary filter to apply to filter microseisms from low amp signals
bwFilterType=BANDPASS
#
#Set bwFilterOrder =0 to skip secondary BW bandpass filtering
bwFilterOrder=6
bwFilterLoFreq=0.5
bwFilterHiFreq=20.0
#When below true, bwFilterHiFreq reduced by distance =1./max(minPeriod, 0.05*exp(.0015*slantDist))
bwFilterScalePassBandByDistance=true
#Do not apply secondary filter if channel mag > bwMicroseismMaxMag
bwMicroseismMaxMag=4.5
#apply secondary filter if WAS amp < bwMicroseismMinAmp cm
bwMicroseismMinAmp=0.025
#
#WAS amp rejected or its derivation from waveform skipped when channel distance exceeds cutoff km
#Used to be if (magValue <= 2.0) max(minDistance(), (66.0*magValue - 2.0))
#else (magValue > 2.) cutoffKm = max(minDistance(), (cutoffDistSlope * magValue) - cutoffDistIntercept)
#cutoffDistIntercept=270.0
#cutoffDistSlope=200.0
#Now code states:
#Pivot at ML=2 (600+200)/200=4.0 (600+174/182)=4.25 (600+222/191)=4.3 (600+270)/200=4.35 (600+205)/170=4.74
# if (magValue <= cutoffPivotMag) {
#  cutoffKm = max(minDistance, ((cutoffPivotMagKm-cutoffMag0Km)*magValue/cutoffPivotMag + cutoffMag0Km))
# }
# else {
#  cutoffDistSlope = (cutoffMaxMagKm-cutoffPivotMagKm)/(cutoffMaxMag-cutoffPivotMag) 
#  cutoffDistIntercept = cutoffPivotMagKm - (cutoffDistSlope * cutoffPivotMag) 
#  cutoffKm = max(minDistance, (cutoffDistSlope * magValue + cutoffDistIntercept))
# }
cutoffMag0Km=25
cutoffPivotMag=2
cutoffPivotMagKm=130
cutoffMaxMag=4.35
cutoffMaxMagKm=600
#
#deleteInvalid WAS amp failing clipping, SNR, or period tests not saved to db 
deleteInvalid=true
#delete all rejected WAS for any reason (eg. seedchan type, no correction) not saved to db
deleteRejected=false
#deleteByDistanceTrim WAS amp exceeding maxDistance not saved to db
deleteByDistanceTrim=false
#deleteByResidualTrim WAS amp exceeding trimResidual not saved to db
deleteByResidualTrim=false
#
#disableMagnitudeDistanceCutoff distance rejection not scaled by event magnitude
disableMagnitudeDistanceCutoff=false
#filterCopiesWf true for jiggle, for hypomag false
filterCopiesWf=true
#globalMagCorr not used in Ml
globalMagCorr=0.0
#acceptable WA amp types for Ml
includedReadingTypes=WAS WASF WA WAU WAC 
#
#default HH WAS clipping amp cm, amp > maxVelAmp reject/delete
maxVelAmp=100.0
#default HH noise floor cm, amp < minAccAmp reject/delete 
minVelAmp=0.0010
#default HL WAS clipping amp cm, amp > maxAccAmp reject/delete
maxAccAmp=12000.0
#default HL noise floor cm, amp < minAccAmp reject/delete 
minAccAmp=0.0010
#
#WAS amp slantDist > maxDistance reject/delete 
#Note: a0(r) calibrated only r=500 im, so maxDistance=500 is set in code
#set method delete options to not delete distance trimmed/rejected amps
#to enable saving amps with slantDist > 500 for future A0 calibration
maxDistance=600.0
#
#WAS amp period < minPeriod always reject/delete 
minPeriod=0.040
#WAS amp period > maxPeriod always reject/delete
maxPeriod=5.0
#coefs for WAS period rejection scaled by distance
#do secondary microseism filter if WAS per > maxPeriodC1*exp(maxPeriodC2*wfSlantDist)
#if filtered WAS period > maxPeriodC1*exp(maxPeriodC2*slantDist) reject from summary
maxPeriodC1=0.70
maxPeriodC2=0.0022
#coefs for WAS period rejection scaled by distance, reject per < minPeriodC1*exp(minPeriodC2*slantDist)
minPeriodC1=0.040
minPeriodC2=0.0050
#
#When screening by distance(magnitude), don't reject WAS amp slantDist < minDistance
minDistance=20.0
#
#min SNR=(peakAmp/avgPreEventNoisePeak) acceptable for a channel mag calc
#where minSNR  < = minSummarySNR
minSNR=3.0
#
#min WAS SNR=(peakAmp/avgPreEventNoisePeak) acceptable for event summary mag calc
#where minSummarySNR  >= minSNR
minSummarySNR=5.0
#
#min valid readings needed for event summary mag calculation (channel mag count)
minValidReadings=5
#
#reject/delete WAS whose channel mag resid > trimResidual
trimResidual=1.0
#useAssignedWts=true,WAS amp inWgt = ampQuality*ml_magparms_view.summary_wgt,  else inWgt=1.
useAssignedWts=true
#simple_response gain required for calculation
requireGain=false
#filterChan not used for Ml only Md 
filterChan=
#maxChannels to use for event summary mag
maxChannels=999999
#lookup any missing channel metadata (corrections etc.) in db
channelDbLookUp=true
#use default value 0.5
chauvenetTrimValue=0.5
#require channel magcorr for summary mag calculation, set false if network has no corrections
requireCorrection=true
#scanAllWaveformsForMag=true always true for Ml, usually false for Md (only picked)
scanAllWaveformsForMag=true
#Don't scan entire trace for peak amplitude, only predicted energy window
scanEnergyWindow=true
#reject/delete WAS amp using summary magnitude statistics (ie. Chauvenet)
sumMagStatTrim=true
#velocity model used for determining energy window
velocityModelClassName=org.trinet.util.velocitymodel.HK_SoCalVelocityModel
#
#output for debugging 
debug=false
#detailed output of calculations
verbose=true