Programs for Stokes resistance functions.

A zip file containing all programs and functions can be downloaded from the following link: File "all_resist.ZIP"

Programs to compute X^A functions

  • File "rxa300.dat". The coefficients for the X^A functions. This file corresponds to 300 terms of the series for X^A, or A^X. The generating programs are given below. The file is read by the demonstration programs, which sum the coefficients and tabulate X^A and A^X.
  • File "axitr1.f". Demonstration program to compute coefficients for X^A functions.This version has most comments and the output allows direct comparison with paper.
  • File "axitr1.res". The output you should get from compiling and running above program.
  • File "axitr2.f". Main program to compute coefficients for X^A functions. This version writes the coefficients in the format that is expected by the evaluation subroutines.
  • File "axitr3.f". Third program for coefficients. This was run in 16-byte precision as a crude way of checking rounding errors.
  • File "rxatst.f". A program to test the coefficients against the asymptotic behaviours predicted by the lubrication analysis. It uses the data rxa300.dat
  • File "rxatst.res". The results from running rxatst.f.
  • File "rxa.f". Demonstration program using function subroutine for X^A.
  • File "rxa.res". The results from running rxa.f.
  • File "rax.f". A program that reads rxa300.dat and prints a table of values for RAX and W1.
  • File "rax.res". The results from running rax.f.
  • File "rax-200.res". The results from calculating RAX using 200 terms.
  • File "w2.f". A program to evaluate the function W2 (lambda) defined in Jeffrey & Onishi (1984). The version here was not used for the paper.
  • File "w2.res". The results from running w2.f.
  • Programs to compute Y^A functions.

  • File "rya300.dat". An ASCII file of 45451 numbers to be read by the other programs.
  • File "trntr1.f". The first of two programs to generate rya300.dat.
  • File "trntr1.res". The output from trntr1.f.
  • File "trntr2.f". The second of these programs.
  • File "trntr3.f". The third of these programs. WARNING - This program uses quadruple precision (REAL*16) which is only supported by a few compilers.
  • File "ryatst.f". A program to test the coefficients against the asymptotic behaviour predicted by the lubrication analysis. It uses the data in rya300.dat.
  • File "ryatst.res". The results from running ryatst.f.
  • File "rya.f". A program that reads rya300.dat and prints a table of values for RYA.
  • File "rya.res". The results from running rya.f.
  • File "ray.f". A program that reads rya300.dat and prints a table of values for RAY.
  • File "ray.res". The results from running ray.f.
  • Programs to compute Y^B functions.

  • File "ryb300.dat". An ASCII file of 45451 numbers to be read by the other programs.
  • File "trntr1.f". The first of two programs to generate ryb300.dat.
  • File "trntr1.res". The output from trntr1.f.
  • File "trntr2.f". The second of these programs.
  • File "trnro1.f". An alternative program to generate ryb300.dat.
  • File "trnro1.res". The results from running trnro1.f.
  • File "ryb.f". A program that reads ryb300.dat and prints a table of values for RYB.
  • File "ryb.res". The results from running ryb.f.
  • File "rby.f". A program that reads ryb300.dat and prints a table of values for RBY.
  • File "rby.res". The results from running rby.f.
  • File "rybtst.f". A program to test the coefficients against the asymptotic behaviour predicted by lubrication analysis.
  • File "rybtst.res". The results from running rybtst.f.
  • Programs to compute X^C functions.

  • File "rxc300.dat". A file of 45451 numbers to be read by the other programs. They are the coefficients for 300 terms of the series.
  • File "axiro1.f". The first of two programs to generate rxc300.dat.
  • File "axiro1.res". The output from axiro1.f.
  • File "axiro2.f". The second of these programs.
  • File "rxc.f". A program that reads rxc300.dat and prints a table of values for RXC.
  • File "rxc.res". The results from running rxc.f .
  • File "rcx.f". A program that reads rxc300.dat and prints a table of values for RCX.
  • File "rcx.res". The results from running rcx.f .
  • File "rxctst.f". A program to test the coefficients against the asymptotic behaviour predicted by lubrication analysis. The program uses rxc300.dat.
  • File "rxctst.res". The results from running rxctst.f.
  • Programs to compute Y^C functions.

  • File "ryc300.dat". An ASCII file of 45451 numbers to be read by the other programs. They are the coefficients for 300 terms of the series.
  • File "trnro1.f". The first of two programs to generate ryc300.dat.
  • File "trnro2.f". The second of these programs.
  • File "trnro3.f". The third of these programs. WARNING - this program uses quadruple precision (REAL*16) which is only supported by a few compilers.
  • File "trnro1.res". The output from trnro1.f.
  • File "ryc.f". A program that reads ryc300.dat and prints a table of values for RYC.
  • File "ryc.res". results from running ryc.f.
  • File "rcy.f". A program that reads ryc300.dat and prints a table of values for RCY.
  • File "rcy.res". The results from running rcy.f.
  • File "ryctst.f". A program to test the coefficients against the behaviour predicted by lubrication analysis. It uses the data ryc300.dat.
  • File "ryctst.res". The results from running ryctst.f.
  • Programs to compute X^G functions.

  • File "rxg300.dat". A file of 45451 numbers to be read by the other programs. They are the coefficients for 300 terms of the series.
  • File "rgx.f". A program that reads rxg300.dat and prints a table of values for RGX.
  • File "rgx.res". The results from running rgx.f.
  • File "rxg.f". A program that reads rxg300.dat and prints a table of values for RXG.
  • File "rxg.res". The results from running rxg.f.
  • File "rxgtst.f". A program to test the coefficients against the asymptotic behaviour predicted by the analysis. It uses the data rxg300.dat.
  • File "rxgtst.res". The results from running rxgtst.f.
  • Programs to compute Y^G functions.

  • File "ryg300.dat". An ASCII file of 45451 numbers to be read by the other programs. They are the coefficients for 300 terms of the series.
  • File "trntr1.f". The first of two programs to generate ryg300.dat.
  • File "trntr1.res". The output from trntr1.f.
  • File "trntr2.f". The second of these programs.
  • File "trnst1.f". An alternative program to generate the coefficients.
  • File "trnst1.res". The output from trnst1.f.
  • File "rygtst.f". A program to test the coefficients against the asymptotic behaviour predicted by the lubrication analysis.
  • File "rygtst.res". The results from running rygtst.f using 300 terms.
  • File "ryg.f". A program that reads ryg300.dat and evaluates RYG.
  • File "ryg.res". The results from running ryg.f.
  • File "rgy.f". A program that reads ryg300.dat and prints a table of values for RGY.
  • File "rgy.res". The results from running rgy.f.
  • Programs to compute Y^H functions.

  • File "ryh300.dat". An ASCII file of 45451 numbers to be read by the other programs. They are the coefficients for 300 terms of the series.
  • File "trnro1.f". The first of two programs to generate ryh300.dat.
  • File "trnro1.res". The output from trnro1.f.
  • File "trnro2.f". The second of these programs.
  • File "trnst1.f". Alternative program to generate coefficients .
  • File "trnst1.res". Output from trnst1.
  • File "ryh.f". A program that reads ryh300.dat and prints a table of values for RYH.
  • File "ryh.res". The results from running ryh.f.
  • File "rhy.f". A program that reads rhy300.dat and prints a table of values for RHY.
  • File "rhy.res". The results from running rhy.f.
  • File "rhytst.f". A program to test the coefficients against the asymptotic behaviour predicted by lubrication analysis.
  • File "rhytst.res". The results for running rhytst.f.
  • Programs to compute X^M functions.

  • File "rxm300.dat". An ASCII file of 45451 numbers to be read by the other programs. They are the coefficients for 200 terms of the series.
  • File "axist1.f". The first of two programs that generate rmx300.dat.
  • File "axist1.res". The results from running axist1.f.
  • File "axist2.f". The second of these programs. Its results are in a form suitable for the subroutines RXG, RGX, RXM and RMX.
  • File "rxmtst.f". A program to test the coefficients against the asymptotic behaviour predicted by the lubrication analysis. It uses the data in rxm300.dat.
  • File "rxmtst.res". The results of running rxmtst.f.
  • File "rxm.f". A program that reads rxm300.dat and prints out a table of values for RXM.
  • File "rxm.res". The results from running rxm.f.
  • File "rmx.f". A program that reads rxm300.dat and prints out a table of values for RMX.
  • File "rmx.res". The results from calculating RMX using 300 terms.
  • Programs to compute Y^M functions.

  • File "rym300.dat". An ASCII file of 45451 numbers to be read by the other programs. They are the coefficients for 300 terms of the series.
  • File "trnst1.f". The first of two programs used to generate rym300.dat.
  • File "trnst1.res". The output from trnst1.f.
  • File "trnst2.f". The second of these programs, which opens output files for the results. The results of trnst2.f are in a form suitable for the subroutines RYM, RMY, RYG, RGY, RYH and RHY.
  • File "rymtst.f". A program to test the coefficients against the asymptotic behaviour predicted by the lubrication analysis.
  • File "rymtst.res". The results from running rymtst.f.
  • File "rym.f". A program that reads rym300.dat and prints a table of values for RYM.
  • File "rym.res". The results from running rym.f.
  • File "rmy.f". A program that reads rym300.dat and prints a table of values for RMY.
  • File "rmy.res". The results from running rmy.f.
  • Programs to compute Z^M functions.

  • File "rzm300.dat". An ASCII file of 45451 numbers to be read by the other programs. They are the coefficients for 300 terms of the series.
  • File "prpst1.f". The first of two programs used to generate rzm300.dat.
  • File "prpst2.f". The second of these programs, which opens output files for the results. The results of prpst2.f are in a form suitable for the subroutines RZM and RMZ.
  • File "prpst1.res". The output from prpst1.f.
  • File "rzmtst.f". A program to test the coefficients against the asymptotic behaviour predicted by the lubrication analysis.
  • File "rzmtst.res". The results from running rzmtst.f.
  • File "rzm.f". A program that reads rzm300.dat and prints a table of values for RZM.
  • File "rzm.res". The results from running rzm.f.
  • File "rmz.f". A program that reads rzm300.dat and prints a table of values for RMZ.
  • File "rmz.res". The results from running rmz.f.
  • Programs to compute X^P functions.

  • File "rxp300.dat". A file of 45451 numbers to be read by the other programs. They are the coefficients for 300 terms of the series.
  • File "axitr3.f". The first of two programs used to generate rxp300.dat.
  • File "axitr3.res". The output from axitr3.f.
  • File "axitr4.f". The second of these programs. Its results are in a form suitable for the subroutines RXP,RPX.
  • File "rpx.f". A program that reads rxp300.dat and prints a table of values for RPX.
  • File "rpx.res". The results from running rpx.f.
  • File "rxp.f". A program that reads rxp300.dat and prints a table of values for RXP.
  • File "rxp.res". The results from running rxp.f.
  • File "rxptst.f". A program to test the coefficients against the asymptotic behaviour predicted. It uses the data rxp300.dat.
  • File "rxptst.res". The results from running rxptst.f.
  • Programs to compute X^Q functions.

  • File "rxq300.dat". A file of 45451 numbers to be read by the other programs. They are the coefficients for 300 terms of the series.
  • File "axist3.f". The first of two programs used to generate rxq300.dat.
  • File "axist3.res". The output from axist3.f.
  • File "axist4.f". The second of these programs. Its results are in a form suitable for the subroutines RXQ,RQX.
  • File "rqx.f". A program that reads rxq300.dat and prints a table of values for RQX.
  • File "rqx.res". The results from running rqx.f.
  • File "rxq.f". A program that reads rxq300.dat and prints a table of values for RXQ.
  • File "rxq.res". The results from running rxq.f.
  • File "rxqtst.f". A program to test the coefficients against the asymptotic behaviour predicted. It uses the data rxq300.dat.
  • File "rxqtst.res". The results from running rxqtst.f.
  • free counters