The focus point is a variable set in the pywrfplotParams.py:
#This defines the grid box for which skewT-plots etc are based lat_focuspoint = 60.2 lon_focuspoint = 11.08
Unfortunately, I don't know what focus point is referring to, so when I use the values set by the author, obviously I get an error, and that error is - IndexError: index out of bounds.
I have found this link useful in providing the definition for the attributes asked in that question, but I'm not entirely sure if I should also get the cell indices from latitude and longitude in the WRF model grids.
Your thoughts are very much appreciated!
def getXY(lon,lat): x_nr = 0 y_nr = 0 while (lon[x_nr] < lon_focuspoint): x_nr += 1 while (lat[y_nr] < lat_focuspoint): y_nr += 1
finds the closest grid point west and south of your focus point by setting the x and y gridpoints to 0 and incrementing the index as long as the latitude/longitude values at that grid point are less than the focus point. This is what the focus point is for -- the lat/lon location you want the skew-t/log-p plot at.
The exception IndexError: index out of bounds.
means that the above code tried to access a value lon[x_nr]
or lat[y_nr]
that was out of bounds. This means that the latitude/longitude coordinate of your focus point is not within your domain.
To run the code successfully, edit the focus point to be a point within your model domain, specifically the value of latitude and longitude that you want the skew-t profile for.