Default for autocorrelation was not working correctly. Fixed!

1 parent 105dc561
Showing with 5 additions and 3 deletions
......@@ -75,6 +75,7 @@
N_Back = 1
N_auto = 0
OPEN(UNIT=5,FILE='parameters',STATUS='old',ACTION='read',IOSTAT=ierr)
IF (ierr /= 0) THEN
WRITE(*,*) 'unable to open <parameters>',ierr
......@@ -140,7 +141,7 @@
Write(6,*) "# of bins: ", Nbins
nbins = Nbins - n_skip
Write(6,*) "Effective # of bins: ", Nbins
N_auto=min(N_auto,Nbins/3)
! Allocate space
Allocate ( bins(Nunit,Nbins), bins_r(Nunit,Nbins), Phase(Nbins), V_help(Nbins), V_help_R(Nbins), Bins0(Nbins,Norb))
......
......@@ -46,7 +46,7 @@
Implicit none
Integer :: Nunit, Norb
Integer :: Nunit, Norb, N_auto
Integer :: no, no1, n, nbins, n_skip, nb, NT, NT1, Lt, N_rebin, N_cov, ierr, N_Back
real (Kind=Kind(0.d0)):: X, Y, dtau, X_diag
Complex (Kind=Kind(0.d0)), allocatable :: Xmean(:), Xcov(:,:)
......@@ -59,12 +59,13 @@
Real (Kind=Kind(0.d0)), allocatable :: Xk_p(:,:)
Character (len=64) :: File_out
NAMELIST /VAR_errors/ n_skip, N_rebin, N_Cov, N_Back
NAMELIST /VAR_errors/ n_skip, N_rebin, N_Cov, N_Back, N_auto
N_Back = 1
N_auto = 0
OPEN(UNIT=5,FILE='parameters',STATUS='old',ACTION='read',IOSTAT=ierr)
IF (ierr /= 0) THEN
WRITE(*,*) 'unable to open <parameters>',ierr
......
Markdown is supported
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!