C ______________________________ Profiles for Galaxy Simulations
C		             
C		                     
      INCLUDE 'PMparameters.h'
      INCLUDE 'PMgalaxy.h'
      REAL         INPUT,moment
      Character    FileASCII*100,Analysis*3
      Character  File1*100,File2*100
      EQUIVALENCE (Ndisk,extras(90)),
     &               (Nbulge,extras(97)),(aMassDisk,extras(93)),
     &               (Box,extras(100))


      Analysis ='GD2'
      ScaleM =2.146e5  !1.e6/2. ! /2.
      ScaleT =9.8e8
      Box    = 1.
C...................................................................
C			Read data and open files
       write(*,*)' enter 3 digits of the snapshot'
       read(*,*) imoment
         write(File1,'("FILES/snapshot_",i3.3)')imoment
         write(*,*) File1

      CALL  ReadGadget(File1,N,ScaleM,ScaleT)
      ISTEP =imoment
      CALL BLOCK_DATA

         write(FileASCII,'("RESULTS/Result",a,".DM.",f6.4,".dat")')
     &    Analysis,AEXPN
      OPEN(37,FILE=FileASCII,STATUS='UNKNOWN') 
        write(FileASCII,'("RESULTS/Result",a,".SPH.",f6.4,".dat")')
     &    Analysis,AEXPN
      OPEN(39,FILE=FileASCII,STATUS='UNKNOWN') 

         write(FileASCII,'("RESULTS/Result",a,".CYL.",f6.4,".dat")')
     &    Analysis,AEXPN
      OPEN(41,FILE=FileASCII,STATUS='UNKNOWN') 

         write(FileASCII,'("RESULTS/Result",a,".BAR.",f6.4,".dat")')
     &    Analysis,AEXPN
      OPEN(43,FILE=FileASCII,STATUS='UNKNOWN') 
      OPEN(18,FILE='Stat.dat',STATUS='UNKNOWN',access='append')
         write(FileASCII,'("RESULTS/Result",a,".PRT.",f6.4,".dat")')
     &    Analysis,AEXPN
      OPEN(20,FILE=FileASCII,STATUS='UNKNOWN')
c      OPEN(30,FILE='particlesu.dat',STATUS='UNKNOWN',form='UNFORMATTED')


      aMassDisk = Ndisk*ScaleM
      time      = AEXPN*1.e9
      Nbulge    = 0

      write (*,*) ' mass of disk=',aMassDisk
      write (*,*) ' RDTAPE is done Time=',time,' a=',AEXPN
          WRITE (*,100) time/1.e9,
     +                  Ndisk,Nbulge,aMassDisk,
     +                  bar_max*1.e3,dZ_slice*1.e3
      Do i=37,43,2
          WRITE (i,100) time/1.e9,
     +                  Ndisk,Nbulge,aMassDisk,
     +                  bar_max*1.e3,dZ_slice*1.e3

      EndDo 
100   FORMAT(T2,'Time(Gyrs)',
     &    T15,'Ndisk',T25,'Nbulge',T35,'Mdisk(Msun)',/
     &    T2,f8.5,T15,i8,T25,i7,T35,g11.4,/
     +            ' Max Bar radius(kpc)=',
     &             f7.3,' Slice in Z for cylindr shells=',f8.3)

c      Box  =INPUT(' Enter box size in comoving Mpc/h =')
c      Box =0.5 ! NG3109
c      Box = 2   ! 'bulge'

      BoxV   = Box*100.    ! Box size in km/s
      ScaleV = BoxV/AEXPN/NGRID      ! scale factor for Velocities (km/s)

      ScaleH  = 0.
      if(Ndisk.le.0)write (*,*) ' wrong number of disk particles:',Ndisk
      if(Ndisk.le.0)Stop



      Call ChangePnts(N,Box,ScaleH)
      Call DumpPnts(N,Box,ScaleM)

c      Call AngMom(N,30000,3.,Box,ScaleM)
c      CALL PV(N,Box,ScaleM)
c      CALL RandV(N,Box,ScaleM)
      Call Profile(N,Box,ScaleM)
      Call Statist(N,Box,ScaleM)
c      Call Projection(N,Box,0.)
c      Call Projection(N,Box,5.)

      END