Changeset 502

Timestamp:

10/11/16 20:40:46 (8 years ago)

Author:

smasson

Message:

bugfix in extrapolate for 1d array andwhen using smwin

File:

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

trunk/SRC/Interpolation/extrapolate.pro

@@ -71 +71 @@
 ; check the number of iteration used in the extrapolation.
   szin = size(zinput)
+  IF szin[0] EQ 1 THEN BEGIN
+    zinput = reform(zinput, szin[1], 1)
+    maskinput = reform(maskinput, szin[1], 1)
+    szin = size(zinput)
+    fillxdir = 1
+  ENDIF
+  
   IF szin[0] NE 2 THEN return, -1. ELSE szin = szin[1:2]
   nx = szin[0]
   ny = szin[1]
   IF n_elements(nb_iteration) EQ 0 THEN nb_iteration = max(szin)
+  IF n_elements(smwin) EQ 0 THEN smwin = 0
+  smwin = smwin < min(szin)
+  IF n_elements(nsmooth) EQ 0 THEN nsmooth = 25
+  IF smwin GT 1 AND nsmooth NE 0 THEN nb_iteration = max(szin)
   IF nb_iteration EQ 0 THEN return, zinput
 ; take care of the boundary conditions...
@@ -231 +242 @@
 ; smooth the filled values
 ;---------------------------------------------------------------
-  IF n_elements(smwin) EQ 0 THEN smwin = 0
-  IF smwin EQ 0 OR smwin EQ 1 THEN return, z
-  IF n_elements(nsmooth) EQ 0 THEN nsmooth = 25
-  IF nsmooth EQ 0 THEN return, z
+  IF smwin LE 1 OR nsmooth EQ 0 THEN return, z
 
 ; add extra bands to avoid edge errors...