Changed calculation of receive power from squared I+Q
authormleech <mleech@221aa14e-8319-0410-a670-987f0aec2ac5>
Tue, 5 Sep 2006 04:13:44 +0000 (04:13 +0000)
committermleech <mleech@221aa14e-8319-0410-a670-987f0aec2ac5>
Tue, 5 Sep 2006 04:13:44 +0000 (04:13 +0000)
git-svn-id: http://gnuradio.org/svn/gnuradio/trunk@3490 221aa14e-8319-0410-a670-987f0aec2ac5

gr-radio-astronomy/src/python/local_calibrator.py

index 6c5fb4f1f809456ac98c0033039a947bcb487839..578d07e9adeabb9c94ff14808fb647fe889efa08 100755 (executable)
@@ -65,10 +65,11 @@ def calib_numogate_ridge_observatory_total_power(data):
 
     rainbow_file = open (filenamestr+".tpdat","a")
   
-    r = (math.sqrt(data) / 2048) * 1000.0
+    r = (data / 4096.0)
+    flt = "%6.3f" % r
     #r = calib_default_total_power(data)
     inter = globals()["calib_decln"]
-    rainbow_file.write(str(ephem.hours(sidtime))+" "+str(r)+" "+str(inter)+"\n")
+    rainbow_file.write(str(ephem.hours(sidtime))+" "+flt+" "+str(inter)+"\n")
     rainbow_file.close()
     return(r)
 
@@ -101,7 +102,7 @@ def calib_numogate_ridge_observatory_fft(data,l):
     if not "calib_then" in globals():
         globals()["calib_then"] = now
 
-    delta = 5
+    delta = (l/1024)*5
                
     if (now - globals()["calib_then"]) >= delta: