f_table=open('run_MECA.tbl','w')
f_table.write('Cluster# lon lat depth num_polarities STR DIP RAKE error_F error_A pol_mis STDR\n')
##### Go through each index_file in a loop
for index_file in list_index:
clus_number=int(index_file[-8:-4])
f_log.write('Processing Cluster %i associated with file %s\n' %(clus_number,index_file))
if clus_number>22:
continue
##### Run HASH 1st time to get number of polarity pspolar.inp
print 'running HASH for cluster number %i' %clus_number
run_HASH(data_path+index_file,parameter=None)
os.system('read_FOCMEC.py focmec.dat')
f_scratch=open('foc_polar.out','r')
total_pol=len(f_scratch.readlines())
f_scratch.close()
f_log.write('Total number of polarities: %i \n' %total_pol)
print total_pol
if total_pol!=0:
#### RUN HASH 2nd time with correct parameter
#### Define parameters for HASH
grid_angle=5
max_pol_error=np.floor(0.1*total_pol)
os.system(cmd)
list_index=glob.glob('index*[0-9][0-9][0-9][0-9].out')
print list_index
####### Start loop in index files
for file in list_index:
event_number=int(file[-8:-4])
###### Compute Focal Mechanism
##### Run HASH 1st time to get number of polarity pspolar.inp
print 'running HASH 1st time to get number of polarities'
run_HASH('./'+file,parameter=None)
##### Read hyp.out to get lon/lat/depth of event
cmd='nor2xyz hyp.out > tmp/temp_hyp.xyz'
os.system(cmd)
fat=open('tmp/temp_hyp.xyz','r')
for line in fat:
B=line.split()
lon=float(B[0])
lat=float(B[1])
depth=float(B[2])
print lon,lat,depth
fat.close()
os.remove('tmp/temp_hyp.xyz')
