Hold graph to plot many data sets

I'm having trouble plotting my data points. Below is the Subroutine I use to plot coastal lines for a world map. The default %pl connects every data point but I don't want a line connecting each Continent. Can anyone tell me how to keep the plot window open and hold it so I can add 1 Continent at a time?

        SUBROUTINE SIMPLE(X,Y)                                                         !
c--SimplePlot.for                                                                      !
        WINAPP                                                                         !
        INTEGER N,winio@,i                                                             !
        PARAMETER(N=5116)                                                              !
        DOUBLE PRECISION x(N),y(N)                                                     ! 
c--Create the data to be plotted.                                                      !
C        CCR=1                                                                          !
C        ENDPLOT=0                                                                      !
9911    DO i=1,N                                                                       !
C        CCR=2                                                                          !
          READ(4,5555) x(i),y(i)                                                       !
5555    FORMAT(9X,I4,8X,I4)                                                            !
C          IF(x(i).EQ.88) GO TO 9119                                                    !
C          IF(x(i).EQ.22) ENDPLOT=1
C          IF(x(i).EQ.22) GO TO 9119                                                     !                                                                               !
        ENDDO                                                                          !
9119    CONTINUE                                                                       !
c--Plot the data.                                                                      !
C        write(6,*) CCR                                                                 !
C        IF(CCR.EQ.1) GO TO 3310                                                        !
C        IF(CCR.EQ.2) GO TO 3311                                                        !            
3310    i=winio@('%ca[SIMPLEPLOT]%bg[white]&')                                         !
                                                                                       !
3311    i=winio@('%pl[x_array]&',950,650,N,x,y)                                        !
                                                                                       !
C        IF(ENDPLOT.EQ.0) GO TO 9911                                                    !
        i=winio@('%ff%nl%cn%tt[OK]')                                                   !
        RETURN                                                                         !
        END

[img]https://docs.google.com/document/d/1c-KQif-3xD_Wi9Rx6UskoHNpIDU_7JglVwduD1NAMbM/edit[/img]

This is my latest attempt at rewriting my subroutine. Perhaps it would be easier for everyone if I just posted both my datafile and subroutine as a program so others can attempt to run it and see what is actually plotting, yes?

Here is my data file: [

Here is version 1 where it collects all datapoints then plots:

        Program SIMPLE1
C                                                                                      C
C......................................Program1........................................C
                                                                        !
c--SimplePlot                                                           !
        WINAPP                                                          !
        INTEGER N,winio@,i                                              !
        PARAMETER(N=4890)                                               !
        DOUBLE PRECISION x(N),y(N)                                      !
c--Create the data to be plotted.                                       !
                                                                        !
        OPEN(UNIT=4, FILE='PRINTER(4).TXT', STATUS='OLD')               !                                                           !
9911    DO i=1,N                                                        !
          READ(4,*) x(i),y(i)                                           !
5555    FORMAT(9X,I4,8X,I4)                                             !                                                                            !
        ENDDO                                                           !                                                        !
                                                                        !
c--Plot the data.                                                       !
                                                                        ! 
        i=winio@('%ca[SIMPLEPLOT]%bg[white]&')                          !                                                                                     !
        i=winio@('%pl[x_array]&',970,621,N,x,y)                         !  
        i=winio@('%ff%nl%cn%tt[OK (NOT)]')                                 !                                                 !
        END                                                             !
                                                                        !
C......................................................................................C
C                                                                                      C

And here, I altered to try and and plots of each dataset one on top of the next but if you run this you'll see what happens!

C                                                                                      C
C......................................PROGRAM2........................................C
                                                                        !
        PROGRAM SIMPLE2                                                 !
c--SimplePlot.for                                                       !
        WINAPP                                                          !
        INTEGER N,winio@,i                                              !
        PARAMETER(N=4867)                                               !
        DOUBLE PRECISION x(N),y(N)                                      !
                                                                        ! 
c--Create the data to be plotted.                                       !
                                                                        !
        OPEN(UNIT=4, FILE='PRINTER(4).TXT', STATUS='OLD')               !
        ENDPLOT=0                                                       !
        j=0                                                             !
        GO TO 9911                                                      !
3333    CONTINUE                                                        !
        j=j+1                                                           !
9911    DO i=1,N                                                        !                                                !
          READ(4,*) x(i),y(i)                                           !
5555    FORMAT(9X,I4,8X,I4)                                             !
          IF(x(i).EQ.88) GO TO 9119                                     !
          IF(x(i).EQ.22) ENDPLOT=1                                      !
          IF(x(i).EQ.22) GO TO 9119                                     !                                                                               !
        ENDDO                                                           !
9119    CONTINUE                                                        !
                                                                        !
c--Plot the data.                                                       !
                                                                        !
        IF (i.EQ.1) GO TO 9911                                          !
        x(i)= x(i-1) !IF(x(i).EQ.88)                                    !
        y(i)= y(i-1) !IF(y(i).EQ.88)                                    !
        IF (j.GE.1) GO TO 3311                                          !
        i=winio@('%ca[SIMPLEPLOT]%bg[white]&')                          !                                                                                     !
3311    i=winio@('%pl[x_array]&',925,630,N,x,y)                         !  
        i=winio@('%ff%nl%cn%tt[NEXT!]')                                 !
C        WRITE(6,*) 'TYPE 1 TO CONTINUE :'                              !
C        READ(5,*) R                                                    !
C        IF(R.NE.1) GO TO 2222                                          !
        IF(ENDPLOT.EQ.0) GO TO 3333                                     !      
2222    CONTINUE                                                        !                                                          !
        END                                                             !
                                                                        !
C......................................................................................C
C                                                                                      C

I use an old Dell running windows XP pro (SP2) and Plato3 to run this code. If you can help please post anything at all, thanks.](https://docs.google.com/open?id=0Bzy5FZxLELbaVTNOUWZydzFSZS14NHV5aVlJOHBHdw[/url)

The second code again:

C                                                                                      C
C......................................PROGRAM2........................................C
                                                                        !
        PROGRAM SIMPLE2                                                 !
c--SimplePlot.for                                                       !
        WINAPP                                                          !
        INTEGER N,winio@,i                                              !
        PARAMETER(N=4867)                                               !
        DOUBLE PRECISION x(N),y(N)                                      !
                                                                        ! 
c--Create the data to be plotted.                                       !
                                                                        !
        OPEN(UNIT=4, FILE='PRINTER(4).TXT', STATUS='OLD')               !
        ENDPLOT=0                                                       !
        j=0                                                             !
        GO TO 9911                                                      !
3333    CONTINUE                                                        !
        j=j+1                                                           !
9911    DO i=1,N                                                        !                                                !
          READ(4,*) x(i),y(i)                                           !
5555    FORMAT(9X,I4,8X,I4)                                             !
          IF(x(i).EQ.88) GO TO 9119                                     !
          IF(x(i).EQ.22) ENDPLOT=1                                      !
          IF(x(i).EQ.22) GO TO 9119                                     !                                                                               !
        ENDDO                                                           !
9119    CONTINUE                                                        !
                                                                        !
c--Plot the data.                                                       !
                                                                        !
        IF (i.EQ.1) GO TO 9911                                          !
        x(i)= x(i-1) !IF(x(i).EQ.88)                                    !
        y(i)= y(i-1) !IF(y(i).EQ.88)                                    !
        IF (j.GE.1) GO TO 3311                                          !
        i=winio@('%ca[SIMPLEPLOT]%bg[white]&')                          !                                                                                     !
3311    i=winio@('%pl[x_array]&',925,630,N,x,y)                         !  
        i=winio@('%ff%nl%cn%tt[NEXT!]')                                 !
C        WRITE(6,*) 'TYPE 1 TO CONTINUE :'                              !
C        READ(5,*) R                                                    !
C        IF(R.NE.1) GO TO 2222                                          !
        IF(ENDPLOT.EQ.0) GO TO 3333                                     !      
2222    CONTINUE                                                        !                                                          !
        END                                                             !
                                                                        !
C......................................................................................C
C                                                                                      C

I suspect that you cannot do what you want via %pl or, at least the solution will be non-standard (you might be able to access the simpleplot dll directly with the aid of the simpleplot manuals).

More likely you will need to use a different approach based on %gr or maybe even %og.

This is just an opinion. I don't have the time to look into the details for you.

It looks like I could possibly use %gr inside my loop that grabs the x/y coordinates from my text file called PRINTER(4).TXT. Can anyone help me make a code that looks like: ???

 x(1) = 0
 y(1) = 0
 DO i=2,3
   read(4,*) x(i), y(i)
   i=winio@('%ww%gr[white,rgb_colours]&',900,600)
   i=winio@('bt[OK]%lw',ctrl)
   CALL draw_line_between@(x(i-1),y(i-1),x(i),y(i),RGB@(255,0,0))
 END DO

If I can control which lines are plotted then my problem is solved. Does anyone have an example code like the one above they can share? If so, can you keep plotting new lines on the existing plot of old lines? Any examples will help me a lot! Thanks :roll:

You have been mixing integer and real values. I'd also recommend changing file unit 4 to 14. Also a format statement with error testing for the reads would help. I've assumed 2 columns of 10 characters. Comma delimeter is also supported. Hopefully the following changes may give you some ideas of things to consider: include <clearwin.ins> ! integer4, parameter :: n = 6000 integer4 ix(n), iy(n), ni, i, ctrl, iostat real8 x(n), y(n), xmin, xmax, ymin, ymax, factor ! ! create window i = winio@ ('%ww&') i = winio@ ('%gr[white,rgb_colours]&',900,600) i = winio@ ('bt[OK]&') i = winio@ ('%lw',ctrl) ! ! read the values xmin = 0 xmax = 0 ymin = 0 ymax = 0 x(1) = 0 y(1) = 0 do i = 2,n read(14,fmt='(bn,2f10.0)',iostat=iostat) x(i), y(i) if (iostat < 0) exit if (iostat > 0) write (,*) 'Error reading line',i,' of data' if (x(i) < 0) exit if (xmin > x(i)) xmin = x(i) if (xmax < x(i)) xmax = x(i) if (ymin > y(i)) ymin = y(i) if (ymax < y(i)) ymax = y(i) end do ni = i-1 ! ! scale the values factor = min ( 800. / (xmax-xmin), 500. / (ymax-ymin) ) do i = 1,ni ix(i) = (x(i)-xmin) * factor + 50 iy(i) = (y(i)-ymin) * factor + 50 end do ! ! plot values X,y where X must be in the range 0:900 and y 0:600 DO i=2,ni CALL draw_line_between@ (ix(i-1),iy(i-1), ix(i),iy(i), RGB@(255,0,0)) END DO ! end

[quote='Natetran77']I'm having trouble plotting my data points. Below is the Subroutine I use to plot coastal lines for a world map. The default %pl connects every data point but I don't want a line connecting each Continent. Can anyone tell me how to keep the plot window open and hold it so I can add 1 Continent at a time?

Use the technique demonstrated in simdem46 which has many afvantages.

1. Areas can be outlined or filled with colour
2. Additional structures can be added (see the K-means example plot)
3. High resolution eps files can be archived
4. You can compile the main program in 64-bit and still use the code

Bill

Okay! I've made much progress since my last post. I implemented JohnCampbell's example code into my main program and it's running great 😄 These are the links to the files needed to run if anyone wants to try it out: (again, I'm running XP SP2 32 bit... so I don't need .txt extentions)

NUPLOTz.FOR (The Main Program) https://docs.google.com/open?id=0Bzy5FZxLELbaZHk2V1Q1X29SbVc3NGJTa0x1Y3RlZw

TAH (The first file my program asks for) https://docs.google.com/open?id=0Bzy5FZxLELbaUnlfeEJpTlRUeTI5OTh5dHZBYk9Hdw

MAP1 (The second file asked for containing map coordinates) https://docs.google.com/open?id=0Bzy5FZxLELbaY2FkLWg5MVpTaXE5aGhNMW9oWEdNUQ

Next, I'd like to be able to interact with the map, adding lines using the mouse pointer, plotting points and maybe even adding a [PRINT] button if possible. I've tried using the windows print commands but I get an Access error. Thanks John and everyone else that has commented.