Made improvements to spec_avg so it prints out about 5/16 Hz around 850 Hz.

runspecS6H1feedforward

 #/bin/csh
 
-/home/keithr/master/opt/lscsoft/lalapps/bin/lalapps_spec_avg --SFTs /archive/frames/S6/pulsar/sfts/feedforward/hann/H-1_H1_1800SFT_AMPSC02L2-931/\*.sft\;/archive/frames/S6/pulsar/sfts/feedforward/hann/H-1_H1_1800SFT_AMPSC02L2-932/\*.sft  --startGPS 931000000 --endGPS 932810000 --fMin 40. --fMax 1999.999 --freqRes 0.10 --timeBaseline 1800 -v 3 --IFO H1 --outputBname S6H1_Hann_40_2000feedforward
+/home/gmeadors/master/opt/lscsoft/lalapps/bin/lalapps_spec_avg --SFTs /archive/frames/S6/pulsar/sfts/feedforward/hann/H-1_H1_1800SFT_AMPSC02L2-931/\*.sft\;/archive/frames/S6/pulsar/sfts/feedforward/hann/H-1_H1_1800SFT_AMPSC02L2-932/\*.sft  --startGPS 931000000 --endGPS 932810000 --fMin 40. --fMax 1999.999 --freqRes 0.10 --timeBaseline 1800 -v 3 --IFO H1 --outputBname S6H1_Hann_40_2000feedforward

runspecS6H1test

 #/bin/csh
 
-/home/keithr/master/opt/lscsoft/lalapps/bin/lalapps_spec_avg --SFTs /archive/frames/S6/pulsar/sfts/feedforward/beforeHann/H-1_H1_1800SFT_AMPSC02L2-931/\*.sft\;/archive/frames/S6/pulsar/sfts/feedforward/beforeHann/H-1_H1_1800SFT_AMPSC02L2-932/\*.sft  --startGPS 931000000 --endGPS 932810000 --fMin 40. --fMax 1999.999 --freqRes 0.10 --timeBaseline 1800 -v 3 --IFO H1 --outputBname S6H1_Hann_40_2000test
+/home/gmeadors/master/opt/lscsoft/lalapps/bin/lalapps_spec_avg --SFTs /archive/frames/S6/pulsar/sfts/feedforward/beforeHann/H-1_H1_1800SFT_AMPSC02L2-931/\*.sft\;/archive/frames/S6/pulsar/sfts/feedforward/beforeHann/H-1_H1_1800SFT_AMPSC02L2-932/\*.sft  --startGPS 931000000 --endGPS 932810000 --fMin 40. --fMax 1999.999 --freqRes 0.10 --timeBaseline 1800 -v 3 --IFO H1 --outputBname S6H1_Hann_40_2000test
 

spec_avg.c

@@ -253,19 +253,22 @@ int main(int argc, char **argv)
 		}
 	      }
 	    thisavepower /= count;
-            if (i >= 1457995 && i <= 1458005)
+            // Print out the calibrated ASD bin value over about 5/16 Hz,
+            // centered on 850 Hz.
+            if (i >= 1457719 && i <= 1458281)
               {
                 f = f0 + ((REAL4)i)*deltaF;
-                printf("SFT number=%d\n",j);
-                printf("Frequency=%f\n",f);
-                printf("Bin=%d\n",i);
-                printf("Real component=%f\n", sft_vect->data[0].data->data[i].re);
-                printf("Imaginary component=%f\n", sft_vect->data[0].data->data[i].im);
-                printf("Bin power=%f\n", thispower);
-                printf("Bin-neighbor average power=%f\n", thisavepower);
-                printf("Scale-factor=%f\n", scalefactor);
-                printf("Time baseline=%f\n", timebaseline);
-                fprintf(fp5,"%d %16.8f %g %g %g %g\n",j, f, thispower, thisavepower, scalefactor, timebaseline);
+                //printf("SFT number=%d\n",j);
+                //printf("Frequency=%f\n",f);
+                //printf("Bin=%d\n",i);
+                //printf("Real component=%f\n", sft_vect->data[0].data->data[i].re);
+                //printf("Imaginary component=%f\n", sft_vect->data[0].data->data[i].im);
+                //printf("Bin power=%f\n", thispower);
+                //printf("Bin-neighbor average power=%f\n", thisavepower);
+                //printf("Scale-factor=%f\n", scalefactor);
+                //printf("Time baseline=%f\n", timebaseline);
+                printf("%f\n", pow(2*thispower/scalefactor/scalefactor/timebaseline, 0.5));
+                fprintf(fp5,"%d %16.8f %g %g %g %g\n",j, f, thispower, pow(2*thispower/scalefactor/scalefactor/timebaseline,0.5), scalefactor, timebaseline);
               }
 	    weight = 1./thisavepower;
 	    //	    weight = 1.;