artjomp, код этот написан очень дотошно. Построчно, можно сказать. Если вы собираетесь его профессионально использовать, то лучше переписать под современный диалект. Например, смотрите как компактно записываются последние функции:
| Код |
subroutine D10(KP)
common /W3/ C(6,20), NP(20)
KP = 5
NP(1:KP) = [(i + 1, i = 1, KP)]
C(:,1:KP) = Reshape([&
11800., 8083., 6138., 4950., 4150.,&
11900., 8083., 6125., 4926., 4133.,&
12090., 8180., 6170., 4956., 4138.,&
12225., 8260., 6212., 4976., 4150.,&
12345., 8280., 6228., 4992., 4155.,&
12400., 8307., 6238., 4996., 4165.],&
[6, KP], [0, 0], [2, 1])
end subroutine D10
subroutine D12(KP)
common /W3/ C(6,20), NP(20)
KP = 8
NP(1:KP) = [(i + 1, i = 1, KP)]
C(:,1:KP) = Reshape([&
...
...],&
[6, KP], [0, 0], [2, 1])
end subroutine D12
subroutine D14(KP)
common /W3/ C(6,20), NP(20)
KP = 16
NP(1:KP) = [3, 4, 6, (i + 4, i = 4, KP)]
C(2:,1:KP) = Reshape([&
...
...],&
[5, KP], [0, 0], [2, 1])
end subroutine D14
subroutine D16(KP)
common /W3/ C(6,20), NP(20)
KP = 7
NP(1:KP) = [(i * 2 + 2, i = 1, KP)]
C(2:,1:KP) = Reshape([&
...
...],&
[5, KP], [0, 0], [2, 1])
end subroutine D16
subroutine D18(KP)
common /W3/ C(6,20), NP(20)
KP = 9
NP(1:KP) = [5, 6, (i * 2 + 2, i = 3, KP)]
C(2:,1:KP) = Reshape([&
...
...],&
[5, KP], [0, 0], [2, 1])
end subroutine D18
subroutine D20(KP)
common /W3/ C(6,20), NP(20)
KP = 10
NP(1:KP) = [(i * 2 + 4, i = 1, KP)]
C(2:,1:KP) = Reshape([&
...
...],&
[5, KP], [0, 0], [2, 1])
end subroutine D20
subroutine D22(KP)
common /W3/ C(6,20), NP(20)
KP = 9
NP(1:KP) = [(i * 2 + 4, i = 1, 4), (i * 2 + 6, i = 5, KP)]
C(2:,1:KP) = Reshape([&
...
...],&
[5, KP], [0, 0], [2, 1])
end subroutine D22
subroutine D25(KP)
common /W3/ C(6,20), NP(20)
KP = 8
NP(1:KP) = [10, 12, (i * 4 + 4, i = 3, 6), 30, 32]
C(2:,1:KP) = Reshape([&
...
...],&
[5, KP], [0, 0], [2, 1])
end subroutine D25
subroutine D28(KP)
common /W3/ C(6,20), NP(20)
KP = 8
NP(1:KP) = [(i * 4 + 8, i = 1, KP)]
C(2:,1:KP) = Reshape([&
...
...],&
[5, KP], [0, 0], [2, 1])
end subroutine D28
subroutine D32(KP)
common /W3/ C(6,20), NP(20)
KP = 9
NP(1:KP) = [(i * 4 + 12, i = 1, KP)]
C(2:,1:KP) = Reshape([&
...
...],&
[5, KP], [0, 0], [2, 1])
end subroutine D32
subroutine D6(KP)
common /W3/ C(6,20), NP(20)
KP = 9
end subroutine D6
subroutine D15(KP)
common /W3/ C(6,20), NP(20)
KP = 9
NP(1:KP) = [(i + 9, i = 1, KP)]
C(2,1:KP) = [1445., 1174., 990., 857., 756., 677., 613., 560., 515.]
end subroutine D15
|
Я вовсю использовал конструкторы массивов. Троеточия можете дополнить по образу и подобию. Первую часть программы в переделанном виде тоже показывал. Работать будет на любом современном компиляторе или на компиляторе, поддерживающем Fortran 90.
А теперь по делу.
Цитата(artjomp @ 3.2.2009, 16:05  ) | | Cr@$h она определена само собой |
- Если честно, похоже на массив OK, но вот только индексы у него через деление определяются.
- Возможно, эта функция была определена где-то ещё. Посмотрите файлы *.mod, например.
- Она точно не входит ни в один стандарт, поэтому является пользовательской.
- Уверен, нужно попробовать ещё раз изучить историю вопроса (возникновения этой программы).
|
Вы уверены, что никакой кусок программы при сращивании в файл не потеряли?
Мы все тут держимся. Заходите, будет интересно.