wrf_helicity, segmentation fault

Hi everyone,

Whenever I try to use the function wrf_helicity I get a segmentation fault
and I'm not exactly sure why. I'm using ncl version 6, but have also tried
the earlier 5.2.1. Below is the script to see if I have overlooked anything.

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;
;
; Script to calculate storm relative helicity from WRF output
data ;
;
;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

begin

;load in WRF file
wrffile =
addfile("/gpfs/env/nhm06sfu/WRFattempt3/WRFV3/test/em_real/bow_echo/wrfout_bow44_2006-12-02_12:00:
00.nc","r")
ua = wrf_user_getvar(wrffile,"ua",15)
va = wrf_user_getvar(wrffile,"va",15)
gpt = wrf_user_getvar(wrffile,"geopt",15)
ter = wrf_user_getvar(wrffile,"ter",15)

;define workstn, colourmap and resources
wks = gsn_open_wks("x11","plt_helicity03")
gsn_define_colormap(wks,"hotcolr_19lev")
pltres = True
mpres = True
res = True
mpres@mpDataBaseVersion = "highres"
mpres@mpGeophysicalLineThicknessF = 2
mpres@mpGeophysicalLineColor = "black"

;calculate helicity
helicity = wrf_helicity(ua,va,gpt,ter,3000)

;create plot
contour = wrf_contour(wrffile,wks,helicity,res)
plot = wrf_map_overlays(wrffile,wks,contour,pltres,mpres)

end

cheers

Chris

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Nov 9 06:42:33 2011