@parameters
default,k0,10
pc_read_ts,obj=ts
plot_io,ts.t,sqrt(ts.divu2m)/(k0*ts.urms),yr=[.1,10]
oplot,ts.t,sqrt(ts.o2m)/(k0*ts.urms),col=122
;
pc_read_param,/param2,obj=param
nu=param.nu
Re=ts.urms/(k0*nu)
plot,ts.t,Re
print,max(Re)
END