iread=0
it=1364
it=120
;
xr_br1=[0.0,1.6]
xr_br2=[1.3,11.1]
xr_br3=[10.0,21.3]
yr_br1=[-1,1]*1.8e-14
yr_br2=[-1,1]*3.2e-14
yr_br3=[-1,1]*1.7e-15