質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.48%
FORTRAN

FORTRAN(フォートラン)は科学時術計算に向いた手続き型プログラミング言語です。 並列計算の最適化が行いやすい特性上、数値予報および気候モデルなどの大規模な計算を行う分野のスーパーコンピュータで使われています。

Q&A

解決済

1回答

4475閲覧

fortran:モジュールとメインプログラムで同じtxtをopenしたいができなかった。

Okuhira007

総合スコア13

FORTRAN

FORTRAN(フォートラン)は科学時術計算に向いた手続き型プログラミング言語です。 並列計算の最適化が行いやすい特性上、数値予報および気候モデルなどの大規模な計算を行う分野のスーパーコンピュータで使われています。

0グッド

0クリップ

投稿2021/01/23 02:36

編集2021/01/23 02:37

At line 13 of file metropolis2.f90 (unit = 11, file = 'data14-1.txt')
Fortran runtime error: End of file

Error termination. Backtrace:

Could not print backtrace: libbacktrace could not find executable to open
というエラーが出ています。
モヂュールとメインプログラムで1つのtxtファイル(data14-1)を読み込んだところ、エラーが出ました。

fortran

1コードmodule hh 2 implicit none 3contains 4 function h(kupper, x_c) result(value) 5 implicit none 6 integer i 7 real(8) kupper, x_c, value 8 real(8) :: x(0:199) 9 10 open(11, file='data14-1.txt', status='old') 11 value = 0.5*(x_c-2.0)*(x_c-2.0) 12 do i = 0, 199 13 read(11, *) x(i) 14 value = value + 0.5*kupper*(x(i)-x_c)*(x(i)-x_c) 15 enddo 16 close(11) 17 end function h 18end module hh 19 20program metropolis2 21 use hh 22 implicit none 23 integer i, j 24 real(8) :: epsi1(0:99), epsi2(0:99), x(0:199) 25 real(8) :: r(0:99) 26 real(8) :: h_prime, x_c = 1.0, kupper = 5.0 27 28 open(11, file='data14-1.txt', status='old') 29 open(12, file='data14-2.txt', status='old') 30 open(13, file='data14-3.txt', status='old') 31 open(14, file='data14-4.txt', status='old') 32 open(15, file='metropolis2.dat', status='replace') 33 do i = 0, 99 34 read(12,*) epsi1(i) 35 read(13,*) epsi2(i) 36 read(11,*) x(i) 37 x_c = x_c + epsi1(i) 38 kupper = kupper + epsi2(i) 39 h_prime = h(kupper, x_c) 40 read(14,*) r(i) 41 if (r(i) <= exp(h(kupper-epsi2(i), x_c-epsi1(i))-h_prime)) then 42 continue 43 else 44 x_c = x_c - epsi1(i) 45 kupper = kupper - epsi2(i) 46 endif 47 write(15, *) i, kupper, x_c 48 enddo 49 close(11) 50 close(12) 51 close(13) 52 close(14) 53 close(15) 54end program metropolis2

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

自己解決

x(i)をメインプログラムの方で消して、open(11)の11の装置関連を消すと動いたので、暫定ですが解決しました。

投稿2021/01/23 02:49

Okuhira007

総合スコア13

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.48%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問