diff options
Diffstat (limited to 'dfsynthe/dfsortp.for')
| -rw-r--r-- | dfsynthe/dfsortp.for | 37 |
1 files changed, 37 insertions, 0 deletions
diff --git a/dfsynthe/dfsortp.for b/dfsynthe/dfsortp.for new file mode 100644 index 0000000..ec00bb3 --- /dev/null +++ b/dfsynthe/dfsortp.for @@ -0,0 +1,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 |