clc
clear all
filename = 's_meter_log_2008_06_23-11_55_31_64_437432000.txt '; % day 1.5, 0.21 Hz, 0.42 Hz
interval = 2500:2700;
filename = 's_meter_log_2008_07_06-22_12_02_70_437432000.txt'; % day 13.9, 0.387 Hz, 0.7742 Hz, 1.16 Hz
interval = 3098:3170;
filename = 's_meter_log_2008_07_18-11_40_00_86_437432000.txt'; % day 25.5,  0.51455 Hz  = 30.873 RPM, 0.99333 Hz = 59.59 RPM
interval = 2300:3300;
filename = 's_meter_log_2008_07_27-11_15_09_0_437432000.txt'; % day 34.5, 0.0348 Hz, 0.592 Hz, 1.173 Hz
interval = 3250:3900;
filename = 's_meter_log_2008_07_30-10_34_53_27_437432000.txt'; % day 37.4, 1.195Hz 71.7 RPM
interval = 2234:2330;
filename = 's_meter_log_2008_08_04-21_34_50_57_437432000.txt'; % day 42.8, 0.646 Hz, 1.292 Hz
interval = 2343:2365;
filename = 's_meter_log_2008_08_04-23_11_44_24_437432000.txt'; % day 42.9, 1.304 Hz
interval = 2968:3006;
filename = 's_meter_log_2008_08_05-12_26_38_37_437432000.txt'; % day 43.5, 0.65 Hz, 1.315 Hz
interval = 3060:4170;
filename = 's_meter_log_2008_08_06-11_08_54_52_437432000.txt';% day 44.5, 0.6697 Hz, 1.3394 Hz
interval = 3820:4950;
filename = 's_meter_log_2008_08_12-11_24_23_70_437432000.txt';% day 50.5, 0.72 Hz, 1.4 Hz
interval = 2910:4000;
filename = 's_meter_log_2008_08_18-11_39_52_88_437432000.txt';% day 56.5, 0.7185 Hz, 1.44Hz
interval = 3600:4600;
filename = 's_meter_log_2008_08_21-10_59_29_43_437432000.txt';% day 59.5, 0.715, 1.395
interval = 2000:3000;
filename = 's_meter_log_2008_08_22-11_18_05_61_437432000.txt';% day 60.5, 0.7(0.0), 1.395 
interval = 1400:2400;
filename = 's_meter_log_2008_08_25-10_44_42_28_437432000.txt';% day 63.5, 0.7, 1.395 
interval = 200:1200;
filename = 's_meter_log_2008_08_25-21_40_43_64_437432000.txt';% day 63.9, 0.6815, 1.395
interval = 3100:4100;
filename = 's_meter_log_2008_08_26-21_59_17_88_437432000.txt';% day 64.9, 0.65, 1.35
interval = 2650:3650;
filename = 's_meter_log_2008_08_28-11_33_32_80_437432000.txt'; % day 66.5, 0.61, 1.225
interval = 2700:3200;
filename = 's_meter_log_2008_08_29-10_16_04_20_437432000.txt';% day 67.4, 0.48, 0.959
interval = 2300:3300;
filename = 's_meter_log_2008_08_29-11_52_09_72_437432000.txt';% day 67.5, 0.46, 0.92
interval = 3000:4000;
filename = 's_meter_log_2008_08_29-21_25_50_42_437432000.txt';% day 67.9, 0.431 0.855
interval = 200:1200;
filename = 's_meter_log_2008_08_31-21_56_04_83_437432000.txt';% day 69.9, 0.292, 0.584
interval = 2400:3400;
filename = 's_meter_log_2008_09_01-11_11_43_54_437432000.txt';% day 70.5, 0.28, 0.566
interval = 2400:3400;
filename = 's_meter_log_2008_09_01-22_14_41_70_437432000.txt';% day 70.9, 0.264, 0.57
interval = 2300:3300;
filename = 's_meter_log_2008_09_03-21_15_57_40_437432000.txt';% day 72.9, 0.258, 0.53
interval = 2900:3900;
filename = 's_meter_log_2008_09_05-10_49_56_35_437432000.txt';% day 74.5, 0.233, 0.474
interval = 2600:3600;
filename = 's_meter_log_2008_09_05-21_52_52_78_437432000.txt';% day 74.9, 0.23, 0.462
interval = 2700:3700;
filename = 's_meter_log_2008_09_07-11_29_18_72_437432000.txt';% day 76.5, 0.22, 0.464
interval = 3600:4400;
filename = 's_meter_log_2008_09_09-12_04_31_58_437432000.txt';% day 78.5, 0.25, 0.58
interval = 3300:4300;
filename = 's_meter_log_2008_09_09-21_31_18_53_437432000.txt';% day 78.9, 0.256, 0.512
interval = 2850:3350;
filename = 's_meter_log_2008_09_10-21_49_46_74_437432000.txt';% day 79.9, 0.24, 0.489
interval = 2100:2700;
filename = 's_meter_log_2008_09_12-22_27_02_56_437432000.txt';% day 81.9, 0.17, 0.38
interval = 2060:2600;
filename = 's_meter_log_2008_09_13-21_09_45_35_437432000.txt';% day 82.9, 0.17, 0.4
interval = 1590:2060;
filename = 's_meter_log_2008_09_16-22_05_09_83_437432000.txt';% day 85.9, 0.27, 0.54
interval = 2390:2530;
interval = 3450:3950;
filename = 's_meter_log_2008_09_17-22_23_48_60_437432000.txt';% day 86.9, 0.255, 0.53
interval = 2400:3400;
filename = 's_meter_log_2008_09_18-21_06_38_33_437432000.txt';% day 87.9, 0.245, 0.478
interval = 3100:3700;
filename = 's_meter_log_2008_09_19-10_21_59_21_437431000.txt';% day 90.4, 0.22, 0.454
intereval = 3000:4000;
filename = 's_meter_log_2008_09_19-21_24_58_46_437432000.txt';% day 90,9, 0.212, 0.417
interval = 2866:3800; %frequencies =   0.0160    0.2120    0.4170
filename = 's_meter_log_2008_09_20-10_40_29_28_437432000.txt';% day 91.4
interval = 1402:2487; % frequencies =  0.0140    0.1830    0.3730
filename = 's_meter_log_2008_09_20-12_16_44_46_437432000.txt';% day 91.5
interval = 2320:3377;% frequencies =   0.0140    0.1720    0.3570
filename = 's_meter_log_2008_09_20-21_43_24_66_437432000.txt';% day 91.9
interval = 1945:2635;% frequencies =    0.0210    0.1660    0.3320
filename = 's_meter_log_2008_09_21-10_59_03_40_437432000.txt';% day 92.5
interval = 1684:2768;% frequencies =    0.0130    0.1670    0.3420
filename = 's_meter_log_2008_09_21-22_01_56_88_437432000.txt';% day 92.9
interval = 1676:2743; %frequencies =    0.0130    0.1740    0.3810
filename = 's_meter_log_2008_09_22-11_17_38_59_437432000.txt';% day 93.45
interval = 1789:2917;% frequencies =    0.0130    0.1810    0.3420
filename = 's_meter_log_2008_09_22-12_54_03_25_437432000.txt';% day 93.55
interval = 2242:3281;% frequencies =    0.0140    0.1750    0.3560
filename = 's_meter_log_2008_09_22-20_45_13_25_437432000.txt';% day 93.7
interval = 2166:2522;% frequencies =    0.0140    0.1630    0.3250   
filename = 's_meter_log_2008_09_22-22_20_37_64_437432000.txt';% day 93.85
interval = 3604:3815;% frequencies =    0.0140    0.1770    0.3180    
filename = 's_meter_log_2008_09_23-19_30_05_0_437432000.txt';% day 94.8
interval = 1458:1856; % frequencies =     0.0370    0.1490    0.2970
filename = 's_meter_log_2008_09_23-22_39_21_44_437432000.txt';% day 94.9
interval = 1527:2621;% frequencies =    0.0200    0.1470    0.3070
filename = 's_meter_log_2008_09_24-10_18_47_20_437432000.txt';% day 95.4
interval = 1179:2293;% frequencies =    0.0200    0.1520    0.2980
filename = 's_meter_log_2008_09_24-11_54_51_70_437432000.txt';% day 95.5
interval = 2415:3485;% frequencies =    0.0200    0.1450    0.2970
filename = 's_meter_log_2008_09_24-21_21_46_43_437432000.txt';% day 95.8
interval = 2065:2343;% frequencies =    0.1570    0.2880    
filename = 's_meter_log_2008_09_24-22_58_12_30_437432000.txt';% day 95.9
interval = 1906:2975;% frequencies =    0.0140    0.1430    0.2920
filename = 's_meter_log_2008_09_30-12_10_14_52_437432000.txt';% day 101.5
interval = 2330:3370;% frequencies =    0.0200    0.1690    0.3310
filename = 's_meter_log_2008_09_30-21_37_01_57_437432000.txt';% day 101.8
interval = 2532:2748;% frequencies =    0.170    0.3430   
filename = 's_meter_log_2008_09_30-23_13_56_24_437432000.txt';% day 101.9
interval = 2728:3242;% frequencies =    0.0690    0.1800    0.3470
filename = 's_meter_log_2008_10_01-20_20_49_17_437432000.txt';% day 102.75
interval = 2198:3265;% frequencies =    0.0130    0.1950    0.3970
filename = 's_meter_log_2008_10_02-11_11_10_51_437432000.txt';% day 103.45
interval = 2893:3983;% frequencies =    0.0130    0.2150    0.4300
filename = 's_meter_log_2008_10_02-20_38_56_22_437432000.txt';% day 103.75
interval = 2073:2945;% frequencies =    0.0160    0.2120    0.4560
filename = 's_meter_log_2008_10_03-11_29_46_73_437432000.txt';% day 104.5
interval = 2763:3793;% frequencies =    0.0140    0.2340    0.4540
filename = 's_meter_log_2008_10_03-20_57_08_28_437432000.txt'; % day 104.8
interval = 1997:2998;% frequencies =    0.0130    0.2360    0.4720
filename = 's_meter_log_2008_10_05-21_33_50_54_437432000.txt';% day 105.8
interval = 2070:2704;% frequencies =    0.0200    0.2600    0.5490
filename = 's_meter_log_2008_10_05-23_10_37_25_437432000.txt';% day 105.9
interval = 1663:1992;% frequencies =    0.0610    0.2660    0.5530
filename = 's_meter_log_2008_10_06-21_52_18_76_437432000.txt';% day 106.8
interval = 2091:2946;% frequencies =    0.0150    0.2610    0.4720
filename = 's_meter_log_2008_10_07-22_10_52_77_437432000.txt';% day 107.8
interval = 1607:2618;% frequencies =    0.0130    0.2570    0.4890
filename = 's_meter_log_2008_10_15-21_27_17_0_437432000.txt';% day 115.8
interval = 2563:2815;% frequencies =    0.0770    0.2320    0.4380
filename = 's_meter_log_2008_10_16-23_22_54_21_437432000.txt';% day 116.9
interval = 1334:2249;% frequencies =    0.0280    0.2410    0.4600
filename = 's_meter_log_2008_10_19-21_05_42_31_437432000.txt';% day 119.8 
interval = 2014:3019;% frequencies =    0.0190    0.1810    0.2330 (0.223 0.46)
filename = 's_meter_log_2008_10_20-10_20_57_0_437432000.txt';% day 120.4
interval = 2074:3051;% frequencies =    0.0200    0.2220    0.4630
filename = 's_meter_log_2008_10_23-22_19_34_66_437432000.txt';% day 123.9
interval = 1972:3046;% frequencies =    0.0250    0.2150    0.4540
filename = 's_meter_log_2008_10_27-20_57_37_0_437432000.txt';% day 127.8
interval = 2143:3162;% frequencies =    0.0120    0.1630    0.3800
filename = 's_meter_log_2008_10_31-11_08_54_55_437432000.txt';% day 131.5
interval = 1958:2943;%frequencies =    0.0180    0.2000    0.4180
filename = 's_meter_log_2008_11_01-20_54_15_0_437432000.txt';% day 132.9
interval = 2076:2655;% frequencies =    0.0210    0.2420    0.5050
filename = 's_meter_log_2008_11_02-21_12_49_76_437432000.txt';% day 133.8
interval = 2379:3420;% frequencies =    0.0120    0.2670    0.5520
filename = 's_meter_log_2008_11_04-10_46_55_82_437432000.txt';% day 135.5
interval = 2318:3416;% frequencies =    0.0180    0.2720    0.5390
filename = 's_meter_log_2008_11_05-20_32_27_50_437432000.txt';% day 136.9
interval = 1713:2801;% frequencies =    0.0170   0.272 0.5210 
filename = 's_meter_log_2008_11_07-21_09_25_80_437432000.txt';% day 138.8 
interval = 2154:3163;% frequencies =    0.0180   0.272  0.4880  
filename = 's_meter_log_2008_11_09-10_43_29_88_437432000.txt';% day 140.5
interval = 2581:3644;% frequencies =    0.0180    0.2640    0.5150
filename = 's_meter_log_2008_11_11-11_20_40_0_437432000.txt';% day 142.5
interval = 2050:3035;% frequencies =    0.0240    0.2560    0.5180
filename = 's_meter_log_2008_11_12-21_06_01_85_437432000.txt';% day 143.9
interval = 2008:3037;%frequencies =    0.0130    0.2840    0.5610
filename = 's_meter_log_2008_11_14-10_40_08_86_437432000.txt';% day 145.4
interval = 2463:3359;%frequencies =    0.0260    0.2780    0.5700
filename = 's_meter_log_2008_11_16-11_17_18_0_437432000.txt';% day 147.5
interval = 2170:3116;%frequencies =    0.0190    0.3010    0.6080
filename = 's_meter_log_2008_11_17-21_02_35_88_437432000.txt';% day 148.85
interval = 2151:2949;%frequencies =    0.0280    0.301   0.5910
filename = 's_meter_log_2008_11_18-10_18_06_57_437432000.txt';% day 149.45
interval = 2079:3097;% frequencies =    0.0180    0.3250    0.6930
filename = 's_meter_log_2008_11_20-10_55_16_72_437432000.txt'; % day 151.45
interval = 1866:2786;% frequencies =    0.0190    0.3260    0.6650
filename = 's_meter_log_2008_11_22-20_59_09_0_437432000.txt';% day 153.9
interval = 2269:3193;% frequencies =    0.0200    0.314    0.648
filename = 's_meter_log_2008_11_23-21_17_43_0_437432000.txt';% day 154.9
interval = 1944:2923;% frequencies =    0.0180    0.3040    0.6490
filename = 's_meter_log_2008_11_28-10_11_14_50_437432000.txt';% day 159.5
interval = 1307:1835;% frequencies =    0.0380    0.3040    0.6590
filename = 's_meter_log_2008_11_28-21_14_13_75_437432000.txt';% day 159.8
interval = 1420:2140;% frequencies =     0.0270    0.3260    0.6430
filename = 's_meter_log_2008_11_30-10_51_23_81_437432000.txt';% day 161.45
interval = 528:1521;% frequencies =    0.0200    0.3450    0.7100
filename = 's_meter_log_2008_12_01-11_06_57_59_437432000.txt';% day 162.5
interval = 1199:1875;%frequencies =    0.0200    0.1570    0.3520 %pfft this cant be right
filename = 's_meter_log_2008_12_03-21_10_46_80_437431000.txt';% day 164.9
interval = 2593:3594;%frequencies =    0.0130    0.3810    0.7500
filename = 's_meter_log_2008_12_06-20_30_22_48_437431000.txt';% day 167.9
interval = 1552:2546;%frequencies =    0.0200    0.3690    0.7510
filename = 's_meter_log_2008_12_22-20_38_22_56_437432000.txt';% day 183.8
interval = 1721:2714;%frequencies =    0.0190    0.3260    0.6850
filename = 's_meter_log_2008_12_31-12_22_28_0_437432000.txt';% day 192.5
interval = 2395:2914;%frequencies =    0.0380    0.3270    0.6540
filename = 's_meter_log_2009_01_02-20_49_44_69_437432000.txt';% day 194.8
interval = 898:1505;% frequencies =    0.0340    0.3090    0.6170
filename = 's_meter_log_2009_01_05-20_09_33_32_437432000.txt';% day 197.75
interval = 1200:1667;% frequencies =    0.0390    0.308    0.602
filename = 's_meter_log_2009_01_05-21_45_30_0_437432000.txt'; % day 197.85
interval = 1400:2000;% frequencies =    0.0290    0.5720    1.1440
filename = 's_meter_log_2009_01_09-10_20_18_0_437432000.txt';% day 201.45
interval = 1719:2260;% frequencies =    0.0580    0.6540    1.3070
filename = 's_meter_log_2009_01_13-21_01_05_82_437432000.txt';% day 205.9
interval = 1519:1894;% frequencies =    0.3310    0.6790    1.0280
filename = 's_meter_log_2009_01_19-11_49_30_29_437432000.txt';% day 211.5
interval = 1355:1877;% frequencies =    0.0390    0.3340    0.6810
filename = 's_meter_log_2009_01_22-11_08_45_58_437432000.txt';% day 214.5
interval = 1780:2761;%frequencies =    0.0130    0.1000    0.7270
filename = 's_meter_log_2009_02_02-11_20_02_47_437432000.txt';% day 225.5
interval = 1238:2220;% frequencies =    0.0210    0.7650    1.5650
filename = 's_meter_log_2009_02_03-21_05_13_87_437432000.txt';% day  226.85
interval = 1671:2609;%frequencies =    0.0140    0.2850    0.5770
filename = 's_meter_log_2009_02_05-12_15_45_0_437432000.txt';% day 228.5 
interval = 1378:2452;%frequencies =    0.2790    0.5710    1.1490
filename = 's_meter_log_2009_02_09-21_20_07_71_437432000.txt';% day 232.8
interval = 1686:1987;%frequencies =    0.0660    0.2440    0.5530
filename = 's_meter_log_2009_02_11-20_21_13_38_437432000.txt';% day 234.8
interval = 1552:2233;%frequencies =    0.2840    0.5780    0.8630
filename = 's_meter_log_2009_02_14-21_22_35_76_437432000.txt';%day 237.8
interval = 182:518;%frequencies =    0.0570    0.2290    0.5530
filename = 's_meter_log_2009_02_15-19_59_23_0_437432000.txt';% day 238.8
interval = 1621:1986;%frequencies =    0.0740    0.2950    0.5720
filename = 's_meter_log_2009_02_15-21_35_00_53_437432000.txt';% day 238.9
interval = 1786:2551;%frequencies =    0.0510    0.5490    1.1060
filename = 's_meter_log_2009_02_17-20_35_51_0_437432000.txt';% day 240.9
interval = 1917:2369;%frequencies =    0.0560    0.2820    0.5790
filename = 's_meter_log_2009_02_18-20_54_14_0_437432000.txt';% 241.9
interval = 2452:2828;% frequencies =    0.0350    0.1390    0.2950  Ahh kan det nu også passe? 
filename = 's_meter_log_2009_02_19-11_46_04_0_437432000.txt'; % day 242.5, BDOT turned off
interval = 1492:2019;%frequencies =    0.0380    0.2940    0.6010
filename = 's_meter_log_2009_02_20-12_04_38_0_437432000.txt';% day 243.5
interval = 1402:1900;%  frequencies =      0.0260    0.307    0.601
filename = 's_meter_log_2009_02_20-19_55_46_0_437432000.txt';% day 243.9
interval = 2653:3654;%frequencies =    0.0210    0.3160    0.6330
filename = 's_meter_log_2009_02_21-10_46_44_88_437432000.txt';% day 244.45
interval = 2304:3302;%frequencies =    0.0190    0.3090    0.6380
filename = 's_meter_log_2009_02_22-20_32_12_0_437432000.txt';% day 245.85
interval = 1973:2510;%frequencies =    0.0360    0.3200    0.6390
filename = 's_meter_log_2009_02_23-20_50_34_67_437432000.txt';% day 246.85
interval = 1303:2341;%frequencies =    0.0130    0.3290    0.6770
filename = 's_meter_log_2009_02_25-10_24_30_60_437432000.txt';% day 248.45 Bdot turned on after log was taken down
interval = 1622:2795;%frequencies =    0.0110    0.6340    0.6660
filename = 's_meter_log_2009_02_25-12_00_55_25_437433000.txt';% day 248.5
interval = 1369:2017;% frequencies =    0.0200    0.0500    0.6420
filename = 's_meter_log_2009_02_26-10_50_29_84_437432000.txt';%day 249.45
interval = 54:1071;%frequencies =    0.0200    0.3260    0.6580
filename = 's_meter_log_2009_02_27-11_01_30_69_437432000.txt';% day 250.45
interval = 1722:2449;%frequencies =    0.0260    0.2720    0.5780
filename = 's_meter_log_2009_02_27-12_38_04_0_437432000.txt'; %day 250.55
interval = 1650:2349;%frequencies =    0.0820    0.2820    0.5820
filename = 's_meter_log_2009_03_01-22_44_28_16_437432000.txt';% day 252.9
interval = 733:1516;%frequencies =    0.3310    0.6620    1.3160
filename = 's_meter_log_2009_03_02-21_23_48_67_437432000.txt'; % day 253.85
interval = 1101:2110;%frequencies =   0.0320    0.2910    0.5830
filename = 's_meter_log_2009_03_05-09_40_07_0_437432000.txt';% day 256.4 Bdot deactivated
interval = 1729:2715;% frequencies =    0.2950    0.5900    1.1730
filename = 's_meter_log_2009_03_07-11_53_24_28_437432000.txt';% day 258.5
interval = 1142:2104;%frequencies =    0.0560    0.3000    0.6060
filename = 's_meter_log_2009_03_08-21_38_39_0_437432000.txt';% day 259.85 Bdot active, Gain2:1 is inverted to -0.1
interval = 1274:2142;%frequencies =    0.0300    0.3    0.5860
filename = 's_meter_log_2009_03_08-23_17_35_8_437432000.txt';% day 259.95
interval = 1019:2002;%frequencies =    0.2940    0.5880    1.1760
filename = 's_meter_log_2009_03_09-10_57_19_79_437432000.txt';% day 260.45
interval = 126:1102;% frequencies =    0.0260    0.2820    0.5630
filename = 's_meter_log_2009_03_09-12_30_31_17_437432000.txt';% day 260.55
interval = 1709:2710;% frequencies =    0.0130    0.5730    1.1390
filename = 's_meter_log_2009_03_09-20_21_07_37_437432000.txt';%day 260.8
interval = 903:1808;%frequencies =    0.2770    0.5540    0.8310
filename = 's_meter_log_2009_03_09-21_57_22_0_437432000.txt';% day 260.9
interval = 1675:2169;% frequencies =    0.0350    0.567    1.0630
filename = 's_meter_log_2009_03_10-09_36_23_25_437432000.txt';%day 261.4
interval = 1619:2610;% frequencies =    0.0130    0.2700    0.5340
filename = 's_meter_log_2009_03_10-11_12_31_56_437432000.txt';% day 261.5
interval = 1860:2853;%frequencies =    0.0250    0.2660    0.5390
filename = 's_meter_log_2009_03_10-20_39_24_0_437432000.txt';% day 261.8
interval = 1701:2656;%frequencies =    0.0200    0.2600    0.5140
filename = 's_meter_log_2009_03_11-09_54_47_33_437432000.txt';% day 262.4 g1 = -9.5, g2 = -0.1 (grønne gains)
interval = 1480:1590;
filename = 's_meter_log_2009_03_11-11_31_04_40_437432000.txt';% day 262.5
interval = 1150:1388;%frequencies =    0.0540    0.2710    0.5410
filename = 's_meter_log_2009_03_11-13_07_41_10_437432000.txt';% day 262.6
interval = 1080:1150;% frequencies =    0.0450    0.1340    0.5040
filename = 's_meter_log_2009_03_11-20_57_47_76_437432000.txt';% day 262.85
interval = [1765:1824 1865:1945 1976:2032 2064:2138 2181:2236];%frequencies =   0.29 0.559
filename = 's_meter_log_2009_03_12-10_13_15_0_437432000.txt';% day 263.4
interval = [1626:1658 1687:1799 1809:1906 1918:2027 2076:2175 2191:2280 ...
            2292:2389 2408:2516];%frequencies =    0.336 0.65
filename = 's_meter_log_2009_03_12-11_49_37_30_437432000.txt'; %day 263.5 
interval = [1026:1105 1115:1271 1287:1355 1389:1464 1490:1555 1677:1744 ...
  1751:1828 1836:1898 1910:1991 2043:2106 2171:2236 2245:2307 2286:2532 2543:2507];% Freq = 0.24  0.6
filename = 's_meter_log_2009_03_12-19_41_09_21_437432000.txt';%day 263.85
interval = [1094:1162 1193:1273 1292:1360 1404:1483 1615:1685 1702:1782];% Freq = 0.317 0.678 , blue gains activated
filename = 's_meter_log_2009_03_13-10_31_43_69_437432000.txt'; % day 264.4
interval = [1664:1717 1855:1953 1975:2264 2277:2578 2598:2699];%frequencies =    0.0370    0.0990    0.6140
filename = 's_meter_log_2009_03_15-22_12_21_27_437432000.txt';% day 266.9
interval = [1669:1823 1846:2013 2051:2092 2979:3956];% freq = 0.28 0.559
filename = 's_meter_log_2009_03_16-11_27_21_43_437432000.txt';% day 267.5
interval = [1375:1545 1575:1758 1853:2031];% freq = 0.281 0.575
filename = 's_meter_log_2009_03_19-20_13_45_33_437432000.txt';%day 270.8
interval = 1685:1990;%frequencies =    0.2650    0.5300    0.8170
filename = 's_meter_log_2009_03_22-19_33_53_19_437432000.txt';% day 273.8
interval = 1704:2102; %frequencies =    0.2520    0.5050 
filename = 's_meter_log_2009_03_28-11_56_53_27_437432000.txt';% day 279.5
interval = 1435:2081;%frequencies =    0.0300    0.2510    0.5220
filename = 's_meter_log_2009_03_29-11_38_59_78_437432000.txt';% day 280.5
interval = 2263:3250;%frequencies =    0.2790    0.5510    0.5830
filename = 's_meter_log_2009_03_30-11_57_30_75_437432000.txt';% day 281.5
interval = 2689:3438;%frequencies =    0.0260    0.2620    0.5410
filename = 's_meter_log_2009_04_06-12_30_45_41_437432000.txt'; % day 288.5
interval = 1283:2337;%frequencies =    0.0190    0.2390    0.4910
filename = 's_meter_log_2009_04_06-21_57_21_74_437432000.txt';% day 288.9
interval = 2494:3508;%frequencies =    0.0200    0.5030    0.5230
filename = 's_meter_log_2009_04_08-21_00_50_27_437432000.txt';% day 290.8
interval = 790:1792;%frequencies =    0.2670    0.5340    0.8000
filename = 's_meter_log_2009_04_12-22_11_58_84_437432000.txt';% day 294.8
interval = 2269:3295;%frequencies =    0.3550    0.7040    0.7240
filename = 's_meter_log_2009_04_14-21_13_11_33_437432000.txt';% day 296.8 Gains set back to 0.1 for all g2
interval = 3402:4199;%frequencies =    0.0180    0.3860    0.7630
filename = 's_meter_log_2009_04_15-21_31_24_45_437432000.txt';% day 297.8
interval = 2446:3061;%frequencies =    0.3630    0.7150 
filename = 's_meter_log_2009_04_26-21_42_09_56_437432000.txt';% day 308.8
interval = 1941:2163;%frequencies =    0.0870    0.4330    0.8670
filename = 's_meter_log_2009_06_23-12_27_16_44_437432000.txt';% day 369.5
interval = 2250:3015;% frequencies =    0.0250    1.0000    2.0490
filename = 's_meter_log_2009_06_24-11_09_26_0_437432000.txt';% day 370.5
interval = 1825:2530;%frequencies =    0.0370    1.5220    2.1130
filename = 's_meter_log_2009_06_29-11_05_23_0_437432000.txt';% day 375.5
interval = 1881:2677;%frequencies =    0.0570    1.4890    2.1680
filename = 's_meter_log_2009_09_24-12_08_12_66_437432000.txt';% day 462.5
interval = 2574:3584;% frequencies =    0.0190    1.3060    1.4020
filename = 's_meter_log_2009_11_20-12_33_46_18_437432000.txt';% day 
interval = 1457:2407;%frequencies =    0.0350    1.9780    2.3820

day = 55 + [1.5 13.9 25.5 34.5 37.4 42.8 42.9 43.5 44.5 50.5 56.5 59.5 60.5 63.5 63.9 64.9...
  66.5 67.4 67.5 67.9 69.9 70.5 70.9 72.9 74.5 74.9 76.5 78.5 78.9 79.9 81.9 82.9 85.9 86.9 87.9...
  90.4 90.9 91.4 91.5 91.9 92.5 92.9 93.45 93.55 93.7 93.85 94.8 94.9 95.4 95.5 95.8 95.9 101.5...
  101.8 101.9 102.75 103.45 103.75 104.5 104.8 105.8 105.9 106.8 107.8 115.8 116.9 119.8 120.4 ...
  123.9 127.8 131.5 132.9 133.8 135.5 136.9 138.8 140.5 142.5 143.9 145.4 147.5 148.85 149.45...
  151.45 153.9 154.9 159.5 159.8 161.45 162.5 164.9 167.9 183.8 192.5 194.8 197.75 197.85 201.45...
  205.9 211.5 214.5 225.5 226.85 228.5 232.8 234.8 237.8 238.8 238.9 240.9 241.9 242.5 243.5 243.9...
  244.45 245.85 246.85 248.45 248.5 249.45 250.45 250.55 252.9 253.85 256.4 258.5 259.85 259.95...
  260.45 260.55 260.8 260.9 261.4 261.5 261.8 262.5 262.6 262.85 263.4 263.5 263.85 264.4 266.9 ...
  267.5 270.8 273.8 279.5 280.5 281.5 288.5 288.9 290.8 294.8 296.8 297.8 308.8 369.5 370.5 375.5... 
  462.5];
freq1 = [0.21 0.387 0.51455 0.592 0.592 0.646 0.646 0.650 0.6697 0.72 0.7185 0.715 0.7 0.7 0.6815...
  0.65 0.61 0.48 0.46 0.431 0.292 0.28 0.264 0.258 0.233 0.23 0.22 0.25 0.25 0.24 0.17 0.17 0.27...
  0.255 0.245 0.22 0.212 0.183 0.172 0.166 0.167 0.174 0.181 0.175 0.163 0.177 0.149 0.147 0.152...
  0.145 0.15 0.143 0.169 0.170 0.18 0.195 0.215 0.212 0.234 0.236 0.26 0.266 0.261 0.257 0.232...
  0.241 0.223 0.222 0.215 0.163 0.2 0.242 0.267 0.272 0.272 0.272 0.264 0.256 0.284 0.278 0.301...
  0.301 0.325 0.326 0.314 0.304 0.304 0.326 0.345 0.352 0.381 0.369 0.326 0.327 0.309 0.308 0.292...
  0.292 0.331 0.334 0.334 0.334 0.285 0.279 0.244 0.284 0.229 0.295 0.255 0.282 0.292 0.294 0.307...
  0.316 0.309 0.32 0.329 0.329 0.329 0.326 0.272 0.282 0.331 0.291 0.295 0.3 0.3 0.295 0.282 0.282...
  0.277 0.271 0.270 0.266 0.26 0.271 0.271 0.29 0.336 0.24 0.317 0.35 0.28 0.281 0.265 0.252 0.251...
  0.279 0.262 0.239 0.239 0.267 0.355 0.386 0.363 0.433 1 1.552 1.489 1.402];
freq2 = [0.42 0.7742 0.9933 1.173 1.195 1.292 1.304 1.315 1.3394 1.4 1.44 1.395 1.395 1.395...
  1.395 1.35 1.225 0.959 0.92 0.855 0.584 0.566 0.57 0.53 0.474 0.462 0.464 0.58 0.512 0.489...
  0.38 0.4 0.54 0.53 0.478 0.454 0.417 0.373 0.357 0.332 0.342 0.381 0.342 0.356 0.325 0.318 0.297...
  0.307 0.298 0.297 0.3 0.292 0.331 0.343 0.347 0.397 0.43 0.456 0.454 0.472 0.549 0.553 0.472 0.489...
  0.438 0.46 0.46 0.463 0.454 0.38 0.418 0.505 0.552 0.539 0.521 0.488 0.515 0.518 0.561 0.570 0.608...
  0.591 0.693 0.665 0.648 0.649 0.659 0.643 0.71 0.71 0.75 0.751 0.685 0.654 0.617 0.602 0.595 0.654...
  0.679 0.681 0.727 0.765 0.577 0.571 0.553 0.578 0.553 0.572 0.549 0.579 0.579 0.601 0.601 0.633...
  0.638 0.639 0.677 0.666 0.642 0.658 0.578 0.582 0.662 0.583 0.590 0.606 0.586 0.588 0.563 0.573...
  0.554 0.567 0.534 0.539 0.514 0.541 0.568 0.559 0.65 0.6 0.678 0.7 0.559 0.575 0.53 0.505 0.522 ...
  0.551 0.541 0.491 0.513 0.534 0.724 0.763 0.715 0.867 2.049 2.113 2.168 2.83];
%%
d = importdata(filename);

for(k=1:size(d.colheaders,2))
  evalstr = [d.colheaders{k} ' = ' 'd.data(:,k)'];
  evalc(evalstr);
end

rxok = rxok - rxok(1);
rxtotal = rxtotal - rxtotal(1);

%
dfreq1 = diff(freq1)./diff(day);
dfreq2 = diff(freq2)./diff(day);
dfreq1(end+1) = dfreq1(end);
dfreq2(end+1) = dfreq2(end);

figure1 = figure(1)
clf
axes1 = axes(...
  'XColor',[0.502 0.502 0.502],...
  'XGrid','on',...
  'YColor',[0.502 0.502 0.502],...
  'YGrid','on',...
  'ZColor',[0.502 0.502 0.502],...
  'Parent',figure1);
whitebg('w')
title('test')
subplot(311)
plot(smeter,'b')
hold on
plot(rxtotal,'k')
plot(rxok,'r')
plot(abs(abs(elevation-90)-90),'m')
ylim([0 200])
grid on
handle = legend('S-meter','rx-total','rx-ok','elevation','Location','BestOutSide');
set(handle, 'TextColor', 'k')
xlabel('Sample Number [-]')
title(strrep(filename,'_','.'))

subplot(312)
y = LimitSignal(smeter(interval), [0 200]);
rep = 1;
interval1 = interval(1):interval(1)+length(y)*rep-1;
y1 = repmat(y,rep,1);
%plot(time(interval1)-time(interval1(1)),y1,'b')
plot((time(interval)-time(interval(1))),y,'b')
% plot(y,'b')
hold on
plot((time(interval)-time(interval(1))),rxtotal(interval),'k')
plot((time(interval)-time(interval(1))),rxok(interval),'r')
plot(time(interval)-time(interval(1)),abs(abs(elevation(interval)-90)-90),'m')
grid on
handle = legend('S-meter','rx-total','rx-ok','elevation','Location','BestOutSide');
set(handle, 'TextColor', 'k')
xlabel('Time [s]')
ylim([0 200])

subplot(313)
[frq,data] = kfft(1/mean(diff(time(interval))), LimitSignal(smeter(interval),[0 200]));
[b,a] = butter(1,0.2);
data_filt = filtfilt(b,a,data);
plot(frq,data,'b');
hold on
[maxtab mintab] = peakdet(data_filt, 0.00001);
[data_sorted indices] = sort(maxtab(:,2),1,'descend');
plot(frq,data_filt,'r');
plot(frq(maxtab(indices(1:3),1)),data_sorted(1:3),'ks')
plot(frq,data,'b');
plot(frq,data_filt,'r');
frequencies = sort(round(frq(maxtab(indices(1:3),1))*1000)/1000,'ascend');
for(i = 1:3)
  text(frq(maxtab(indices(i))),data_sorted(i)+5,num2str(round(frq(maxtab(indices(i)))*1000)/1000),'FontSize',10)
end
handle = legend('FFT','FFT-filtered','Sign. Freq.','Location','NorthEast');
set(handle, 'TextColor', 'k')
xlabel('Frequency [Hz]')
grid on
title('FFT');
xlabel('Frequency [Hz]');
xlim([0 max(frequencies)*4])
ylim([0 15])
print('-dpng',['tumbling_' filename(1:end-4) '.png'])

figure2 = figure(2)
clf
axes1 = axes(...
  'XColor',[0.502 0.502 0.502],...
  'XGrid','on',...
  'YColor',[0.502 0.502 0.502],...
  'YGrid','on',...
  'ZColor',[0.502 0.502 0.502],...
  'Parent',figure2);
whitebg('w')
plot(day,freq1,'b-x')
hold on
stairs(day,dfreq1,'m')
plot(day,freq2,'r-x')
stairs(day,dfreq2,'k')
grid on
handle = legend('freq1','dfreq1/dt','freq2','dfreq2/dt', 'Location', 'NorthWest');
set(handle, 'TextColor', 'k')
xlabel('Day')
ylabel('Freqency [Hz]')
%xlim([0 140])
ylim([-0.5 3])
print('-dpng','tumbling_freq_over_time.png')


figure2 = figure(3)
clf
axes1 = axes(...
  'XColor',[0.502 0.502 0.502],...
  'XGrid','on',...
  'YColor',[0.502 0.502 0.502],...
  'YGrid','on',...
  'ZColor',[0.502 0.502 0.502],...
  'Parent',figure2);
whitebg('w')
len = 30; 
plot(day(end-len:end),freq1(end-len:end),'b-x')
hold on
stairs(day(end-len:end),dfreq1(end-len:end),'m')
plot(day(end-len:end),freq2(end-len:end),'r-x')
stairs(day(end-len:end),dfreq2(end-len:end),'k')
grid on
handle = legend('freq1','dfreq1/dt','freq2','dfreq2/dt', 'Location', 'NorthWest');
set(handle, 'TextColor', 'k')
xlabel('Day')
ylabel('Freqency [Hz]')
%xlim([0 140])
ylim([-0.5 3])
print('-dpng','tumbling_freq_over_time_latest.png')