<div dir="ltr"><div>It would be nice if you provided mor information on the file.</div><div><br></div><div>[1] The binary file is "LittleEndian" not "BigEndian"</div><div>[2] I think there is a header record because nlat*mlon=424*1327=562648 <br></div><div> is less than the total file length.</div><div><br></div><div>Attached is a guess at how to read the file.</div><div><br></div><div>%> ncl VIENTO_MUNDO.ncl</div><div><br></div><div>Good luck<br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 9, 2020 at 11:47 AM Paimazumder, Debasish via ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@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 lang="EN-US">
<div class="gmail-m_5781353221990621946WordSection1">
<p class="MsoNormal"><span style="font-family:"Georgia","serif"">Hi NCL users<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-family:"Georgia","serif"">I’m trying to open/read/plot a grd file in NCL. Frist I read the file in R to identify the dimensions<u></u><u></u></span></p>
<p class="MsoNormal" style="line-height:11.25pt;background:white none repeat scroll 0% 0%;word-break:break-all">
<span style="font-size:9pt;font-family:"Lucida Console";color:rgb(0,112,192)">> r<u></u><u></u></span></p>
<p class="MsoNormal" style="line-height:11.25pt;background:white none repeat scroll 0% 0%;word-break:break-all">
<span style="font-size:9pt;font-family:"Lucida Console";color:rgb(0,112,192)">class : RasterLayer
<u></u><u></u></span></p>
<p class="MsoNormal" style="line-height:11.25pt;background:white none repeat scroll 0% 0%;word-break:break-all">
<span style="font-size:9pt;font-family:"Lucida Console";color:rgb(0,112,192)">dimensions : 424, 1327, 562648 (nrow, ncol, ncell)<u></u><u></u></span></p>
<p class="MsoNormal" style="line-height:11.25pt;background:white none repeat scroll 0% 0%;word-break:break-all">
<span style="font-size:9pt;font-family:"Lucida Console";color:rgb(0,112,192)">resolution : 0.2713213, 0.2713213 (x, y)<u></u><u></u></span></p>
<p class="MsoNormal" style="line-height:11.25pt;background:white none repeat scroll 0% 0%;word-break:break-all">
<span style="font-size:9pt;font-family:"Lucida Console";color:rgb(0,112,192)">extent : -180.1357, 179.9078, -55.13566, 59.90459 (xmin, xmax, ymin, ymax)<u></u><u></u></span></p>
<p class="MsoNormal" style="line-height:11.25pt;background:white none repeat scroll 0% 0%;word-break:break-all">
<span style="font-size:9pt;font-family:"Lucida Console";color:rgb(0,112,192)">crs : NA
<u></u><u></u></span></p>
<p class="MsoNormal" style="line-height:11.25pt;background:white none repeat scroll 0% 0%;word-break:break-all">
<span style="font-size:9pt;font-family:"Lucida Console";color:rgb(0,112,192)">source : H:/MunichRe_Wind/VIENTO_MUNDO_TR100_INT1.grd
<u></u><u></u></span></p>
<p class="MsoNormal" style="line-height:11.25pt;background:white none repeat scroll 0% 0%;word-break:break-all">
<span style="font-size:9pt;font-family:"Lucida Console";color:rgb(0,112,192)">names : VIENTO_MUNDO_TR100_INT1
<u></u><u></u></span></p>
<p class="MsoNormal" style="line-height:11.25pt;background:white none repeat scroll 0% 0%;word-break:break-all">
<span style="font-size:9pt;font-family:"Lucida Console";color:rgb(0,112,192)">values : 0, 444.2024 (min, max)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-family:"Georgia","serif"">When I read it in NCL, the max and min number are different
<u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:"Courier New";color:rgb(0,112,192)">Variable: x<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)">Type: float<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)">Total Size: 2250592 bytes<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> 562648 values<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)">Number of Dimensions: 2<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)">Dimensions and sizes: [lat | 424] x [lon | 1327]<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)">Coordinates:<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> lat: [-55.13566..59.90459]<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> lon: [-180.1357..179.9078]<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)">Number Of Attributes: 3<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> units : ??<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> long_name : 100 year windspeed<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> _FillValue : -999<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)">(0) min(x)=-3.39951e+38<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)">(0) max(x)=3.40086e+38<u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><span style="font-family:"Georgia","serif"">Max/Min number in R looks more reasonable. Here is file –<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-family:"Georgia","serif""><a href="https://data.humdata.org/dataset/cyclone-wind-100-years-return-period" target="_blank"><span style="color:blue">https://data.humdata.org/dataset/cyclone-wind-100-years-return-period</span></a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-family:"Georgia","serif""><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-family:"Georgia","serif"">here is my NCL code-<u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><span style="font-size:9pt;color:rgb(0,112,192)"> </span><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> diri = "./"<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> fili = "VIENTO_MUNDO_TR100_INT1.grd" ; DSET<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> fName = diri+fili<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> nlat = 424 ; YDEF<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> mlon = 1327 ; XDEF<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> ; generate lat/lon<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> lon = fspan(-180.1357,179.9078,1327)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> lon!0 = "lon"<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> lon@long_name = "longitude"<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> lon@units = "degrees_east"<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> lat = fspan(-55.13566,59.90459,424)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> lat!0 = "lat"<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> lat@long_name = "latitude"<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> lat@units = "degrees_north"<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> ; create an array to contain data<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> UNDEF = -999. ; UNDEF<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> x = new ( (/nlat,mlon/), float, UNDEF)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> x!0 = "lat"<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> x!1 = "lon"<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> x&lat = lat<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> x&lon = lon<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> x@long_name = "100 year windspeed" ; VARS<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> x@units = "??"<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)">;setfileoption("bin","ReadByteOrder","Native") ; a033 default<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)">;setfileoption("bin","ReadByteOrder","LittleEndian")<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> setfileoption("bin","ReadByteOrder","BigEndian")<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> ; read each record: store in x<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> ; the ::-1 changes the latitude order<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> x = fbinrecread(fName,0, (/nlat,mlon/), "float")<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> printVarSummary(x)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> print ("min(x)="+min(x))<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9pt;font-family:"Courier New";color:rgb(0,112,192)"> print ("max(x)="+max(x))</span><span style="font-size:9pt;font-family:"Courier New""><u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><span style="font-family:"Georgia","serif"">with regards<u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><b><span style="font-size:12pt;color:rgb(0,176,240)">Debasish PaiMazumder, PhD |</span></b><span style="font-size:12pt;color:rgb(102,102,102)">Tel
<a href="tel:%2B1%20215%20255%206371" target="_blank"><span style="color:blue">+1 215 255 6052</span></a> | Cell
<a href="tel:%2B1%20718%20913%206946" target="_blank"><span style="color:blue">+1 267 449 5226</span></a> |</span><u><span style="font-size:12pt;color:blue"><a href="mailto:Debasish.Paimazumder@aig.com" target="_blank">Debasish.Paimazumder@aig.com</a></span></u><span style="font-size:12pt;color:blue">
</span><b><span style="font-size:12pt;color:rgb(0,176,240)"><u></u><u></u></span></b></p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
_______________________________________________<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/mailman/listinfo/ncl-talk</a></blockquote></div>