テーマ:Fortran90

やっとできた

 gnuplotのコマンドが間違っていました。ubuntu日本語フォーラムに投稿してようやく理解しました。 計算結果 case1 Fr=0.5 case2 Fr=1.5 case3 Fr=1.2
トラックバック:0
コメント:0

続きを読むread more

gnuplotのバグではないの

 gnuplotのWindows版をインストールして同じファイルを図化してみた。画面上にグラフが表示されたが、保存はできなかった。そうなると読み込んだファイルがどこか変ということになる。計算条件を変えて、結果をグラフにしてみたが同じだ。 直接クリップボードにコピーしてWordに貼りつけて、グラフの比較を行った。
トラックバック:0
コメント:0

続きを読むread more

プログラム完成

 何とか完成にこぎつけた。中間データをファイルに出力して読み込むという工程を入れたのが不満だが、それで計算ができるようになった。  しかしgnuplotでグラフを保存できないというエラーは致命的だ。
トラックバック:0
コメント:0

続きを読むread more

もう一歩のところで進まない

 一次元浅水流方程式(波動方程式)の数値計算がうまくいかない。実行可能ファイルは作成できるが、計算結果がおかしい。途中の変数を出力してみるのだが、プログラムの改善点を見つけるには至らない。  以下がプログラムだ。 module globals !implicit none integer,save::…
トラックバック:0
コメント:0

続きを読むread more

デバッグ終わらず

 波動方程式を計算するプログラムを作成しましたが、デバッグが終了しません。 !module globals module globals implicit none integer,save:: n,itrmax,pintv,i real(8),save:: dt,xl,fr,g…
トラックバック:0
コメント:0

続きを読むread more

HPに写真を追加

 HPのプログラミングの項目にルンゲクッタ法を載せてあったが、昨日完成させたプログラムリストをグラフとともに載せた。
トラックバック:0
コメント:0

続きを読むread more

Scilabで答えを確認した

 Scilabでルンゲクッタ法の答えが正しいかどうか確認した。正解だった。 1.00 1.000 1.10 0.995 1.20 0.979 1.30 0.951 1.40 0.913 1.50 0.863 1.60 0.803 1.70 0.735 1.80 0.662 1.90 0.586…
トラックバック:0
コメント:0

続きを読むread more

ルンゲクッタ法修正版

 昨日のルンゲクッタ法プログラムは未完成だった。 program main implicit none real ystart,xstart,xend,h real y,x,f real k0,k1,k2,k3 real f0,f1,f2,f3 charact…
トラックバック:0
コメント:0

続きを読むread more

ルンゲクッタ法

 4次のルンゲクッタ法の計算をFortran90で作成したのだが、正解かどうかなんだか自信が持てない。ソースコードを画像ファイルとして添付しようとしたが、コピーできなかった。
トラックバック:0
コメント:0

続きを読むread more

研究室のプロジェクトに格上げ

 北海道の寒地土木研究所寒地河川チームのプロジェクトに、現場のための水理学を書き直すというプロジェクトが加わった。私の村上さんへの働きかけが大きかったらしい。研究室として、ソースを書き直すという作業に取り組むという。完成はいつ頃になるだろう、楽しみだ。
トラックバック:0
コメント:0

続きを読むread more

同じような現象に出会う

 Fortranで作成したプログラムと同じ処理をするプログラムをC言語で作成した。デバッグしてコンパイルエラーがなくなったのだが、出力が変だ。Fortranと同じ箇所の関数で処理がストップしてしまう。原因不明なので、今日は一日寝ていた。
トラックバック:0
コメント:0

続きを読むread more

今朝ベッドが壊れまして

 今朝7時前にベッドが壊れ、起きざるを得なくなった。出勤前の父親に話すと、今から修理するから  午前中会社を休み、溶接してさび止めまで塗って修理してくれた。  今日は、北海道のM氏が作成したプログラムのデバッグをおこなった。  gdbというデバッガーを使うとエラー箇所の特定が簡単にできる。しかし、計算結果がおかしいとい…
トラックバック:0
コメント:0

続きを読むread more

RIC-Naysで浮遊砂の河床変動計算

 Fortran90/95で作成したプログラムはどういう訳か、正しく計算できない。そのことは北海道のM氏にも伝えた。同氏も正しく計算できないことを確認。  そこでいったんFortranを離れて、RIC-Naysで計算することにした。7時過ぎに計算開始したので、終了は何時間後になるかな。 数値計算&可視化ソフトYorick楽…
トラックバック:0
コメント:0

続きを読むread more

河床変動計算(浮遊砂)ではまる

 浮遊砂の河床変動計算で、はまってしまいました。にっちもさっちも行かず、おかしな出力になってしまいました。北海道の研究者に質問をぶつけてみたところです。自力でのデバッグも限界です。
トラックバック:0
コメント:0

続きを読むread more

順列と組み合わせの計算

 2つの数(n,r)を入力し、順列nPrと組み合わせnCrを計算するプログラム。 program junretsukumiawase integer n,r,w1,w2,w3,i,j,k write(*,*)'input n(n<=10):' read(*,*)n write(*,*)'input r…
トラックバック:0
コメント:0

続きを読むread more

Fortran入門書を購入しました

 森北出版 数値計算のためのFortran90/95プログラミング入門 牛島 省著をAmazonで購入し、今日届きました。 数値計算のためのFortran 90/95プログラミング入門楽天ブックス牛島省 森北出版発行年月:2007年07月 ページ数:244p サイズ:単行本 ISBN:97846楽天市場 by JIS FO…
トラックバック:0
コメント:0

続きを読むread more