00001
00002 c
00003 c Copyright (c) 1986,1987,1988,1989,1990,1991,1992,1993,
00004 c by Steve McMillan, Drexel University, Philadelphia, PA.
00005 c
00006 c All rights reserved.
00007 c
00008 c Redistribution and use in source and binary forms are permitted
00009 c provided that the above copyright notice and this paragraph are
00010 c duplicated in all such forms and that any documentation,
00011 c advertising materials, and other materials related to such
00012 c distribution and use acknowledge that the software was developed
00013 c by the author named above.
00014 c
00015 c THIS SOFTWARE IS PROVIDED ``AS IS'' AND WITHOUT ANY EXPRESS OR
00016 c IMPLIED WARRANTIES, INCLUDING, WITHOUT LIMITATION, THE IMPLIED
00017 c WARRANTIES OF MERCHANTIBILITY AND FITNESS FOR A PARTICULAR PURPOSE.
00018 c
00019
00020
00021 subroutine aerrbars(xc,yc,dx,dy)
00022 save
00023 c
00024 c Draw one-sided error bars from (xc,yc), of length(s)
00025 c dx and/or dy, all measured in user units.
00026 c
00027 common/scales/xl,xr,dinchx,ybot,ytop,dinchy,rlen,slen
00028 data cap/.07/
00029 c
00030 iboth=1
00031 go to 10
00032 c
00033 entry axerr(xc,yc,dx)
00034 c
00035 iboth=-1
00036 10 rc=(xc-xl)*dinchx
00037 sc=(yc-ybot)*dinchy
00038 dr=dx*dinchx
00039 c
00040 c *** horizontal bar ***
00041 c
00042 call plotin(rc,sc,3)
00043 r=rc+dr
00044 call plotin(r,sc,2)
00045 call plotin(r,sc-cap,3)
00046 call plotin(r,sc+cap,2)
00047 if(iboth.lt.0)return
00048 go to 20
00049 c
00050 entry ayerr(xc,yc,dy)
00051 c
00052 rc=(xc-xl)*dinchx
00053 sc=(yc-ybot)*dinchy
00054 20 ds=dy*dinchy
00055 c
00056 c *** vertical bar ***
00057 c
00058 call plotin(rc,sc,3)
00059 s=sc+ds
00060 call plotin(rc,s,2)
00061 call plotin(rc-cap,s,3)
00062 call plotin(rc+cap,s,2)
00063 c
00064 end
