<div dir="ltr">The data has been closely examined. The problem was related to downloading ERA5 data, for some reason a dimension "expver" (presumably experiment version?) appeared in the variable. This dimension did not appear after a fresh download of the data (minor changes in the script but the exact reason remains obscure to me).<div>In any case, if someone happens to have problems plotting ERA5 data, pls check if the dimension "expver" appears, and if necessary re-download the data to end up with e.g. z(time, lev, lat, lon).</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Feb 8, 2021 at 8:38 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>The 'golden rule' of data processing is to look at the data.</div><div><br></div><div>netcdf geop_500_975_89_20 {<br>dimensions:<br> longitude = 1440 ;<br> latitude = 721 ; <b> ; has pole points</b><br> level = 2 ;<br> expver = 2 ; <b>; not sure what this is BUT, possibly, two different physics packages</b></div><div> <b>; OR, possibly, an ensemble of size two OR, ???</b><br></div><div> time = 160 ;<br>variables:<br> float longitude(longitude) ;<br> longitude:units = "degrees_east" ;<br> longitude:long_name = "longitude" ;<br> float latitude(latitude) ;<br> latitude:units = "degrees_north" ;<br> latitude:long_name = "latitude" ;<br> int level(level) ;<br> level:units = "millibars" ;<br> level:long_name = "pressure_level" ;<br> int expver(expver) ;<br> expver:long_name = "expver" ;<br> int time(time) ;</div><div><div> time:units = "hours since 1900-01-01 00:00:00.0" ;<br> time:long_name = "time" ;<br> time:calendar = "gregorian" ;<br> <b><span style="color:rgb(0,0,255)">short</span> z(time, expver, level, latitude, longitude)</b> ; <====== rank 5<br><b> z:scale_factor = 0.898208868468462 ; netCDF CF [and COARDS] conventions<br></b> <b>z:add_offset = 28656.6907881439 </b>; <b>; scale_factor and add_offet must be present</b><br> <b>z:_FillValue = -32767s ; ; CF </b><br> z:missing_value = -32767s ; ; COARDS <br> z:units = "m**2 s**-2" ; ; COARDS and CF conventions<br> z:long_name = "Geopotential" ; ; COARDS convention<br> z:standard_name = "geopotential" ; ; CF convention<br><br>======================================================</div><div><br></div><div><div><div><div><div><div id="m_-1914730850972550663gmail-m_6867493540619869841m_949011946592702635gmail-m_-4094324489795839166gmail-:2j2"><div><div><div><div><div id="m_-1914730850972550663gmail-m_6867493540619869841m_949011946592702635gmail-m_-4094324489795839166gmail-:2zl"><div id="m_-1914730850972550663gmail-m_6867493540619869841m_949011946592702635gmail-m_-4094324489795839166gmail-:2j4"><div dir="ltr"><span><div>dir="/Users/nuovoaccount/Documents/PYTHON/"<br>fil ="<a href="http://geop_500_975_89_20.nc" target="_blank">geop_500_975_89_20.nc</a>"<br>f = addfile(dir+fil,"r")<br>z = <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml" target="_blank"><b>short2flt</b></a>(f->z)<br></div></span><div>printVarSummary(z) ; <b>(time, expver, level, latitude, longitude)</b> </div><div>printMinMax(z,0)</div><div><br></div><div><div><div>---<br></div><div>or, <br></div><div> [a] noting that NCL<i> eliminates degenerate dimensions</i> [size 1]</div><div> [b] no information is lost, since NCL will add variable attributes for the degenerate dimensions: <b>expver, level</b></div><div><br></div><div>z = <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml" target="_blank"><b>short2flt</b></a>( f->z(:,1,1,:,:) ) ; second ensemble member and second level<br></div><div><div>printVarSummary(z) ; <b>(time, latitude, longitude) </b>with attributes <b>level </b>and<b> expver</b> </div><div>printMinMax(z,0)</div></div><div><div id="m_-1914730850972550663gmail-m_6867493540619869841gmail-m_-5170570634263361908m_-4094324489795839166gmail-:2zw"><img src="https://ssl.gstatic.com/ui/v1/icons/mail/images/cleardot.gif"></div><div><br></div><div>HTH</div><div><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Feb 8, 2021 at 8:46 AM Anne Seidenglanz via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>I am not new to NCL but encounter a rather new problem when dealing with ERA5 data (in .nc format, see script attached). I am trying to produce a simple test plot for geopotential height, and somehow all my data seems to be missing? The error message appears: </div><div>warning:ContourPlotInitialize: no valid values in scalar field; ContourPlot not possible:[errno=1101] <br></div><div>When I printed a subset of the dataset, all the values seemed to be missing values. </div><div><br></div><div>I have:</div><div>- tested that not all data is missing in Z with "print("Is Z all missing ? " + all(ismissing(Z)))". It gave me False.</div><div>- I have tried both, reading in the data using "short2flt" to account for scale factor and offset value, and without short2flt. In both cases appears the missing value problem and the same error message. </div><div><br></div><div>Where have I gone wrong?</div><div><br></div><div>Thx for any help</div><div>Anne </div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div></div>
</blockquote></div></div>