[ncl-talk] Missing value in data variable

[Nishtha Agrawal]

Hi,

I have recently started using NCL, here I intend to use
mjo_wavenum_freq_season to obtain the frequency-wavelet spectrum of my
variable. But I am getting the missing value error while accessing the
function.

I am using the following code for my 4D variable 'var'

date=(/20170315, 20170330, 20170415, 20170430, 20170515, 20170530,
20170615, 20170630/)
arr=ndtooned(var)
indi=ind_resolve(ind(ismissing(arr)),dimsizes(var))
ind_size=dimsizes(indi)
npts     = ind_size(0)       ; number of missing elements
ndim     = ind_size(1)       ; rank of var
print(ind_size)
print(npts)
print(ndim)
var at _FillValue = -999.9
do n=0,npts-1
  var(indi(n,0),indi(n,1),indi(n,2),indi(n,3)) = var at _FillValue
end do
wave=mjo_wavenum_freq_season(var(:,0,180,:),date, "summer",False)

And I am getting

(0) mjo_wavenum_freq_season: currently, missing data not allowed:
nMsg=308355

I don't understand what else should I keep in place of -999.9. Please
provide any solution to this problem.

Thanks and Regards,
Nishtha
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190524/65aba2b5/attachment.html>