In [1]:

    

%run basics


    

In [2]:

    

ncname=qcio.get_filename_dialog(path="../../Sites")


    

In [3]:

    

ds=qcio.nc_read_series(ncname)
Fc_SOLO,f,a=qcutils.GetSeriesasMA(ds,"Fc_SOLO")
Fc_SOLO_mask=numpy.ma.getmaskarray(Fc_SOLO)
idx_array = qcutils.contiguous_regions(Fc_SOLO_mask)


    

In [16]:

    

print idx_array.shape,idx_array.shape[0],idx_array.shape[1]


    

    




(6, 2) 6 2

In [5]:

    

print idx_array


    

    




[[  7297  11568]
 [ 40849  42336]
 [ 59857  61248]
 [ 62737  64176]
 [ 86113  89040]
 [103633 105120]]

In [6]:

    

si_gap=idx_array[0,0]
ei_gap=idx_array[0,1]


    

In [14]:

    

ldt = ds.series["DateTime"]["Data"]
Fc,f,a=qcutils.GetSeriesasMA(ds,"Fc")


    

In [15]:

    

si_gfperiod = si_gap
ei_gfperiod = ei_gap
min_points = int((ei_gfperiod-si_gfperiod)/2)
num_good_points = numpy.ma.count(Fc[si_gfperiod:ei_gfperiod])
print min_points,num_good_points
while num_good_points<min_points:
    si_gfperiod = si_gfperiod - 48
    ei_gfperiod = ei_gfperiod + 48
    min_points = int((ei_gfperiod-si_gfperiod)/2)
    num_good_points = numpy.ma.count(Fc[si_gfperiod:ei_gfperiod])
    print ldt[si_gfperiod],ldt[ei_gfperiod],num_good_points,min_points


    

    




2135 0
2009-02-27 00:30:00 2009-05-29 00:00:00 96 2183
2009-02-26 00:30:00 2009-05-30 00:00:00 192 2231
2009-02-25 00:30:00 2009-05-31 00:00:00 288 2279
2009-02-24 00:30:00 2009-06-01 00:00:00 384 2327
2009-02-23 00:30:00 2009-06-02 00:00:00 480 2375
2009-02-22 00:30:00 2009-06-03 00:00:00 576 2423
2009-02-21 00:30:00 2009-06-04 00:00:00 672 2471
2009-02-20 00:30:00 2009-06-05 00:00:00 768 2519
2009-02-19 00:30:00 2009-06-06 00:00:00 864 2567
2009-02-18 00:30:00 2009-06-07 00:00:00 960 2615
2009-02-17 00:30:00 2009-06-08 00:00:00 1056 2663
2009-02-16 00:30:00 2009-06-09 00:00:00 1152 2711
2009-02-15 00:30:00 2009-06-10 00:00:00 1248 2759
2009-02-14 00:30:00 2009-06-11 00:00:00 1344 2807
2009-02-13 00:30:00 2009-06-12 00:00:00 1440 2855
2009-02-12 00:30:00 2009-06-13 00:00:00 1536 2903
2009-02-11 00:30:00 2009-06-14 00:00:00 1632 2951
2009-02-10 00:30:00 2009-06-15 00:00:00 1728 2999
2009-02-09 00:30:00 2009-06-16 00:00:00 1824 3047
2009-02-08 00:30:00 2009-06-17 00:00:00 1920 3095
2009-02-07 00:30:00 2009-06-18 00:00:00 2016 3143
2009-02-06 00:30:00 2009-06-19 00:00:00 2112 3191
2009-02-05 00:30:00 2009-06-20 00:00:00 2208 3239
2009-02-04 00:30:00 2009-06-21 00:00:00 2304 3287
2009-02-03 00:30:00 2009-06-22 00:00:00 2400 3335
2009-02-02 00:30:00 2009-06-23 00:00:00 2496 3383
2009-02-01 00:30:00 2009-06-24 00:00:00 2592 3431
2009-01-31 00:30:00 2009-06-25 00:00:00 2688 3479
2009-01-30 00:30:00 2009-06-26 00:00:00 2784 3527
2009-01-29 00:30:00 2009-06-27 00:00:00 2880 3575
2009-01-28 00:30:00 2009-06-28 00:00:00 2976 3623
2009-01-27 00:30:00 2009-06-29 00:00:00 3072 3671
2009-01-26 00:30:00 2009-06-30 00:00:00 3168 3719
2009-01-25 00:30:00 2009-07-01 00:00:00 3264 3767
2009-01-24 00:30:00 2009-07-02 00:00:00 3360 3815
2009-01-23 00:30:00 2009-07-03 00:00:00 3456 3863
2009-01-22 00:30:00 2009-07-04 00:00:00 3552 3911
2009-01-21 00:30:00 2009-07-05 00:00:00 3648 3959
2009-01-20 00:30:00 2009-07-06 00:00:00 3744 4007
2009-01-19 00:30:00 2009-07-07 00:00:00 3840 4055
2009-01-18 00:30:00 2009-07-08 00:00:00 3936 4103
2009-01-17 00:30:00 2009-07-09 00:00:00 4032 4151
2009-01-16 00:30:00 2009-07-10 00:00:00 4128 4199
2009-01-15 00:30:00 2009-07-11 00:00:00 4224 4247
2009-01-14 00:30:00 2009-07-12 00:00:00 4320 4295

In [ ]: