Re: Segmentation Fault

Hi Ibraheem,
Can you send us the complete error message? And if the line number isn't
specified in the error message, can you tell us exactly what line the
error message is occurring? (You can find that out via print statements)..
Adam

Ibraheem Khan wrote:
> Hi,
>
> I am using an NCL script to calculate the zonal mean and precipitation anomaly. The number of "*.nc" files are 30. It crashes with an error "segmentation fault". Please help me with to rectify this error.
> The script is given below:
>
>
> ;************************************************
> 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/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
> ;************************************************
> begin
> ;************************************************
> ; variable and file handling
> ;************************************************
> fils = systemfunc("ls 19*.nc") ; define filename
> f = addfiles(fils,"r") ; open netcdf file
> p = f[:]->precip ; select variable to ave
> ;************************************************
> ; calculate deviation from zonal mean
> ;************************************************
> anom = p ; trick to copy cv's and atts
> anom = dim_rmvmean(p)
> ;************************************************
> ; plot parameters
> ;************************************************
> wks = gsn_open_wks("ps","dev_zonal") ; Opens a ps file
> plot = new(2,graphic) ; create graphical array
> gsn_define_colormap(wks,"radar")
> res = True
> ;************************************************
> ; original data
> ;************************************************
> res_at_gsnDraw = False ; don't draw yet
> res_at_gsnFrame = False ; don't advance frame yet
> res_at_gsnPaperOrientation = "portrait"
> res_at_gsnMaximize = True
>
> res_at_lbLabelBarOn = False
> res_at_cnFillOn = True ; color plot
> res_at_gsnSpreadColors = True ; subset color map
> res_at_cnLinesOn = False ; no cn lines
> res_at_cnLevelSelectionMode = "ManualLevels" ; set manual contour level
>
> res_at_cnMinLevelValF = -10. ; set min contour level
> res_at_cnMaxLevelValF = 50 ; set max contour level
> res_at_cnLevelSpacingF = 5 ; set contour spacing
> res_at_gsnSpreadColors = True ; full colors
>
> res_at_gsnAddCyclic = False
> res_at_gsnZonalMean = True ; add zonal plot
> res_at_mpFillOn = False ; no grey continent
> res_at_cnInfoLabelOrthogonalPosF = -0.17 ; move info label up
> res_at_cnLineLabelFontHeightF = .012 ; increase font size
> res_at_gsnLeftString = " " ; left string name
> res_at_gsnRightString = "mm" ; units
> res_at_gsnCenterString = "Original Data"
>
> ; panel expects plots to be of the same size. Since these two plots are
> ; very different, we are forcing them to be the same size by setting the
> ; width.
>
> res_at_vpWidthF = 0.7
> res_at_vpYF = 0.5
> plot(0) = gsn_csm_contour_map(wks,p(0,:,:),res)
>
> zres=True
> zres_at_amOrthogonalPosF = 0.2
> zonal_id=gsn_csm_attach_zonal_means(wks,plot(0),p(0,:,:),zres)
> ;delete(res_at_gsnZonalMean) ; delete zonal plot
> ;************************************************
> ; anomaly data
> ;************************************************
> gsn_define_colormap(wks,"radar")
>
> res_at_cnFillOn = True ; color plot
> res_at_gsnSpreadColors = True ; subset color map
> res_at_gsnAddCyclic = False
> plot(1) = gsn_csm_contour_map_ce(wks,anom(0,:,:),res)
>
> ; panel the two plots together
> resP = True ; modify the panel plot
>
> resP_at_txPosYF = 0.88
> resP_at_txString = "Precipitation Climatology (1980-2002)"
>
> resP_at_gsnPanelYF = (/0.8 , 0.5/)
>
> resP_at_gsnPanelDebug = True ; Maybe helpful
> resP_at_gsnPanelYWhiteSpacePercent = 10 ; it can make plots closer or farther
> ; between each other.
>
>
> resP_at_lbLabelFontHeightF = 0.005
> resP_at_gsnPanelLabelBar = True ; add common colorbar
> resP_at_gsnPanelBottom = 0.03
> gsn_panel(wks,plot,(/2,1/),resP)
> end
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
--------------------------------------------------------------
Adam Phillips			             asphilli_at_ucar.edu
National Center for Atmospheric Research   tel: (303) 497-1726
ESSL/CGD/CAS                               fax: (303) 497-1333
P.O. Box 3000				
Boulder, CO 80307-3000	  http://www.cgd.ucar.edu/cas/asphilli
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk