core
/
eccodes
mirror of https://github.com/ecmwf/eccodes.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Improve comments

This commit is contained in:
Shahram Najm 2014-10-28 13:13:57 +00:00
parent b0ad650d4a
commit 1d13e3eb2b
2 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
src/grib_api.h
View File

@ -1305,6 +1305,7 @@ grib_handle* grib_handle_new_from_partial_message(grib_context* c,void* data, si
int grib_is_missing(grib_handle* h, const char* key, int* err);
int grib_is_defined(grib_handle* h, const char* key);
int grib_set_missing(grib_handle* h, const char* key);
/* The truncation is the Gaussian number (or order) */
int grib_get_gaussian_latitudes(long truncation,double* latitudes);
int grib_julian_to_datetime(double jd, long *year, long *month, long *day, long *hour, long *minute, long *second);

4
src/grib_geography.c
View File

@ -60,6 +60,8 @@ static void gauss_first_guess(long trunc, double* vals)
}
}
/* 'trunc' is the Gaussian number (or order) */
/* i.e. Number of parallels between a pole and the equator */
int grib_get_gaussian_latitudes(long trunc, double *lats)
{
long jlat, iter, legi;
@ -102,7 +104,7 @@ int grib_get_gaussian_latitudes(long trunc, double *lats)
conv = legfonc / ((((double)nlat) * (mem2 - root * legfonc) ) / (1.0 - (root *root)));
root -= conv;
/* Routine fails if no convergence after JPMAXITER iterations */
/* Routine fails if no convergence after MAXITER iterations */
if( iter++ > MAXITER )
{
return GRIB_GEOCALCULUS_PROBLEM;