
      subroutine samples(fr,id1,id2,id3,j)
      dimension fr(id1,id2,id3)
      imxp1 = id1
      kmx = id2
c
      write(6,"(' ')")
      write(6,"('j=',i3,' NT (1) fr(i=1)=',/(6e12.4))")
     +  j,(fr(1,k,1),k=1,kmx)
      write(6,"('fr(i=imxp1)=',/(6e12.4))") (fr(imxp1,k,1),k=1,kmx)
c
      write(6,"(' ')")
      write(6,"('j=',i3,' NTNM (28) fr(i=1)=',/(6e12.4))")
     +  j,(fr(1,k,28),k=1,kmx)
      write(6,"('fr(i=imxp1)=',/(6e12.4))") (fr(imxp1,k,28),k=1,kmx)

      write(6,"(' ')")
      write(6,"('j=',i3,' NOP (13) fr(i=1)=',/(6e12.4))")
     +  j,(fr(1,k,13),k=1,kmx)
      write(6,"('fr(i=imxp1)=',/(6e12.4))") (fr(imxp1,k,13),k=1,kmx)

      write(6,"(' ')")
      write(6,"('j=',i3,' NOPNM (40) fr(i=1)=',/(6e12.4))")
     +  j,(fr(1,k,40),k=1,kmx)
      write(6,"('fr(i=imxp1)=',/(6e12.4))") (fr(imxp1,k,40),k=1,kmx)
c
      return 
      end
