Changeset 417

Timestamp:

12/03/09 16:27:44 (15 years ago)

Author:

smasson

Message:

add new keywords to file_interp

File:

• trunk/SRC/Interpolation/file_interp.pro (modified) (7 diffs)

trunk/SRC/Interpolation/file_interp.pro

@@ -5 +5 @@
 ;-
 FUNCTION call_interp2d, data, inlon, inlat, inmask, outlon, outlat $
-                      , INIRR=inirr, METHOD=method, SMOOTH=smooth $
-                      , WEIG=weig, ADDR=addr, MISSING_VALUE=missing_value, _EXTRA=ex
+                        , INIRR = inirr, METHOD = method, SMOOTH = smooth $
+                        , WEIG = weig, ADDR = addr, MISSING_VALUE = missing_value $
+                        , GETHAN = gethan, LETHAN = lethan, _EXTRA = ex
 ;
   compile_opt idl2, strictarrsubs
@@ -44 +45 @@
   ENDELSE
 
+  IF n_elements(gethan) EQ 1 THEN data = gethan > temporary(data)
+  IF n_elements(lethan) EQ 1 THEN data = temporary(data) < lethan
+
   if intype LE 3 THEN data = round(temporary(data))
 
@@ -217 +221 @@
 ; same as inyaxisname but for output grid file gridout
 ;
+; @keyword GETHAN
+; to force interpolated data to be always > value defined by gethan (for example gethan = 0.)
+;
+; @keyword LETHAN
+; to force interpolated data to be always < value defined by lethan (for example lethan = 0.)
+;
 ; @keyword  _EXTRA
 ; to use <pro>extrapolate</pro>, <pro>extrapsmooth</pro> and <pro>fromreg</pro> keywords
@@ -246 +256 @@
 ;
 ;-
-PRO file_interp, filein, fileout, gridout, GRIDIN=gridin, MASKIN=maskin, MASKOUT=maskout $
-               , KEEP=keep, REMOVE=remove, METHOD=method, SMOOTH=smooth $
-               , SET_XDIMNAME=set_xdimname, SET_YDIMNAME=set_ydimname $
-               , SET_XAXISNAME=set_xaxisname, SET_YAXISNAME=set_yaxisname $
-               , INMASKNAME=inmaskname, ININVMASK=ininvmask $
-               , INUSEASMASK=inuseasmask, INMISSING_VALUE=inmissing_value $
-               , INADDSCL_BEFORE=inaddscl_before, INTESTOP = intestop $
-               , OUTMASKNAME=outmaskname, OUTINVMASK=outinvmask $
-               , OUTUSEASMASK=outuseasmask, OUTMISSING_VALUE=outmissing_value $
-               , OUTADDSCL_BEFORE=outaddscl_before, OUTTESTOP = outtestop $
-               , MISSING_VALUE=MISSING_VALUE, WEIG=weig, ADDR=addr $
-               , INXAXISNAME=inxaxisname, INYAXISNAME=inyaxisname $
-               , OUTXAXISNAME=outxaxisname, OUTYAXISNAME=outyaxisname $
-               , _EXTRA=ex
+PRO file_interp, filein, fileout, gridout, GRIDIN = gridin, MASKIN = maskin, MASKOUT = maskout $
+                 , KEEP = keep, REMOVE = remove, METHOD = method, SMOOTH = smooth $
+                 , SET_XDIMNAME = set_xdimname, SET_YDIMNAME = set_ydimname $
+                 , SET_XAXISNAME = set_xaxisname, SET_YAXISNAME = set_yaxisname $
+                 , INMASKNAME = inmaskname, ININVMASK = ininvmask $
+                 , INUSEASMASK = inuseasmask, INMISSING_VALUE = inmissing_value $
+                 , INADDSCL_BEFORE = inaddscl_before, INTESTOP = intestop $
+                 , OUTMASKNAME = outmaskname, OUTINVMASK = outinvmask $
+                 , OUTUSEASMASK = outuseasmask, OUTMISSING_VALUE = outmissing_value $
+                 , OUTADDSCL_BEFORE = outaddscl_before, OUTTESTOP = outtestop $
+                 , MISSING_VALUE = MISSING_VALUE, WEIG = weig, ADDR = addr $
+                 , INXAXISNAME = inxaxisname, INYAXISNAME = inyaxisname $
+                 , OUTXAXISNAME = outxaxisname, OUTYAXISNAME = outyaxisname $
+                 , GETHAN = gethan, LETHAN = lethan $
+                 , _EXTRA = ex
 ;
   compile_opt idl2, strictarrsubs
@@ -439 +450 @@
               IF interp THEN data = call_interp2d(temporary(data), inlon, inlat, inmask[*, *, 0], outlon, outlat $
                                                   , INIRR = inirr, METHOD = method, SMOOTH = smooth $
-                                                  , WEIG = weig, ADDR = addr, MISSING_VALUE = var_missing_value, _extra = ex)
+                                                  , WEIG = weig, ADDR = addr, MISSING_VALUE = var_missing_value $
+                                                  , GETHAN = gethan, LETHAN = lethan, _extra = ex)
               IF interp AND n_elements(data) EQ 1 THEN data = replicate(data, jpiout, jpjout)
               ncdf_varput, outid, outvarid[i], temporary(data)
@@ -455 +467 @@
                 IF interp THEN data = call_interp2d(temporary(data), inlon, inlat, temporary(tmpmsk), outlon, outlat $
                                                     , INIRR = inirr, METHOD = method, SMOOTH = smooth $
-                                                    , WEIG = weig, ADDR = addr, MISSING_VALUE = var_missing_value, _extra = ex)
+                                                    , WEIG = weig, ADDR = addr, MISSING_VALUE = var_missing_value $
+                                                    , GETHAN = gethan, LETHAN = lethan, _extra = ex)
                 IF interp AND n_elements(data) EQ 1 THEN data = replicate(data, jpiout, jpjout)
                 ncdf_varput, outid, outvarid[i], temporary(data), offset = off, count = outcnt
@@ -472 +485 @@
                   IF interp THEN data = call_interp2d(temporary(data), inlon, inlat, temporary(tmpmsk), outlon, outlat $
                                                       , INIRR = inirr, METHOD = method, SMOOTH = smooth $
-                                                      , WEIG = weig, ADDR = addr, MISSING_VALUE = var_missing_value, _extra = ex)
+                                                      , WEIG = weig, ADDR = addr, MISSING_VALUE = var_missing_value $
+                                                      , GETHAN = gethan, LETHAN = lethan, _extra = ex)
                   IF interp AND n_elements(data) EQ 1 THEN data = replicate(data, jpiout, jpjout)
                   ncdf_varput, outid, outvarid[i], temporary(data), offset = off, count = outcnt