--- XFOIL694.orig/plotlib/ps_subs.f	2006-09-10 22:04:24.000000000 -0700
+++ XFOIL694/plotlib/ps_subs.f	2006-09-10 22:07:47.000000000 -0700
@@ -559,8 +559,8 @@
 c                   in the iout array (3,1,3,1,3,1,3,1)
 c     Note: the bit mask is shifted to always start counting on a '1' bit.
 c
-c---  Uses the library routines and() and rshift() for bit manipulation 
-c     which are present in most fortran libraries as extensions to f77.
+c---  Uses the library routines and() and ishft() for bit manipulation 
+c---  which are present in gfortran.
 c
       dimension iout(*)
 c
@@ -570,7 +570,7 @@
         ibitold = and(1,imask)
         if(ibitold.NE.0) go to 5
           nshft = n
-          imask = rshift(imask,1)
+          imask = ishft(imask,-1)
       end do
 c
  5    nout  = 0
@@ -594,7 +594,7 @@
 c
         ibitold = ibit
         nbits = nbits + 1
-        imask = rshift(imask,1)
+        imask = ishft(imask,-1)
       end do
 c--- Add final bit(s) to end of 16 bits checked
 c--- Now append any zero bits shifted out originally