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
29
30
31
32
33
34
35
36
37
|
PROGRAM dfsortp
c sorts ODFs: for each nu there are 25 ODfs(nu,Pgas)
INTEGER*4 IDFOUT(7)
c INTEGER*2 IDF(12),ISTEPS(12,25,57)
INTEGER*2 IDF(12),ISTEPS(57,25,12)
EQUIVALENCE (IDFOUT(2),IDF(1))
dimension istepnu(1540,25,12),idfoutnu(1540,25)
open(unit=1,type='old',form='unformatted',
1 recordtype='fixed',blocksize=32000,recl=8)
do 5 it=1,1
DO 6 IP=1,25
DO 6 nu=1,1540
READ(1)IDFOUT
DO 3 ISTEP=1,12
ISTEPS(IT,IP,istep)=IDF(ISTEP)
istepnu(nu,ip,istep)=isteps(it,ip,istep)
3 continue
idfoutnu(nu,ip)=idfout(1)
c write(6,4),IDFOUT(1),(isteps(it,IP,istep),istep=1,12)
c write(6,66)(IDFOUT(I),I=1,7)
66 FORMAT(1X,7I10)
4 format(1x,i10,1X,12i8)
6 continue
5 continue
it=1
do 600 nu=1,1540
do 601 ip=1,25
do 602 istep=1,12
isteps(it,ip,istep)=istepnu(nu,ip,istep)
602 continue
write(2,22)idfoutnu(nu,ip),(isteps(it,ip,istep),istep=1,12)
c write(2)idfoutnu(nu,ip),(isteps(it,ip,istep),istep=1,12)
22 format(1x,i10,1x,12i8)
601 continue
600 continue
CALL EXIT
END
|
