[MSNoise] Question

[Sammarco, Carmelo]

Dear users,


I want find a way to insert into the CC header's  files  the pair stations distance for each CC.


In the file 'database_tools.py' I found in 'MISCS' this code:



############ MISCS ############


def azimuth(coordinates, x0, y0, x1, y1):
    if coordinates == "DEG":
        dist, azim, bazim = gps2DistAzimuth(y0, x0, y1, x1)
        # print dist, azim, bazi
        return azim
    elif coordinates == 'UTM':
        azim = 90. - np.arctan2((y1 - y0), (x1 - x0)) * 180. / np.pi
        # print azim
        return azim
    else:
        print "woooooow, please consider having a single coordinate system for all stations"
        return 0


So my idea was insert this code inside the 'correlation' section to insert these information in the header's files  like it happens with other parameters below:


 if sac_format == "doublets":
        tr.SetHvalue('A', 120)
    else:
        tr.SetHvalue('B', -maxlag)
        tr.SetHvalue('DEPMIN', np.min(corr))
        tr.SetHvalue('DEPMAX', np.max(corr))
        tr.SetHvalue('DEPMEN', np.mean(corr))
        tr.SetHvalue('SCALE', 1)
        tr.SetHvalue('NPTS', len(corr))
    tr.WriteSacBinary(filename)
    del st, tr
    return


I think also that insert this type of implementation in the code could be very helpful. I hope that someone can help me.


cheers,


Carmelo


-----------
Mr Carmelo Sammarco MSc  BSc  FGS
PhD student in Geology and Petroleum Geology
University of Aberdeen


The University of Aberdeen is a charity registered in Scotland, No SC013683.
Tha Oilthigh Obar Dheathain na charthannas clàraichte ann an Alba, Àir. SC013683.