;iread=0 igtheta5=0 lam=.8e11 mu0=40. it2=120 it3=15 it1=0 iit=[200,500,1180,4000,13716] tmin=5900. umax=4e-7 bmax=8e-6 pspec_select_noy=1 xr=[1,144]/40.