anchor.diff:  			ind=index(nameconfig,' ') -1
   			nameconfig=nameconfig(1:ind)//'.datbin'
   			! write out using 2DECOMP&FFT MPI-IO routines
  -			CALL decomp_2d_write_one(1,u,nameconfig)
  +			CALL decomp_2d_write_one(1,uabs,nameconfig)
   		END IF
   	END DO		
   	IF (myid.eq.0) THEN
  
positive.code: 			END IF
  			CALL decomp_2d_fft_3d(v,u,DECOMP_2D_FFT_BACKWARD)
  			u=u*modescalereal
  			nameconfig='./data/u'
  			plotnum=plotnum+1
  			WRITE(numberfile,'(i0)') 10000000+plotnum
  
negative.code: 	nameconfig=nameconfig(1:ind)//numberfile
  	ind=index(nameconfig,' ') -1
  	nameconfig=nameconfig(1:ind)//'.datbin'
  	CALL decomp_2d_write_one(1,uabs,nameconfig)
    
  	CALL decomp_2d_fft_3d(u,v,DECOMP_2D_FFT_FORWARD)
  	IF (myid.eq.0) THEN
  