t2m.en[[i]] <- get.var.ncdf(t2m,paste(“TMP__”,i,”_2maboveground”,sep=”"))

I tested the script with that change and it runs correctly.

Why does the following line from your code return a “variable not found” error?

t2m.en[[i]] <- get.var.ncdf(t2m,"TMP_2maboveground")

When I look at the variable names for the t2m variable (i.e. print(t2m) ) it suggest (I think) that somewhere in the process the variable name has been modified in the follwoing way:

from TMP_2maboveground to TMP__1_2maboveground
from TMP_2maboveground to TMP__2_2maboveground
from TMP_2maboveground to TMP__3_2maboveground
…
from TMP_2maboveground to TMP__20_2maboveground

Will look at the raster doccumetnation but any assistance you could proveivde for this new user would be very much appreciated.

as a sanity check, try making a map of the ensemble mean. it should be obvious what the problem is e.g. the map is inverted.

by the way, i recommend R’s raster package. it gives a nice concise way of handling weather ensembles.

TMP__1_2maboveground
TMP__2_2maboveground
TMP__3_2maboveground
…
…
TMP__19_2maboveground

for each of the ensebme runs?

I tested my wgrib2 statement mannually in the command window and it seems to produced the same prefixed variable names listed above which prevent the following line from working:

t2m.en[[i]] <- get.var.ncdf(t2m,"TMP_2maboveground")

I tried adding and modifying the above line to the followig to dynamicaly create the variable name:

nm=paste("TMP__",i,"_2maboveground",sep="");
t2m.en[[i]] <- get.var.ncdf(t2m,nm)

While this does seem to produce the charts the results dont seem reasonable the histograms produce results that show temps in the -30 range (Kelvin?)

I'm so close but cant seem to crack it, and I have become obsessed with figguring this out. Any assistance would be greatly appreciated.

t2m.en[[i]] <- get.var.ncdf(t2m,"TMP_2maboveground")

…should be changed to:

t2m.en[[i]] <- get.var.ncdf(t2m,"TMP__1_2maboveground")

Is this correct?

One item to point out as well, I had to add a specific reference to the “RColorBrewer” package as well, which your version did not include.

Thanks for making this available. There is some R based processing I’ve been thinking of doing with the GFS ensemble for some time. This gets me started.

Jan

I used wgrib2 executables available here:

ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib2/Windows_XP/

shell(paste(“wgrib2 -s “,t,” | grep \”TMP:2 m\” | wgrib2 -i “,t,” -netcdf “,tmp,sep=”"),intern=T)

how did you get wgrib2 to run in windows? I can find where I might download a windows ready version.

Yes. The graphs were produced running R (2.10 I think) on windows.