namihuzi - @zigen 's note

c =====================================================================
c　波裏富士のシミュレーション
c yr:density　 yu(or yv):x(or y)velocity　 yp:pressure
c grx(or gry):differential of x(or y) density
c
c purpose　 : Output Graphic on File
c
c date　　　 : 2006 Oct. 19
c ---------------------------------------------------------------------
　　 subroutine graph

　　 include "dim.h"
　　 include "namihuzi.h"

c　　 dimension f(0:im,0:jm)
　　 character*1 num(0:9)
　 data num /'0','1','2','3','4','5','6','7','8','9'/
　 character*8 filename

　　 filename='nami'//num(mod(istep/1000,10))//num(mod(istep/100,10))
　　 &　　　　　　　 //num(mod(istep/10,10))//num(mod(istep,10))
　　 open(10,file=filename,status='unknown')

　　 do 100 i=1, nx-1

　　 do 110 k=1,3

　　 do 120 j=1, ny-1
120 write(10,1000) yrn(i,j)
1000 format(f5.2,$)
110 continue

　　 write(10,1010)
1010 format(1h )

100 continue
　　 close(10)
　　 return
　　 end