Error: glibc detected *** ncl: double free or corruption (out)

[sahidul@cc3 sahidul]$ ncl max-min.ncl
 Copyright (C) 1995-2009 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 5.1.1
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.
*** glibc detected *** ncl: double free or corruption (out):
0x00002b07deadb010 ***
======= Backtrace: =========
/lib64/libc.so.6[0x34fde71634]
/lib64/libc.so.6(cfree+0x8c)[
0x34fde74c5c]
ncl(NclFree+0x25)[0x5681de]
ncl[0x5495bb]
ncl(_NclDestroyObj+0x56)[0x553e09]
ncl(CallLIST_READ_FILEVAR_OP+0x2900)[0x5e0398]
ncl(_NclExecute+0x442)[0x5e09c5]
ncl(yyparse+0x792)[0x51fc69]
ncl(main+0xa3e)[0x51b4ce]
/lib64/libc.so.6(__libc_start_main+0xf4)[0x34fde1d8b4]
ncl(__gxx_personality_v0+0x14a)[0x51a9fa]
======= Memory map: ========
00400000-01437000 r-xp 00000000 b5:84fa2 9213089
/home/sahidul/ncl/bin/ncl
01537000-01722000 rw-p 01037000 b5:84fa2 9213089
/home/sahidul/ncl/bin/ncl
01722000-0252a000 rw-p 01722000 00:00 0
0a70b000-122c0000 rw-p 0a70b000 00:00 0
34fce00000-34fce1a000 r-xp 00000000 68:02 160307
/lib64/ld-2.5.so
34fd01a000-34fd01b000 r--p 0001a000 68:02 160307
/lib64/ld-2.5.so
34fd01b000-34fd01c000 rw-p 0001b000 68:02 160307
/lib64/ld-2.5.so

Following script, I am trying to run
====================================
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin

      fils = systemfunc("ls -1 /WRF/output-wrf/2010053000/wrfout_d02_*")
      a = addfiles(fils+".nc","r")

      ;temp = a[:]->T ; T = temperature or theta? If theta use wrf_tk to
convert if desired
      wks = gsn_open_wks("ps","test")

      P = a[:]->P ; perturbation
      Pb = a[:]->PB ; base state pressure
      P = P + Pb ; total pressure
      T = a[:]->T ; perturbation potential temperature
(theta+t0)
      ;T = T + 300.

      TK = wrf_tk (P, T)
      temp = wrf_tk (P, T)

      dailymax = temp(0::24,:,:,:) ; copy metadata
      dailymax = dailymax@_FillValue
      dailymin = dailymax

    cntr = 0
    do gg = 0,dimsizes(temp&Time)-1
       dailymax(cntr,:,:,:) = dim_max_n(temp(cntr:cntr+23,:,:,:),0)
       dailymin(cntr,:,:,:) = dim_min_n(temp(cntr:cntr+23,:,:,:),0)
       cntr = cntr+24
    end do
end
====================================

I am having 72 WRF out files which is reading /WRF/output/2010053000
directory, (1 hr interval files)
I would like to calculate Maximum & minimun temperature for a day(24
hourly).

Kindly help me to sort out this error.

Thanking you,

Sahidul

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jun 17 22:05:14 2010