diff options
| author | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-02-16 12:40:45 -0500 |
|---|---|---|
| committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-02-16 12:40:45 -0500 |
| commit | 0373ffdfaaaa3845306ca71243d535fdffd941d4 (patch) | |
| tree | 194c3c278d7e352e39d555d31aae93c0be2dfc03 /dfsynthe/dfsortp.for | |
| parent | 01b51f73bd06b2d6eabb776ba6cc69e4abfaa0b3 (diff) | |
| download | kasym-0373ffdfaaaa3845306ca71243d535fdffd941d4.tar.gz | |
Initial commit
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 |