[ncl-talk] Function in NCL

[Md. Jalal Uddin]

Thank you. It works now.

On Tue, 27 Aug 2019 at 21:24, Rick Brownrigg <brownrig at ucar.edu> wrote:

> Hi,
>
> 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:
>
>   a = new((/p/), float)
>
> 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
>
>   a = new((/ p /), float)
>   a = 0    ; assigns all p-values of the array to zero
>
> Rick
>
>
> On Tue, Aug 27, 2019 at 3:56 AM Rashed Mahmood via ncl-talk <
> ncl-talk at ucar.edu> wrote:
>
>> There are obvious errors in the code...e.g.
>> 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.
>> 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(,,,)
>>
>> Cheers,
>> Rashed
>>
>> On Tue, Aug 27, 2019 at 10:54 AM Md. Jalal Uddin via ncl-talk <
>> ncl-talk at ucar.edu> wrote:
>>
>>> Hi all,
>>>
>>> I made a function in NCL that showed an error like *fatal:Undefined
>>> identifier: (a) is undefined, can't continue.  The error line highlighted
>>> in green color. Any suggestions would be appreciated.*
>>>
>>> *Code:*
>>>
>>> undef ("fftwavenumber")
>>> function fftwavenumber (in,k1,k2)
>>>
>>> ; K1, k2 are high and low-frequency wave numbers.
>>> ; You can determine bandpass, lowpass, highpass filtering,
>>> ; and in is the input data time series according to your needs.
>>>
>>> local in,p,a,b,f,k1,k2
>>>
>>> begin
>>>     N=sizeof(in)
>>>
>>>    if mod(N,2).eq.0 then
>>>          p=N/2     ; N even
>>>    else
>>>          p=(N-1)/2 ; N odd
>>>    end if
>>>
>>> pi=3.1416
>>>
>>> do k=0,p-1
>>>       a(k)=0
>>>       b(k)=0
>>>
>>>   do t=0,N-1
>>>       a(k)=a(k)+(2/N)*( in(t)*cos(2*pi*k*t/N) )
>>>       b(k)=b(k)+(2/N)*( in(t)*sin(2*pi*k*t/N) )
>>>   end do
>>> end do
>>>
>>> do t=0,N-1
>>>    do k=0,p-1
>>>         f(k,t)=a(k)*cos(2*pi*k*t/N)+b(k)*sin(2*pi*k*t/N)
>>>    end do
>>>    return (sum(f(k1:k2,t)))
>>> end do
>>>
>>> end
>>>
>>>
>>> Jalal
>>> --
>>> *Md. Jalal Uddin*
>>> MSc in Applied Meteorology
>>> Nanjing University of Information, Science and Technology, China
>>> B.Sc. in Disaster Management (Hons.)
>>> Patuakhali Science and Technology University, Bangladesh.
>>> Cell: +8613260859092, +8801792052662
>>> Web: www.dmjalal90.weebly.com
>>> Facebook: jalal.hossen.39
>>> LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>

-- 
*Md. Jalal Uddin*
MSc in Applied Meteorology
Nanjing University of Information, Science and Technology, China
B.Sc. in Disaster Management (Hons.)
Patuakhali Science and Technology University, Bangladesh.
Cell: +8613260859092, +8801792052662
Web: www.dmjalal90.weebly.com
Facebook: jalal.hossen.39
LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190828/0d57fcdb/attachment.html>