anchor.diff:         InputFileData%kp_coordinate(i,1) =  TmpReAry(3)
          InputFileData%kp_coordinate(i,4) = -TmpReAry(4)
      ENDDO
  +   
      
      !---------------------- MESH PARAMETER -----------------------------------------
      CALL ReadCom(UnIn,InputFile,'Section Header: Mesh Parameter',ErrStat2,ErrMsg2,UnEc)
  
positive.code:    ErrStat = ErrID_None
     ErrMsg  = ""
           
     IF(InputFileData%analysis_type .NE. 1 .AND. InputFileData%analysis_type .NE. 2) &
         CALL SetErrStat ( ErrID_Fatal, 'Analysis type must be 1 (static) or 2 (dynamic)', ErrStat, ErrMsg, RoutineName )
     IF(InputFileData%rhoinf .LT. 0.0_BDKi .OR. InputFileData%rhoinf .GT. 1.0_BDKi) &
  
negative.code:    CALL ReadCom(UnIn,InputFile,'key point and initial twist units',ErrStat2,ErrMsg2,UnEc)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
        
     !InputFileData%kp_coordinate = 0.0_BDKi
     DO i=1,InputFileData%kp_total
         CALL ReadAry( UnIn, InputFile, TmpReAry, 4, 'kp_coordinate', 'Key point coordinates and initial twist', ErrStat2, ErrMsg2, UnEc )       
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  