このスレッドは、他のスレッドでは書き込めない超低レベル、
もしくは質問者自身何が何だか分からない質問を勇気を持って書き込むスレッドです。
FORTRAN使いが優しくコメントを返しますが、
お礼はFORTRANの布教と初心者の救済と次期Fortran2008規格でのCOMEFROM文採用をお願いします。
●注意事項
・質問する前にGoogle等の検索サイトで検索しましょう。
・回答者にわかりやすい様に、質問内容はできる限り詳しく書きましょう。
・エラーの場合は起きた状況、環境（OS・コンパイラ・バージョン）、エラーメッセージも詳しく書きましょう。
●前スレ
くだすれFORTRAN（超初心者用）その5
http://toro.2ch.net/test/read.cgi/tech/1269704830/
●過去スレ
くだすれFORTRAN（超初心者用）その4
http://pc12.2ch.net/test/read.cgi/tech/1232789521/
くだすれFORTRAN（超初心者用）その3
http://pc11.2ch.net/test/read.cgi/tech/1196384126/
くだすれFORTRAN（超初心者用）その2
http://pc11.2ch.net/test/read.cgi/tech/1164121236/
くだすれFORTRAN（超初心者用）
http://pc8.2ch.net/test/read.cgi/tech/1138063703/
探検
くだすれFORTRAN（超初心者用）その6
■ このスレッドは過去ログ倉庫に格納されています
2012/08/16(木) 20:58:00.31
2012/08/16(木) 21:32:02.53
, r‐‐ ､＿
, ´￣ ´￣ ￣ミ ､ ＿＿＿＿／: : : : ヽ、
__.／ ｀¨ ＼／: : : : : : : : :: : : : : : : : : : : }
/ ／// ／／ ／ ,.､ / / : : : : : : : : : : : : : : : ／
_,ｧ≫イ/ / / ／ <､＼>′ ヽ｛ : :＼ フ: : : : :／
ｲ/| //ｰ-､ />/｀Y ∧＼.:: : : :ヽ ／: : : :／ ） （
彡' ﾉ,ｲ{/{__／/／「/¨ヽ.! ! | .}} /: : :.i /: : : : :/ （ ）
彡' ﾉ }ｲ:::} ´´ r=ミ､ }ﾉ 1 小 /: : : : :/ （ ノ こ、これは>>1乙じゃなくて
片ｿ /:::しﾍ从 j/从.ゝ 彡/: : : : :/ |＼ レンコンしそ餃子の湯気なんだからね
｛:::::､ 廴::／}{/ ,'{/ ／ ,': : : : / （ （. |: : ＼ 今後は、気をつけなさい！！
ヽ ｰ ::::: ∨{彡'.:.:.:.:. :/i: : : : :{ （ ） |: : : :.ヽ
＼＿＿ イ≦_¨＼ ＼＿ {: : : : :丶＿＿＿＿＿__.ノ: : : : : }
､＿___xｰ|. |, -､｀ ヽ ＿.: '.: : : : : : : : : : : : : : : : : : : : : : ／＿
/ // / ー ／ ／》ｰ--､ l-‐ ＼: : : : : : : : : : : : : : : : : : イ ﾞ ヽ
ノ // /ヘ／ /// ＼ ／ ／,,r"i／ ,r"i／,,r"i／,,r"i ､ ヽ
/ // ,'⌒/. /// }/ ./ /#; / /#.; / //#; / ,/#; / .," i |
V // i / /// ｀ヽi. /#; / ,/#; ./ /#;; / ,/#; / .,' / /
片厶≦j/ヘ //〈 ／. |〃/ ,|〃/ ,|〃/ , |〃/ .,"..;;／ /
｀ー'⌒Y⌒ヽ // ヽ-=＝＝=< ヽ ヽゝ' `ヽゝ~ ヽゝ ~ヽゝ' " ／
ﾞ ー-- ､､... _ ＿＿,,,,, .... --―"
, ´￣ ´￣ ￣ミ ､ ＿＿＿＿／: : : : ヽ、
__.／ ｀¨ ＼／: : : : : : : : :: : : : : : : : : : : }
/ ／// ／／ ／ ,.､ / / : : : : : : : : : : : : : : : ／
_,ｧ≫イ/ / / ／ <､＼>′ ヽ｛ : :＼ フ: : : : :／
ｲ/| //ｰ-､ />/｀Y ∧＼.:: : : :ヽ ／: : : :／ ） （
彡' ﾉ,ｲ{/{__／/／「/¨ヽ.! ! | .}} /: : :.i /: : : : :/ （ ）
彡' ﾉ }ｲ:::} ´´ r=ミ､ }ﾉ 1 小 /: : : : :/ （ ノ こ、これは>>1乙じゃなくて
片ｿ /:::しﾍ从 j/从.ゝ 彡/: : : : :/ |＼ レンコンしそ餃子の湯気なんだからね
｛:::::､ 廴::／}{/ ,'{/ ／ ,': : : : / （ （. |: : ＼ 今後は、気をつけなさい！！
ヽ ｰ ::::: ∨{彡'.:.:.:.:. :/i: : : : :{ （ ） |: : : :.ヽ
＼＿＿ イ≦_¨＼ ＼＿ {: : : : :丶＿＿＿＿＿__.ノ: : : : : }
､＿___xｰ|. |, -､｀ ヽ ＿.: '.: : : : : : : : : : : : : : : : : : : : : : ／＿
/ // / ー ／ ／》ｰ--､ l-‐ ＼: : : : : : : : : : : : : : : : : : イ ﾞ ヽ
ノ // /ヘ／ /// ＼ ／ ／,,r"i／ ,r"i／,,r"i／,,r"i ､ ヽ
/ // ,'⌒/. /// }/ ./ /#; / /#.; / //#; / ,/#; / .," i |
V // i / /// ｀ヽi. /#; / ,/#; ./ /#;; / ,/#; / .,' / /
片厶≦j/ヘ //〈 ／. |〃/ ,|〃/ ,|〃/ , |〃/ .,"..;;／ /
｀ー'⌒Y⌒ヽ // ヽ-=＝＝=< ヽ ヽゝ' `ヽゝ~ ヽゝ ~ヽゝ' " ／
ﾞ ー-- ､､... _ ＿＿,,,,, .... --―"
2012/08/17(金) 02:14:17.73
> 初心者に適切な指導もなくプログラム書かせているのはアカハラの域に達してるな。
これって割と多いような気がするな。
まあ分野によっちゃいずれ必要になることが前前から分かっていそうな気もするのだが
それも指導のうちと言えばそうだ
これって割と多いような気がするな。
まあ分野によっちゃいずれ必要になることが前前から分かっていそうな気もするのだが
それも指導のうちと言えばそうだ
4デフォルトの名無しさん
2012/08/17(金) 04:02:00.57 do i=2,1000
write(*,*) '1otsu'
enddo
write(*,*) '1otsu'
enddo
2012/08/17(金) 19:12:19.36
>>3
昔は先輩が教えることが前提となっていたから成立したが、
最近、先輩がいなかったり、知らなかったりで、伝承が途切れている。
時代の変化に、サヨクかぶれの爺さんたちの現状認識が追いついていない。
連中は鳩山みたいに夢を見て、管みたいに怒鳴り散らしてるからｗ
昔は先輩が教えることが前提となっていたから成立したが、
最近、先輩がいなかったり、知らなかったりで、伝承が途切れている。
時代の変化に、サヨクかぶれの爺さんたちの現状認識が追いついていない。
連中は鳩山みたいに夢を見て、管みたいに怒鳴り散らしてるからｗ
2012/08/17(金) 20:00:57.70
前スレくだすれFORTRAN（超初心者用）その5の>>999さん
アドバイスありがとうございます。引き続きよろしくお願いします。
変更点の件ですが、
入力文
character (len = 20) :: filenm
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm
open(11,file= TRIM(filenm) )
出力文
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm
open(14,file=filenm)
close(14)
という修正をくわえればよろしいでしょうか？
アドバイスありがとうございます。引き続きよろしくお願いします。
変更点の件ですが、
入力文
character (len = 20) :: filenm
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm
open(11,file= TRIM(filenm) )
出力文
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm
open(14,file=filenm)
close(14)
という修正をくわえればよろしいでしょうか？
2012/08/17(金) 20:51:34.50
>>6
それだと出力してない。
OPEN と CLOSE の間で出力する。
ただ一度に欲張らず、まず固定された入力ファイルと出力ファイルで
ちゃんと計算できるプログラムをつくって確認しておくべき。
その後、色々ループさせたり、ファイル名を変えたり出来るようにすべき。
ものには順番があり、近くから遠くへ及ぼすべしとは、四書のはじめ大学にも書いてあるｗ
それだと出力してない。
OPEN と CLOSE の間で出力する。
ただ一度に欲張らず、まず固定された入力ファイルと出力ファイルで
ちゃんと計算できるプログラムをつくって確認しておくべき。
その後、色々ループさせたり、ファイル名を変えたり出来るようにすべき。
ものには順番があり、近くから遠くへ及ぼすべしとは、四書のはじめ大学にも書いてあるｗ
2012/08/17(金) 21:14:43.76
これまでの教えていただいたプログラムが下記になります。
program test
implicit none
integer, parameter :: ndata = 400
real :: a(18, ndata), x(3, ndata), y(3, ndata)
real :: s, sx, sy, tx(ndata), ty(ndata)
integer :: i, j, icol
character (len = 20) :: filenm
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm
open(11,file= TRIM(filenm) )
open(12, file = 'input2.txt')
open(13, file = 'input3.txt')
!
do i = 1, 6 ! skip 6 lines
read(11, *)
end do
do i = 1, ndata
read(11, *) (a(j, i), j = 1, 18)
end do
! exel data x
do i = 1, ndata
read(12, *) (x(j, i), j = 1, 3)
end do
program test
implicit none
integer, parameter :: ndata = 400
real :: a(18, ndata), x(3, ndata), y(3, ndata)
real :: s, sx, sy, tx(ndata), ty(ndata)
integer :: i, j, icol
character (len = 20) :: filenm
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm
open(11,file= TRIM(filenm) )
open(12, file = 'input2.txt')
open(13, file = 'input3.txt')
!
do i = 1, 6 ! skip 6 lines
read(11, *)
end do
do i = 1, ndata
read(11, *) (a(j, i), j = 1, 18)
end do
! exel data x
do i = 1, ndata
read(12, *) (x(j, i), j = 1, 3)
end do
2012/08/17(金) 21:16:18.67
! exel data y
do i = 1, ndata
read(12, *) (y(j, i), j = 1, 3)
end do
!
do i = 1, ndata
tx(i) = (i * 50.0 - 25.0) * x(3, i)
ty(i) = (i * 50.0 - 25.0) * y(3, i)
end do
print *, 'input column '
read *, icol
sx = sum( tx * a(icol, :) )
sy = sum( ty * a(icol, :) )
s = sx / sy
write(*,'(a,$)') 'File Name : '
stop
end program test
do i = 1, ndata
read(12, *) (y(j, i), j = 1, 3)
end do
!
do i = 1, ndata
tx(i) = (i * 50.0 - 25.0) * x(3, i)
ty(i) = (i * 50.0 - 25.0) * y(3, i)
end do
print *, 'input column '
read *, icol
sx = sum( tx * a(icol, :) )
sy = sum( ty * a(icol, :) )
s = sx / sy
write(*,'(a,$)') 'File Name : '
stop
end program test
2012/08/17(金) 21:38:35.93
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm
open(14,file=filenm)
write(*,*) 'result=', s ! <--だから結果を書きだせってｗ
close(14)
read(*,'(a)') filenm
open(14,file=filenm)
write(*,*) 'result=', s ! <--だから結果を書きだせってｗ
close(14)
2012/08/17(金) 21:39:27.86
write(14,*) 'result=', s ! <--だから結果を書きだせってｗ
思わず間違えたｗｗｗｗ
思わず間違えたｗｗｗｗ
2012/08/17(金) 21:52:24.90
2012/08/17(金) 22:03:15.16
そうそう。欲しい結果を書きだす。
本来は、入力ファイル名や計算したカラム（列）番号のような
付加的な情報も書いておくべきだろうが・・・
open(14,file = trim(filenm) ) ! <- こっちもTRIMで尻の空白を削っておくべし。
本来は、入力ファイル名や計算したカラム（列）番号のような
付加的な情報も書いておくべきだろうが・・・
open(14,file = trim(filenm) ) ! <- こっちもTRIMで尻の空白を削っておくべし。
1412
2012/08/17(金) 22:16:41.91 >>13
最後の方は、こういう感じでプログラムすればうまく動きますかね。
print *, 'input column '
read *, icol
sx = sum( tx * a(icol, :) )
sy = sum( ty * a(icol, :) )
s = sx / sy
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm
open(14,file = trim(filenm) ) ! <- TRIMで尻の空白を削る。
write(14,*) 'result=', s ! <--結果を書きだす。
close(14)
stop
end program test
こちらの掲示板に来て良かったです。
色々こんなデキの悪い自分に付き合ってくださってありがとうございました。
最後の方は、こういう感じでプログラムすればうまく動きますかね。
print *, 'input column '
read *, icol
sx = sum( tx * a(icol, :) )
sy = sum( ty * a(icol, :) )
s = sx / sy
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm
open(14,file = trim(filenm) ) ! <- TRIMで尻の空白を削る。
write(14,*) 'result=', s ! <--結果を書きだす。
close(14)
stop
end program test
こちらの掲示板に来て良かったです。
色々こんなデキの悪い自分に付き合ってくださってありがとうございました。
2012/08/17(金) 23:15:49.56
character (len = 20) :: filenm_in, filenm_out
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm_in
open(11,file= TRIM(filenm_in) )
・・・・
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm_out
open(14,file=trim(filenm_out))
write(14, *) 'input file name', filenm_in
write(14, *) 'column no.', icol
write(14,*) 'result=', s ! <--結果を書き出すプログラム。
close(14)
として、読み込んだファイルと、使ったコラムも出力ファイルに書きだしておくと、
あとで何の計算したのかヒントになっていいぞ。
とろくさいと、どのファイルに何を書きだしたかすぐ忘れるだろうから、こうしておけｗｗｗ
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm_in
open(11,file= TRIM(filenm_in) )
・・・・
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm_out
open(14,file=trim(filenm_out))
write(14, *) 'input file name', filenm_in
write(14, *) 'column no.', icol
write(14,*) 'result=', s ! <--結果を書き出すプログラム。
close(14)
として、読み込んだファイルと、使ったコラムも出力ファイルに書きだしておくと、
あとで何の計算したのかヒントになっていいぞ。
とろくさいと、どのファイルに何を書きだしたかすぐ忘れるだろうから、こうしておけｗｗｗ
1612
2012/08/17(金) 23:31:33.54 >>15
度々すみません。親切にありがとうございます。
プログラムに組み込ませていただきます！
今日もかれこれ5時間くらいネカフェにいますが、
そろそろおいとまします。自転車こいで橋を越えた甲斐がありました。
また色々教えてください。
度々すみません。親切にありがとうございます。
プログラムに組み込ませていただきます！
今日もかれこれ5時間くらいネカフェにいますが、
そろそろおいとまします。自転車こいで橋を越えた甲斐がありました。
また色々教えてください。
1712
2012/08/17(金) 23:46:42.05 最後にプログラム全体を書き込ませてください。
program test
implicit none
integer, parameter :: ndata = 400
real :: a(18, ndata), x(3, ndata), y(3, ndata)
real :: s, sx, sy, tx(ndata), ty(ndata)
integer :: i, j, icol
character (len = 20) :: filenm_in, filenm_out
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm_in
open(11,file= TRIM(filenm_in) )
open(12, file = 'input2.txt')
open(13, file = 'input3.txt')
!
do i = 1, 6 ! skip 6 lines
read(11, *)
end do
do i = 1, ndata
read(11, *) (a(j, i), j = 1, 18)
end do
! exel data x
do i = 1, ndata
read(12, *) (x(j, i), j = 1, 3)
end do
! exel data y
do i = 1, ndata
read(12, *) (y(j, i), j = 1, 3)
end do
program test
implicit none
integer, parameter :: ndata = 400
real :: a(18, ndata), x(3, ndata), y(3, ndata)
real :: s, sx, sy, tx(ndata), ty(ndata)
integer :: i, j, icol
character (len = 20) :: filenm_in, filenm_out
write(*,'(a,$)') 'File Name : '
read(*,'(a)') filenm_in
open(11,file= TRIM(filenm_in) )
open(12, file = 'input2.txt')
open(13, file = 'input3.txt')
!
do i = 1, 6 ! skip 6 lines
read(11, *)
end do
do i = 1, ndata
read(11, *) (a(j, i), j = 1, 18)
end do
! exel data x
do i = 1, ndata
read(12, *) (x(j, i), j = 1, 3)
end do
! exel data y
do i = 1, ndata
read(12, *) (y(j, i), j = 1, 3)
end do
■ このスレッドは過去ログ倉庫に格納されています
ニュース
- 中国「国連安保理の許可なしに日本攻撃可能」 Xで旧敵国条項に言及… ★5 [BFU★]
- 中国「国連安保理の許可なしに日本攻撃可能」 Xで旧敵国条項に言及… ★6 [BFU★]
- 【サッカー】U-17W杯 日本代表、無念のベスト8敗退… チャンスは多く作ったが仕留め切れず オーストリアに0－1で敗戦 [冬月記者★]
- 【千葉】コンビニに尿入りペットボトル並べた疑い、２６歳男「むしゃくしゃして」…購入した客が飲もうとしたところ臭いに違和感 [ぐれ★]
- 植田日銀総裁 「円安進行が物価高を起こしている」 ★4 [お断り★]
- 中国官製報道「日本経済はもう持たない」にネット民ツッコミ「ニュースだけ見てたら日本はもう百回くらい爆発してる」 [１ゲットロボ★]
- 現役JKのお茶会スレ( ¨̮ )︎︎𖠚ᐝ157
- 中国大使館、「高市早苗の正体」を完璧に絵にしてしまう。こら才能あるでぇ！ [592058334]
- 自分、馬鹿なのでおっしゃる意味がわかりません
- WTO世界のコメ🌾価格は記録的な豊作により1年で35%下落(5キロで200円程度)と発表※日本は1年で3倍値上がり [709039863]
- 【ぺこ専🐰】なんG 兎田ぺこら突発配信実況スレ🏡【ホロライブ▶】
- NHKニュースウオッチ9「日本側は対話にオープンな姿勢で安定した日中関係を築きたい考えなのに中国が意固地で糸口が見いだせない」 [904151406]