Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- open(unit=46,file=fich)
- write(46,'(1A26)') '# vtk DataFile Version 2.0'
- write(46,*) 'Flow'
- write(46,*) 'ASCII'
- write(46,*) 'DATASET STRUCTURED_GRID'
- write(46,*) 'FIELD FieldData 1'
- write(46,*) 'TIME 1 1 double'
- write(46,*) time
- write(46,*) 'DIMENSIONS' ,nc+1,mc+1,1
- write(46,*) 'POINTS', (nc+1)*(mc+1), 'double'
- write(46,*) (( nd(n,m)%center(1), nd(n,m)%center(2) , 0.0d0, n=0,nc), m=0,mc)
- write(46,*) 'CELL_DATA', nc*mc
- write(46,*) 'SCALARS rho double'
- write(46,*) 'LOOKUP_TABLE default'
- write(46,*) (( sngl( cl(n,m)%w(1) ), n=1,nc), m=1,mc)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement