anchor.diff:  
      logical, parameter                      :: OutputFields = .FALSE. ! due to confusion about what fields mean on a surface, we are going to just output the basic meshes if people ask for fields
      INTEGER(IntKi)                          :: NumBl, k, Twidth
  +   CHARACTER(1024)                         :: VTK_path
      INTEGER(IntKi)                          :: ErrStat2
      CHARACTER(ErrMsgLen)                    :: ErrMSg2
      CHARACTER(*), PARAMETER                 :: RoutineName = 'WrVTK_Surfaces'
  
positive.code: !! displays the copyright notice, and displays some version information (including addressing scheme and precision).
  SUBROUTINE FAST_ProgStart(ThisProgVer)
     TYPE(ProgDesc), INTENT(IN) :: ThisProgVer     !< program name/date/version description
     character(200) :: name, version, date
     character(200) :: git_commit, architecture, precision
     character(200) :: execution_date, execution_time, execution_zone
     
  
negative.code:    INTEGER(IntKi),               INTENT(  OUT) :: ErrStat          !< Error status of the operation
     CHARACTER(*),                 INTENT(  OUT) :: ErrMsg           !< Error message if ErrStat /= ErrID_None
  
     REAL(SiKi)                              :: RefPoint(3), RefLengths(2)               
     REAL(SiKi)                              :: x, y                
     REAL(SiKi)                              :: TwrDiam_top, TwrDiam_base, TwrRatio, TwrLength
     INTEGER(IntKi)                          :: topNode, baseNode
     INTEGER(IntKi)                          :: tipNode, rootNode, cylNode
     INTEGER(IntKi)                          :: NumBl, k
     CHARACTER(1024)                         :: VTK_path
     INTEGER(IntKi)                          :: ErrStat2
     CHARACTER(ErrMsgLen)                    :: ErrMsg2
     CHARACTER(*), PARAMETER                 :: RoutineName = 'SetVTKParameters'
  
     ErrStat = ErrID_None
     ErrMsg  = ""
     
     ! create the VTK directory if it does not exist
     call MKDIR( get_vtkdir_path(p_FAST%OutFileRoot) )
  
     ! initialize the vtk data
     p_FAST%VTK_Surface%NumSectors = 18
     p_FAST%VTK_Surface%HubRad     = InitOutData_ED%HubRad
     ! NOTE: we set p_FAST%VTK_Surface%GroundRad in SetVTKParameters_B4HD
  
     ! write the ground or seabed reference polygon:
     VTK_path = get_vtkroot_path( p_FAST%OutFileRoot )
     RefPoint = p_FAST%TurbinePos
     if (p_FAST%CompHydro == MODULE_HD) then
        RefLengths = p_FAST%VTK_Surface%GroundRad*VTK_GroundFactor/2.0_SiKi
        
        RefPoint(3) = p_FAST%TurbinePos(3) - InitOutData_HD%WtrDpth      
        call WrVTK_Ground ( RefPoint, RefLengths, trim(VTK_path) // '.SeabedSurface', ErrStat2, ErrMsg2 )   
        
        RefPoint(3) = p_FAST%TurbinePos(3) - InitOutData_HD%MSL2SWL    
        call WrVTK_Ground ( RefPoint, RefLengths, trim(VTK_path) // '.StillWaterSurface', ErrStat2, ErrMsg2 )       
     else
        RefLengths = p_FAST%VTK_Surface%GroundRad !array = scalar
        call WrVTK_Ground ( RefPoint, RefLengths, trim(VTK_path) // '.GroundSurface', ErrStat2, ErrMsg2 )         
     end if
     
     
  