anchor.diff:        
         print *, "read ", ntIn, " time steps from input file."
         print *, PtfmMotIn
  -   
  -   ELSE
  -      ntIn = 0   ! flag to indicate no motion input file
      END IF
   
     
  
positive.code:    Integer(IntKi)                        :: iIn
     integer(intKi)                        :: Un
     
     TYPE(ProgDesc), PARAMETER             :: progdesc = ProgDesc( 'MoorDyn Driver' )
    
     CALL NWTC_Init()
        ! Display the copyright notice
     CALL DispCopyrightLicense( progdesc%Name )
        ! Obtain OpenFAST git commit hash
     git_commit = QueryGitVersion()
        ! Tell our users what they're running
     CALL WrScr( ' Running '//TRIM( progdesc%Name )//' a part of OpenFAST - '//TRIM(git_Commit)//NewLine//' linked with '//TRIM( NWTC_Ver%Name )//NewLine )
  
     ! -------------------------------------------------------------------------
     ! Read in prescribed motions from text file if available 
     ! (single 6DOF platform for now, plus one active tensioning command)
     ! (to be updated for versatile coupling in future)
     ! -------------------------------------------------------------------------
     
     CALL GetNewUnit( UnPtfmMotIn ) 
     
     CALL OpenFInpFile ( UnPtfmMotIn, "PtfmMotion.txt", ErrStat, ErrMsg ) 
     
     IF (ErrStat == 0 ) THEN
     
        print *, "Reading platform motion input data from PtfmMotion.txt"
           
        ! Read through length of file to find its length
        i = 1  ! start counter
        DO
  
negative.code:    Allocate(MD_Input(MD_interp_order + 1))
       
     ! set the input file name and other environment terms.
     !MD_InitInput%NStepWave   = 1        ! an arbitrary number > 0 (to set the size of the wave data, which currently contains all zero values)     
     MD_InitInput%g           = 9.81     ! This need to be according to g used in ElastoDyn 
     MD_InitInput%rhoW        = 1025     ! This needs to be set according to seawater density in HydroDyn      
  