anchor.diff:  			END IF
   			CALL decomp_2d_fft_3d(v,u,DECOMP_2D_FFT_BACKWARD)
   			u=u*modescalereal
  +        		DO k=decomp%xst(3),decomp%xen(3)
  +                		DO j=decomp%xst(2),decomp%xen(2)
  +                        		DO i=decomp%xst(1),decomp%xen(1)
  +                                		uabs(i,j,k)=abs(u(i,j,k))
  +                        		END DO
  +                		END DO
  +        		END DO
   			nameconfig='./data/u'
   			plotnum=plotnum+1
   			WRITE(numberfile,'(i0)') 10000000+plotnum
  
positive.code: 	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_fft_3d(u,v,DECOMP_2D_FFT_FORWARD)
  	IF (myid.eq.0) THEN
  
negative.code: 	! clean up 
    	CALL decomp_2d_fft_finalize
    	CALL decomp_2d_finalize
   	DEALLOCATE(u,v,pot,uabs,&
     		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
  