Changeset 326 for trunk/SRC/ToBeReviewed/PLOTS/DESSINE/pltz.pro

Timestamp:

12/13/07 13:55:52 (17 years ago)

Author:

smasson

Message:

bugfix with max values of strictfill keyword

File:

• trunk/SRC/ToBeReviewed/PLOTS/DESSINE/pltz.pro (modified) (3 diffs)

trunk/SRC/ToBeReviewed/PLOTS/DESSINE/pltz.pro

@@ -129 +129 @@
 ; @keyword STRICTFILL
 ; Activate this keyword to that the filling of contours be
-; precisely done between the min and the max specified letting values inferior at the
-; specified min and values superior at the specified max in white.
+; precisely done between the min and the max specified. Values inferior/superior at the
+; specified min/max are not filled. Note that max values a considered
+; as missing values so cells containing one or more corners with
+; values above max will have no contours drawn through them.  
 ;
 ; @keyword STYLE {default=style=0}
@@ -305 +307 @@
     IF NOT keyword_set(nofill) THEN z2d[where(masknan EQ 0)] = max
   ENDIF ELSE masknan = 1
+  IF keyword_set(strictfill) THEN BEGIN
+    tmp = z2d ge max
+    tmp = 1b - byte(extrapolate(tmp, tmp, 1))
+    key_save = key_periodic & key_periodic = 0
+    trifield = triangule(temporary(tmp)*mask, coinmonte = coinmontemask $
+                         , coindescend = coindescendmask, keep_cont = 0)
+    key_periodic = key_save
+  ENDIF
   IF n_elements(fillxdir) EQ 0  THEN fillxdir = 1
   IF keyword_set(fillxdir) THEN BEGIN
@@ -335 +345 @@
 ;   if keyword_set(notri) then trifield = -1 $
 ;   ELSE trifield = triangule(mask,/basic)
-  if (usetri GE 1 AND (vargrid EQ 'T' OR vargrid EQ 'W')) $
-    OR (usetri EQ 2 AND (vargrid NE 'T' AND vargrid NE 'W')) THEN $
-    trifield = triangule(mask, /basic)
+  if NOT keyword_set(strictfill) AND ((usetri GE 1 AND (vargrid EQ 'T' OR vargrid EQ 'W')) $
+                                      OR (usetri EQ 2 AND (vargrid NE 'T' AND vargrid NE 'W'))) THEN $
+                                         trifield = triangule(mask, /basic)
 ;
   IF NOT keyword_set(endpoints)  THEN BEGIN