program contdata
      parameter(ixx=501,jxx=501,ijxx=ixx*jxx)
      real x(ixx),y(jxx)
      real cp(ixx,jxx)
      data cp/ijxx*0.0/
c*****files
      open(unit=28,file='tsd.cpcon',status='old')
      open(unit=29,file='cp.dat',status='replace')
c*****constants
      read(28,*)ix,jx
      read(28,*)(x(i),i=1,ix)
      read(28,*)(y(j),j=1,jx)
      do 1 j=1,jx
         read(28,*)(cp(i,j),i=1,ix)
         write(29,100,advance='yes')(cp(i,j),i=1,ix)
 1    continue
      write(*,*)' cp.dat is generated from tsd.cpcon'
	write(*,*)' To plot Cp countor, put 3 input files:'
	write(*,*)'       cp.dat '
	write(*,*)'       tsd.ix'
	write(*,*)'       tsd.jy'
	write(*,*)'   in the same folder and execute the MATLAB file:'
	write(*,*)'        cpcont.m'
	write(*,*)
	write(*,*)'Press Enter key to close this window.'
	read(*,*)
  100 FORMAT(500(1x,E14.7))
      end