Changeset 448 for trunk

Timestamp:

01/26/11 14:44:08 (13 years ago)

Author:

smasson

Message:

add SET_OUTMSKVAL in file_interp

File:

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

trunk/SRC/Interpolation/file_interp.pro

@@ -27 +27 @@
                         , INIRR = inirr, METHOD = method, SMOOTH = smooth $
                         , WEIG = weig, ADDR = addr, MISSING_VALUE = missing_value $
-                        , DIVX = divx, DIVY = divy, OUTMASK = outmask $
-                        , GETHAN = gethan, LETHAN = lethan, _EXTRA = ex
+                        , DIVX = divx, DIVY = divy, OUTMASK_IND = outmask_ind $
+                        , GETHAN = gethan, LETHAN = lethan, SET_OUTMSKVAL = set_outmskval $
+                        , _EXTRA = ex
 ;
   compile_opt idl2, strictarrsubs
@@ -73 +74 @@
   IF n_elements(gethan) EQ 1 THEN data = gethan > temporary(data)
   IF n_elements(lethan) EQ 1 THEN data = temporary(data) < lethan
-
+  IF outmask_ind[0] NE -1 THEN data[outmask_ind] = set_outmskval
+  
   if intype LE 3 THEN data = round(temporary(data))
-  IF outmask[0] NE -1 THEN data = temporary(data) * outmask
 
   RETURN, data
@@ -226 +227 @@
 ; used (in association with (IN/OUT)USEASMASK) to built the mask (that
 ; does not change from one record to another).
-; Note: this value will be apply to all interpolated variables
+; Note: this value will be applied to all interpolated variables
+;
+; @keyword SET_OUTMSKVAL {type=scalar}{default=defined by attribute missing_value or _fillvalue of input file}
+; (Re)define the masked (over land) value in output data.
+; Note: output mask as to be defined through the keyword maskout and
+; its associated keywords...
+; Note: do not mismatch with OUTMISSING_VALUE which are missing value
+; used (in association with OUTUSEASMASK) to built the mask (that
+; does not change from one record to another).
 ;
 ; @keyword ADDR {type=2d array or variable name}
@@ -303 +312 @@
                  , INXAXISNAME = inxaxisname, INYAXISNAME = inyaxisname $
                  , OUTXAXISNAME = outxaxisname, OUTYAXISNAME = outyaxisname $
-                 , GETHAN = gethan, LETHAN = lethan $
+                 , GETHAN = gethan, LETHAN = lethan, SET_OUTMSKVAL = set_outmskval $
                  , _EXTRA = ex
 ;
@@ -431 +440 @@
 ;
   outvarid = lonarr(ininq.nvars)
+  outmiss  = fltarr(ininq.nvars)
   FOR i = 0, ininq.nvars-1 DO BEGIN
     varinq = ncdf_varinq(inid, i)
@@ -469 +479 @@
 ; Variables attributes
       IF outvarid[i] NE - 1 THEN BEGIN
+
+        IF varinq.ndims GE 2 THEN BEGIN
+          interp = varinq.dim[0] EQ indimx AND varinq.dim[1] EQ indimy
+        ENDIF ELSE interp = 0b
+        CASE 1 OF
+          strlowcase(varinq.name) EQ strlowcase(inxaxisname):interp = 0b
+          strlowcase(varinq.name) EQ strlowcase(inyaxisname):interp = 0b
+          strlowcase(varinq.name) EQ strlowcase(inmaskname):interp = 0b
+        ENDCASE
+        
         FOR j = 0, varinq.natts-1 DO BEGIN
           name = ncdf_attname(inid, i, j)
-          dummy = ncdf_attcopy(inid, i, name, outid, outvarid[i])
+          CASE 1 OF
+            keyword_set(interp) AND strlowase(name) EQ '_fillvalue':BEGIN
+              ncdf_attget, inid, i, '_fillvalue', tmp & outmiss[i] = tmp
+            END
+            keyword_set(interp) AND strlowase(name) EQ 'missing_value':BEGIN
+              ncdf_attget, inid, i, 'missing_value', tmp & outmiss[i] = tmp
+            END
+            ELSE:dummy = ncdf_attcopy(inid, i, name, outid, outvarid[i])
+          ENDCASE
         ENDFOR
+ 
+        IF keyword_set(interp) AND outmask[0] NE -1 THEN BEGIN
+          IF n_elements(set_outmskval) NE 0 THEN outmiss[i] = set_outmskval
+          ncdf_attput, outid, outvarid[i], '_fillvalue', outmiss[i]
+          ncdf_attput, outid, outvarid[i], 'missing_value', outmiss[i] 
+        ENDIF
+
       ENDIF
     ENDIF ELSE outvarid[i] = -1
@@ -478 +513 @@
 ;
   ncdf_control, outid, /endef
+;
+  IF outmask[0] NE -1 THEN outmask_ind = where(outmask EQ 0) ELSE outmask_ind = -1
 ;
   FOR i = 0, ininq.nvars-1 DO BEGIN
@@ -514 +551 @@
                                                   , INIRR = inirr, METHOD = method, SMOOTH = smooth $
                                                   , WEIG = weig, ADDR = addr, MISSING_VALUE = var_missing_value $
-                                                  , DIVX = divx, DIVY = divy, OUTMASK = outmask $
+                                                  , SET_OUTMSKVAL = outmiss[i] $
+                                                  , DIVX = divx, DIVY = divy, OUTMASK_IND = outmask_ind $
                                                   , GETHAN = gethan, LETHAN = lethan, _extra = ex)
               IF interp AND n_elements(data) EQ 1 THEN data = replicate(data, jpiout, jpjout)
@@ -533 +571 @@
                                                     , INIRR = inirr, METHOD = method, SMOOTH = smooth $
                                                     , WEIG = weig, ADDR = addr, MISSING_VALUE = var_missing_value $
-                                                    , DIVX = divx, DIVY = divy, OUTMASK = outmask $
+                                                    , SET_OUTMSKVAL = outmiss[i] $
+                                                    , DIVX = divx, DIVY = divy, OUTMASK_IND = outmask_ind $
                                                     , GETHAN = gethan, LETHAN = lethan, _extra = ex)
                 IF interp AND n_elements(data) EQ 1 THEN data = replicate(data, jpiout, jpjout)
@@ -553 +592 @@
                                                       , INIRR = inirr, METHOD = method, SMOOTH = smooth $
                                                       , WEIG = weig, ADDR = addr, MISSING_VALUE = var_missing_value $
-                                                      , DIVX = divx, DIVY = divy, OUTMASK = outmask $
+                                                      , SET_OUTMSKVAL = outmiss[i] $
+                                                      , DIVX = divx, DIVY = divy, OUTMASK_IND = outmask_ind $
                                                       , GETHAN = gethan, LETHAN = lethan, _extra = ex)
                   IF interp AND n_elements(data) EQ 1 THEN data = replicate(data, jpiout, jpjout)