Hi there,
I am using NCL to read a couple of satellite data files (HDF format) and to write into NETCDF format.The original data is monthly and has a dimension of 360x720x9x4x6. I read 48 files, and store them as vart(48,360,720,9,4,6). Till now, the reading process is OK. When I try to get a segment of "vart" (i.e. var = vart(:,:,:,0,:,0)), I got a "Segmentation fault (core dumped)". But, if I delete this line (i.e. var = vart(:,:,:,0,:,0)), the script can go on. Does someone know why is that? Thanks!
P.S. My platform is 64-bit Linux, with 32GB memory.
Part of my ncl script is as follow:
=====================================
......
vname = "data_name"
;;;; Read the varialbe ;;;;
;;;; switch to "join" mode ;;;;
ListSetType (f, "join")
vart = addfiles_GetVar (f, fils, vname) ; fils is a string array with filenames
printVarSummary(vart) ; The information of vart can be printed successfully.
var = vart(:,:,:,0,:,0) * 10.0 ; If this line is removed, the script can go on ...
printVarSummary(var) ; Information about var is NOT printed.
......
=====================================
Best,
-Eddy
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Dec 17 16:07:41 2009
This archive was generated by hypermail 2.1.8 : Thu Dec 17 2009 - 17:15:52 MST