<div dir="ltr">In addition to what Dennis has pointed out, the thing that seems problematic to me is that you have used 'addfiles' in a loop to open individual files.<div><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
do j = 0,numfils-1,1<br>a = addfiles(fils(j)+".nc","r")<br>time=a[:]->time</div><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
u1=a[:]->vrect_x(:,:,:)<br>u1@_FillValue =-99</div><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">v1=a[:]->vrect_y(:,:,:)<br>v1@_FillValue =-99<br>
</div><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">ntim = dimsizes(time)<br> end do</div></div><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
<br></div><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">This is not how addfiles is used normally. It can open all the files at once. Eliminate the do loop and use:</div>
<div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px"><br></div><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
a = addfiles(fils + ".nc","r")</div><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px"><br></div><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
Otherwise if you want to keep the loop use "addfile" and eliminate the list-type syntax (the brackets [:]).</div><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
-dave</div><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px"><br></div><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
<br></div><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px"><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Aug 20, 2014 at 7:17 AM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@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><div><div><div><div><div><div><div><div><div><div><div><div>The email text you wrote does not agree with the script you included.<br>
<br></div>[1] Please clean the script coding so it is nicely formatted. People do not have time to unravel code.<br>
<br></div>[2] As noted frequently on ncl-talk, users should include output from printVarSummary(...)<br></div> for selected variables. This is not just for ncl-talk to examine but also you. Please<br></div><div> examine the output.<br>
</div><div><br></div>[3] The email text states "The netcdf files are from the instrument but are not in the regular <br><div class=""> netcdf format meaning that the time is not included in the dimensions of
the variables <br></div> to be extracted". This is not clear.<br><br></div> Yet your code has<br><br> u1 = a[:]->vrect_x ; (:,:,:)<br><br></div> and subsequently<br><br> data(npt) = data(npt) + sprintf("%7.1f ",u1(nt,nl,ml))<br>
<br></div> where u1 references a time dimension via subscrip 'nt' . Very confusing.<br><br></div>[4] You use<br><br> x = (/ntim,nlat,mlon/)<br><br></div> If you were to use<br><br></div> printVarSummary(x)<br>
<br></div> you would see that it is a one dimensional array of tyope integer and length 3.<br><br></div> Later you use<br><br> data(npt) = data(npt) + sprintf("%10.3f ", x(nt))<br><br></div> why x(nt)????<br>
<br><br><div><div><div><div><div><div><div><div><div><div><div><br><br><br></div></div></div></div></div></div></div></div></div></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote"><div><div class="h5">
On Tue, Aug 19, 2014 at 11:52 PM, Muhammad Omer Mughal <span dir="ltr"><<a href="mailto:m.mughal1@postgrad.curtin.edu.au" target="_blank">m.mughal1@postgrad.curtin.edu.au</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" style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">
<p><br>
</p>
<div dir="ltr" style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">
<div><div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
Hi everyone</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
I am running NCL version 6.1.0. I have to write output from multiple netcdf file into the ascii fromat containing wind speed comparison at a particular latitude and longitude.This aim is to get the comparison from the wrf ouput and the instrument at a particular
location. The netcdf files are from the instrument but are not in the regular netcdf format meaning that the time is not included in the dimensions of the variables to be extracted and each file contains one value of time.The first problem that I am facing
is that the fill value is not getting deleted in the ascii output.I use the following to generated the ascii output .<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
b=addfile("/20090711/<a href="http://20090711_000239_01_g.nc" target="_blank">20090711_000239_01_g.nc</a>","r")<br>
lon=b->yp(0,0,:)<br>
lat=b->xp(0,:,0)</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
DATADir = "./"<br>
fils = systemfunc (" ls -1 " + DATADir + "/20090711* ") <br>
numfils = dimsizes(fils) <br>
do j = 0,numfils-1,1<br>
a = addfiles(fils(j)+".nc","r")<br>
time=a[:]->time</div>
</div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
::::::::::::::</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
::::::::::<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
npt = npt + 1 <br>
data(npt) = sprinti("%0.5i", (npt+1) ) <br>
data(npt) = data(npt) + sprintf("%7.1f ",lat(nl))<br>
data(npt) = data(npt) + sprintf("%7.1f ",lon(ml))</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
do nt=0,ntim-1<br>
data(npt) = data(npt) + sprintf("%10.3f ", x(nt))</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
data(npt) = data(npt) + sprintf("%7.1f ",u1(nt,nl,ml))</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
data(npt) = data(npt) + sprintf("%7.1f ",v1(nt,nl,ml))<br>
end do<br>
end do<br>
end do<br>
asciiwrite (fName , data)</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">
In this way I get the output that is attached .<span style="font-family:Tahoma;font-size:13px">Any help to correct this code will be appreciated</span></div>
<p><br>
</p>
<p><br>
</p>
<div>
<div style="font-family:Tahoma;font-size:13px">
<div style="font-family:Tahoma;font-size:13px"><b><br>
</b></div>
<div style="font-family:Tahoma;font-size:13px"><b>Regards</b></div><span><font color="#888888">
<div style="font-family:Tahoma;font-size:13px"><b>Muhammad Omer Mughal</b>
<div><b>Department of Applied Physics</b></div>
<div><br>
</div>
</div>
</font></span></div>
</div>
</div>
</div>
</div>
<br></div></div>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>