UO-Geophysics / MudPy

Kinematic and static rupture forward modeling and inversion code
Other
116 stars 63 forks source link

Running insar and gps statics runs into issues with dataframe #89

Open crowellbw opened 1 year ago

crowellbw commented 1 year ago

within parallel.py, if you are creating static synthetics for both SAR and GPS, the code ends up overwriting the statics file, i.e. the following .to_csv lines. I added a mode='a' tag within the write function, but this is not ideal since you need to delete all GFs then recompute since you could end up with duplicates and I don't think the code is indexed correctly to handle this downstream.

    if write_df==True and static ==1: #Note to self: stop using 0,1 and swithc to True/False

        #Strike slip
        SSdf = df(data=None, index=None, columns=['staname','n','e','u','beta'])
        SSdf.staname=staname
        SSdf.n=staticsSS[:,0]
        SSdf.e=staticsSS[:,1]
        SSdf.u=staticsSS[:,2]
        SSdf.beta=staticsSS[:,3]
        SSdf.to_csv(green_path+'subfault'+num+'.SS.static.neu',sep='\t',index=False,header=False)

        DSdf = df(data=None, index=None, columns=['staname','n','e','u','beta'])     
        DSdf.staname=staname
        DSdf.n=staticsDS[:,0]
        DSdf.e=staticsDS[:,1]
        DSdf.u=staticsDS[:,2]
        DSdf.beta=staticsDS[:,3]
        DSdf.to_csv(green_path+'subfault'+num+'.DS.static.neu',sep='\t',index=False,header=False)