diff --git a/Makefile.in b/Makefile.in index 91f9450c2..16bdf431b 100644 --- a/Makefile.in +++ b/Makefile.in @@ -65,7 +65,7 @@ SRCS2F77 = indexx.f gen65.f chkmsg.f \ packmsg.f packtext.f setup65.f unpackcall.f unpackgrid.f \ unpackmsg.f unpacktext.f astro.f azdist.f coord.f dcoord.f \ deg2grid.f dot.f ftsky.f geocentric.f GeoDist.f grid2deg.f \ - moon2.f MoonDop.f sun.f toxyz.f pfxdump.f \ + moon2.f MoonDop.f sun.f toxyz.f pfxdump.f f77_wisdom.f \ ccf65.f trimlist.f chkhist.f decode1a.f \ filbig.f fil6521.f twkfreq.f decode65b.f \ afc65b.f fchisq.f ccf2.f rfile2.f diff --git a/f77_wisdom.f b/f77_wisdom.f new file mode 100755 index 000000000..b0a72585e --- /dev/null +++ b/f77_wisdom.f @@ -0,0 +1,45 @@ + subroutine write_char(c, iunit) + character c + integer iunit + write(iunit,1000) c + 1000 format(a,$) + end + + subroutine export_wisdom_to_file(iunit) + integer iunit + external write_char +c call dfftw_export_wisdom(write_char, iunit) + call sfftw_export_wisdom(write_char, iunit) + end + + subroutine read_char(ic, iunit) + integer ic + integer iunit + character*256 buf + save buf + integer ibuf + data ibuf/257/ + save ibuf + if (ibuf .lt. 257) then + ic = ichar(buf(ibuf:ibuf)) + ibuf = ibuf + 1 + return + endif + read(iunit,1000,end=10) buf + 1000 format(a256) + ic = ichar(buf(1:1)) + ibuf = 2 + return + 10 ic = -1 + ibuf = 257 + rewind iunit + return + end + + subroutine import_wisdom_from_file(isuccess, iunit) + integer isuccess + integer iunit + external read_char +c call dfftw_import_wisdom(isuccess, read_char, iunit) + call sfftw_import_wisdom(isuccess, read_char, iunit) + end diff --git a/runqqq.F90 b/runqqq.F90 index 07ac5dd8c..943a3921a 100644 --- a/runqqq.F90 +++ b/runqqq.F90 @@ -10,7 +10,7 @@ subroutine runqqq(fname,cmnd,iret) #ifdef CVF iret=runqq(fname,cmnd) #else - iret=system('KVASD_g95 -q > dev_null') + iret=system('kvasd -q > dev_null') #endif return