Changeset 381 for trunk

Timestamp:

10/03/08 09:26:07 (16 years ago)

Author:

smasson

Message:

bugfix whith NaN values in extrapolate

File:

• trunk/SRC/Interpolation/extrapolate.pro (modified) (3 diffs)

trunk/SRC/Interpolation/extrapolate.pro

@@ -78 +78 @@
 ; the x and y periodicity.
 ;
-  ztmp = bytarr(nx+2, ny+2)
+  msk = bytarr(nx+2, ny+2)
   IF n_elements(maskinput) EQ 1 AND maskinput[0] EQ -1 THEN maskinput = replicate(1b, nx, ny)
   IF n_elements(maskinput) NE nx*ny THEN BEGIN
@@ -84 +84 @@
     return, -1
   ENDIF
-  ztmp[1:nx, 1:ny] = byte(maskinput)
-  msk = temporary(ztmp)
+  msk[1:nx, 1:ny] = byte(maskinput)
 ;
   ztmp = replicate(1.e20, nx+2, ny+2)
@@ -94 +93 @@
   ENDIF
 ; remove NaN points if there is some...
-  nan = where(finite(ztmp) EQ 0, cnt_nan)
-  IF cnt_nan NE 0 THEN ztmp[temporary(nan)] = 1.e20
+  finztmp = finite(ztmp)
+  nan = where(finztmp EQ 0, cnt_nan)
+  IF cnt_nan NE 0 THEN BEGIN 
+    ztmp[temporary(nan)] = 1.e20
+    msk = temporary(msk) * temporary(finztmp)
+  ENDIF ELSE finztmp = -1 ; free memory
   z = temporary(ztmp)
   nx2 = nx+2