1 subroutine check_vecgrad
2 implicit real*8 (a-h,o-z)
4 include 'COMMON.IOUNITS'
9 include 'COMMON.VECTORS'
10 dimension uygradt(3,3,2,maxres),uzgradt(3,3,2,maxres)
11 dimension uyt(3,maxres),uzt(3,maxres)
12 dimension uygradn(3,3,2),uzgradn(3,3,2),erij(3)
13 double precision delta /1.0d-7/
16 crc write(iout,'(2i5,2(3f10.5,5x))') i,1,dc_norm(:,i)
17 crc write(iout,'(2i5,2(3f10.5,5x))') i,2,uy(:,i)
18 crc write(iout,'(2i5,2(3f10.5,5x)/)')i,3,uz(:,i)
19 cd write(iout,'(2i5,2(3f10.5,5x))') i,1,
20 cd & (dc_norm(if90,i),if90=1,3)
21 cd write(iout,'(2i5,2(3f10.5,5x))') i,2,(uy(if90,i),if90=1,3)
22 cd write(iout,'(2i5,2(3f10.5,5x)/)')i,3,(uz(if90,i),if90=1,3)
29 uygradt(l,k,j,i)=uygrad(l,k,j,i)
30 uzgradt(l,k,j,i)=uzgrad(l,k,j,i)
43 cd write (iout,*) 'i=',i
51 dc_norm(j,i)=dc_norm(j,i)+delta
52 c fac=dsqrt(scalar(dc_norm(1,i),dc_norm(1,i)))
54 c dc_norm(k,i)=dc_norm(k,i)/fac
56 c write (iout,*) (dc_norm(k,i),k=1,3)
57 c write (iout,*) (erij(k),k=1,3)
60 uygradn(k,j,1)=(uy(k,i)-uyt(k,i))/delta
61 uygradn(k,j,2)=(uy(k,i-1)-uyt(k,i-1))/delta
62 uzgradn(k,j,1)=(uz(k,i)-uzt(k,i))/delta
63 uzgradn(k,j,2)=(uz(k,i-1)-uzt(k,i-1))/delta
65 c write (iout,'(i5,3f8.5,3x,3f8.5,5x,3f8.5,3x,3f8.5)')
66 c & j,(uzgradt(k,j,1,i),k=1,3),(uzgradn(k,j,1),k=1,3),
67 c & (uzgradt(k,j,2,i-1),k=1,3),(uzgradn(k,j,2),k=1,3)
73 cd write (iout,'(i5,3f8.5,3x,3f8.5,5x,3f8.5,3x,3f8.5)')
74 cd & k,(uygradt(k,l,1,i),l=1,3),(uygradn(k,l,1),l=1,3),
75 cd & (uygradt(k,l,2,i-1),l=1,3),(uygradn(k,l,2),l=1,3)
76 cd write (iout,'(i5,3f8.5,3x,3f8.5,5x,3f8.5,3x,3f8.5)')
77 cd & k,(uzgradt(k,l,1,i),l=1,3),(uzgradn(k,l,1),l=1,3),
78 cd & (uzgradt(k,l,2,i-1),l=1,3),(uzgradn(k,l,2),l=1,3)