From 0373ffdfaaaa3845306ca71243d535fdffd941d4 Mon Sep 17 00:00:00 2001
From: Joseph Hunkeler <jhunkeler@gmail.com>
Date: Mon, 16 Feb 2015 12:40:45 -0500
Subject: Initial commit

---
 dfsynthe/dfsortp.for | 37 +++++++++++++++++++++++++++++++++++++
 1 file changed, 37 insertions(+)
 create mode 100644 dfsynthe/dfsortp.for

(limited to 'dfsynthe/dfsortp.for')

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
-- 
cgit