anchor.diff:                 CALL Cleanup()
                  RETURN
               END IF
  -
  +         
         do iBody = 1, p%NBody
   
            theta = (/ 0.0_R8Ki, 0.0_R8Ki, InitInp%PtfmRefztRot(iBody)/)
  
positive.code: 
        CLOSE ( UnWh ) ! Close file.
  
  
  
           ! Linear, frequency-dependent hydrodynamic added mass and damping from the
           !   radiation problem:
  
negative.code:                                      , iBody              &
                                                   )
              CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, 'WAMIT_Init')
              
        end do
  
        CALL MeshCommit ( u%Mesh              &
  