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/separatedf.for | |
| parent | 01b51f73bd06b2d6eabb776ba6cc69e4abfaa0b3 (diff) | |
| download | kasym-0373ffdfaaaa3845306ca71243d535fdffd941d4.tar.gz | |
Initial commit
Diffstat (limited to 'dfsynthe/separatedf.for')
| -rw-r--r-- | dfsynthe/separatedf.for | 22 |
1 files changed, 22 insertions, 0 deletions
diff --git a/dfsynthe/separatedf.for b/dfsynthe/separatedf.for new file mode 100644 index 0000000..53ad80c --- /dev/null +++ b/dfsynthe/separatedf.for @@ -0,0 +1,22 @@ + program separatedf +C sums the ODFs computed for each T +C separates ODFs in BIG and LITTLE ODFs + integer*2 isteps(25,57,12) + open(unit=2,form='unformatted',status='new') + open(unit=3,form='unformatted',status='new') + do 1500 inu=1,1540 + do 150 it=1,57 + itape=it+10-1 + do 151 ip=1,25 + read(itape,*)idfout,(isteps(ip,it,istep),istep=1,12) +c write(6,22)idfout,(isteps(ip,it,istep),istep=1,12) + 22 format(1x,i10,1x,12i8) + 151 continue + if(inu.le.328)write(2)((isteps(ip,it,istep),istep=1,12),ip=1,25) + if(inu.gt.328)write(3)((isteps(ip,it,istep),istep=1,12),ip=1,25) + 150 continue + 1500 continue + close(unit=2) + close(unit=3) + stop + end |