<div dir="ltr"><div>No progress this side !<br></div><div><br></div><div>Included the suggested changes:<br><br>  nlat = dimsizes(xlat)<br>  nlon = dimsizes(xlon)<br><br>;;*************************<br>  do it = 0,ntimes-1               ; TIME LOOP<br><br>  td = wrf_user_getvar(a,&quot;td&quot;,it)<br>  p  = wrf_user_getvar(a,&quot;pressure&quot;,it)<br><br>  td_plane = new((/nlevels,nlat,nlon/),typeof(td))<br><br>  do level = 0,nlevels-1                 ; LOOP OVER LEVELS<br>  pressure = plevs(level)<br><br>  td_plane(level,:,:) = wrf_user_intrp3d(td,p,&quot;h&quot;,pressure,0.,False)<br><br> end do ; end of level loop<br> end do ; end of time loop<br> end do ; end of files loop<br> end<br><br></div><div class="gmail_extra">Error: <br>fatal:_NclBuildArray: each element of a literal array must have the same dimension sizes, at least one item doesn&#39;t<br>^Mfatal:[&quot;Execute.c&quot;:8128]:Execute: Error occurred at or near line 49 in file wrf_td.ncl<br><br></div><div class="gmail_extra">Thanks,<br></div><div class="gmail_extra">Ancy<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Sep 10, 2015 at 5:47 PM, Alessandra Giannini <span dir="ltr">&lt;<a href="mailto:alesall@iri.columbia.edu" target="_blank">alesall@iri.columbia.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
All right… Making progress!<br>
<br>
What you need to do is define the variable td_plane BEFORE the loop, using the new command:<br>
<br>
&lt;<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/new.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/new.shtml</a>&gt;<br>
<br>
Something like:<br>
<br>
td_plane = new((/nlevels,nlat, nlon/),typeof(td))<br>
<div class=""><div class="h5"><br>
<br>
<br>
<br>
<br>
<br>
<br>
On Sep 10, 2015, at 8:02 AM, Ancy Thomas &lt;<a href="mailto:ancyt.thomas@gmail.com">ancyt.thomas@gmail.com</a>&gt; wrote:<br>
<br>
&gt; Hello Madam !<br>
&gt;<br>
&gt; Done the changes as you suggested :<br>
&gt;<br>
&gt;  do it = 0,ntimes-1               ; TIME LOOP<br>
&gt;   td = wrf_user_getvar(a,&quot;td&quot;,it)<br>
&gt;   p  = wrf_user_getvar(a,&quot;pressure&quot;,it)<br>
&gt;<br>
&gt;   do level = 0,nlevels-1                 ; LOOP OVER LEVELS<br>
&gt;   pressure = plevs(level)<br>
&gt;<br>
&gt;   td_plane(level,:,:) = wrf_user_intrp3d(td,p,&quot;h&quot;,pressure,0.,False)<br>
&gt;<br>
&gt;  end do ; end of level loop<br>
&gt;  end do ; end of time loop<br>
&gt;  end do ; end of files loop<br>
&gt;  end<br>
&gt;<br>
&gt; Still I am getting the same error:<br>
&gt; fatal:Undefined identifier: (td_plane) is undefined, can&#39;t continue<br>
&gt; ^Mfatal:[&quot;Execute.c&quot;:8128]:Execute: Error occurred at or near line 49 in file wrf_td.ncl<br>
&gt;<br>
&gt;<br>
&gt; Thanks,<br>
&gt; regards,<br>
&gt; Ancy<br>
&gt;<br>
&gt;<br>
&gt; On Thu, Sep 10, 2015 at 5:10 PM, Alessandra Giannini &lt;<a href="mailto:alesall@iri.columbia.edu">alesall@iri.columbia.edu</a>&gt; wrote:<br>
&gt;<br>
&gt; Hi there,<br>
&gt;<br>
&gt; [it&#39;s a madam speaking, actually!]<br>
&gt; You should always reply to ncl-talk, not just to the person who responds…<br>
&gt;<br>
&gt; It would help if you explained what is not working…<br>
&gt; At any rate, what is the structure of the variable you call here:<br>
&gt;<br>
&gt; &gt;   td = wrf_user_getvar(a,&quot;td&quot;,it)<br>
&gt;<br>
&gt;<br>
&gt; if it has all the levels you need, just write it out.<br>
&gt;<br>
&gt; It may be that these two lines that define your variables need to sit outside/before the loops:<br>
&gt;<br>
&gt; &gt;   td = wrf_user_getvar(a,&quot;td&quot;,it)<br>
&gt; &gt;   p  = wrf_user_getvar(a,&quot;pressure&quot;,it)<br>
&gt;<br>
&gt;<br>
&gt; Otherwise you redefine them each time.<br>
&gt; Please think through what you want to accomplish, and the correct order of the steps to get there.<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; --<br>
&gt; Alessandra Giannini<br>
&gt; IRI for Climate and Society - The Earth Institute at Columbia University<br>
&gt; P.O. Box 1000, Palisades NY 10964-8000<br>
&gt; phone/fax: +1 845 680-4473/4864 - email: <a href="mailto:alesall@iri.columbia.edu">alesall@iri.columbia.edu</a><br>
&gt;<br>
&gt; On Sep 10, 2015, at 6:54 AM, Ancy Thomas &lt;<a href="mailto:ancyt.thomas@gmail.com">ancyt.thomas@gmail.com</a>&gt; wrote:<br>
&gt;<br>
&gt; &gt; Hello Sir,<br>
&gt; &gt;<br>
&gt; &gt; Thanks for the response.<br>
&gt; &gt;<br>
&gt; &gt; Yes, I want to write the data variable for each level in a loop. The program is given below: How should I modify .<br>
&gt; &gt;<br>
&gt; &gt; ************************************************<br>
&gt; &gt; begin<br>
&gt; &gt; datadir = systemfunc(&quot;tail -1 fil&quot;)<br>
&gt; &gt; FILES = systemfunc (&quot;ls -1 &quot; + datadir + &quot;wrfout_d02* &quot; )<br>
&gt; &gt; numFILES = dimsizes(FILES)<br>
&gt; &gt; ;;***************Pressure levels **********<br>
&gt; &gt;  plevs = (/ 925., 850., 700., 600., 500., 300., 200./)   ; pressure levels to plot<br>
&gt; &gt;  nlevels         = dimsizes(plevs)     ; number of pressure levels<br>
&gt; &gt;<br>
&gt; &gt; do ifil = 0,numFILES-1<br>
&gt; &gt;    a = addfile(FILES(ifil)+&quot;.nc&quot;,&quot;r&quot;)<br>
&gt; &gt;<br>
&gt; &gt; ; What times and how many time steps are in the data set<br>
&gt; &gt; times  = wrf_user_list_times(a)  ; get times in the file<br>
&gt; &gt;   ntimes = dimsizes(times)       ; number of times in the file<br>
&gt; &gt; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br>
&gt; &gt;<br>
&gt; &gt; ;Read in constant variables<br>
&gt; &gt;       xlat  = wrf_user_getvar(a, &quot;XLAT&quot;,0)<br>
&gt; &gt;       xlon  = wrf_user_getvar(a, &quot;XLONG&quot;,0)<br>
&gt; &gt;       hght   = wrf_user_getvar(a, &quot;HGT&quot;,0)<br>
&gt; &gt;<br>
&gt; &gt; ; writing data into netcdf format<br>
&gt; &gt;     system(&quot;rm -f <a href="http://wrfout.nc" rel="noreferrer" target="_blank">wrfout.nc</a>&quot;)<br>
&gt; &gt;     fw = addfile(&quot;<a href="http://wrfout.nc" rel="noreferrer" target="_blank">wrfout.nc</a>&quot;, &quot;c&quot;)<br>
&gt; &gt;<br>
&gt; &gt; ; writing the dimensions time, level, lat, lon<br>
&gt; &gt;   fw-&gt;xlat = xlat<br>
&gt; &gt;   fw-&gt;xlon = xlon<br>
&gt; &gt;<br>
&gt; &gt; ;;*************************<br>
&gt; &gt;   do it = 0,ntimes-1               ; TIME LOOP<br>
&gt; &gt;<br>
&gt; &gt;   do level = 0,nlevels-1                 ; LOOP OVER LEVELS<br>
&gt; &gt;<br>
&gt; &gt;   pressure = plevs(level)<br>
&gt; &gt;<br>
&gt; &gt;   td = wrf_user_getvar(a,&quot;td&quot;,it)<br>
&gt; &gt;   p  = wrf_user_getvar(a,&quot;pressure&quot;,it)<br>
&gt; &gt;<br>
&gt; &gt;   td_plane(level,:,:) = wrf_user_intrp3d(td,p,&quot;h&quot;,pressure,0.,False)<br>
&gt; &gt;<br>
&gt; &gt;  end do ; end of level loop<br>
&gt; &gt;  end do ; end of time loop<br>
&gt; &gt;  end do ; end of files loop<br>
&gt; &gt;  end<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; Thanks,<br>
&gt; &gt; Best regards,<br>
&gt; &gt; Ancy<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; On Thu, Sep 10, 2015 at 4:04 PM, Alessandra Giannini &lt;<a href="mailto:alesall@iri.columbia.edu">alesall@iri.columbia.edu</a>&gt; wrote:<br>
&gt; &gt;<br>
&gt; &gt; Hi there,<br>
&gt; &gt;<br>
&gt; &gt; I am not familiar with the &quot;wrf&quot; commands you use, but you are asking it to overwrite levels with this line:<br>
&gt; &gt;<br>
&gt; &gt; &gt; fw-&gt;u_plane = u_plane<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; Do you need to write the levels out separately, meaning do you need a do-loop to begin with?<br>
&gt; &gt; If you have all the data in a variable, it may be preferable to select the time and levels you want to write out with coordinate selection, and write out all in one go.<br>
&gt; &gt;<br>
&gt; &gt; warm regards<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; --<br>
&gt; &gt; Alessandra Giannini<br>
&gt; &gt; IRI for Climate and Society - The Earth Institute at Columbia University<br>
&gt; &gt; P.O. Box 1000, Palisades NY 10964-8000<br>
&gt; &gt; phone/fax: +1 845 680-4473/4864 - email: <a href="mailto:alesall@iri.columbia.edu">alesall@iri.columbia.edu</a><br>
&gt; &gt;<br>
&gt; &gt; On Sep 10, 2015, at 12:26 AM, Ancy Thomas &lt;<a href="mailto:ancyt.thomas@gmail.com">ancyt.thomas@gmail.com</a>&gt; wrote:<br>
&gt; &gt;<br>
&gt; &gt; &gt; Could anyone please help for this query.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; On Wed, Sep 9, 2015 at 6:04 PM, Ancy Thomas &lt;<a href="mailto:ancyt.thomas@gmail.com">ancyt.thomas@gmail.com</a>&gt; wrote:<br>
&gt; &gt; &gt; Hello all,<br>
&gt; &gt; &gt; Greetings!<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; I want to write a variable in do loop. But it is writing the data only for the last level(200hpa), ie; data is getting over written for each level. How to save the u_plane data as array?<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; ;;***************Pressure levels **********<br>
&gt; &gt; &gt;  plevs = (/ 925., 850., 700., 600., 500., 300., 200./)   ; pressure levels to plot<br>
&gt; &gt; &gt;  nlevels         = dimsizes(plevs)     ; number of pressure levels<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; ; writing data into netcdf format<br>
&gt; &gt; &gt;     system(&quot;rm -f <a href="http://wrfout.nc" rel="noreferrer" target="_blank">wrfout.nc</a>&quot;)<br>
&gt; &gt; &gt;     fw = addfile(&quot;<a href="http://wrfout.nc" rel="noreferrer" target="_blank">wrfout.nc</a>&quot;, &quot;c&quot;)<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; ; writing the dimensions time, level, lat, lon<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; ;  fw-&gt;xlev = xlev<br>
&gt; &gt; &gt;   fw-&gt;xlat = xlat<br>
&gt; &gt; &gt;   fw-&gt;xlon = xlon<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;  do it = 0,ntimes-1               ; TIME LOOP<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;  do level = 0,nlevels-1                 ; LOOP OVER LEVELS<br>
&gt; &gt; &gt; u    = wrf_user_getvar(a,&quot;U&quot;,it)<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; u_plane    = wrf_user_intrp3d(u,p,&quot;h&quot;,pressure,0.,False)<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; fw-&gt;u_plane = u_plane<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; end do ; level loop<br>
&gt; &gt; &gt; end do ; time loop<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Kindly help to resolve the issue.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Thank you,<br>
&gt; &gt; &gt; With regards,<br>
&gt; &gt; &gt; --<br>
&gt; &gt; &gt; Ancy Thomas,<br>
&gt; &gt; &gt; Computational Earth Science Group / HPC application Group<br>
&gt; &gt; &gt; CDAC<br>
&gt; &gt; &gt; Pune<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; --<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; _______________________________________________<br>
&gt; &gt; &gt; ncl-talk mailing list<br>
&gt; &gt; &gt; <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
&gt; &gt; &gt; List instructions, subscriber options, unsubscribe:<br>
&gt; &gt; &gt; <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
<br>
</div></div></blockquote></div><br><div class="gmail_signature"><div dir="ltr"><br></div></div>
</div></div>