/// <summary>
        /// Get Y coordinate array
        /// </summary>
        /// <returns>Y array</returns>
        private double[] GetYArray()
        {
            double[] Y = new double[NY];

            switch (GridType)
            {
            case 4:        //Gaussian Lat/Lon grid
                Y = (double[])DataMath.Gauss2Lats(NY)[0];
                break;
            //default:
            //    double dy = Math.Abs(DY);
            //    //double sLat = Math.Min(Lat1, Lat2);
            //    double sLat = Lat1;

            //    for (int i = 0; i < NY; i++)
            //    {
            //        Y[i] = sLat + dy * i;
            //    }
            //    break;
            default:
                for (int i = 0; i < NY; i++)
                {
                    Y[i] = Lat1 + DY * i;
                }
                if (DY < 0)
                {
                    Array.Reverse(Y);
                }
                break;
            }

            return(Y);
        }
        /// <summary>
        /// Get Y coordinate array of Gaussian grid
        /// </summary>
        /// <returns>Y coordinate array</returns>
        public double[] GetGaussYArray()
        {
            double[] Y = new double[NY];
            Y = (double[])DataMath.Gauss2Lats(NY)[0];

            //double ymin = Y[0];
            //double ymax = Y[Y.Length - 1];
            //double delta = (ymax - ymin) / (NY - 1);
            //for (int i = 0; i < NY; i++)
            //    Y[i] = ymin + i * delta;

            return(Y);
        }