Ejemplo n.º 1
0
                                     paz,
                                     db_bins=(-80, 20, 1.0),
                                     period_limits=[2, 1200],
                                     special_handling="ringlaser")
                        Zppsd = PPSD(z[0].stats,
                                     paz,
                                     db_bins=(-80, 20, 1.0),
                                     period_limits=[2, 1200],
                                     special_handling="ringlaser")

                    #Done, update variable
                    ppsd_init = True

                #Now check if the start time of the current trace is already in the
                #object, if it is not, tne add the trace, otherwise skip
                if Eppsd._PPSD__check_time_present(
                        e[0].stats.starttime) == False:

                    #Did you find some new data?
                    updated_ppsd = True

                    #process the traces to zero out gaps
                    e = analysis.prepare_for_ppsd(e)
                    n = analysis.prepare_for_ppsd(n)
                    z = analysis.prepare_for_ppsd(z)

                    Eppsd.add(e)
                    Nppsd.add(n)
                    Zppsd.add(z)

        #Done with that site, save the god danged ppsd
        if updated_ppsd == True: