Fortranとgnuplotについて困っています

質問

Fortranとgnuplotについて困っています

こんにちは。

現在、Fortranプログラミングでgnuplotを使用した座標出力プログラムを考えているのですが、いまいちわかりません。

中心(0,0,0),半径1の級の球面上の座標を出力するプログラムを考えています。その際
・z>=0のときだけでよい
・gnuplotで確認できる
・Fortranのプログラムファイルでまとめる。

以上を満たすものを考えています。内容自体は基本的なことだとは思いますが、いまいち勝手がわからず困っています。

拡張子.f90でどのようなファイルを作ればいいのか教えていただければ幸いです。

回答No.3

neKo_deux · ベストアンサー

#2です。

途中で送信しちゃいました。
最後のサンプルプログラムに、x, y, zを計算する処理が抜けてます。
printする前に、

x=sin(th)*cos(phi)
y=sin(th)*sin(phi)
z=cos(th)

って計算が必要かも。

回答No.2

neKo_deux

> ・gnuplotで確認できる

のためには、x, y, zが空白区切りで出力されていればＯＫです。

座標の表現を、

同志社大学 工学部電気系 応用数学研究室 - 2005年度講義 - 解析学II - 2.23 3 次元空間の極座標
http://gandalf.doshisha.ac.jp/~kon/lectures/2005.calculus-II/html.dir/node31.html

に倣い、球面上の１点って事なら、

x = sinθcosφ
y = sinθsinφ
z = cosθ

で計算され、

> ・z>=0のときだけでよい

の条件ならば、

θは0～π/2までループ
φは0～2πまでループ
で、それぞれの状況でのx, y, zの値を出力します。

FORTRANはしばらく見ていないのでうろ覚えですが、ざっくり書くと、

double pi, th, phi, x, y, z
pi=3.14159265358979
do th=0.0, pi/4, pi/60
　do phi=0.0, pi*2, pi/30
　　print *, x, y, z
　end do
end do

とか？

回答No.1

Tacosan

「中心(0,0,0),半径1の級の球面上の座標」といわれてもそれがどういうものか全くわからないし「gnuplotで確認できる」が何を意味するのかさっぱり分からない.