anchor.diff:     IF ( p_FAST%CompHydro == Module_HD ) THEN 
         
         if (p_FAST%CompSub == Module_NONE) then
  -         call MeshWrVTK(p_FAST%TurbinePos, HD%y%AllHdroOrigin, trim(p_FAST%VTK_OutFileRoot)//'.HD_AllHdroOrigin', y_FAST%VTK_count, p_FAST%VTK_fields, ErrStat2, ErrMsg2, p_FAST%VTK_tWidth, HD%Input(1)%Mesh )
  +         call MeshWrVTK(p_FAST%TurbinePos, HD%y%WAMITMesh, trim(p_FAST%VTK_OutFileRoot)//'.HD_AllHdroOrigin', y_FAST%VTK_count, p_FAST%VTK_fields, ErrStat2, ErrMsg2, p_FAST%VTK_tWidth, HD%Input(1)%WAMITMesh )
            outputFields = .false.
         else         
            OutputFields = p_FAST%VTK_fields
         end if
  -      
  -      call MeshWrVTK(p_FAST%TurbinePos, HD%Input(1)%Morison%DistribMesh, trim(p_FAST%VTK_OutFileRoot)//'.HD_MorisonDistrib', &
  -                     y_FAST%VTK_count, OutputFields, ErrStat2, ErrMsg2, p_FAST%VTK_tWidth, Sib=HD%y%Morison%DistribMesh )
  +     !TODO: Fix for Visualization GJH 4/23/20 
  +     ! call MeshWrVTK(p_FAST%TurbinePos, HD%Input(1)%Morison%DistribMesh, trim(p_FAST%VTK_OutFileRoot)//'.HD_MorisonDistrib', &
  +     !                y_FAST%VTK_count, OutputFields, ErrStat2, ErrMsg2, p_FAST%VTK_tWidth, Sib=HD%y%Morison%DistribMesh )
      END IF
      
      
  
positive.code: !   IF ( p_FAST%CompSub == Module_SD ) THEN
  !     call MeshWrVTK(p_FAST%TurbinePos, SD%Input(1)%TPMesh, trim(p_FAST%VTK_OutFileRoot)//'.SD_TPMesh_motion', y_FAST%VTK_count, OutputFields, ErrStat2, ErrMsg2 )     
  !      call MeshWrVTK(p_FAST%TurbinePos, SD%y%y2Mesh, trim(p_FAST%VTK_OutFileRoot)//'.SD_y2Mesh_motion', y_FAST%VTK_count, OutputFields, ErrStat2, ErrMsg2 )        
  !   END IF     
        
     IF ( HD%Input(1)%Morison%DistribMesh%Committed ) THEN 
        !if ( p_FAST%CompSub == Module_NONE ) then ! floating
        !   OutputFields = .false.
        !else
        !   OutputFields = p_FAST%VTK_fields
        !end if
           
        call MeshWrVTK_Ln2Surface (p_FAST%TurbinePos, HD%Input(1)%Morison%DistribMesh, trim(p_FAST%VTK_OutFileRoot)//'.MorisonSurface', &
                                   y_FAST%VTK_count, OutputFields, ErrStat2, ErrMsg2, p_FAST%VTK_tWidth, p_FAST%VTK_Surface%NumSectors, &
                                   p_FAST%VTK_Surface%MorisonRad, Sib=HD%y%Morison%DistribMesh )
     END IF
     
     
  ! Mooring Lines?            
  
negative.code:    ! morison surfaces
     !.......................
     
     IF ( HD%Input(1)%Morison%Mesh%Committed ) THEN      
      !TODO: FIX for visualization GJH 4/23/20  
      !  call move_alloc(InitOutData_HD%Morison%Morison_Rad, p_FAST%VTK_Surface%MorisonRad)
        
     END IF
     
  