lccor.py 文件源码

def lccor(rc,bs=0,fs=1,step=1,kind='int'):
   import numpy as np
   from AnalysisFunctions import fcorr 
   ie=1/np.e
   rc.vars2load(['bx','by','bz'])
   tt  = np.zeros((fs-bs)/step)
   lxc = np.zeros((fs-bs)/step)
   lyc = np.zeros((fs-bs)/step)
   lc  = np.zeros((fs-bs)/step)
   for i in range(bs,fs,step):
      print i; idx = (i-bs)/step
      rc.loadslice(i); 
      tt[idx] = rc.time
      rx,bxcor=fcorr(rc.bx,rc.bx,ax=0,dx=rc.dx)
      ry,bycor=fcorr(rc.by,rc.by,ax=1,dx=rc.dy)
      if kind == "ie":
         lxc[idx]=rx[abs(bxcor-ie).argmin()]
         lyc[idx]=ry[abs(bycor-ie).argmin()]
      elif kind == "int":
         lxc[idx]=np.sum(bxcor)*rc.dx
         lyc[idx]=np.sum(bycor)*rc.dy
      lc[idx] = 0.5*(lxc[idx]+lyc[idx])
      print tt[idx],lxc[idx],lyc[idx],lc[idx]
   return tt,lxc,lyc,lc