<div dir="ltr">Hello, world!<div><br></div><div>I want to draw several plots with data which is a binary file contains 6 variables(3-dimensional) and made like below:</div><div><br></div><div>open(2,file=ofile,status=&#39;<wbr>unknown&#39;,form=&#39;unformatted&#39;)<br></div><div><div>     write(2) (((hgtm(i,j,k),i=1,imx),j=1,<wbr>jmx),k=1,lx)</div><div>     write(2) (((ugrdm(i,j,k),i=1,imx),j=1,<wbr>jmx),k=1,lx)</div><div>     write(2) (((vgrdm(i,j,k),i=1,imx),j=1,<wbr>jmx),k=1,lx)</div><div>     write(2) (((vvelm(i,j,k),i=1,imx),j=1,<wbr>jmx),k=1,lx)</div><div>     write(2) (((tmpm(i,j,k),i=1,imx),j=1,<wbr>jmx),k=1,lx)</div><div>     write(2) (((rhm(i,j,k),i=1,imx),j=1,<wbr>jmx),k=1,lx)</div><div>  close(2)</div></div><div><br></div><div>And then I wrote a NCL script to read and plot them. </div><div>Here is the prob: <br>fatal:No coordinate variable exists for dimension (level) in variable (hgt)</div><div><br></div><div>My script:</div><div><br></div><div><div>  nlon = 659</div><div>  nlat = 539</div><div>  nlev = 17</div></div><div><br></div><div><div>  dims  = (/nlat,nlon,nlev/)</div><div>  lat      = fspan(slat,elat,nlat)</div><div>  lon      = fspan(slon,elon,nlon)</div><div> level = (/1000.,925.,850.,700.,600.,<wbr>500.,400.,300.,\</div><div>           250.,200.,150.,100.,70.,50.,<wbr>30.,20.,10./)</div><div><br></div><div>  hgt       = fbindirread(path_prs, 0, dims, &quot;float&quot;)</div><div>  u       = fbindirread(path_prs, 1, dims, &quot;float&quot;)</div><div>  v       = fbindirread(path_prs, 2, dims, &quot;float&quot;)</div><div>  vvel       = fbindirread(path_prs, 3, dims, &quot;float&quot;)</div><div>  temp       = fbindirread(path_prs, 4, dims, &quot;float&quot;)</div><div>  rh      = fbindirread(path_prs, 5, dims, &quot;float&quot;)</div></div><div><br></div><div><div>  lon!0         = &quot;lon&quot;</div><div>  lon@long_name = &quot;lon&quot;</div><div>  lon@units     = &quot;degrees_east&quot;</div><div>  lon&amp;lon       = lon</div><div><br></div><div>  lat!0         = &quot;lat&quot;</div><div>  lat@long_name = &quot;lat&quot;</div><div>  lat@units     = &quot;degrees_north&quot;</div><div>  lat&amp;lat       = lat</div><div><br></div><div>  level!0                 = &quot;level&quot;</div><div>  level@units             = &quot;hPa&quot;</div><div>  level@long_name         = &quot;isobaric_surface&quot;</div><div>  level&amp;level           = level</div><div><br></div><div>hgt!0                   = &quot;level&quot;</div><div>hgt!1                   = &quot;lat&quot;</div><div>hgt!2                   = &quot;lon&quot;</div><div>hgt&amp;level               = level &lt;--- this one makes error!</div><div>hgt&amp;lat                 = lat</div><div>hgt&amp;lon                 = lon</div><div>hgt&amp;level@units         = &quot;hPa&quot;</div><div>hgt&amp;lat@units           = &quot;degrees_north&quot;</div><div>hgt&amp;lon@units           = &quot;degrees_east&quot;</div><div>hgt@units               = &quot;gpm&quot;</div><div>hgt@long_name           = &quot;geopotential height&quot;</div></div><div><br></div><div><br></div><div><br></div></div>