From 9fe84e07b3a37ed017b04a816f0a95d9f322e137 Mon Sep 17 00:00:00 2001 From: caes Date: Mon, 30 Jan 2017 03:09:59 -0500 Subject: [PATCH] new frequency bin paradigm, should be more robust# Please enter the commit message for your changes. Lines starting --- scripts/analyze_lightcurve.py | 9 +++------ 1 file changed, 3 insertions(+), 6 deletions(-) diff --git a/scripts/analyze_lightcurve.py b/scripts/analyze_lightcurve.py index 0f4feb7..7512b94 100644 --- a/scripts/analyze_lightcurve.py +++ b/scripts/analyze_lightcurve.py @@ -18,7 +18,7 @@ except getopt.GetoptError: # Time resolution determined from inspection and testing. This script # does not expect evenly spaced data in time. -dt = 0.01 +dt = 0.1 #### Get the psd for the first light curve #### @@ -54,12 +54,9 @@ fqL = np.array([0.0049999999, 0.018619375, 0.044733049, 0.069336227, 0.10747115, # fqL = np.concatenate(([0.5/seg_length], fqL)) # -f1 = 1/270. +f1 = 1/175. f2 = 0.5/dt - -fqL = np.hstack((np.array(0.5*f1),np.logspace(np.log10(0.9*f1),np.log10(.5*f2),7),np.array(2*f2))) - - +fqL = np.hstack((np.array(0.5*f1),np.logspace(np.log10(0.9*f1),np.log10(0.3*f2),9),np.array(2*f2))) nfq = len(fqL) - 1 fqd = 10**(np.log10( (fqL[:-1]*fqL[1:]) )/2.)