zaxisCompare: set epsilon from 0 to 1e-9

b3e30be5 · Uwe Schulzweida · 617208b1 · b3e30be5 · b3e30be5
Commit b3e30be5 authored Oct 11, 2011 by Uwe Schulzweida
Hide whitespace changes
Inline Side-by-side

Showing with 8 additions and 4 deletions

ChangeLog ChangeLog +4 -0

src/varscan.c src/varscan.c +4 -4

No files found.
--- a/ChangeLog
+++ b/ChangeLog
+2011-10-11  Uwe Schulzweida  <Uwe.Schulzweida@zmaw.de>
+
+       * zaxisCompare: set epsilon from 0 to 1e-9 [request: Felicia Brisc]
+	
 2011-10-06  Uwe Schulzweida  <Uwe.Schulzweida@zmaw.de>

 	* added level type ZAXIS_TOA, ZAXIS_SEA_BOTTOM, ZAXIS_ATMOSPHERE [request: Drte Liermann]

--- a/src/varscan.c
+++ b/src/varscan.c
@@ -701,11 +701,11 @@ int zaxisCompare(int zaxisID, int zaxistype, int nlevels, int lbounds, double *l
 	  const double *dlevels;
 	  char zlongname[256];
 	  char zunits[256];
-	  
+
 	  dlevels = zaxisInqLevelsPtr(zaxisID);
 	  for ( levelID = 0; levelID < nlevels; levelID++ )
 	    {
-	      if ( fabs(dlevels[levelID] - levels[levelID]) > 0 )
+	      if ( fabs(dlevels[levelID] - levels[levelID]) > 1.e-9 )
 		break;
 	    }

@@ -731,8 +731,8 @@ int zaxisCompare(int zaxisID, int zaxistype, int nlevels, int lbounds, double *l
 }


-int varDefZaxis(int vlistID, int zaxistype, int nlevels, double *levels, int lbounds, 
-		double *levels1, double *levels2, int vctsize, double *vct, char *name, 
+int varDefZaxis(int vlistID, int zaxistype, int nlevels, double *levels, int lbounds,
+		double *levels1, double *levels2, int vctsize, double *vct, char *name,
 		char *longname, char *units, int prec, int mode, int ltype)
 {
  /*