現在(&lastmod;)作成中です。
既に書いている内容も&color(#ff0000){大幅に変わる};可能性が高いので注意。

-------

神戸大学 大学院システム情報学研究科 計算科学専攻 陰山 聡

-------

-------
【目次】
#contents
-------

前回の宿題は、「thermal_diffusion.f90を理解しておくこと」であったので、
その中身がわかっているものと仮定して今回は話を進める。


* 2次元可視化 [#k39ed7ce]
今回は2次元データの可視化を行う。
前回と同様gnuplotを利用する。

gnuplotで2次元データをプロットする場合、そのデータのフォーマットは
 x00 y00 関数値
 x01 y00 関数値
 x02 y00 関数値
  ・
  ・
  ・
 x09 y00 関数値
 (空行)
 x00 y01 関数値
 x01 y01 関数値
 x02 y01 関数値
  ・
  ・
  ・
 x09 y01 関数値
 (空行)
  ・
  ・
  ・
 x09 y09 関数値
である。


* 2次元データ出力ルーチン [#a2717789]

thermal_diffusion.f90による2次元データ出力ルーチン
temperature__output_2d_profile
は、正方形上(x,y平面上)に分布する温度をすべて書き出すものである。
中身を見てみよう。

  subroutine temperature__output_2d_profile
    real(DP), dimension(0:MESH_SIZE+1,    &
                        0:MESH_SIZE+1) :: prof
    integer                            :: counter = 0   ! saved
    integer                            :: ierr          ! use for MPI
    integer                            :: jstart_, jend_
    character(len=4)                   :: serial_num    ! put on file name
    character(len=*), parameter        :: base = "../data/temp.2d."
    integer :: i, j
    call set_jstart_and_jend(jstart_,jend_)
    write(serial_num,'(i4.4)') counter
    prof(:,:) = parallel__set_prof_2d(jstart_,jend_,temp(:,jstart_:jend_))
    if ( myrank==0 ) then
       open(10,file=base//serial_num)
       do j = 0 , MESH_SIZE+1
          do i = 0 , MESH_SIZE+1
             write(10,*) i, j, prof(i,j)
          end do
          write(10,*)' ' ! gnuplot requires a blank line here.
       end do
       close(10)
    end if
    counter = counter + 1
  end subroutine temperature__output_2d_profile



*  &color(#0000ff){【演習】}; [#n213fe1f] 

- thermal_diffusion.f90をmpif90コマンドでコンパイルし、
- ジョブスクリプトthermal_diffusion.jsで並列計算せよ。
- 連番つきファイルtemp.2d.????が出力されるはずである。
その中身は以下のようになっているはず。確認せよ。

           55           35   0.1165598588705999     
           56           35   9.9624877672293416E-002
           57           35   8.1734108631726782E-002
           58           35   6.2857224006520482E-002
           59           35   4.2963409431420671E-002
           60           35   2.2021479795155254E-002
           61           35    0.000000000000000     
  
            0           36    0.000000000000000     
            1           36   2.1867122785152873E-002
            2           36   4.2655284590767971E-002
            3           36   6.2396502500601705E-002
            4           36   8.1122529495226178E-002
           ・
           ・
           ・
           57           36   8.1122529495226178E-002
           58           36   6.2396502500601705E-002
           59           36   4.2655284590767971E-002
           60           36   2.1867122785152873E-002
           61           36    0.000000000000000     
  
            0           37    0.000000000000000     
            1           37   2.1680615643577435E-002
            2           37   4.2282992534786040E-002
            3           37   6.1839860798781864E-002
            4           37   8.0383668447763873E-002
            5           37   9.7946446090882752E-002
            6           37   0.1145596782417825     
           ・
           ・
           ・

* gnuplotによる2次元等高線の表示 [#o6645045]
静止画による説明

* &color(#0000ff){【演習】アニメーション}; [#n213fe1f] 

* gnuplotによる2次元カラー表示 [#n1e01e1c]
静止画による説明

* &color(#0000ff){【演習】アニメーション}; [#n213fe1f] 

* gnuplotによるsurface表示 [#v677383d]
静止画による説明

* &color(#0000ff){【演習】アニメーション}; [#n213fe1f] 

#ref()

* 熱源分布が一様でない場合の平衡温度分布 [#w71b98ed]

ソースコードの変更。 thermal_diffusion2.f90

* レポート課題 01 [#u26817b8]
thermal_diffusion02.f90で求めたチェッカーボード型熱源分布の平衡解(最終状態解)を可視化せよ。
温度の最大値と最小値も記せ。


* 授業アンケート [#qa65e3f2]
今回の演習内容はどうでしたか?(どれか一つ、一度だけ押してください。)
#vote(簡単すぎた[0], 難しすぎた[0], ちょうどよかった[0])

* 2次元領域分割による並列化 [#x20338f4]

* 質問、コメントなど自由にどうぞ [#k16ca111]

「お名前」欄は空欄で可。

#pcomment


------------------------------

as of &_now; (&counter;)