Posting my issue here as well, because I didn't get any help from the Google Group. I've got a bit of an oddball WRF dataset that isn't organized like a normal wrfout file. Instead of all of the variables being in one single netCDF file for some time x, all of the variables are split into their own netCDF files. This means that anything that calls for a singular input wrf file doesn't work, so I have to use the raw computational algorithms in wrf python; getvar won't work.
That said, I've used wrf.rh just fine with this dataset in the past. However, when I load in files from the exact same dataset, I can't get wrf.slp to work. It needs geopotential height, tk, total pressure, and qv as inputs. I have given it all 4 of these, and I double checked the dimensions of all of them (they're all (232, 282) numpy arrays). The units check out as well.
The code (not much to see):
" slpctrl = wrf.slp(zctrl,tkctrl,pctrl,qvaporctrl) "
The error I get is:
" ValueError: invalid shape for argument 0 (got 232, 282), expected (282,)) "
Not sure what it's trying to tell me, other than it doesn't like the geopotential height input, but I can't figure out why it won't accept it. It's the same as the others. Is this a bug perhaps? In the meantime I suppose I could try NCL or just calculate it manually, but it would be nice to save the time.
Regardless, any help would be appreciated.
I'm using python 3.6 miniconda on CentOS 6.8. I can attach code or elaborate as needed if it'll help.
Hi there,
Posting my issue here as well, because I didn't get any help from the Google Group. I've got a bit of an oddball WRF dataset that isn't organized like a normal wrfout file. Instead of all of the variables being in one single netCDF file for some time x, all of the variables are split into their own netCDF files. This means that anything that calls for a singular input wrf file doesn't work, so I have to use the raw computational algorithms in wrf python; getvar won't work.
That said, I've used wrf.rh just fine with this dataset in the past. However, when I load in files from the exact same dataset, I can't get wrf.slp to work. It needs geopotential height, tk, total pressure, and qv as inputs. I have given it all 4 of these, and I double checked the dimensions of all of them (they're all (232, 282) numpy arrays). The units check out as well.
The code (not much to see): " slpctrl = wrf.slp(zctrl,tkctrl,pctrl,qvaporctrl) "
The error I get is: " ValueError: invalid shape for argument 0 (got 232, 282), expected (282,)) "
Not sure what it's trying to tell me, other than it doesn't like the geopotential height input, but I can't figure out why it won't accept it. It's the same as the others. Is this a bug perhaps? In the meantime I suppose I could try NCL or just calculate it manually, but it would be nice to save the time.
Regardless, any help would be appreciated. I'm using python 3.6 miniconda on CentOS 6.8. I can attach code or elaborate as needed if it'll help.
Thank you!