Subroutine MFDIH1
c    ===================
c
      Implicit None
c
c Fills the LT1 common block (MFLT1) starting from mfrh.
c
#include "partap.inc"
#include "mfrh.inc"
#include "mfgeox.inc"
#include "mflt1.inc"
c
      integer i,q
      logical first /.true./
      save first
      real dr
      save dr
c
      if(first) then
         first=.false.
         dr=acos(-1.)/180.
      endif
c
      call vzero(nclur,8)
      call vzero(nclup,8)
c
      do i=1,coutab(mfrh)
         call fettab(mfrh,ID,i)
         if(mfrh_itype.eq.1.and.mfrh_plane.eq.1) then
            if(nclup(mfrh_oct).lt.maxclu) then
               nclup(mfrh_oct)=nclup(mfrh_oct)+1
               PhiX(nclup(mfrh_oct),mfrh_oct)=mfrh_val1*dr
               q=(mfrh_oct+1)/2
              DPhiX(nclup(mfrh_oct),mfrh_oct)=
     &       (mfrh_nb2-mfrh_nb1+1)*spsel_ltp*dr
            endif
         elseif(mfrh_itype.eq.2.and.mfrh_plane.eq.1) then
            if(nclur(mfrh_oct).lt.maxclu) then
               nclur(mfrh_oct)=nclur(mfrh_oct)+1
               RhoX(nclur(mfrh_oct),mfrh_oct)=mfrh_val1
               q=(mfrh_oct+1)/2
              DRhoX(nclur(mfrh_oct),mfrh_oct)=
     &       (mfrh_nb2-mfrh_nb1+1)*spsel_ltr
            endif
         endif
      enddo
      end
c