-
Notifications
You must be signed in to change notification settings - Fork 2
Expand file tree
/
Copy pathplot3d.f
More file actions
executable file
·32 lines (27 loc) · 1001 Bytes
/
Copy pathplot3d.f
File metadata and controls
executable file
·32 lines (27 loc) · 1001 Bytes
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
! PLOT3D (2D FORMAT)
IXYZ = 20
ITEMP = 21
! WRITE OUT THE SOLUTION FOR POST-PROCESSING IN PLOT3D FORMAT
OPEN(UNIT=IXYZ,FILE='XYZ.DAT',FORM='FORMATTED')
WRITE(IXYZ,10) NBLOCKS
10 FORMAT(I10)
WRITE(IXYZ,20) (IMAX(N),JMAX(N),N=1,NBLOCKS)
20 FORMAT(10I10)
DO N = 1,NBLOCKS
WRITE(IXYZ,30) ((X(I,J,N),I=1,IMAX(N)),J=1,JMAX(N)),
& ((Y(I,J,N),I=1,IMAX(N)),J=1,JMAX(N))
ENDDO
30 FORMAT(10E15.8)
OPEN(UNIT=IFLO,FILE='TEMP.DAT',FORM='FORMATTED')
WRITE(IFLO,10) NBLOCKS
WRITE(IFLO,20) (IMAX(N),JMAX(N),N=1,NBLOCKS)
DO N = 1,NBLOCKS
WRITE(IFLO,30) TREF,DUM,DUM,DUM
WRITE(IFLO,30) ((T(N,I,J),I=1,IMAX(N)),J=1,JMAX(N)),
& ((T(N,I,J),I=1,IMAX(N)),J=1,JMAX(N)),
& ((T(N,I,J),I=1,IMAX(N)),J=1,JMAX(N)),
& ((T(N,I,J),I=1,IMAX(N)),J=1,JMAX(N))
ENDDO
! CLOSE THE SOLUTION FILES
CLOSE(IXYZ)
CLOSE(IFLO)