shape_half_ring_2d.f90 Source File

program shape_half_ring_2d

    use forcad, only: rk, nurbs_surface

    implicit none
    type(nurbs_surface) :: shape


    !> Set up a half ring shape centered at 0,0,0 with inner radius 1 and outer radius 2.
    call shape%set_half_ring([0.0_rk, 0.0_rk, 0.0_rk], 1.0_rk, 2.0_rk)

    !> Export the control points to a VTK file for visualization.
    call shape%export_Xc('vtk/shape_half_ring_2d_Xc.vtk')

    !> Create the shape using the specified number of elements in each direction.
    call shape%create(60, 15)

    !> Export the geometry to a VTK file for visualization.
    call shape%export_Xg('vtk/shape_half_ring_2d_Xg.vtk')

    !-----------------------------------------------------------------------------
    ! Visualization using PyVista
    ! Note: PyVista is required for visualization. Install it using `pip install pyvista`
    !-----------------------------------------------------------------------------

    !> Show the control geometry and geometry using PyVista
    call shape%show('vtk/shape_half_ring_2d_Xc.vtk','vtk/shape_half_ring_2d_Xg.vtk')

    !> Finalize and clean up the shape object.
    call shape%finalize()

end program