SUBROUTINE LGLIN

 LGLIN plots a solid curve through a set of values from x to y. Either x or y or both may in the process be converted to  logarithmic form. Symbols may be used at selected intervals. (see also  SCALG and LINE). The plotted x values are computed according to,

    for logarithmic scaling,

       xplotted= (alog10(abs(x(i))+1.e-38)-xm)/dx

    for linear scaling,

       xplotted= (x(i)-xm)/dx

CALL LGLIN (x,y,n,k,j,l,lg,ix,iy,xm,dx,ym,dy)

x     (R): array containing the x coordinates
y     (R): array containing the y coordinates
n     (I): number of data points in x and y
           (the number of data points in x,y must be equal)
k     (I): take every first (k=1), second (k=2) value etc.
           (normally k=1)
j     (I): plotting symbol spacing flag
           > 0 : symbol is plotted at every jth plotted point
                 connected by lines
           = 0 : no symbols, lines only
           < 0 : symbol is plotted at every abs(j)th plotted point
                 with no connecting lines
l     (I): plot symbol number (see SYMBOL)
lg    (I): log option
           = - 2 : x and y are plotted using logarithmic scaling
           = - 1 : x logarithmic, y linear
           =   1 : x linear, y logarithmic
ix,iy (I): start index if arrays (normally ix,iy=1)
xm    (R): minimum value scale factor for x array
dx    (R): increment scale factors for x array
ym    (R): minimum value scale factor for y array
dy    (R): increment scale factors for y array