def run_pecube_map(faultParamsListFiltered):
[initA, initB, accel, srA1, srA2, srB1, srB2] = faultParamsListFiltered
#netExtension = psc.calcNetExtension(initA, accel, srA1, srA2, dipFaultA)
netExtensionA = psc.calcNetExtension(initA, accel, srA1, srA2, dipFaultA)
netExtensionB = psc.calcNetExtension(initB, accel, srB1, srB2, dipFaultB)
netExtension = netExtensionA + netExtensionB
psc.modifyInputFiles(initA, initB, accel, srA1, srA2, srB1, srB2, inputFile,
outputFile)
psc.renameFaultParams(inputFile, outputFile)
pecube_print = psc.run_pecube_cloud()
out_name = psc.save_output(comparison_file, initA, initB, accel, srA1,
srA2, srB1, srB2, netExtension)
out_array_name = psc.append_info(out_name, initA, initB, accel, srA1, srA2,
srB1, srB2, netExtension)
return out_array_name
t1 = time.time()
print t1-t0, 's to generate fault list'
print len(faultParamsList)
# filter variables
faultParamsListFiltered = []
nlrT4_run_list = []
for faultParams in faultParamsList:
[initA, initB, accel, srA1, srA2, srB1, srB2] = faultParams
netExtensionA = psc.calcNetExtension(initA, accel, srA1, srA2, dipFaultA)
netExtensionB = psc.calcNetExtension(initB, accel, srB1, srB2, dipFaultB)
netExtension = netExtensionA + netExtensionB
#print netExtension
if 5 < netExtension < 7.5 and initA > accel and initB > accel and srB1 == srB2:
faultParamsListFiltered.append(faultParams)
run_name = 'nlrT4_{}_{}_{}_{}_{}_{}_{}.npy'.format(initA, initB, accel, srA1, srA2, srB1, srB2)
nlrT4_run_list.append(run_name)
t2 = time.time()
