.
SUBROUTINE PLTOUT(FILENAME)
CHARACTER *18 FILENAME
CHARACTER *70 HEAD
INTEGER I,J,K
HEAD='VARIABLES="X","Y","U","V"'
OPEN(15,FILE=filename,STATUS='UNKNOWN')
WRITE(15,'(A70)') HEAD
WRITE(15,102) M,N
102 FORMAT('zone I=', I2, ' J=', I2, ' f=point')
DO 10 J=2,N+1
DO 10 I=2,M+1
WRITE(15,30) x,y,U,V !对应各网格点的x,y坐标及u,v速度分量
10 CONTINUE
30 FORMAT(2(f8.4,1x),2(1X,e11.4))
CLOSE(15)
RETURN
END |