program undump

C      This program reads the saved or rho file to make
C      a density contour plot. To compile use f77 -o Xcont Xcont.f

      integer  jmax,jmax1,jmax2
      integer  lmax
      integer  kmax,kmax1,kmax2
      parameter (jmax=512, jmax1=jmax+1, jmax2=jmax+2)
      parameter (kmax=32, kmax1=kmax+1, kmax2=kmax+2)
      parameter (lmax=128)
      integer jreq, j, k, i, l ,n, m, lp, lm
      real*8 rof3n,zof3n
      real*8 delt,den,time,elost,sound,omcen,tmassini,tmass,dtheta,pi
      real*8 tmassadd,tmassout,tmassacc,totcool,totheat
      real*8 rho(jmax2,kmax2,lmax),s(jmax2,kmax2,lmax),t(jmax2,kmax2
     &     ,lmax),a(jmax2,kmax2,lmax),eps(jmax2,kmax2,lmax)
      character infile*24

      
      pi=acos(-1.0)
      dtheta=2.0*pi/float(lmax)

      print*,'Input file name?'
      read(*,10)infile
 10   format(a24)
      
      open(unit=11,file=infile,form='unformatted')
      read(11)s
      read(11)t
      read(11)a
      read(11)rho
      read(11)eps
      read(11)ROF3N,ZOF3N,DELT,TIME,ELOST,DEN,SOUND,JREQ,OMCEN
      read(11)tmassini,tmass,tmassadd,tmassout,tmassacc,totcool,totheat
      close(11)
      print*,'check format',rho(208,3,23),rho(235,3,12)
      print*,time,ROF3N,tmassini,tmass,tmassadd,tmassout
     &     ,tmassacc,totcool,totheat

      do j=1,100
         write(*,*) j,rho(j,2,1)
      enddo

      end