I have two variables say x and y and both have around 60 points in them(basically values of the x and y axis of the plot)
x
y
You could use a do loop:
do
do i=1,size(xpts) write(21,"(T2,F10.4: T60,F10.4)"), xpts(i), ypts(i) enddo