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 c
00020 c
00021 c--- subroutine digital
00022 c
00023 c--- creates a digital map of an array
00024 c
00025 c--- created 12.aug.80
00026 c 27.mar.81 for xplot
00027 c
00028 c--- parameters
00029 c a - data array
00030 c m,n - dimensions of a to plot
00031 c isc - icode for nomber
00032 c x0,x1,y0,y1 - delimeters of plot
00033 c msiz,nsiz - dimensions of array a
00034 c
00035 subroutine digital (a,m,n,isc,x0,x1,y0,y1,msiz,nsiz)
00036 save
00037 real*4 a(msiz,nsiz)
00038 common/fontc1/offx,offy,lastinc,xp,yp,xmax,xmin,ymax,ymin
00039 dx=(x1-x0)/m
00040 dy=(y1-y0)/n
00041 ht=dy/2.
00042 do 10 j=1,n
00043 do 10 i=1,m
00044 il=alog10(abs(a(i,j)))+isc+1
00045 if (a(i,j).lt.0.) il=il+1
00046 h1=.5*dx/il
00047 ht=min(h1,ht)
00048 10 continue
00049 do 20 j=1,n
00050 do 20 i=1,m
00051 call nomber (0.,0.,-ht,a(i,j),0.,isc)
00052 call nomber (x0+(i-.5)*dx-.5*(xmax-xmin),
00053 + y0+(j-.5)*dy-.5*(ymax-ymin)+.05,
00054 + ht,a(i,j),0.,isc)
00055 20 continue
00056 end
