;$Id: p3.pro,v 1.2 2007/12/31 07:12:47 brandenb Exp $
;
;  plot alpa and eta results from time series
;
pc_read_ts,o=o
;
;  set minimum time after which averaging begins
;
t2=1e30
spawn,'touch parameters.pro'
@parameters
@data/index.pro
@data/testfield_info.dat
default,delmax,.05
default,etamax,.01
default,epsmax,.002
default,iwrite,0
default,run,''
;
;  introduce abbreviations
;
tt=o.t
urms=o.urms
;
;  del tensor
;
alp11=o.alp11
alp21=o.alp21
;
;  eta tensor as in EMF_i = ... -eta_ij J_j
;  so these new etas are also referred to as eta^*
;  writing this furthermore as EMF = delta x J - etat*J,
;  with delta=(0,0,delta), so we have delta x J = (-Jy, Jx, 0)*delta,
;  eta*_12 = -eta*_21 = delta
;
eta12=-o.eta11
eta22=-o.eta21
;
;  read extra fields
;
if itestfield eq 'B11-B22' then begin
  alp12=o.alp12
  alp22=o.alp22
  eta11=+o.eta12
  eta21=+o.eta22
endif
;
;  range of time where to do the analysis
;
tmin=min(tt)
tmax=max(tt)
default,t1,(tmin+tmax)/2.
good=where(tt gt t1 and tt lt t2)
default,kf,5.
;
;  give modified alpmax values in parameters.pro file
;
!p.charsize=2.4
!p.multi=[0,2,4]
;
alp=.5*(alp11+alp22) & eps_alp=.5*(alp11-alp22)
gam=.5*(alp21-alp12) & eps_gam=.5*(alp21+alp12)
del=.5*(eta12-eta21) & eps_del=.5*(eta12+eta21)
eta=.5*(eta11+eta22) & eps_eta=.5*(eta11-eta22)
;
f=alp & plot,tt,f & pc_error_range,tt(good),f(good),/oplot,/accum,mean=alpm,err=alperr
f=gam & plot,tt,f & pc_error_range,tt(good),f(good),/oplot,/accum,mean=gamm,err=gamerr
f=eta & plot,tt,f & pc_error_range,tt(good),f(good),/oplot,/accum,mean=etam,err=etaerr
f=del & plot,tt,f & pc_error_range,tt(good),f(good),/oplot,/accum,mean=delm,err=delerr
;
f=eps_alp & plot,tt,f & pc_error_range,tt(good),f(good),/oplot,/accum,mean=eps_alpm,err=eps_alperr
f=eps_gam & plot,tt,f & pc_error_range,tt(good),f(good),/oplot,/accum,mean=eps_gamm,err=eps_gamerr
f=eps_eta & plot,tt,f & pc_error_range,tt(good),f(good),/oplot,/accum,mean=eps_etam,err=eps_etaerr
f=eps_del & plot,tt,f & pc_error_range,tt(good),f(good),/oplot,/accum,mean=eps_delm,err=eps_delerr
;
fo='(e9.1,8f8.4,3x,8f8.4,3x,a)'
pc_read_param,/param2,obj=param
print,param.etatest,alpm,alperr,gamm,gamerr,etam,etaerr,delm,delerr,eps_alpm,eps_alperr,eps_gamm,eps_gamerr,eps_etam,eps_etaerr,eps_delm,eps_delerr,run,fo=fo
!p.multi=0
END