anchor.diff:     ! -------------------------------------------------------------------------
      ! Initialize data for VTK output
      ! -------------------------------------------------------------------------
  -            
      if ( p_FAST%WrVTK > VTK_None ) then
         call SetVTKParameters(p_FAST, InitOutData_ED, InitOutData_AD, InitInData_HD, InitOutData_HD, ED, BD, AD, HD, ErrStat2, ErrMsg2)      
            call SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName)
  
positive.code:       ! Add all of the input meshes:
     DO K_local = 1,NumBl
        CALL MeshWrBin( unOut, u_ED%BladePtLoads(K_local), ErrStat, ErrMsg )
        
        !!!write(unOut2,'(A)') '_____________________________________________________________________________'
        !!!write(unOut2,*)     'ElastoDyn blade ', k_local, ' mesh'
        !!!write(unOut2,'(A)') '_____________________________________________________________________________'      
        !!!CALL MeshPrintInfo(unOut2, u_ED%BladePtLoads(K_local))
        !!!write(unOut2,'(A)') '_____________________________________________________________________________'
     END DO            
     CALL MeshWrBin( unOut, u_ED%TowerPtLoads,            ErrStat, ErrMsg )
     CALL MeshWrBin( unOut, u_ED%PlatformPtMesh,          ErrStat, ErrMsg )
  
negative.code: 
     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'
  