[ncl-talk] segmentation fault caused by large variable

[Cristiana Stan]

Hi,

I would like to compute daily anomalies of a large 4D variable that exceeds the NCL size limit and causes segmentation fault. I’m looking for a suggestion to work around this issue.

My code looks like:

 x       =  f->$var$(iStrt:iLast,:,:,:)

ntim = 12419
nlev = 22
nlat = 145
mlon = 512


xClmDay = clmDayTLLL((/x/), yyyyddd)


xAnom = new((/ntim,nlev,nlat,mlon/),float) ; This line causes the segmentation fault.
 do n = 0,nlev-1
    xAnom(:,n,:,:)  = calcDayAnomTLL (x(:,n,:,:), yyyyddd, xClmDay(:,n,:,:))
  end do

;***********************************************************
; Create netCDF using 'simple' method
;***********************************************************


Thanks,

Cristiana



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170125/5f1b36cf/attachment.html