Re: segmentation fault

Hi Ahmad,
I would put in a few print statements to try to figure out what is going
on. For instance, I would definitely put in a print(itime) statement. I
would take a look at your YMDH array. I would also try reading in var129
without switching it to float or subsetting it (via latitude) to see if
that works.

Try the above and if you are still stuck please respond to the ncl-talk
email list with additional information from the print statements.
Adam

On 03/31/2014 12:55 PM, Ahmad Farsyud wrote:
> Hallo NCL folks, would you please to help me, why did I get this error:
>
> ncl 0> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> ncl 1> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> ncl 2> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> ncl 3> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl"
>
> ncl 6>
> ncl 7> ;begin
> ncl 8>
> ncl 9> f = addfile ("ERA_daily.nc" , "r")
> ncl 10>
> ncl 11> TIME = f->time
> ncl 12> YMDH = cd_calendar(TIME, 0)
> ncl 13> itime = ind(YMDH(:,1).gt.1 .and. YMDH(:,1).le.3)
> ncl 14>
> ncl 15> data1 = dble2flt(f->var129(itime,:,{10:90},:))
>
> *Segmentation fault*
> *
> *
> *
> *
> *My input data:*
>
> dimensions:
> lon = 144 ;
> lat = 73 ;
> lev = 23 ;
> time = UNLIMITED ; // (20270 currently)
> variables:
> double lon(lon) ;
> lon:long_name = "longitude" ;
> lon:units = "degrees_east" ;
> lon:standard_name = "longitude" ;
> lon:axis = "X" ;
> double lat(lat) ;
> lat:long_name = "latitude" ;
> lat:units = "degrees_north" ;
> lat:standard_name = "latitude" ;
> lat:axis = "Y" ;
> double lev(lev) ;
> lev:long_name = "pressure" ;
> lev:units = "Pa" ;
> lev:axis = "Z" ;
> double time(time) ;
> time:units = "hours since 1958-01-01 00:00:00" ;
> *time:calendar = "proleptic_gregorian" ;*
> double var129(time, lev, lat, lon) ;
> var129:table = 128 ;
> *
> *
> Thank you
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
______________________________________________________________
Adam Phillips                                asphilli@ucar.edu
NCAR/Climate and Global Dynamics Division       (303) 497-1726
P.O. Box 3000				
Boulder, CO 80307-3000    http://www.cgd.ucar.edu/staff/asphilli

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Apr 1 12:02:13 2014