|
| |
|
Главная Журналы COME HERE TO CREATE DATA FOR PLOTTING С KFL0W=2 WRITE(8,300) HEADER 300 FORMAT(A64) WRITE(8,310) KFLOW,Ll,Ml,NFMAX,MODE,(KPLOT(I),1=1,NFMAX) 310 FORMAT(1815) IBLOK=0 DO 320 J=2,M2 DO 320 1=2,L2 IF(IBLOCK(I,J).EQ.l) THEN IBL0K=1 GO TO 330 ENDIF 320 CONTINUE 330 CONTINUE WRITE(8,310) IBLCK WRITE(8,34 0) (TITLE(N),N=1,NFMAX) 34 0 FORMAT(4A18) WRITE(8,350) (X(I),1=1,Ll),(Y(J),J=l,Ml),(XU(I),1=2,Ll) 1,(YV(J),J=2,M1),(R(J),J=1,M1) 350 FORMAT(5E12.6> DO 360 N=1,NFMAX IF(KPLOT(N).NE.O) WRITE(8,350) ((F(I,J,N),1=1,Ll),J=l,Ml) 360 CONTINUE IF(IBLOK.EQ.l) THEN WRITE(8,310) ((IBLOCK(I,J),1=1,Ll),J=1,M1) ENDIF CLOSE(8) RETURN CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC SUBROUTINE VALUES CREATE A FACILITY TO ASSIGN VALUES TO REAL VARIABLES С ENTRY DATA9(Al,CI,A2,C2,A3,C3,A4,C4,A5,C5,Аб,C6,Al,07, A8,C8,A9,C9) A9=C9 ENTRY DATA8(Al,CI,A2,C2,A3,C3,A4,C4,A5,C5,Аб,C6, A7,C7,A8,C8) A8=C8 ENTRY DATA7(Al, CI,A2,C2, A3,C3,A4,C4,A5,C5,Аб,Сб,A7,C7) A7=C7 ENTRY DATA6(A1,C1,A2,C2,A3,C3,A4,C4,A5,C5,A6,C6) A6=C6 ENTRY DATA5(Al,CI,A2,02,A3,C3,A4,C4,A5,05) А5=С5 ENTRY DATA4(AI,Cl,А2,С2,A3,СЗ,А4,С4) А4 = С4 ENTRY DATA3(AI,Cl,А2,С2,A3,СЗ) АЗ = СЗ ENTRY DATA2(AI,Cl,А2,02) A2=C2 ENTRY DATAl(AI,Cl) A1 = C1 RETURN CREATE A FACILITY TO ASSIGN VALUES TO INTEGER VARIABLES С ENTRY IMTA9 (II, Jl,I2,J2,I3,J3,I4,J4,I5,J5,I6,J6,n, J7, 18, J8, 19,J9) I9 = J9 ENTRY INTA8(I1,J1,I2,J2,I3,J3,I4,J4,I5,J5,I6,J6,I7, J7, IB, J8) 18 = J8 ENTRY INTA7(I1,J1,I2,J2,I3,J3,I4,J4,I5,J5,I6,J6,I7,J7) I7 = J7 ENTRY INTA6(I1,J1,I2,J2,I3,J3,I4,J4,I5,J5,I6,J6) I6 = J6 ENTRY INTA5(II,J1,I2,J2,I3,J3,I4,J4,I5,J5) I5=J5 ENTRY INTA4 (11,Jl,12,J2,I 3,J3,14,J4) I4 = J4 ENTRY INTA3(II,Jl,12,J2,13,J3) I3 = J3 ENTRY INTA2(II,Jl,12,J2) I2 = J2 ENTRY INTAl(II,Jl) I1 = J1 RETURN CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC СПИСОК ИМЕН НА ФОРТРАНе Приведем список, куда включены имена только тех переменных, значение которых существенно. Другие переменные, которые используются для временного хранения величин, в список не входят, но их смысл легко понять из контекста. а1м(1,j) ainr aip(i,j) ajm(i,j) ajp(i,j) alam (i, j) an в ap(i,j) apt area arx(i,j) beta con (i, j) crit(nf) - коэффициент o; - инерция / [см. (5.64)]; - коэффициент О;.; - коэффициент а; - коэффициент Од,; - коэффициент X [см. (3.4)]; - сумма соседних коэффициентов; - коэффициент ар, массив, эквивалентный Sj,; - коэффициент Ор для нестационарного члена [см. (5.20)]; - площадь грани контрольного объема; - площадь грани контрольного объема, перпендикулярной оси х; - коэффициент в выражениях для аппроксимации на границах; - очень большое число, например 1. е20; - коэффициент, используемый в блочной коррекции; - то же; - » - » - постоянный член Ь; массив, эквивалентный S.; - константа для критерия сходимости, используемого в подпрограмме solve; □ iff - диффузионная проводимость D; dt - шаг по времени At; f(i,j,nf) - различные функции ф; fluxil (j,nf) - плотность диффузионного потока через левую границу; fluxjl (i,nf) - то же через нижнюю границу; 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 [ 96 ] 97 98 99 |