Fortran: вызов массива из подпрограммы

@Puma · Регистрация: 11.04.2016

Author24 — интернет-сервис помощи студентам

Доброе время суток! Мне нужно дописать старую программу для расчета некоторых величин. Программа состоит из 4х подпрограмм. В одной из них определяются члены массива PB(1000,1000), после определения его членов мне нужно приравнять его члены к членам массива РВ1(1000,1000) и использовать их в цикле другой подпрограммы в расчетной формуле. Проблема в том, что имея норм значения (0, 10, 20... и т.д.) в первой подпрограмме, во второй все члены массива равны нулю!!! Я неправильно задаю что-то в коде, или упускаю что-то важное? (честно, с фортраном знакома всего две недели)
первая подпрограмма я дописываю массив

Fortran

DIMENSION PB1(1000,1000)  
!далее, внутри цикла, где определаются значения массива РВ(I,J) пишу
PB1(I,J)=PB(I,J)  
!для проверки вывожу их значения  
WRITE(*,*)PB1(I,J)             
WRITE(*,*)'PB1(I,J)'            
 на этом с первой подпрограммой заканчиваю.
Вторая подпрограмма:
в COMMON /POL/PB1(1000,1000)  !прописывается последовательно на своем месте среди запрашиваемых массивов
!и далее в программе вставка формулы для расчета:
    DO 906 I=1,NS
    DO 907 J=1,NS
        PIZ1=RG/3.14 
    PIZ22(I)=PRR(I)/PIZ 
    PIZ2(I)=SQRT(PIZ22(I))
    IF (I.NE.J) THEN
    PIZ3(I,J)=PIZ2(I)*PB1(I,J) !Здесь PB1(I,J)=0 на всем цикле и все решение уходит в бесконечность((
    PIZ4(I,J)=1.12/PIZ3(I,J)
    PIZ5(I,J)=SQRT(PIZ4(I,J))
    PB(I,J)=PIZ1*PIZ5(I,J)
    WRITE(*,*) PB(I,J)
    END IF
907   CONTINUE 
906   CONTINUE

@Krasme · 11.04.2016, 19:49

Сообщение от Puma

COMMON /POL/PB1(1000,1000)

это должно быть в обеих подпрограммах

@Puma · 12.04.2016, 04:38 **[ТС]**

Сообщение от Krasme

это должно быть в обеих подпрограммах

Я тоже придерживалась такого мнения, но при отладке программа ругается.
Грешу на задание типа массива, но вроде бы выбран верно членами являются действительные числа. Задавала следующим образом:

Fortran

1 2	DIMENSION PB1(1000,1000) COMMON /POL/PB1(1000,1000)

@Krasme · 12.04.2016, 06:23

ваша конструкция дважды объявляет, что PB1 массив...

либо

Fortran

1	COMMON /POL/PB1(1000,1000)

либо

Fortran

1 2	DIMENSION PB1(1000,1000) COMMON /POL/PB1

@Puma · 12.04.2016, 18:33 **[ТС]**

Про второй вариант записи не знала, спасибо! попробую завтра проверить!

@Krasme 6930 / 4969 / 2098 Регистрация: 02.02.2014 Сообщений: 13,210
	11.04.2016, 19:49	2
	Сообщение от Puma COMMON /POL/PB1(1000,1000) это должно быть в обеих подпрограммах 1

@Puma 0 / 0 / 0 Регистрация: 11.04.2016 Сообщений: 7
	12.04.2016, 18:33 [ТС]	5
	Про второй вариант записи не знала, спасибо! попробую завтра проверить! 0

Fortran: вызов массива из подпрограммы

Решение