anchor.diff:  
   !IF SOUTION IS NOT YET CONVERGED TRY TAKING LOAD IN 2 STEPS
   	F_total = BD_Parameter%F_ext	!DEFINES DUMMY VARIABLE TO CUT LOAD IN HALF, NJ 3/18/2014
  -!	DO WHILE (n .NE. 0)  
  +	DO WHILE (n .NE. 0)  
       k=n
   		DO i=1,k
             BD_Parameter%F_ext=F_total/n*i
  
positive.code:    Integer(IntKi)                     :: j               ! counter for various loops
     Integer(IntKi)                     :: k               ! counter for various loops
     Integer(IntKi)                     :: n               ! counter load reduction
     Integer(IntKi)                     :: m               ! counter load reduction
  
     OPEN(unit = OutUnit, file = 'BeamDyn.out', status = 'REPLACE',ACTION = 'WRITE')
     OPEN(unit = OutQiUnit, file = 'QiStaticDisp.out', status = 'REPLACE',ACTION = 'WRITE')
  
negative.code:                       &BD_Parameter%elem_total, BD_Parameter%dof_total,BD_Parameter%node_total,&
                        &BD_Parameter%ngp,BD_Parameter%niter,BD_Parameter%piter) 
          
          IF (BD_Parameter%piter .LT. BD_Parameter%niter) THEN
            n=0  
            EXIT
          ELSE
              n=n+1
              WRITE(*,*) "Warning: Load may be too large, BeamDyn will attempt to solve with addition steps"
              WRITE(*,*) "n=",n
              !BD_OtherState%uuNf=0
          ENDIF!CONDITION IF NR ITERATIONS HAVE BEEN REACHED ZERO OUT uunF, NJ 3/18/2014 
  		     
          ENDDO
  	END DO
  	
  
  !IF SOUTION IS NOT YET CONVERGED TRY TAKING LOAD IN 3 STEPS
  