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
Received on Fri Jun 05 2009 - 15:07:07 MDT
This archive was generated by hypermail 2.2.0 : Mon Jun 08 2009 - 09:30:31 MDT