;$Id: panal.pro,v 1.1 2011/06/27 21:54:59 brandenb Exp $ ; ; procedure for computing alpha and etat in units ; of the expected values alp0 and etat0. ; St=.07 k1=1./8. orms=mean(ts.orms) urms=mean(ts.urms) oum=mean(ts.oum) kf=orms/urms epsf=oum/urms^2/kf tau=St/(urms*kf) alp=.5*mean(ts.alp11+ts.alp22) etat=.5*mean(ts.eta11+ts.eta22) alp0=-tau*oum/3. etat0=tau*urms^2/3. print,kf/k1,epsf,tau print,alp/alp0,etat/etat0 ; ; plot as function of time ; plot,ts.t,ts.orms/ts.urms,yr=[0,.01] END