def runPlot_theta(): #(x4_t, x4_kf, x4_cfsr) = dataOutput() (x1kf, x1t, x4kf, x4t, cfsr) = dataOutputJuly() dataPlot = [None]*4 dataPlot[0] = x1kf-cfsr dataPlot[1] = x1t-cfsr dataPlot[2] = x4kf-cfsr dataPlot[3] = x4t-cfsr compare.plot_save_lists(dataPlot, 4, ['x1_kf','x1_t', 'x4_kf', 'x4_t'], ' ', 'dtheta_2007-24-12')
def runPlot_ll():
#ll_x4t, ll_x4kf, ll_x4cfsr = latLonAugust()
(ll_x1t,ll_x1kf,ll_x4t, ll_x4kf, ll_cfsr) = latLonJuly()
#print ll_x4t
#print ll_x4kf
#print ll_x4cfsr
nVals = ll_x4t.shape[0]
dVal = np.empty((4,nVals))
for i in xrange(nVals):
dVal[0,i] = conn.calc_distSphere(6371000., ll_x1t[i,:], ll_cfsr[i,:])
dVal[1,i] = conn.calc_distSphere(6371000., ll_x1kf[i,:], ll_cfsr[i,:])
dVal[2,i] = conn.calc_distSphere(6371000., ll_x4t[i,:], ll_cfsr[i,:])
dVal[3,i] = conn.calc_distSphere(6371000., ll_x4kf[i,:], ll_cfsr[i,:])
dataPlot = [None]*4
for i in xrange(4):
dataPlot[i] = dVal[i,:]/1000.
compare.plot_save_lists(dataPlot, 4, ['x1_t', 'x1_kf','x4_t', 'x4_kf'], ' ', 'dloc_2006-07-24_12')