c      for VRML frame program
       subroutine fram31(x0,y0,z0,xm,ym,zm,xp,yp,zp,ic)
       write(10,10)
   10  format("Shape {"/"geometry IndexedLineSet {")
       write(10,12)
   12  format("color Color {"/"color [ 1.0 1.0 1.0,  1.0 0.0 0.0,
     &  0.0 1.0 0.0,  0.0 0.0 1.0 ]")
       write(10,14)
   14  format("}")
c
       write(10,20)
   20  format("coord Coordinate {"/"point [")
   22  format(3f6.1,",")
       write(10,22) xm,ym,z0
       write(10,22) xm,ym,zp
       write(10,22) xm,yp,zp
       write(10,22) xm,yp,z0
       write(10,22) xm,ym,z0
c
       write(10,22) x0,ym,z0
       write(10,22) x0,ym,zp
       write(10,22) x0,yp,zp
       write(10,22) x0,yp,z0
       write(10,22) x0,ym,z0
c
       write(10,22) xp,ym,z0
       write(10,22) xp,ym,zp
       write(10,22) xp,yp,zp
       write(10,22) xp,yp,z0
       write(10,22) xp,ym,z0
c       
       write(10,22) xm,ym,z0
       write(10,22) xm,ym,zp
       write(10,22) xp,ym,zp
       write(10,22) xp,ym,z0
       write(10,22) xm,ym,z0
c
       write(10,22) xm,y0,z0
       write(10,22) xm,y0,zp
       write(10,22) xp,y0,zp
       write(10,22) xp,y0,z0
       write(10,22) xm,y0,z0
c
       write(10,22) xm,yp,z0
       write(10,22) xm,yp,zp
       write(10,22) xp,yp,zp
       write(10,22) xp,yp,z0
       write(10,22) xm,yp,z0
c
       write(10,22) x0,y0,z0
       write(10,23) x0,y0,zp
   23  format(3f6.1,/"]"/"}")
       write(10,24) ic,ic,ic,ic,ic,ic,ic
   24  format("colorIndex [",7(1x,i1)," ]")
       write(10,26)
   26  format("coordIndex [ 0 1 2 3 4 -1 5 6 7 8 9 -1",
     & " 10 11 12 13 14 -1"/" 15 16 17 18 19 -1 20 21 22 23 24 -1",
     & " 25 26 27 28 29 -1 30 31 ]")
       write(10,28)
   28  format("colorPerVertex FALSE"/"}"/"}")
       return
       end