1 subroutine prev(v,m,n,ndim) 2c 3c ----- print out matrices 4c vectors density etc. 5c 6 implicit double precision (a-h,o-z) 7 common /rdwr/ iun1,iun2,iun3,iun4,iun5 8 dimension v(ndim,*) 9 10 max = 12 11 imax = 0 12 13100 imin = imax + 1 14 imax = imax + max 15 if (imax .gt. m) imax = m 16 write (iun3,9008) 17 write (iun3,8028) (i,i = imin,imax) 18 write (iun3,9008) 19 do j = 1,n 20 write (iun3,8048) j,(v(j,i),i = imin,imax) 21 end do 22 if (imax .lt. m) go to 100 23 24 return 25 9008 format(/) 26 8028 format(7x,12(3x,i3,3x)) 27 8048 format(i5,2x,12f9.5) 28 29 end 30