change cos_lat limits and use clamp

Author: aaronjridley

src/grid_sphere.cpp

@@ -154,7 +154,7 @@ void Grid::create_sphere_grid(Quadtree quadtree) {
   }
 
   arma_cube cos_lat = cos(geoLat_scgc);
-  cos_lat.elem( find(cos_lat < 0.00001) ).fill(0.00001);
+  cos_lat.clamp(0.0001, 1.0);
 
   y_Center = geoLat_scgc;
   x_Center = geoLon_scgc % cos_lat;
@@ -180,7 +180,7 @@ void Grid::create_sphere_grid(Quadtree quadtree) {
   }
 
   arma_cube cos_lat_L = cos(geoLat_Left);
-  cos_lat_L.elem( find(cos_lat_L < 0.00001) ).fill(0.00001);
+  cos_lat_L.clamp(0.0001, 1.0);
 
   x_Left = geoLon_Left % cos_lat_L;
   dy_Left.set_size(nLons, nLats, nAlts);
@@ -206,7 +206,7 @@ void Grid::create_sphere_grid(Quadtree quadtree) {
   }
 
   arma_cube cos_lat_D = cos(geoLat_Down);
-  cos_lat_D.elem( find(cos_lat_D < 0.00001) ).fill(0.00001);
+  cos_lat_D.clamp(0.0001, 1.0);
 
   y_Down = geoLat_Down;
   dx_Down = dlon * cos_lat_D;