Changeset 10

Timestamp:

11/27/07 17:23:47 (16 years ago)

Author:

kolasinski

Message:

Ensure valmask is positive in nc_read.pro

File:

• trunk/procs/nc_read.pro (modified) (2 diffs)

trunk/procs/nc_read.pro

@@ -21 +21 @@
 ; decide which subdomain of data to read
 ;
+
    IF debug_w THEN print, 'keyword_set(ALL_DATA) : ', keyword_set(ALL_DATA) 
 
@@ -537 +538 @@
    ENDFOR 
 
+; set valmask to 1.e20
+
+; ensure valmask is positive
+
+   IF valmask LT 0 THEN BEGIN
+      print, '      *** Warning valmask is negative - changing sign: ', valmask
+      idx_t = where (field.data EQ valmask)
+      IF idx_t(0) NE -1 THEN field.data(idx_t) = -valmask
+      valmask = -valmask
+      idx_t=0 ; free memory
+   ENDIF
+
 ; min/max