Changeset 326 for trunk/SRC/ToBeReviewed/PLOTS/DESSINE/plt.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/plt.pro (modified) (2 diffs)

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

@@ -147 +147 @@
 ; @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}
@@ -342 +344 @@
     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, x_periodic = keyword_set(key_periodic)*(nx EQ jpi) ))
+    trifield = triangule(temporary(tmp)*mask, coinmonte = coinmontemask $
+               , coindescend = coindescendmask, keep_cont = 0)
+    indicezoommask = (lindgen(jpi, jpj))[firstx:lastx, firsty:lasty]
+  ENDIF
 ; filling the mask values
 ; we fill only masknan or we fill mask*masknan?