def pwat(self):
self.var = wrf.get_pwat(self.qvapor, self.rho, self.height)
self.var2 = wrf.get_mslp(self.height, self.press, self.temp, self.qvapor)
self.varTitle = "Total Precipitable Water (in)\n" + self.dataSet.getTime()
#Set short variable title for time series
self.sTitle = "Total Precipitable Water (in)"
def mslp(self):
self.var = wrf.get_mslp(self.height,self.press,self.temp,self.qvapor)
self.var2 = self.var
self.varTitle = "Mean Sea Level Pressure (hPa)\n" + self.dataSet.getTime()
#Set short variable title for time series
self.sTitle = "Mean Sea Level Pressure (hPa)"
def temp_850mb(self):
u = self.dataSet.readNCVariable('U')
v = self.dataSet.readNCVariable('V')
u_corr = wrf.unstaggerX(u)
v_corr = wrf.unstaggerY(v)
ref_val = 85000.
#Switched to Cython
#self.u10 = wrf.loglinear_interpolate(u_corr, self.press, ref_val)
#self.v10 = wrf.loglinear_interpolate(v_corr, self.press, ref_val)
self.u10 = np.array(wrf_cython.loglinear_interpolate(u_corr, self.press, ref_val))
self.v10 = np.array(wrf_cython.loglinear_interpolate(v_corr, self.press, ref_val))
#var1 = wrf.loglinear_interpolate(self.temp, self.press, ref_val)
var1 = np.array(wrf_cython.loglinear_interpolate(self.temp, self.press, ref_val))
self.var = wrf.convertT_KtoC(var1)
self.var2 = wrf.get_mslp(self.height, self.press, self.temp, self.qvapor)
self.varTitle = "850-mb Temperature\n" + self.dataSet.getTime()
#Set short variable title for time series
self.sTitle = "850-mb Temperature"
def vort_850mb(self):
u = self.dataSet.readNCVariable('U')
v = self.dataSet.readNCVariable('V')
u_corr = wrf.unstaggerX(u)
v_corr = wrf.unstaggerY(v)
ref_val = 85000.
#Switched to Cython
#self.u10 = wrf.loglinear_interpolate(u_corr, self.press, ref_val)
#self.v10 = wrf.loglinear_interpolate(v_corr, self.press, ref_val)
self.u10 = np.array(wrf_cython.loglinear_interpolate(u_corr, self.press, ref_val))
self.v10 = np.array(wrf_cython.loglinear_interpolate(v_corr, self.press, ref_val))
self.var = wrf.rel_vort(self.u10,self.v10,
self.dataSet.dx[self.dataSet.currentGrid-1],
self.dataSet.dy[self.dataSet.currentGrid-1])
self.var2 = wrf.get_mslp(self.height, self.press, self.temp, self.qvapor)
self.varTitle = "850-mb Relative Voriticity ($10^{-5}$ $s^{-1}$)\n" +\
self.dataSet.getTime()
#Set short variable title for time series
self.sTitle = "850-mb Relative Voriticity ($10^{-5}$ $s^{-1}$)"