<html>
  <head>
    <meta content="text/html; charset=windows-1252"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    Have you tried a <tt>printVarSummary </tt>to see what your
    metadata actually is, both before and after <tt>where</tt>?<br>
    <br>
    - Walter<br>
    <br>
    <div class="moz-cite-prefix">On 24-Jun-15 12:18, Erik J.
      Schaffernicht wrote:<br>
    </div>
    <blockquote cite="mid:558AD839.5020207@uni-koeln.de" type="cite">Now,
      I rewrote my ncl-script because on first sight it seemed you are
      right.  But when it comes to gsn_csm_contour_map  (or somewhere
      else, there is no line number given with the warning messages)
      <br>
      I get these warnings:
      <br>
      <br>
      (0)    check_for_y_lat_coord: Warning: Data either does not
      contain a valid latitude coordinate array or doesn't contain one
      at all.
      <br>
      (0)    A valid latitude coordinate array should have a 'units'
      attribute equal to one of the following values:
      <br>
      (0)        'degrees_north' 'degrees-north' 'degree_north' 'degrees
      north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN'
      'deg north'
      <br>
      <br>
      and the only thing I changed was this:
      <br>
      <br>
      w2=sediProbab
      <br>
      w2=where(landAreaFrV .lt. 0.5, w2@_FillValue, w2)
      <br>
      copy_VarMeta(sediProbab,w2)
      <br>
      sediProbab=w2
      <br>
      <br>
      to
      <br>
      <br>
      sediProbab=where(landAreaFrV .lt. 0.5, sediProbab@_FillValue,
      sediProbab)
      <br>
      <br>
      <br>
      So, what is happening there with the where()
      <br>
      <br>
      <br>
      <br>
      <br>
      <br>
      On 23.06.2015 14:58, Carl Schreck wrote:
      <br>
      <blockquote type="cite">You should be able to simply do:
        <br>
        <br>
        sediProbab = where(landAreaFrV .lt. 0.5, sediProbab@_FillValue,
        sediProbab)
        <br>
        <br>
        sediProbab's metadata will be unchanged.
        <br>
        <br>
             Carl
        <br>
        <br>
        On Tue, Jun 23, 2015 at 7:14 AM, Erik J. Schaffernicht &lt;
        <br>
        <a class="moz-txt-link-abbreviated" href="mailto:eschaffe@uni-koeln.de">eschaffe@uni-koeln.de</a>&gt; wrote:
        <br>
        <br>
        <blockquote type="cite">Is there an easier/shorter way(?) than
          the following to keep the metaData:
          <br>
          <br>
                                       (sediProbab is a matrix with
          metaData)
          <br>
          <br>
          w2=sediProbab
          <br>
          w2=where(landAreaFrV .lt. 0.5, w2@_FillValue, w2)
          <br>
          copy_VarMeta(sediProbab,w2)
          <br>
          sediProbab=w2
          <br>
          <br>
          <br>
          <br>
          _______________________________________________
          <br>
          ncl-talk mailing list
          <br>
          <a class="moz-txt-link-abbreviated" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>
          <br>
          List instructions, subscriber options, unsubscribe:
          <br>
          <a class="moz-txt-link-freetext" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a>
          <br>
          <br>
          <br>
        </blockquote>
        <br>
      </blockquote>
      <br>
      <br>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <br>
      <pre wrap="">_______________________________________________
ncl-talk mailing list
<a class="moz-txt-link-abbreviated" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>
List instructions, subscriber options, unsubscribe:
<a class="moz-txt-link-freetext" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a>
</pre>
    </blockquote>
    <br>
    <pre class="moz-signature" cols="72">-- 
Walter Kolczynski, Jr.
Global Ensemble Team
NOAA/NWS/NCEP/EMC (via I.M. Systems Group)
(301) 683-3781</pre>
  </body>
</html>