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

Timestamp:

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

Author:

smasson

Message:

bugfix with max values of strictfill keyword

File:

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

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

@@ -150 +150 @@
 ; @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}
@@ -467 +469 @@
          IF NOT keyword_set(nofill) THEN z2d[where(masknan EQ 0)] = max
          mask = temporary(mask)*temporary(masknan)
+       ENDIF
+      usetri = 0 ; default definition
+      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 = coinmontemaskaa $
+                             , coindescend = coindescendmaskaa, keep_cont = 0)
+        key_periodic = key_save
+        usetri = 2
       ENDIF
       z2d = remplit(z2d,nite=2+keyword_set(nan), mask = mask, /basique, _extra=ex)
       if NOT keyword_set(strictfill) then z2d = min > z2d <  max
-      usetri = 0 ; default definition
       if keyword_set(nan) then BEGIN
         triangulation = triangule(mask, /basic, coinmonte = coinmontemask $
@@ -482 +493 @@
                                      , coindescend = coindescendmask)
       ENDIF
+      IF NOT keyword_set(strictfill) THEN trifield = triangulation
 ;----------------------------------------------------------------------
       pltbase, z2d, xx, yy, mask, xx, yy, level_z2d, colnumb, contour = contour, /noerase $
        , c_linestyle=linestyle,c_labels=1-(indgen(n_elements(level_z2d)) MOD 2) $
-       , trichamp = triangulation, trimsk = triangulation, overplot = overplot $
+       , trichamp = trifield, trimsk = triangulation, overplot = overplot $
        , c_thick=thick, performance = key_performance, nofill = nofill, usetri = usetri $
        , coinmontemask=coinmontemask, coindescendmask=coindescendmask, _extra = ex