<div dir="ltr">Hello, Mary, <div><br></div><div>I modified my script following your first suggestion. The problem is still there. The warning msg shows:</div><div>-------------------------------------------------</div><div><div>Variable: o3_p</div><div>Type: float</div><div>Total Size: 831744 bytes</div><div>            207936 values</div><div>Number of Dimensions: 2</div><div>Dimensions and sizes:   [lat | 361] x [lon | 576]</div><div>Coordinates:</div><div>            lat: [ -90..  90]</div><div>            lon: [-180..179.375]</div><div>Number Of Attributes: 13</div><div>  lev :  900</div><div>  time :        540</div><div>  long_name :   ozone mixing ratio</div><div>  units :       ppbv</div><div>  _FillValue :  1e+15</div><div>  missing_value :       1e+15</div><div>  fmissing_value :      1e+15</div><div>  scale_factor :         1</div><div>  add_offset :   0</div><div>  standard_name :       ozone_mass_mixing_ratio</div><div>  vmax :        1e+15</div><div>  vmin :        -1e+15</div><div>  valid_range : ( -1e+15, 1e+15 )</div><div>warning:ContourPlotDraw: out of range coordinates encountered; standard AreaFill rendering method may be unreliable;</div><div> consider setting the resource trGridType to "TriangularMesh" if coordinates contain missing values</div></div><div>-------------------------------------</div><div>And I attached the updated script. I wounder if I need to add triangular meshes resource in my script?</div><div><br></div><div>Thanks! </div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Sep 7, 2017 at 8:46 AM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-size:small">I think the issue is that "o3_p", which is the variable you are trying to plot, doesn't have any metadata attached to it.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">When you do calculations like this in NCL:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div class="gmail_default">   o3= a->O3(3,4,:,:) ; 18GTM = 11AM PST                                        </div><div class="gmail_default">  printVarSummary(o3)</div><div class="gmail_default">   o3_p =(o3*1e9)*(28.97/48)<br></div><div class="gmail_default"><br></div><div class="gmail_default">The "o3" will have metadata because no computations were done when you read the variable off the file.</div><div class="gmail_default"><br></div><div class="gmail_default">However, since o3_p involves calculations, NCL will strip all the metadata and the only thing that o3_p will retain from o3 is the _FillValue attribute, if any.</div><div class="gmail_default"><br></div><div class="gmail_default">As a "trick", you can force o3_p to have the same metadata by doing the following:</div><div class="gmail_default"><br></div><div class="gmail_default"><div class="gmail_default"><font face="monospace, monospace">   o3_p = a->O3(3,4,:,:) ; 18GTM = 11AM PST                                        </font></div><div class="gmail_default"><font face="monospace, monospace">   o3_p =(o3_p*1e9)*(28.97/48)</font></div><div class="gmail_default"><br></div><div class="gmail_default">The first line will read the data into o3_p, and the metadata will all be there. Then, when you do the calculation, since o3_p already exists, no metadata will be removed.</div><div class="gmail_default"><br></div><div class="gmail_default">Note that I removed the use of o3, just to save memory. If you need to keep o3 and o3_p around, then you can do:</div><div class="gmail_default"><br></div><div class="gmail_default"><div class="gmail_default"><font face="monospace, monospace">   o3 = a->O3(3,4,:,:) ; 18GTM = 11AM PST <br>   o3_p = o3           ; trick to copy metadata</font></div><div class="gmail_default"><font face="monospace, monospace">   o3_p =(o3*1e9)*(28.97/48)</font></div></div><div class="gmail_default"><br></div><div class="gmail_default">You should however, update the long_name and units attributes of o3_p, to reflect what the new values represent.</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Wed, Sep 6, 2017 at 5:05 PM, Ying Song <span dir="ltr"><<a href="mailto:ying.song@sjsu.edu" target="_blank">ying.song@sjsu.edu</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Hello, ncl helper,<div><br></div><div>I trying to plot MERRA-2 reanalysis data in NCL using gsn_contour_map. I am only interested in west coast of US.  But I can not get the zoom in plot correctly by using resources:</div><div><div><br></div><div>  res@mpMinLatF               = 30</div><div>  res@mpMaxLatF               = 45</div><div>  res@mpMinLonF               = -130</div><div>  res@mpMaxLonF               = -110</div></div><div><br></div><div>The warning shows:</div><div><div>warning:ContourPlotDraw: out of range coordinates encountered; standard AreaFill rendering method may be unreliable;</div><div> consider setting the resource trGridType to "TriangularMesh" if coordinates contain missing values</div></div><div><br></div><div>Attachment is my script. The info of MERRA-2 O3 data from printVarSummary is:</div><div> -----------------------------<wbr>-</div><div><div>Variable: o3</div><div>Type: float</div><div>Total Size: 279465984 bytes</div><div>            69866496 values</div><div>Number of Dimensions: 4</div><div>Dimensions and sizes:   [time | 8] x [lev | 42] x [lat | 361] x [lon | 576]</div><div>Coordinates:</div><div>            time: [0..1260]</div><div>            lev: [1000..0.1000000014901161]</div><div>            lat: [ -90..  90]</div><div>            lon: [-180..179.375]</div><div>Number Of Attributes: 11</div><div>  long_name :   ozone_mass_mixing_ratio</div><div>  units :       kg kg-1</div><div>  _FillValue :  1e+15</div><div>  missing_value :       1e+15</div><div>  fmissing_value :      1e+15</div><div>  scale_factor :         1</div><div>  add_offset :   0</div><div>  standard_name :       ozone_mass_mixing_ratio</div><div>  vmax :        1e+15</div><div>  vmin :        -1e+15</div><div>  valid_range : ( -1e+15, 1e+15 )</div></div><div>------------------------------<wbr>---------------------------</div></div>
<br></div></div>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>