parser.add_option("--check",
action="store_true",
dest="isCheck",
default=False,
help="Random used for debugging")
(opt, args) = parser.parse_args()
# read in from stdin assuming psmeca format
Ndata = 0
while 1:
thisline = sys.stdin.readline()
if thisline != '':
# Read line as SymMT object
MT, extra = readpsmecaSm(thisline)
Ndata += 1
# get vals and vecs in E-N-Up coord system
(tval, t, bval, b, pval, p) = MT.getPTB()
# debugging option
if opt.isCheck:
# if( sqrt2*MT.mhat[2] > 1 ):
# print NP.sum( MT.mhat[0:3]*MT.mhat[0:3] ) + 2*NP.sum( MT.mhat[3:6]*MT.mhat[3:6] )
print MT.mhat[0:3]
# print thisline.rstrip('\n')
continue
# optional output the position
if opt.isPos:
parser.add_option("--mw",action="store_true", dest="isMw", default=False,
help="Get moment magnitude")
parser.add_option("--check",action="store_true", dest="isCheck", default=False,
help="Random used for debugging")
(opt, args)=parser.parse_args()
# read in from stdin assuming psmeca format
Ndata = 0
while 1:
thisline = sys.stdin.readline()
if thisline != '':
# Read line as SymMT object
MT, extra = readpsmecaSm( thisline )
Ndata += 1
# get vals and vecs in E-N-Up coord system
(tval, t, bval, b, pval, p) = MT.getPTB()
# debugging option
if opt.isCheck:
# if( sqrt2*MT.mhat[2] > 1 ):
# print NP.sum( MT.mhat[0:3]*MT.mhat[0:3] ) + 2*NP.sum( MT.mhat[3:6]*MT.mhat[3:6] )
print MT.mhat[0:3]
# print thisline.rstrip('\n')
continue
# optional output the position
if opt.isPos:
help="Random used for debugging")
(opt, args) = parser.parse_args()
# read in from stdin assuming psmeca format
Ndata = 0
lcount = 0
while 1:
thisline = sys.stdin.readline()
lcount += 1
if thisline != '':
# Read line as SymMT object
try:
MT, extra = readpsmecaSm(thisline, lcount)
Ndata += 1
except:
sys.exit()
continue
# debugging option
if opt.isCheck:
# check decomposition is working
(MTdc, MTclvd, MTiso) = MT.decompose(2)
print "ALL:", MT.getMvec()
print "sum:", MTdc.getMvec() + MTclvd.getMvec() + MTiso.getMvec()
print "DC:", MTdc.getMvec(), MTdc.getFclvd()
print "CLVD:", MTclvd.getMvec(), MTclvd.getFclvd()
print "ISO:", MTiso.getMvec()
(opt, args)=parser.parse_args()
# read in from stdin assuming psmeca format
Ndata = 0
lcount = 0
while 1:
thisline = sys.stdin.readline()
lcount += 1
if thisline != '':
# Read line as SymMT object
try:
MT, extra = readpsmecaSm( thisline, lcount )
Ndata += 1
except:
sys.exit()
continue
# debugging option
if opt.isCheck:
# check decomposition is working
(MTdc, MTclvd, MTiso) = MT.decompose( 2 )
print "ALL:", MT.getMvec()
print "sum:", MTdc.getMvec() + MTclvd.getMvec() + MTiso.getMvec()
print "DC:",MTdc.getMvec(), MTdc.getFclvd()
print "CLVD:",MTclvd.getMvec(), MTclvd.getFclvd()
