!--------------------------------------------------
!
!      Routines to analyze results on PM simulations
!
!--------------------------------------------------
Module Analyze

CONTAINS

!--------------------------------------------
!
subroutine Analysis(mDENSIT)
!
!--------------------------------------------
use Tools
use Density
use LUXURY
use Random
use Power
use LinkerList
        character*120 :: NameF


           !-----  Fill the table for things to analyze
    !iPower         = 1   !-- all particles Power spectrum 
    !iPowerRSD      = 0   !       redshift distortions for DM
    !iDensityDistr  = 1   !-- PDF for DM for different cell-sizes
    ! iBDM  = 0
    Call Timing(6,-1)      ! start reading time

    NGRID_old = NGRID                  ! store old value of NGRID
                                       ! mDENSIT DM density index:
                                       !  = 1 if density was overwritten
                                       !      and needs to be restored
    !moment = 100.*(1./AEXPN-1.)+0.5 ! redshift*100.
    !moment = max(moment,0)
    moment = ISTEP
    Nseed = 13979821
               !------------ Density Distribution in real space
    If(iDensityDistr == 1)Then
       write(NameF,'(2(a,i4.4),3(a,i3.3))')'DensDistrDM.',moment,   &
                                        '.',Nrealization,'.dat'
       OPEN(18,FILE=TRIM(NameF),STATUS='UNKNOWN')
       write(18,'(a)')HEADER
       write(18,'(a,f7.4,a,i4,a,f8.3,a,i4,a,f8.2)')' Aexpn =',AEXPN,' Step=',ISTEP,    &
            ' Redshift= ',1./AEXPN-1.,' Ngrid= ',Ngrid
       If(mDENSIT==1)Then
          CALL DENSIT
          mDENSIT = 0
       End If
      Call DensDistr

    End If
                      !----------------- Power spectrum
    If(iPower == 1)Then
       If(mDENSIT==1)Then
          CALL DENSIT
          mDENSIT = 0
       EndIf       

      Call GetPower(0)  ! pk only real space, Ngrid mesh
    End If
                      !----------------- Save snapshot    

10 continue
 If(iBDM==1)Call BDM(mDENSIT)
     Call Timing(6,1)      ! stop reading time

  end subroutine Analysis



end Module Analyze