anchor.diff:  !   END IF 
   !TODO: Fix below section for new Morison GJH 4/23/20
      !   
  -   !IF ( HD%Input(1)%Morison%DistribMesh%Committed ) THEN 
  +   !IF ( HD%Input(1)%Morison%Mesh%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', &
  +   !   call MeshWrVTK_Ln2Surface (p_FAST%TurbinePos, HD%Input(1)%Morison%Mesh, 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 )
  +   !                              p_FAST%VTK_Surface%MorisonRad, Sib=HD%y%Morison%Mesh )
      !END IF
      
      
  
positive.code:    IF ( p_FAST%CompHydro == Module_HD .and. allocated(HD%Input)) THEN     
        !TODO: Fix for Visualizaton GJH 4/23/20
        !call MeshWrVTK(p_FAST%TurbinePos, HD%Input(1)%Mesh, trim(p_FAST%VTK_OutFileRoot)//'.HD_Mesh_motion', y_FAST%VTK_count, p_FAST%VTK_fields, ErrStat2, ErrMsg2 )     
        !call MeshWrVTK(p_FAST%TurbinePos, HD%Input(1)%Morison%LumpedMesh, trim(p_FAST%VTK_OutFileRoot)//'.HD_MorisonLumped_motion', y_FAST%VTK_count, p_FAST%VTK_fields, ErrStat2, ErrMsg2 )     
        !call MeshWrVTK(p_FAST%TurbinePos, HD%Input(1)%Morison%DistribMesh, trim(p_FAST%VTK_OutFileRoot)//'.HD_MorisonDistrib_motion', y_FAST%VTK_count, p_FAST%VTK_fields, ErrStat2, ErrMsg2 )     
        
        if (p_FAST%CompSub == Module_NONE) then
  !TODO         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 )
           outputFields = .false.
        else         
           call MeshWrVTK(p_FAST%TurbinePos, HD%y%WamitMesh, trim(p_FAST%VTK_OutFileRoot)//'.HD_Mesh', y_FAST%VTK_count, p_FAST%VTK_fields, ErrStat2, ErrMsg2, p_FAST%VTK_tWidth, HD%Input(1)%WAMITMesh )
           outputFields = p_FAST%VTK_fields
        end if
        call MeshWrVTK(p_FAST%TurbinePos, HD%y%Morison%Mesh, trim(p_FAST%VTK_OutFileRoot)//'.HD_Morison', y_FAST%VTK_count, outputFields, ErrStat2, ErrMsg2, p_FAST%VTK_tWidth, HD%Input(1)%Morison%Mesh ) 
                    
     END IF
     
  ! SubDyn   
  
negative.code:          OutputFields = p_FAST%VTK_fields
        end if
       !TODO: Fix for Visualization GJH 4/23/20 
       call MeshWrVTK(p_FAST%TurbinePos, HD%Input(1)%Morison%Mesh, trim(p_FAST%VTK_OutFileRoot)//'.HD_Morison', &
                      y_FAST%VTK_count, OutputFields, ErrStat2, ErrMsg2, p_FAST%VTK_tWidth, Sib=HD%y%Morison%Mesh )
     END IF
     
     
  