Hi Christopher,
I'm not aware of any problems in wrf_helicity that would cause this, but I would like to run your script to see what happens on my system.
Do you get any other errors?
Can you provide the WRF data file on our ftp account?
ftp ftp.cgd.ucar.edu
<log in as "anonymous">
<Use email address as password>
cd incoming
put wrfout_bow44_2006-12-02_12:00:00.nc
quit
Thanks,
--Mary
On Nov 9, 2011, at 6:42 AM, Christopher Steele wrote:
> 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
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Nov 9 08:15:53 2011
This archive was generated by hypermail 2.1.8 : Mon Nov 14 2011 - 10:41:55 MST