drawing an elliptical figure

[ssv123]

hi,
I have written a code for the circular figures within the circle . But I am unable to draw the elliptical figure within the circle .Please help me in doing so...

I am attaching these files for your reference

 
SUBROUTINE Hole_Profile_Calc(dtheta_hole,hole_profile,	&
			   AxialHole,n_slice_net)
	!
	!     Calculation of circular axial hole profile (effective width) via radial slices 
	!			for one hole size
	!     This subroutine is called for each set of different size holes
	!
	!
	!-----------------------------------------------------------------------
		IMPLICIT NONE

	!-- Input Variables
		TYPE(AxialHoleType), INTENT(IN) :: AxialHole
		REAL, INTENT(IN) :: dtheta_hole              ! angular increment per slice, rad
		INTEGER, INTENT(IN) :: n_slice_net           ! net number slices per 360deg detailing hole profile

	!-- Output Variables
		REAL, INTENT(OUT) :: hole_profile(maxRadialSlices) ! chord distance of radial slices through holes over 360 

deg

	!-- Internal Variables
		REAL c                           ! radial (chord) length of each slice, m
		REAL d                           ! circumferential distance from hole center to slice, m
		REAL w_slice_hole                ! width of each hole slice
		REAL tmp                         ! temporary variable
		REAL theta_hole_width            ! one half angular width of each hole, rad
		REAL theta_hole(maxAxialHoles)		! angular location of each hole center, rad

		INTEGER j,jj,k
		INTEGER center_index             ! array index at center of hole
		INTEGER CW_index,CCW_index       ! array indices at position on both sides of hole center
		INTEGER nn


	!-----------------------------------------------------------------------
	!-----------------------------------------------------------------------

	!-- One half angular span of each hole ---
	  theta_hole_width = PI2*(0.5*AxialHole.OD)/(PI2*AxialHole.radius_to_center)
	!-- Number of slices per half hole --
		  nn = NINT(theta_hole_width/dtheta_hole)

	!-- Width of each slice --
		w_slice_hole = AxialHole.OD/(2*nn)

	!-- hole Locations assuming Uniform Distribution ---
		  DO jj = 1,AxialHole.n_holes
			theta_hole(jj) = AxialHole.theta_to_center + (jj-1)*2*PI/AxialHole.n_holes
		  ENDDO

	!-- Initialize hole_profile
		  hole_profile = 0.0

		  DO jj = 1, AxialHole.n_holes
			d = 0

			! index locating hole center
			center_index = NINT(theta_hole(jj)/dtheta_hole)+1 
 
			DO k = 1, nn
				! distance from hole center to slice, m
				d = d + w_slice_hole	
		
				! radial length of slice (chord) inc. bolt, m	
				tmp = 0.25*AxialHole.OD**2 - d*d
			
				IF (tmp < 0 )then
				   c = 0
				ELSE
				   c = 2.0*sqrt(tmp)		
				ENDIF

				CW_index = center_index - k + 1	
				CCW_index = center_index + k
				IF (CW_index < 1)then
				   CW_index = CW_index + n_slice_net
				ENDIF

				! radial slice length in CCW direction from center, m
				hole_profile(CCW_index) = c
		
				! radial slice length in CW direction from center, m
				hole_profile(CW_index) = c	   

			 ENDDO
		  ENDDO
  
	!-----------------------------------------------------------------------
	!-----------------------------------------------------------------------
	END SUBROUTINE Hole_Profile_Calc
	!-----------------------------------------------------------------------
	!-----------------------------------------------------------------------

DO ii = 1,StatorCore.n_axial_hole_sets
			!-- Calc hole profile for each type/size of hole
			
 CALL Hole_Profile_Calc(dtheta_hole,hole_profile,					&
									StatorCore.AxialHole(ii),n_slice_net)
end do

SUBROUTINE MCM_DISPLAY_STATOR_AXIAL_HOLES

USE DFLIB
USE SimulationSettingsModule
USE DefineGlobalVarsModule
IMPLICIT NONE

!-------------------------------------------------------------------------------
! Local Variables
!-------------------------------------------------------------------------------
INTEGER j
INTEGER k
INTEGER(2) status
REAL theta

REAL(8) wx1,wy1
REAL(8) wx2,wy2
REAL(8) wx00,wy00
REAL(8) wx0,wy0
REAL(8) WX3,WX4,WY3,WY4
!-----------------------------------------------------------------------
! Set & Display Stator Axial Circular Holes and Bolts
!-----------------------------------------------------------------------
!status = SETCOLOR(4) ! Red

DO k = 1,StatorCore.n_axial_hole_sets

! Coordinates of first stator hole of type k
wx00 = StatorCore.AxialHole(k).radius_to_center
wy00 = 0.0

DO j = 1,StatorCore.AxialHole(k).n_holes
! Rotate center coordinate for each stator hole
theta = StatorCore.AxialHole(k).theta_to_center &
+ (j-1)*PI2/StatorCore.AxialHole(k).n_holes



wx0 = wx00*cos(theta) - wy00*sin(theta)
wy0 = wx00*sin(theta) + wy00*cos(theta)

! Calculate new coordinates for each stator hole
wx1 = wx0 - 0.5*StatorCore.AxialHole(k).OD
wy1 = wy0 + 0.5*StatorCore.AxialHole(k).OD
wx2 = wx0 + 0.5*StatorCore.AxialHole(k).OD
wy2 = wy0 - 0.5*StatorCore.AxialHole(k).OD




! Display each rotated hole
status = ellipse_W($GBORDER, wx1,wy1,wx2,wy2)


! Display bolt (if present)
!IF (StatorCore.AxialHole(k).Bolt.present) THEN

! Calculate new coordinates for each stator bolt
! wx1 = wx0 - 0.5*StatorCore.AxialHole(k).Bolt.OD
! wy1 = wy0 + 0.5*StatorCore.AxialHole(k).Bolt.OD
! wx2 = wx0 + 0.5*StatorCore.AxialHole(k).Bolt.OD
! wy2 = wy0 - 0.5*StatorCore.AxialHole(k).Bolt.OD

! Display each rotated hole
! status = ELLIPSE_W($GBORDER, wx1,wy1,wx2,wy2)

! ENDIF

ENDDO
ENDDO

!-----------------------------------------------------------------------
END SUBROUTINE MCM_DISPLAY_STATOR_AXIAL_HOLES
!-----------------------------------------------------------------------