fortran 语言中怎么实现两个txt文件中同列的数据加减

Program test
implicit real*8 (a-h,o-z)
parameter (nx=3, ny=120)
dimension c1(120,3), c2(120,3), diff(120,3)
open(1,file='1',status='old',form='formatted')
open(2,file='2',status='old',form='formatted')
open(21,file='D',status='unknown',form='formatted')
do i=1,3
do j=1,120
read(1,*) c1(j,i)
read(2,*) c2(j,i)
enddo
enddo
do i=1,120
diff(i,1)=c2(i,1)
diff(i,2)=c2(i,2)-c1(i,2)
diff(i,3)=c2(i,3)-c1(i,3)
enddo
do i=1,3
do j=1,120
write(21,*) diff(j,i)
enddo
enddo
90 format(3f12.9)
end
求高手找错：文件2第2，3列与文件1第2，3列减，然后输出D（文件2第一列，相减得的2,3列）

do i=1,3
do j=1,120
这两行交换下顺序，输入输出都是

open打开文件 read(文件号,格式串)变量列表 具体的可以看看帮助文件：）