예제 #1
0
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')
예제 #2
0
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')