c------------------------------------------------------------------------------ subroutine chk include 'parameter.par' #include "precision.com" include 'common.com' character*80 fmt c------------------------------------------------------------------------------- c ... writeout open(unit=50,file='in_T.out',status='unknown') write(50,'(''# x, v, T, rho, X_1, X_2, ... X_N '')') write(50,'(''# number of cells in y : '',i9)')n write(50,'(''# N = '',i9)')nsp write(fmt,'(a,i3,a)') '(', nsp+4, 'e22.14)' do j=1,n write(50,fmt)x(j),v(j),temp(j),rho(j), + (xs(j,k),k=1,nsp) end do close(50) open(unit=50,file='vht_T.out',status='unknown') write(50,'(''# x, T, visc_T, htcond_T '')') write(50,'(''# number of cells in y : '',i9)')n write(50,'(''# N = '',i9)')nsp do j=1,n write(50,'(4e22.14)')x(j),temp(j),amu(j),htcond(j) end do close(50) open(unit=50,file='cpsp_T.out',status='unknown') write(50,'(''# x, T, cp_1_T, cp_2_T, ... , cp_nsp_T '')') write(50,'(''# number of cells in y : '',i9)')n write(50,'(''# N = '',i9)')nsp write(fmt,'(a,i3,a)') '(', nsp+2, 'e22.14)' do j=1,n write(50,fmt)x(j),temp(j),(cpsp(j,k),k=1,nsp) end do close(50) open(unit=50,file='D_T.out',status='unknown') write(50,'(''# x, T, D_1_nsp_T, D_2_nsp_T,..., D_nsp_nsp_T'')') write(50,'(''# number of cells in y : '',i9)')n write(50,'(''# N = '',i9)')nsp write(fmt,'(a,i3,a)') '(', nsp+2, 'e22.14)' do j=1,n write(50,fmt)x(j),temp(j),(capd(j,k),k=1,nsp) end do close(50) open(unit=50,file='h_T.out',status='unknown') write(50,'(''# x, T, h_1_T, h_2_T, ... , h_nsp_T '')') write(50,'(''# number of cells in y : '',i9)')n write(50,'(''# N = '',i9)')nsp write(fmt,'(a,i3,a)') '(', nsp+2, 'e22.14)' do j=1,n write(50,fmt)x(j),temp(j),(hsp(j,k),k=1,nsp) end do close(50) c------------------------------------------------------------------------------ return end