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

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

新規登録して質問してみよう
ただいま回答率
85.46%
ファイル

ファイルとは、文字列に基づいた名前又はパスからアクセスすることができる、任意の情報のブロック又は情報を格納するためのリソースです。

データ構造

データ構造とは、データの集まりをコンピュータの中で効果的に扱うために、一定の形式に系統立てて格納する形式を指します。(配列/連想配列/木構造など)

FORTRAN

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

コードレビュー

コードレビューは、ソフトウェア開発の一工程で、 ソースコードの検査を行い、開発工程で見過ごされた誤りを検出する事で、 ソフトウェア品質を高めるためのものです。

C++

C++はC言語をもとにしてつくられた最もよく使われるマルチパラダイムプログラミング言語の1つです。オブジェクト指向、ジェネリック、命令型など広く対応しており、多目的に使用されています。

Q&A

0回答

1072閲覧

fortran read文の仕組み

MAEZAWAMiumi

総合スコア0

ファイル

ファイルとは、文字列に基づいた名前又はパスからアクセスすることができる、任意の情報のブロック又は情報を格納するためのリソースです。

データ構造

データ構造とは、データの集まりをコンピュータの中で効果的に扱うために、一定の形式に系統立てて格納する形式を指します。(配列/連想配列/木構造など)

FORTRAN

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

コードレビュー

コードレビューは、ソフトウェア開発の一工程で、 ソースコードの検査を行い、開発工程で見過ごされた誤りを検出する事で、 ソフトウェア品質を高めるためのものです。

C++

C++はC言語をもとにしてつくられた最もよく使われるマルチパラダイムプログラミング言語の1つです。オブジェクト指向、ジェネリック、命令型など広く対応しており、多目的に使用されています。

0グッド

1クリップ

投稿2020/08/23 07:10

質問内容

Fortran初心者です.
2点わからないことがあったのでお聞かせください.

質問.1
最下部に記載するコード内の後半部分

read(7) ( ni(m), nj(m), nk(m), m = 1, nblocks )

は2行目を読み取りそれぞれに格納だと思うのですが、m=1の部分にも数字の様なものが”3D.x”から読み取られ、入るのでしょうか?

質問.2、質問1でお聞きした行の2行下に記載されている

read(7) & ((( x(i,j,k,m), i=1,ni(m)), j=1,nj(m)), k=1,nk(m)), & ((( y(i,j,k,m), i=1,ni(m)), j=1,nj(m)), k=1,nk(m)), & ((( z(i,j,k,m), i=1,ni(m)), j=1,nj(m)), k=1,nk(m))

は装置番号7に設定したファイルからデータを読み込み、
それぞれのx,y,zに格納していくというものだと思うのですが、
その読み取る”3D.x"の中身はどのx,y,zと対応しているのか、

ni nj nk
x(1,1,1,1), y(2,1,1,1), z(3,1,1,1)
x(4,1,1,1), ..., z(ni,1,1,1)
x(1,2,1,1), y(2,2,1,1), z(3,2,1,1)

の形で読み取っていくのでしょうか、また3D.xもこれと同様の形態をとっているのですよね?
間違っていれば順番も教えていただけるとありがたいです.

以上です、よろしくお願いします.

該当のソースコード

parameter ( imax = 100 )
parameter ( jmax = 100 )
parameter ( kmax = 100 )
parameter ( nbmax = 10 )

integer i integer j integer m integer nblocks integer ni (nbmax) integer nj (nbmax) integer nk (nbmax) real x(imax,jmax,kmax,nbmax) real y(imax,jmax,kmax,nbmax) real z(imax,jmax,kmax,nbmax) open ( unit=7, form='unformatted', file='3D.x' ) read(7) nblocks read(7) ( ni(m), nj(m), nk(m), m = 1, nblocks ) do m = 1, nblocks read(7) & ((( x(i,j,k,m), i=1,ni(m)), j=1,nj(m)), k=1,nk(m)), & ((( y(i,j,k,m), i=1,ni(m)), j=1,nj(m)), k=1,nk(m)), & ((( z(i,j,k,m), i=1,ni(m)), j=1,nj(m)), k=1,nk(m)) enddo

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

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

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

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

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

guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだ回答がついていません

会員登録して回答してみよう

アカウントをお持ちの方は

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

ただいまの回答率
85.46%

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

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

質問する

関連した質問