目前有一笔11x11的资料(二维阵列?) 每一个横行的资料我都会分别做计算所以使用一维阵列去读值 但是我要的资料只有矩阵的一半，对角线是右上到左下 如下： １　１　１　１　１　１　１　１　１　１　１ １　１　１　１　１　１　１　１　１　１　０ １　１　１　１　１　１　１　１　１　０　０ １　１　１　１　１　１　１　１　０　０　０ １　１　１　１　１　１　１　０　０　０　０ １　１　１　１　１　１　０　０　０　０　０ １　１　１　１　１　０　０　０　０　０　０ １　１　１　１　０　０　０　０　０　０　０ １　１　１　０　０　０　０　０　０　０　０ １　１　０　０　０　０　０　０　０　０　０ １　０　０　０　０　０　０　０　０　０　０ １的部分是我要读的值，并做一点计算 该部分的程式码如下： !============================================================= CHARACTER*1000 FMT1,FMT2,FMT3,FMT4,FMT5 INTEGER ND,stat,I,J,K,L,M,N INTEGER, PARAMETER :: LYR=11 REAL U_B,U_A,U_1,U_2,U_3,U_4,U_5,U_6,U_7,U_8,U_9 !============================================================= FMT1='(AA(1X,ES11.4),BB(A0))' FMT2='(A___)' FMT3='(1X,ES11.4,__(12X),1X,ES11.4)' FMT4='(12X,__(1X,ES11.4))' FMT5='(AA(1X,ES11.4))' !============================================================= DO K=1,LYR L=LYR+1-K M=11-L WRITE(FMT1(2:3),'(I2)')L WRITE(FMT1(16:17),'(I2)')M CONTINUE 100 WRITE(109,FMT1)U(K),U(K+11),U(K+22),U(K+33),U(K+44),U(K+55) 1,U(K+66),U(K+77),U(K+88),U(K+99),U(K+110) END DO REWIND 109 DO K=1,LYR L=12*(LYR+1-K) M=LYR-1-K N=LYR+1-K WRITE(FMT2(3:5),'(I3)')L WRITE(FMT3(12:13),'(I2)')M WRITE(FMT4(6:7),'(I2)')M WRITE(FMT5(2:3),'(I2)')N CONTINUE READ(109,FMT2)Rline 101 READ(Rline,FMT3)U_B,U_A 102 READ(Rline,FMT4)U_1,U_2,U_3,U_4,U_5,U_6,U_7,U_8,U_9 103 WRITE(110,FMT5)U_B,U_1,U_2,U_3,U_4,U_5,U_6,U_7,U_8,U_9,U_A END DO !=============================================================== 在这里我没有列出计算的部分 我先把我要的值写成一个类似对角矩阵的形式 再分别以行去读值 101行是去读该一维阵列的两边边界值 102行则是去读取边界内包含的其他值 问题点在於说 当K=2时 102行的U_9是读不到的 以此类推 在往後的情况下102行的项目可以读到值的会越来越少 因此并无法去做计算及输出 想请问对於102/103行所需要读取/输出的项目个数是否可以做调整? 并且想知道说在102行读取到的值有哪些是可以用来做计算的 可以只用阵列的方法去达成我的目标吗? 另外，如果我的阵列大小有改变 在第100行的地方是否有较好的语法可以自动更改我要读取的项目数目? 或是我就只能回到程式码去修改? 感谢各位的教学与帮忙!! --

※ 发信站: 批踢踢实业坊(ptt.cc) ◆ From: 140.115.61.228