%     log space r
      rw=0.4;
      re=1128.38
      red=re/rw;
      x=logspace(0,log10(red),20);
      x=x'
      save rd.dat x -ascii
