Make FoV routine

.gitignore

@@ -2,6 +2,7 @@
 *.so
 *.o
 *~
+#*
 .DS_Store
 
 bin
@@ -14,6 +15,7 @@ bin
 /lib
 /log
 
+codebase/superdarn/src.lib/tk/idl.tar.gz
 codebase/general/src.lib/stdkey.1.5/src/make_palette
 
 *errstr.h

codebase/analysis/src.lib/linearfit.1.0/doc/linearfit.doc.xml

@@ -0,0 +1,20 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<library>
+<project>analysis</project>
+<name>linearfit</name>
+<location>src.lib/linearfit</location>
+
+
+<function>
+<name>linear_regression</name>
+<location>src.lib/linearfit</location>
+<header>analysis/linearfit.h</header>
+<syntax>int linear_regression(float x[], float y[], float sig[], int num, int mwt, float *intercept, float *sig_intercept, float *slope, float *sig_slope, float *chi2, float *q);</syntax>
+<description><p>The <fn href="linear_regression.html">linear_regression</fn> function calculates linear regression for an array of independent and dependent variables, with optional errors.</p>
+</description>
+<returns><p>Returns the fitted slope and intercept (<ar>slope</ar>, <ar>intercept</ar>), their errors (<ar>sig_slope</ar>, <ar>sig_intercept</ar>), as well as the chi-squared and gamma results for the fit (<ar>chi2</ar>, <ar>q</ar>).</p>
+</returns>
+<error>If there are problems assigning data to the subarrays when performing the sorting, the routine will exit with 0.</error>
+</function>
+
+</library>

codebase/analysis/src.lib/linearfit.1.0/include/linearfit.h

@@ -0,0 +1,39 @@
+/* linearfit.h
+   ===========
+   Author: Angeline G. Burrell - NRL - 2021
+   This is a U.S. government work and not under copyright protection in the U.S.
+
+   This file is part of the Radar Software Toolkit (RST).
+
+   Disclaimer: RST is licensed under GPL v3.0. Please visit 
+               <https://www.gnu.org/licenses/> to see the full license
+
+   Modifications:
+
+*/
+
+#ifndef _LINEARFIT_H
+#define _LINEARFIT_H
+
+#ifndef _RTYPES_H
+#include "rtypes.h"
+#endif
+
+#ifndef _STATS_H
+#include "stats.h"
+#endif
+
+static float sqrarg;
+#ifndef SQR
+#define SQR(a) ((sqrarg = (a)) == 0.0 ? 0.0 : sqrarg * sqrarg);
+#endif
+
+#ifndef MINSIG
+#define MINSIG 1.0e-5
+#endif
+
+int linear_regression(float x[], float y[], float sig[], int num, int mwt,
+		      float *intercept, float *sig_intercept, float *slope,
+		      float *sig_slope, float *chi2, float *q);
+
+#endif

codebase/analysis/src.lib/linearfit.1.0/src/linearfit.c

@@ -0,0 +1,110 @@
+/* linearfit.c
+   ===========
+   Author: Angeline G. Burrell - NRL - 2021
+   This is a U.S. government work and not under copyright protection in the U.S.
+
+   This file is part of the Radar Software Toolkit (RST).
+
+   Disclaimer: RST is licensed under GPL v3.0. Please visit 
+               <https://www.gnu.org/licenses/> to see the full license
+
+   Modifications:
+
+*/
+
+#include <stdio.h>
+#include <math.h>
+
+#include "linearfit.h"
+
+/**
+ * @brief perform a linear regression
+ *
+ * @params[in] x   - independent variable array
+ *             y   - dependent variable array
+ *             sig - error of the dependent variable array
+ *             num - number of variables in x, y, and sig
+ *             mwt - boolean, if true sig is used to calculate weights
+ *
+ * @params[out] intercept     - pointer to the fitted intercept
+ *              sig_intercept - pointer to the error of the intercept
+ *              slope         - pointer to the fitted slope
+ *              sig_slope     - pointer to the error of the slope
+ *              chi2          - pointer to the Chi-squared value
+ *              q             - pointer to the gamma function output
+ *
+ * @references Numerical Recipes in C, 2nd Edition
+ **/
+
+int linear_regression(float x[], float y[], float sig[], int num, int mwt,
+		      float *intercept, float *sig_intercept, float *slope,
+		      float *sig_slope, float *chi2, float *q)
+{
+  int n, status = 1;
+
+  float wt, t, sxoss, sx, sy, st2, ss, sigdat;
+
+  sx = sy = st2 = 0.0;
+  *slope  = 0.0;
+
+  if(mwt)
+    {
+      for(ss = 0.0, n = 0; n < num; n++)
+	{
+	  if(sig[n] == 0.0) sig[n] = MINSIG;
+
+	  wt  = 1.0 / (sig[n] * sig[n]);
+	  ss += wt;
+	  sx += x[n] * wt;
+	  sy += y[n] * wt;
+	}
+    }
+  else
+    {
+      for(ss = num, n = 0; n < num; n++)
+	{
+	  sx += x[n];
+	  sy += y[n];
+	}
+    }
+
+  sxoss = sx / ss;
+
+  for(n = 0; n < num; n++)
+    {
+      t       = (mwt) ? (x[n] - sxoss) / sig[n] : x[n] - sxoss;
+      st2    += t * t;
+      *slope += t * y[n];
+    }
+
+  *slope        /= st2;
+  *intercept     = (sy - sx * (*slope)) / ss;
+  *sig_intercept = sqrt((1.0 + sx * sx / (ss * st2)) / ss);
+  *sig_slope     = sqrt(1.0 / st2);
+
+  /* Calculate Chi2 */
+
+  *chi2 = 0.0;
+  *q    = 1.0;
+
+  if(mwt == 0)
+    {
+      for(n = 0; n < num; n++)
+	*chi2 += SQR(y[n] - (*intercept) - (*slope) * x[n]);
+
+      sigdat          = sqrt((*chi2) / (float)(num - 2));
+      *sig_slope     *= sigdat;
+      *sig_intercept *= sigdat;
+      status = 0;
+    }
+  else
+    {
+      for(n = 0; n < num; n++)
+	*chi2 += SQR((y[n] - (*intercept) - (*slope) * x[n]) / sig[n]);
+
+      if(num > 2) *q = (float)gammaq(0.5 * (double)(num - 2),
+				     0.5 * (double)(*chi2));
+    }
+
+  return(status);
+}

codebase/analysis/src.lib/linearfit.1.0/src/makefile

@@ -0,0 +1,24 @@
+#  Makefile for the linearfit library
+#  ==================================
+# Author: Angeline G. Burrell, NRL (2021)
+# This is a U.S. government work and not under copyright protection in the U.S.
+#
+# Disclaimer:  RST is licensed under GPL v3.0. Please visit 
+#              <https://www.gnu.org/licenses/> to see the full license
+#
+# Modifications:
+
+include $(MAKECFG).$(SYSTEM)
+
+INCLUDE	= -I$(IPATH)/base -I$(IPATH)/general -I$(IPATH)/analysis -I../include
+INC	= $(IPATH)/analysis
+OUTPUT	= linearfit
+DSTPATH	= $(LIBPATH)
+
+SRC	= linearfit.c
+OBJS	= $(SRC:.c=.o)
+
+LINK	= "1"
+
+include $(MAKELIB).$(SYSTEM)
+