anchor.diff:  
   	!initial condition
   	CALL MPI_BARRIER(MPI_COMM_WORLD,ierr)
  -  CALL initialdata(Nx,Ny,Nz,x,y,z,uhigh,vhigh,decomp)
  +        CALL initialdata(Nx,Ny,Nz,x,y,z,uhigh,vhigh,decomp)
   	n=1
   	time(n)=0.0
   	plotnum=10000000
   	!Write initial data to disc
  -    CALL savedata(Nx,Ny,Nz,plotnum,name,savefield,uhigh,vhigh,decomp)
  +        CALL savedata(Nx,Ny,Nz,plotnum,name,savefield,uhigh,vhigh,decomp)
   	IF (myid.eq.0) THEN
   	 PRINT *,'Got initial data, starting timestepping'
   	END IF
  
positive.code: 	CALL decomp_info_init(Nx,Ny,Nz,decomp) ! physical domain
  	! initialise FFT library
  	CALL decomp_2d_fft_init
  	
  
  	ALLOCATE(Uhigh(decomp%xst(1):decomp%xen(1),&
     				   decomp%xst(2):decomp%xen(2),&
  
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
  