home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #18 / NN_1992_18.iso / spool / bit / listserv / sasl / 3703 < prev    next >
Encoding:
Text File  |  1992-08-13  |  3.5 KB  |  110 lines

  1. Comments: Gated by NETNEWS@AUVM.AMERICAN.EDU
  2. Path: sparky!uunet!zaphod.mps.ohio-state.edu!darwin.sura.net!paladin.american.edu!auvm!IBM-B.RUTHERFORD.AC.UK!MJF
  3. X-Delivery-Notice: SMTP MAIL FROM does not correspond to sender.
  4. Message-ID: <SAS-L%92081309040151@VTVM2.BITNET>
  5. Newsgroups: bit.listserv.sas-l
  6. Date:         Thu, 13 Aug 1992 14:02:09 BST
  7. Reply-To:     "Mike Froggatt" <MJF@IBM-B.RUTHERFORD.AC.UK>
  8. Sender:       "SAS(r) Discussion" <SAS-L@UGA.BITNET>
  9. From:         "Mike Froggatt" <MJF@IBM-B.RUTHERFORD.AC.UK>
  10. Subject:      GCONTOUR question
  11. Lines: 97
  12.  
  13. A colleague who does not subscribe to the list has a SAS/GRAPH problem.
  14. I wondered if any kind SAS-Ler has any suggestions.
  15. Thanks in advance for any help.
  16.  
  17. Mike Froggatt
  18. Rutherford Appleton Laboratory
  19. Oxford, England
  20.  
  21. ------------------------------------------------------------------------
  22.  
  23. I am producing contour maps using GCONTOUR. I am also trying to overlay a
  24. grid on top of the plot using the HREF and VREF options in the PLOT
  25. statement but to no avail. I have even specified the colour of the lines to
  26. be black to ensure that the grid stands out. I am also using the JOIN and
  27. PATTERN options to produced joined up areas of the same contour level. I
  28. just wonder whether this has any effect on the HREF and VREF lines.
  29.  
  30. ...
  31.  
  32. Here is the SAS file. I have since found that if you do not specify
  33. PATTERN and JOIN the grid lines are drawn. I am therefore wondering
  34. whether the lines are being overwritten and indeed whether you can have
  35. grid lines if the contour pattern is filled in.
  36. FILENAME GSASFILE 'VNC06 POST H';
  37. filename sitedata 'SITES DATA A';
  38. GOPTIONS RESET=GLOBAL DEVICE=APPLELW
  39.          HANDSHAKE=XONXOFF
  40.          ctext=black
  41.          ctitle=black
  42.          BORDER
  43.          FTEXT=SWISS
  44.          GUNIT=PCT
  45.          HTITLE=3
  46.          HTEXT=2
  47.          GACCESS=GSASFILE
  48.          COLORS=(white grayf0 grayd3  graya7
  49.                    gray7b  gray4f  black);
  50. /* GRID POINT DATA SET VALUES TO CREATE CONTOURS  */
  51. DATA cont;
  52.     infile contdata;
  53.     INPUT x rb4. y rb4. z rb4.;
  54.     IF sqrt((x-255/2)**2+(y-255/2)**2) > 255*7/20
  55.                 THEN z=-999.0;
  56.     x=x*20/255+300;
  57.     y=y*20/255+960;
  58. RUN;
  59. /* Read in sites + coordinates */
  60. data sites;
  61.      infile sitedata;
  62.      input projsite $ 1-17 xcoord ycoord;
  63.      put projsite;
  64. run;
  65. /* DEFINE ANNOTATION */
  66. data anno;
  67.     set sites;
  68.     length function style color $ 8 text $ 17;
  69.     retain xsys '2' ysys '2' hsys '1' when 'a';
  70. /* Create observation for a star */
  71.    function='label'; color='black'; size=5; text='M'; position='5';
  72.   style='special'; x=xcoord; y=ycoord;
  73.   output;
  74. /* Create observation for the text */
  75. function='label'; color='black'; size=3;
  76. text=projsite; position='8'; style='swissb';
  77.   output;
  78. run;
  79. /* Define titles */
  80. TITLE1 C=BLACK F=SWISS 'ND06 Normalised Velocity Contours';
  81. TITLE2 C=BLACK F=SWISS 'Mean over all directions';
  82. TITLE3 C=Black f=swiss 'Reference speed=13.81 m/s';
  83. /* Define pattern definitions */
  84. pattern value=msolid;
  85. /* DEFINE ANNOTATION */
  86. LEGEND VALUE=(J=L) LABEL=('Norm. Vel. contours:') ACROSS=3
  87.                     FRAME cborder=black;
  88. /* Define axes */
  89.  axis1 order=(303 to 317 by 2)
  90.        label=(justify=c 'E (km)');
  91.  axis2 order=(963 to 977 by 2)
  92.        label=('N (km)');
  93. PROC Gcontour data=cont;
  94.       plot y*x=z
  95.           /
  96.            levels=0.8 to 1.4 by 0.1
  97.            haxis=axis1
  98.            vaxis=axis2
  99.            cv=black
  100.            ch=black
  101.            ctext=black
  102.            caxis=black
  103.            annotate=anno
  104.            pattern
  105.            join
  106.            href=(303 to 317 by 2)
  107.            vref=(963 to 977 by 2)
  108.                   LEGEND=LEGEND1;
  109. RUN;
  110.