anchor.diff:  	!
   	! 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})
  +        ! v_t = bu - u^2v + Dv (v_{xx} + v_{yy} + v_{zz})
   	!
   	! .. Parameters ..
   	!  Nx				= number of modes in x - power of 2 for FFT
  
positive.code: 	CHARACTER*100		::  name
  	name='./data/brusselator'
  
      plottime=plotgap
  	    ! initialisation of MPI
  	CALL system_clock(starttot,count_rate)	
  	CALL MPI_INIT(ierr)
  
negative.code:  	DEALLOCATE(uhigh,vhigh,&
  		        uhat,vhat,&
     		        kx,ky,x,y,savefield,&	
     		        time,stat=AllocateStatus)	
  	IF (AllocateStatus .ne. 0) STOP 
  
  	IF (myid.eq.0) THEN
  