*  Another plot of some ozone data.  This graph shows the ozone for
   three stations and adds station labels.  ;
OPTIONS LS=80 NODATE;

*  Ozone hourly ozone, 8am to 10pm, for 10 Houston stations.  8/9/88;
DATA one;
  ARRAY statid(10) $ ('0024' '0045' '0047' '0059' '0061' '1003' '1010'
    '1025' '1035' '1037');
  ARRAY all(*) stat1-stat10;
  INPUT stat1-stat10;
  n+1;
    DO i = 1 TO 10;
      hour = n + 7;
      station = statid(i);
      ozone = all(i);
      KEEP hour station ozone;
      OUTPUT;
      END;
  CARDS;
  40   40   10   10   10   10   20   20   20    10
  80   90   30    .   40   10   40   80   90    20
  60    .  110    .  140   60   70   60  140   130
  50   70  170    .  170  100  120   70  110   190
  60   50  150    .  140    .  180   90  100   130
  70   50  130    .  120    .  130   80  100   110
  70   60  110    .  100    .  120   90   90    90
  70   60  100    .  110  140  120   60   70    70
  70   60   70    .  100  100  100   50   70    70
  50   50   60    .   80   60   70   60   70   100
  40   40   50    .   40   40   40   40   30    10
  30   20   50    .   30   40   50   30   20    10
;
PROC PRINT;

DATA two;
  SET one;
  IF (station = '1010') or (station = '1025') or (station = '1035');

PROC PRINT;

FILENAME grafout 'ozone2.ps';
GOPTIONS DEVICE=psll GSFMODE=REPLACE  GSFLEN=64  GSFNAME=grafout
  FTEXT=SWISS FTITLE=SWISS HTEXT=1.5 HTITLE=3;

PROC GPLOT;
SYMBOL1 INTERPOL=JOIN VALUE=STAR COLOR=RED;
SYMBOL2 INTERPOL=JOIN VALUE=DIAMOND COLOR=BLUE;
SYMBOL3 INTERPOL=JOIN VALUE=CIRCLE COLOR=GREEN;
  PLOT ozone*hour=station;
TITLE 'Hourly ozone concentrations for 3 stations';
RUN;