anchor.diff:  	nameconfig=nameconfig(1:ind)//numberfile
   	ind=index(nameconfig,' ') -1
   	nameconfig=nameconfig(1:ind)//'.datbin'
  -	CALL decomp_2d_write_one(1,u,nameconfig)
  +	CALL decomp_2d_write_one(1,uabs,nameconfig)
     
   	CALL decomp_2d_fft_3d(u,v,DECOMP_2D_FFT_FORWARD)
   	IF (myid.eq.0) THEN
  
positive.code: 	! clean up 
    	CALL decomp_2d_fft_finalize
    	CALL decomp_2d_finalize
   	DEALLOCATE(u,v,pot,&
     			kx,ky,kz,x,y,z,&
     			time,stat=AllocateStatus)	
  	IF (AllocateStatus .ne. 0) STOP 
  	IF (myid.eq.0) THEN
  	   	PRINT *,'Program execution complete'
  	END IF
  	CALL MPI_FINALIZE(ierr)		
  	END PROGRAM main
negative.code: 	!  j				= loop counter in y direction
  	!  k				= loop counter in z direction
  	!  n				= loop counter for timesteps direction	
  	!  allocatestatus		= error indicator during allocation
  	!  start			= variable to record start time of program
  	!  finish			= variable to record end time of program
  	!  count_rate			= variable for clock count rate
  	!  dt				= timestep
  	!  modescalereal		= Number to scale after backward FFT 
  	!  myid				= Process id
  	!  ierr				= error code
  	!  p_row			= number of rows for domain decomposition
  