Minor fixes

Author: thwbecker

Makefile.include

@@ -6,10 +6,11 @@
 #
 #
 # quad precision
-ADD_FLAGS = -DHC_PRECISION=32 -O2 
+#ADD_FLAGS = -DHC_PRECISION=32 -O2 
 #
 # double precision
-#ADD_FLAGS = -O2 
+ADD_FLAGS = -O2
+
 
 GGRD_INC_FLAGS = -I$(GMTHOME)/include -I$(NETCDFHOME)/include 
 GGRD_LIBS_LINKLINE = -lggrd -lgmt -lpsl -lnetcdf 

README.TXT

@@ -165,7 +165,9 @@ KNOWN LIMITATIONS
 The propagator matrix approach can be unstable for moderately high
 maximum degrees. This behavior can be addressed partially by compiling
 HC in quadruple precision, and kinematic boundary can be addressed
-with a kludge, see -cbckl. 
+with a kludge, see -cbckl. If you compile with this trick and
+quadruple precision (the default within SEATREE), you should be good
+to up to L=127.
 
 
 SPHERICAL HARMONICS FORMAT

ggrd_grdtrack_util.c

@@ -500,7 +500,7 @@ int ggrd_grdtrack_init(double *west, double *east,
   FILE *din;
   float dz1,dz2;
   struct GRD_HEADER ogrd;
-  int i,one_or_zero,nx,ny,mx,my,nn;
+  int i,one_or_zero,nx,ny,mx,my;
   char filename[BUFSIZ*2],*cdummy;
   static int gmt_init = FALSE;
   /* 
@@ -716,7 +716,7 @@ int ggrd_grdtrack_init(double *west, double *east,
   nx = irint ( (*east - *west) / (*grd)[0].x_inc) + one_or_zero;
   ny = irint ( (*north - *south) / (*grd)[0].y_inc) + one_or_zero;
   /* real size of data */
-  nn = nx * ny;
+  //nn = nx * ny;
 
   /* padded */
   mx = nx + 4;

make_tar

@@ -2,7 +2,7 @@
 #
 # create a tar file 
 #
-ver=${1-1.0.3}
+ver=${1-1.0.4}
 date=`date '+%m%d%y'`
 
 tf=$HOME/tmp/hc-$ver.$date.tgz