Commit 9f6877f6 authored by Uwe Schulzweida's avatar Uwe Schulzweida
Browse files

Renamed gaussaw() to gaussianLatitudes().

parent e374d488
......@@ -1282,7 +1282,7 @@ int vlistInqVarSubtype(int vlistID, int varID);
void gribapiLibraryVersion(int *major_version, int *minor_version, int *revision_version);
void gaussaw(double *pa, double *pw, size_t nlat);
void gaussianLatitudes(double *latitudes, double *weights, size_t nlat);
#ifdef __cplusplus
}
......
......@@ -176,10 +176,10 @@ void gauaw(size_t kn, double *restrict pl, double *restrict pw)
}
void gaussaw(double *pa, double *pw, size_t nlat)
void gaussianLatitudes(double *latitudes, double *weights, size_t nlat)
{
//gauaw_old(pa, pw, nlat);
gauaw(nlat, pa, pw);
//gauaw_old(latitudes, weights, nlat);
gauaw(nlat, latitudes, weights);
}
......@@ -187,12 +187,12 @@ bool isGaussGrid(size_t ysize, double yinc, const double *yvals)
{
bool lgauss = false;
if ( IS_EQUAL(yinc, 0) && ysize > 2 ) /* check if gaussian */
if ( IS_EQUAL(yinc, 0) && ysize > 2 ) // check if gaussian
{
size_t i;
double *yv = (double *) Malloc(ysize*sizeof(double));
double *yw = (double *) Malloc(ysize*sizeof(double));
gaussaw(yv, yw, ysize);
gaussianLatitudes(yv, yw, ysize);
Free(yw);
for ( i = 0; i < ysize; i++ )
yv[i] = asin(yv[i])/M_PI*180.0;
......@@ -203,7 +203,7 @@ bool isGaussGrid(size_t ysize, double yinc, const double *yvals)
if ( i == ysize ) lgauss = true;
/* check S->N */
// check S->N
if ( lgauss == false )
{
for ( i = 0; i < ysize; i++ )
......
......@@ -381,7 +381,7 @@ static
void calc_gaussgrid(double *restrict yvals, size_t ysize, double yfirst, double ylast)
{
double *restrict yw = (double *) Malloc(ysize * sizeof(double));
gaussaw(yvals, yw, ysize);
gaussianLatitudes(yvals, yw, ysize);
Free(yw);
for (size_t i = 0; i < ysize; i++ )
yvals[i] = asin(yvals[i])/M_PI*180.0;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment