anchor.diff:  	! Declare variables
   	IMPLICIT NONE	
   	INCLUDE 'mpif.h'	
  -	INTEGER(kind=4), PARAMETER	::  Nx=64
  -	INTEGER(kind=4), PARAMETER 	::  Ny=64
  -	INTEGER(kind=4), PARAMETER 	::  Nz=64
  -	INTEGER(kind=4), PARAMETER ::  nmax=10000000
  +	INTEGER(kind=4), PARAMETER	 ::  Nx=64
  +	INTEGER(kind=4), PARAMETER 	 ::  Ny=64
  +	INTEGER(kind=4), PARAMETER 	 ::  Nz=64
  +	INTEGER(kind=4), PARAMETER       ::  nmax=10000000
   	REAL(kind=8), PARAMETER		 ::  Tmax=60.0
   	REAL(kind=8), PARAMETER		 ::  plotgap=0.05
  -	REAL(kind=8), PARAMETER	   ::&
  +	REAL(kind=8), PARAMETER	         ::&
   		pi=3.14159265358979323846264338327950288419716939937510d0
   	REAL(kind=8), PARAMETER		 ::  Lx=1.0,Ly=1.0,Lz=1.0
   	!equation specific
   	REAL(kind=8), PARAMETER		 ::  a=2.0d0
  -	REAL(kind=8), PARAMETER    ::  b=18.20d0
  +	REAL(kind=8), PARAMETER          ::  b=18.20d0
   	REAL(kind=8), PARAMETER		 ::  Du=0.050d0
   	REAL(kind=8), PARAMETER		 ::  Dv=0.050d0
  -	REAL(kind=8)		           ::  dt=0.1**3
  -	REAL(kind=8)	             ::  plottime=0.0
  -	INTEGER(kind=4) 	        ::  plotnum=0
  -	COMPLEX(kind=8)		::   utemp, vtemp, uhatemp
  +	REAL(kind=8)		         ::  dt=0.1**3
  +	REAL(kind=8)	                 ::  plottime=0.0
  +	INTEGER(kind=4) 	         ::  plotnum=0
  +	COMPLEX(kind=8)		         ::  utemp, vtemp, uhatemp
   	COMPLEX(kind=8), 	DIMENSION(:), ALLOCATABLE	::  kx,ky,kz
   	REAL(kind=8),  	        DIMENSION(:), ALLOCATABLE	::  x,y,z 
   	COMPLEX(kind=8), 	DIMENSION(:,:,:), ALLOCATABLE	::  uhigh,vhigh
   	REAL(kind=8),       DIMENSION(:,:,:), ALLOCATABLE	::  savefield
  -	COMPLEX(kind=8), DIMENSION(:,:,:), ALLOCATABLE	::  uhat,vhat
  -	REAL(kind=8), 	 DIMENSION(:), ALLOCATABLE	::  time
  +	COMPLEX(kind=8), DIMENSION(:,:,:), ALLOCATABLE	        ::  uhat,vhat
  +	REAL(kind=8), 	 DIMENSION(:), ALLOCATABLE	        ::  time
   	REAL(kind=8) 		::  modescalereal, myerr,allerr,mymaxv,mymaxu,maxu,maxv
   	INTEGER(kind=4)		::  i,j,k,l,n,m,mm,ll,AllocateStatus,stat
   	INTEGER(kind=4)		::  myid,numprocs,ierr
  
positive.code: 	IF ((myid.eq.0).and.(n.eq.nmax)) THEN
  	 PRINT *,'time',time(n),'n',n,'plot',plotnum
  	END IF
      CALL savedata(Nx,Ny,Nz,plotnum,name,savefield,uhigh,vhigh,decomp)
    	CALL system_clock(finish,count_rate)
  	IF (myid.eq.0) THEN
  	 PRINT *,'Finished time stepping steps:', n
  	 PRINT*,'Program took ',REAL(finish-start)/REAL(count_rate),'for time stepping', &
  	 REAL(finish-start)/REAL(count_rate)*numprocs
       CALL saveresults(n,nmax,time,error,tries,name)
  	 PRINT *,'Saved data'
  	END IF
  	CALL MPI_BARRIER(MPI_COMM_WORLD,ierr)
  	 ! clean up 
    	CALL decomp_2d_fft_finalize
  
negative.code: 	!
  	! This program solves Brusselator equations in 3 dimensions
  	! u_t = a + u^2v - (b+1)u + Du (u_{xx} + u_{yy} + u_{zz})
          ! v_t = bu - u^2v + Dv (v_{xx} + v_{yy} + v_{zz})
  	!
  	! .. Parameters ..
  	!  Nx				= number of modes in x - power of 2 for FFT
  