🚀 go-pugleaf

Sorry, I jumped the gun.  The problem is not with boundaries or 
alignment, but instead there is a bizarre problem with double 
precision.  for example, try compiling the following:

      subroutine plot (x, n)
      double precision x(n)
      double precision xmin
      integer w
      j = (x(i) - xmin) * w    <== this works, BUT
                                   j = w * (x(i) - xmin) does not!
      return
      end

changing the minus sign to a division or multication works, as does other
random things, like switching xmin with x(i).