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

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

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

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

Q&A

解決済

1回答

996閲覧

実数と整数が混ざったデータを配列として読み込むときの宣言の方法が分からない(fortran)

tatsu31

総合スコア12

FORTRAN

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

0グッド

0クリップ

投稿2018/08/11 10:18

編集2018/08/12 10:18

海洋系の観測データの解析を行っています。元のデータである”160909000000_refL.out”というファイルには、一行に日付や時刻、観測値などの16のデータが羅列しており、それらが縦に約90万行ほど並んでいます。それらをある条件を満たした時の行だけ抽出したいです。
今回は配列としてデータを読み込み、行数をカウント、条件に当てはまる配列だけを”160909000000_refL.outZ”へ書き込むという方法でデータの抽出を行いたいです。
プログラムは以下になります。

fortran

1c refL 2 Program refL 3 Open (10,file='160909000000_refL.out',status='old') 4 Open (12,file='160909000000_refL.outZ',status='replace') 5 Integer z(5000) 6 ii=0 7 100 Continue 8 ii=ii+1 9 Read (10,*,end=200) z(ii) 10 print *,'hello1' 11 Go to 100 12 200 Continue 13 N=ii-1 14 do 300 i=1,N 15 if (z(25)==3 .and. z(26)==2) Then 16 print *,'hello2' 17 18 write(12,*) z(ii) 19 print *,'hello3' 20 end if 21 300 print *,z(ii) 22 Close (10) 23 Close (12) 24 STOP 25 END

元のデータは実数と整数が混じっているのでおそらく配列の宣言がおかしいのではないかと思うのですが、調べてもわかりません。
おかしいところと解決方法を教えていただけたらと思います。
ちなみにエラーは下のようになります。
イメージ説明

訂正後

fortran

1c refL 2 Program refL 3 Real z(5000) 4 Open (10,file='160909000000_refL.out',status='old') 5 Open (12,file='160909000000_refL.outZ',status='replace') 6 ii=0 7 100 Continue 8 ii=ii+1 9 Read (10,*,end=200) z(ii) 10 print *,z(10) 11 Go to 100 12 200 Continue 13 N=ii-1 14 do 300 i=1,N 15 if (z(10)==3 .and. z(11)==2) Then 16 print *,'hello2' 17 18 write(12,*) z(16) 19 print *,'hello3' 20 end if 21 300 print *,z(16) 22 Close (10) 23 Close (12) 24 STOP 25 END 26

このように訂正してsyntaxエラーは改善したのですが、下記のように同じ数字が羅列してエラーが起きているのですが、宣言の仕方が間違っているのでしょうか。
イメージ説明

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

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

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

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

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

guest

回答1

0

ベストアンサー

宣言文(Integer)の前に実行文(open)があるからですね。

Syntaxエラーなら、初心者でも、文法を確認すれば自分で必ず解決できます。

あと、実数と整数が混じった物を一つの配列に入れるなら、実数型にすると思いますが、なぜ整数型で宣言しているのでしょうか?

投稿2018/08/11 10:48

編集2018/08/11 10:50
otn

総合スコア84531

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

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

tatsu31

2018/08/11 11:07

先日に引き続き回答ありがとうございます。 >宣言文(Integer)の前に実行文(open)があるからですね。 わかりました。訂正します。 >あと、実数と整数が混じった物を一つの配列に入れるなら、実数型にすると思いますが、なぜ整数型で宣言しているのでしょうか? ごもっともなご指摘ありがとうございます。つまりrealで型宣言ということでよろしいですか。 少しとんちんかんな質問が多いとは思いますが、できる限り自分でできるように考えていき、少しずつ減らしていこうと思います。どうかよろしくお願いします。
tatsu31

2018/08/12 10:19

訂正後上のようなエラーになったのですが宣言の仕方がおかしいのでしょうか。
otn

2018/08/12 11:35

同じ変数を表示しているので同じ値が表示されます。それが意図通りでないなら、意図通りにプログラムを書きましょう。 SIGSEGVの原因として考えられるのは、データファイルが5000行以上あったとか。
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

ただいまの回答率
85.48%

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

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

質問する

関連した質問