<div dir="ltr"><div><div><div><div>Dear All,<br><br></div>I am trying to understand how
 to do the wavelet analysis using NCL. I am using ncl v6.4.0 in cygwin/x. I started using an example code and the 
data that are available online. But I got an error message as &quot;segmentation fault
 (core dumped)&quot; when I ran the program. I have given the code that I 
used as below for you to review. Also I have attached the data I used in
 the email. <br><br></div>So, could you please suggest me how I solve this problem? <br><br></div>Thanks in advance.<br><br></div><div>Dipti S.<br></div><div><br>THE CODE is here and the data is attached in the email,<br><br> load &quot;$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_csm.ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/contributed.<wbr>ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/shea_util.ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/shea_util.ncl&quot;<br><br>;-----------------------------<wbr>------------------------------<wbr>---------<br>; CREATE SAME WAVELET FIGURE AS TORRENCE &amp; COMPO using seasonal Nino3<br>; from 1871-1997<br>begin<br>ninoseas = asciiread(&quot;C:\diptic\sst\sst_<wbr>nino3.dat&quot;,-1,&quot;float&quot;)<br>  ninoseas!0 = &quot;time&quot;<br>  ntime = dimsizes(ninoseas)<br>  timeo = fspan(1871.25,1996.,ntime)<br>  print(ntime)<br>  ninoseas&amp;time = timeo<br>  ninomam = dim_avg_Wrap(ninoseas(0::4))<br>  ninojja = dim_avg_Wrap(ninoseas(1::4))<br>  ninoson = dim_avg_Wrap(ninoseas(2::4))<br>  ninodjf = dim_avg_Wrap(ninoseas(3::4))<br>  ninoseas(0::4) = ninoseas(0::4) - ninomam<br>  ninoseas(1::4) = ninoseas(1::4) - ninojja<br>  ninoseas(2::4) = ninoseas(2::4) - ninoson<br>  ninoseas(3::4) = ninoseas(3::4) - ninodjf<br>  time  = timeo<br>  N     = dimsizes(time)<br>;_FillValue=1e+36<br>;printMinMax(ninoseas,False)<br>;return<br>;*****************************<wbr>*******<br>; compute wavelet<br>;*****************************<wbr>*******<br>  mother  = 0<br>  param   = 6.0<br>  dt      = 0.25    ;timesteps in units of years<br>  s0      = dt<br>  dj      = 0.25<br>  jtot    = 1+floattointeger(((log10(N*dt/<wbr>s0))/dj)/log10(2.))<br>  npad    = N<br>  nadof   = 0.0<br>  noise   = 1<br>  siglvl  = 0.05<br>  isigtest= 0<br>print(nadof)<br>w =wavelet(ninoseas,mother,dt,<wbr>param,s0,dj,jtot,npad,noise,<wbr>isigtest,siglvl,nadof)<br>printVarSummary(w)<br>return<br>end<br></div>;;; This is the output I received,<br><br><div>Variable: ntime<br>Type: integer<br>Total Size: 4 bytes<br>            1 values<br>Number of Dimensions: 1<br>Dimensions and sizes:   [1]<br>Coordinates: <br>(0)     504<br>(0)     min=-1.8504   max=2.49921<br><br>Variable: nadof<br>Type: float<br>Total Size: 4 bytes<br>            1 values<br>Number of Dimensions: 1<br>Dimensions and sizes:   [1]<br>Coordinates: <br>(0)      0<br>Segmentation fault (core dumped)<br><br></div><div></div><div><br></div></div>