anchor.diff:      CALL CheckError( ErrStat2, ErrMsg2 )
       IF (ErrStat >= AbortErrLev) RETURN
   
  -    InitOut%WriteOutputHdr = p%OutParam(1:p%NumOuts)%Name
  -    InitOut%WriteOutputUnt = p%OutParam(1:p%NumOuts)%Units     
  +    do i=1,p%NumOuts    
  +       InitOut%WriteOutputHdr(i) = p%OutParam(i)%Name
  +       InitOut%WriteOutputUnt(i) = p%OutParam(i)%Units     
  +    end do
       InitOut%Ver = FEAM_Ver  
   
       !............................................................................................
  
positive.code:     p%Eps      = InputFileData%Eps           ! Static iteration tolerance
  
      DO I = 1, p%NumLines
          p%GSL        (I,:,:)  = InputFileData%GSL     (I,:,:)
          p%LineCI     (I)      = InputFileData%LineCI  (I)
          p%LineCD     (I)      = InputFileData%LineCD  (I)
          p%Elength    (I)      = InputFileData%LUnstrLen(I)/p%NumElems
  
negative.code:     END DO   
  
      CALL FEAM_Solve(t, n, u, p, x, xd, z, OtherState, .true., ErrStat2, ErrMsg2 ) ! @yb: This is the time integration, correct?
        call setErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName)
        
      !Save current fairlead position for next step integration   
      DO I = 1,u%PtFairleadDisplacement%NNodes
          OtherState%FAST_FPA(I,1) = OtherState%FAST_FP(I,1)
  