def run_pecube_map(faultParamsListFiltered): [initA, accel, srA1, srA2] = faultParamsListFiltered netExtension = psc.calcNetExtension(initA, accel, srA1, srA2, dipFaultA) psc.modifyInputFiles(initA, accel, srA1, srA2, inputFile, outputFile) psc.renameFaultParams(inputFile, outputFile) pecube_print = psc.run_pecube_cloud() out_name = psc.save_output(comparison_file, initA, accel, srA1, srA2) out_array_name = psc.append_info(out_name, initA, accel, srA1, srA2, netExtension) #cloud.files.get(out_array_name, '/home/itchy/src/Pecube/nlrT2/picloud/{}'.format(out_array_name)) #cloud.files.get(out_array_name) return out_array_name
slipRate1FaultAs = [0., 0.25, 0.5, 0.75, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5] slipRate2FaultAs = [0.5, 0.75, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5] dipFaultA = [0.4716] faultParamsList = [[j, k, l, m] for j in initFaultAs for k in accelFaults for l in slipRate1FaultAs for m in slipRate2FaultAs] faultParamsListFiltered = [] for faultParams in faultParamsList: [initA, accel, srA1, srA2] = faultParams netExtension = psc.calcNetExtension(initA, accel, srA1, srA2, dipFaultA) if netExtension < 18: if netExtension >13: faultParamsListFiltered.append(faultParams) print len(faultParamsList), 'total params' print len(faultParamsListFiltered), 'filtered params' def run_pecube_map(faultParamsListFiltered): [initA, accel, srA1, srA2] = faultParamsListFiltered netExtension = psc.calcNetExtension(initA, accel, srA1, srA2, dipFaultA) psc.modifyInputFiles(initA, accel, srA1, srA2, inputFile, outputFile)