<div dir="ltr">Hi, <div><br></div><div>I don't see where you declare "a" to be an array before attempting to use it as such. You'll need to do something like:</div><div><br></div><div> a = new((/p/), float)</div><div><br></div><div>Also, as just a hint, its not necessary to set a(k) = 0 in the outer loop; that could be done more simply by writing</div><div><br></div><div> a = new((/ p /), float)</div><div> a = 0 ; assigns all p-values of the array to zero</div><div><br></div><div>Rick</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Aug 27, 2019 at 3:56 AM Rashed Mahmood 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 dir="ltr"><div>There are obvious errors in the code...e.g.</div><div>1) I do not think that you should use function sizeof, i guess what you need is: N=dimsizes(in), please read the documents for details.<br></div><div>2) a and b were never predefined, so a(0) = 0 and b(0) = 0, would NOT work. You need to define both of them using something like: a = new(,,,)</div><div><br></div><div>Cheers,</div><div>Rashed<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Aug 27, 2019 at 10:54 AM Md. Jalal Uddin via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">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 dir="ltr">Hi all,<div><br></div><div>I made a function in NCL that showed an error like <span style="background-color:rgb(243,243,243)"><b>fatal:Undefined identifier: (a) is undefined, can't continue. The error line highlighted in green color. Any suggestions would be appreciated.</b></span></div><div><span style="background-color:rgb(243,243,243)"><b><br></b></span></div><div><b style="background-color:rgb(255,255,255)">Code:</b></div><div><span style="background-color:rgb(243,243,243)"><b><br></b></span></div><div>undef ("fftwavenumber") <br>function fftwavenumber (in,k1,k2) <br><br>; K1, k2 are high and low-frequency wave numbers. <br>; You can determine bandpass, lowpass, highpass filtering, <br>; and in is the input data time series according to your needs.<br><br>local in,p,a,b,f,k1,k2<br><br>begin <br> N=sizeof(in)<br><br> if mod(N,2).eq.0 then<br> p=N/2 ; N even<br> else<br> p=(N-1)/2 ; N odd<br> end if<br><br>pi=3.1416<br><br>do k=0,p-1<br> <font color="#00ff00"> a(k)=0<br> b(k)=0</font><br> <br> do t=0,N-1<br> a(k)=a(k)+(2/N)*( in(t)*cos(2*pi*k*t/N) )<br> b(k)=b(k)+(2/N)*( in(t)*sin(2*pi*k*t/N) )<br> end do<br>end do<br><br>do t=0,N-1<br> do k=0,p-1<br> f(k,t)=a(k)*cos(2*pi*k*t/N)+b(k)*sin(2*pi*k*t/N)<br> end do<br> return (sum(f(k1:k2,t)))<br>end do<br><br>end <span style="background-color:rgb(243,243,243)"><b><br></b></span></div><div><b><br></b></div><div><b><br clear="all"></b><div>Jalal</div>-- <br><div dir="ltr" class="gmail-m_8019298114968827860gmail-m_6330338382275390821gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font size="2" color="#000000"><b>Md. Jalal Uddin</b></font></div><div dir="ltr">MSc in Applied Meteorology<br>Nanjing University of Information, Science and Technology, China</div><div dir="ltr"><span style="letter-spacing:0.2px">B.Sc. in Disaster Management (Hons.) </span><br></div><div dir="ltr"><span style="letter-spacing:0.2px">Patuakhali Science and Technology University, Bangladesh.</span></div><div dir="ltr">Cell: +8613260859092, <span style="letter-spacing:0.2px">+8801792052662</span></div><div dir="ltr"><font size="2">Web: <a href="http://www.dmjalal90.weebly.com" target="_blank">www.dmjalal90.weebly.com</a> <br>Facebook: jalal.hossen.39 <br>LinkedIn: <a href="https://bd.linkedin.com/in/md-jalal-uddin-80a026b0" target="_blank">https://bd.linkedin.com/in/md-jalal-uddin-80a026b0</a> </font><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></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>
_______________________________________________<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>