***********************************************************************

      SUBROUTINE UFLDINT(MAPFILE)

*      READ TOSCA FIELD MAP FROM EXTERNAL FILE
*
*      created 12/12/98 R.Schicker
*
***********************************************************************

       implicit none

#include "fieldt.inc"

       integer iz,ir,ip,j1,j2,j3
       character*(*) mapfile 
       character mode*3 / 'UNF' /   ! data format: unformatted

       if(mode.eq.'FOR')then
          open(unit=31,file=mapfile,form='formatted',status='old')
       elseif(mode.eq.'UNF')then
          open(unit=31,file=mapfile,form='unformatted',status='old')
       endif

       write(6,99)mapfile
 99    format(' mapfile = ' ,a32)

 	if(mode.eq.'FOR')then

        read(31,*)iz,ir,ip

      if(iz.ne.nfz .or. ir.ne.nfr .or. ip.ne.nfp)then

         write(6,101)mapfile
 101     format(' mismatch in mapfile ',a32)
         write(6,102)nfz,nfr,nfp,iz,ir,ip
 102     format(' nfz,nfr,nfp = ',3i6,' iz,ir,ip = ',3i6)
         write(6,103)mapfile
 103     format(' mapfile ',a32,' not opened ')
         go to 998
       else
         write(6,106)mapfile
 106     format(' mapfile opened ',a32)
       endif
       

        read(31,*)zflmin,zflmax,zfldel
        read(31,*)rflmin,rflmax,rfldel
        read(31,*)pflmin,pflmax,pfldel
        do j3=1,ip
        do j2=1,ir
        read(31,*)(tzfl(j1,j2,j3),j1=1,iz)
        read(31,*)(trfl(j1,j2,j3),j1=1,iz)
        read(31,*)(tpfl(j1,j2,j3),j1=1,iz)
        end do
        end do

       write(6,107)zflmin,zflmax,zfldel,nfz
 107   format(' zflmin,zflmax,zfldel = ',3f9.3,' nfz = ',i4)
       write(6,108)rflmin,rflmax,rfldel,nfr
 108   format(' rflmin,rflmax,rfldel = ',3f9.3,' nfr = ',i4)
       write(6,109)pflmin,pflmax,pfldel,nfp
 109   format(' pflmin,pflmax,pfldel = ',3f9.3,' nfp = ',i4)

 998	continue

	elseif(mode.eq.'UNF')then

        read(31)iz,ir,ip

       if(iz.ne.nfz .or. ir.ne.nfr .or. ip.ne.nfp)then
         write(6,101)mapfile
         write(6,102)nfz,nfr,nfp,iz,ir,ip
         write(6,103)mapfile
         go to 999
       else
         write(6,106)mapfile
       endif
         
       read(31)zflmin,zflmax,zfldel
       read(31)rflmin,rflmax,rfldel
       read(31)pflmin,pflmax,pfldel
       read(31)trfl
       read(31)tpfl
       read(31)tzfl

       write(6,107)zflmin,zflmax,zfldel,nfz
       write(6,108)rflmin,rflmax,rfldel,nfr
       write(6,109)pflmin,pflmax,pfldel,nfp

 999    continue

  	end if
       close(31)
         
       write(6,111)mapfile
 111   format(' mapfile closed ',a32)


       RETURN
       END