anchor.diff:  !   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
  +!   END IF 
  +!TODO: Fix below section for new Morison GJH 4/23/20
  +   !   
  +   !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?            
  
positive.code:    ! morison surfaces
     !.......................
     
     IF ( HD%Input(1)%Morison%DistribMesh%Committed ) THEN      
        
        call move_alloc(InitOutData_HD%Morison%Morison_Rad, p_FAST%VTK_Surface%MorisonRad)
        
     END IF
     
  
negative.code:    CALL MeshWrBin( unOut, u_ED%PlatformPtMesh,          ErrStat, ErrMsg )
     CALL MeshWrBin( unOut, u_SD%TPMesh,                  ErrStat, ErrMsg )
     CALL MeshWrBin( unOut, u_SD%LMesh,                   ErrStat, ErrMsg )
     CALL MeshWrBin( unOut, u_HD%Morison%Mesh,      ErrStat, ErrMsg )
     CALL MeshWrBin( unOut, u_HD%WAMITMesh,                    ErrStat, ErrMsg )
     CALL MeshWrBin( unOut, u_MAP%PtFairDisplacement,     ErrStat, ErrMsg )
        ! Add how many BD blade meshes there are:
  