      dimension al(150),be(150),h(150)
      read(5,*)
      read(5,*)
      do 10 j=1,138
        read(5,*) h(j),a,b,rho
        if(h(j).eq.h(j-1)) h(j)=h(j)+0.002
        rs = 6371.0-h(j)
        call emdlv(rs,as,bs)
        al(j) = as
        be(j) = bs
        write(6,fmt='(f12.3,3f12.4)') h(j),as,bs,rho
  10  continue
      stop
      end
