Programs for Stokes resistance functions
A zip file containing all programs and functions can be downloaded from the following link:
Programs to compute X^A functions.
- File "rxa300.dat.txt". 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.txt". 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.txt". The output you should get from compiling and running above program.
- File "axitr2.f.txt". 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.txt". Third program for coefficients. This was run in 16-byte precision as a crude way of checking rounding errors.
- File "rxatst.f.txt". A program to test the coefficients against the asymptotic behaviours predicted by the lubrication analysis. It uses the data rxa300.dat
- File "rxatst.res.txt". The results from running rxatst.f.
- File "rxa.f.txt". Demonstration program using function subroutine for X^A.
- File "rxa.res.txt". The results from running rxa.f.
- File "rax.f.txt". A program that reads rxa300.dat and prints a table of values for RAX and W1.
- File "rax.res.txt". The results from running rax.f.
- File "rax-200.res.txt". The results from calculating RAX using 200 terms.
- File "w2.f.txt". 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.txt". The results from running w2.f.
Programs to compute Y^A functions.
- File "rya300.dat.txt". An ASCII file of 45451 numbers to be read by the other programs.
- File "trntr1.f.txt". The first of two programs to generate rya300.dat.
- File "trntr1.res.txt". The output from trntr1.f.
- File "trntr2.f.txt". The second of these programs.
- File "trntr3.f.txt". The third of these programs. WARNING - This program uses quadruple precision (REAL*16) which is only supported by a few compilers.
- File "ryatst.f.txt". 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.txt". The results from running ryatst.f.
- File "rya.f.txt". A program that reads rya300.dat and prints a table of values for RYA.
- File "rya.res.txt". The results from running rya.f.
- File "ray.f.txt". A program that reads rya300.dat and prints a table of values for RAY.
- File "ray.res.txt". The results from running ray.f.
Programs to compute Y^B functions.
- File "ryb300.dat.txt". An ASCII file of 45451 numbers to be read by the other programs.
- File "trntr1.f.txt". The first of two programs to generate ryb300.dat.
- File "trntr1.res.txt". The output from trntr1.f.
- File "trntr2.f.txt". The second of these programs.
- File "trnro1.f.txt". An alternative program to generate ryb300.dat.
- File "trnro1.res.txt". The results from running trnro1.f.
- File "ryb.f.txt". A program that reads ryb300.dat and prints a table of values for RYB.
- File "ryb.res.txt". The results from running ryb.f.
- File "rby.f.txt". A program that reads ryb300.dat and prints a table of values for RBY.
- File "rby.res.txt". The results from running rby.f.
- File "rybtst.f.txt". A program to test the coefficients against the asymptotic behaviour predicted by lubrication analysis.
- File "rybtst.res.txt". The results from running rybtst.f.
Programs to compute X^C functions.
- File "rxc300.dat.txt". 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.txt". The first of two programs to generate rxc300.dat.
- File "axiro1.res.txt". The output from axiro1.f.
- File "axiro2.f.txt". The second of these programs.
- File "rxc.f.txt". A program that reads rxc300.dat and prints a table of values for RXC.
- File "rxc.res.txt". The results from running rxc.f .
- File "rcx.f.txt". A program that reads rxc300.dat and prints a table of values for RCX.
- File "rcx.res.txt". The results from running rcx.f .
- File "rxctst.f.txt". A program to test the coefficients against the asymptotic behaviour predicted by lubrication analysis. The program uses rxc300.dat.
- File "rxctst.res.txt". The results from running rxctst.f.
Programs to compute Y^C functions.
- File "ryc300.dat.txt". 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.txt". The first of two programs to generate ryc300.dat.
- File "trnro2.f.txt". The second of these programs.
- File "trnro3.f.txt". The third of these programs. WARNING - this program uses quadruple precision (REAL*16) which is only supported by a few compilers.
- File "trnro1.res.txt". The output from trnro1.f.
- File "ryc.f.txt". A program that reads ryc300.dat and prints a table of values for RYC.
- File "ryc.res.txt". results from running ryc.f.
- File "rcy.f.txt". A program that reads ryc300.dat and prints a table of values for RCY.
- File "rcy.res.txt". The results from running rcy.f.
- File "ryctst.f.txt". A program to test the coefficients against the behaviour predicted by lubrication analysis. It uses the data ryc300.dat.
- File "ryctst.res.txt". The results from running ryctst.f.
Programs to compute X^G functions.
- File "rxg300.dat.txt". 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.txt". A program that reads rxg300.dat and prints a table of values for RGX.
- File "rgx.res.txt". The results from running rgx.f.
- File "rxg.f.txt". A program that reads rxg300.dat and prints a table of values for RXG.
- File "rxg.res.txt". The results from running rxg.f.
- File "rxgtst.f.txt". A program to test the coefficients against the asymptotic behaviour predicted by the analysis. It uses the data rxg300.dat.
- File "rxgtst.res.txt". The results from running rxgtst.f.
Programs to compute Y^G functions.
- File "ryg300.dat.txt". 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.txt". The first of two programs to generate ryg300.dat.
- File "trntr1.res.txt". The output from trntr1.f.
- File "trntr2.f.txt". The second of these programs.
- File "trnst1.f.txt". An alternative program to generate the coefficients.
- File "trnst1.res.txt". The output from trnst1.f.
- File "rygtst.f.txt". A program to test the coefficients against the asymptotic behaviour predicted by the lubrication analysis.
- File "rygtst.res.txt". The results from running rygtst.f using 300 terms.
- File "ryg.f.txt". A program that reads ryg300.dat and evaluates RYG.
- File "ryg.res.txt". The results from running ryg.f.
- File "rgy.f.txt". A program that reads ryg300.dat and prints a table of values for RGY.
- File "rgy.res.txt". The results from running rgy.f.
Programs to compute Y^H functions.
- File "ryh300.dat.txt". 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.txt". The first of two programs to generate ryh300.dat.
- File "trnro1.res.txt". The output from trnro1.f.
- File "trnro2.f.txt". The second of these programs.
- File "trnst1.f.txt". Alternative program to generate coefficients .
- File "trnst1.res.txt". Output from trnst1.
- File "ryh.f.txt". A program that reads ryh300.dat and prints a table of values for RYH.
- File "ryh.res.txt". The results from running ryh.f.
- File "rhy.f.txt". A program that reads rhy300.dat and prints a table of values for RHY.
- File "rhy.res.txt". The results from running rhy.f.
- File "rhytst.f". (Available upon request). A program to test the coefficients against the asymptotic behaviour predicted by lubrication analysis.
- File "rhytst.res". (Available upon request). The results for running rhytst.f.
Programs to compute X^M functions.
- File "rxm300.dat.txt". 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.txt". The first of two programs that generate rmx300.dat.
- File "axist1.res.txt". The results from running axist1.f.
- File "axist2.f.txt". The second of these programs. Its results are in a form suitable for the subroutines RXG, RGX, RXM and RMX.
- File "rxmtst.f.txt". 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.txt". The results of running rxmtst.f.
- File "rxm.f.txt". A program that reads rxm300.dat and prints out a table of values for RXM.
- File "rxm.res.txt". The results from running rxm.f.
- File "rmx.f.txt". A program that reads rxm300.dat and prints out a table of values for RMX.
- File "rmx.res.txt". The results from calculating RMX using 300 terms.
Programs to compute Y^M functions.
- File "rym300.dat.txt". 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.txt". The first of two programs used to generate rym300.dat.
- File "trnst1.res.txt". The output from trnst1.f.
- File "trnst2.f.txt". 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.txt". A program to test the coefficients against the asymptotic behaviour predicted by the lubrication analysis.
- File "rymtst.res.txt". The results from running rymtst.f.
- File "rym.f.txt". A program that reads rym300.dat and prints a table of values for RYM.
- File "rym.res.txt". The results from running rym.f.
- File "rmy.f.txt". A program that reads rym300.dat and prints a table of values for RMY.
- File "rmy.res.txt". The results from running rmy.f.
Programs to compute Z^M functions.
- File "rzm300.dat.txt". 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.txt". The first of two programs used to generate rzm300.dat.
- File "prpst2.f.txt". 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.txt". The output from prpst1.f.
- File "rzmtst.f.txt". A program to test the coefficients against the asymptotic behaviour predicted by the lubrication analysis.
- File "rzmtst.res.txt". The results from running rzmtst.f.
- File "rzm.f.txt". A program that reads rzm300.dat and prints a table of values for RZM.
- File "rzm.res.txt". The results from running rzm.f.
- File "rmz.f.txt". A program that reads rzm300.dat and prints a table of values for RMZ.
- File "rmz.res.txt". The results from running rmz.f.
Programs to compute X^P functions.
- File "rxp300.dat.txt". 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.txt". The first of two programs used to generate rxp300.dat.
- File "axitr3.res.txt". The output from axitr3.f.
- File "axitr4.f.txt". The second of these programs. Its results are in a form suitable for the subroutines RXP,RPX.
- File "rpx.f.txt". A program that reads rxp300.dat and prints a table of values for RPX.
- File "rpx.res.txt". The results from running rpx.f.
- File "rxp.f.txt". A program that reads rxp300.dat and prints a table of values for RXP.
- File "rxp.res.txt". The results from running rxp.f.
- File "rxptst.f.txt". A program to test the coefficients against the asymptotic behaviour predicted. It uses the data rxp300.dat.
- File "rxptst.res.txt". The results from running rxptst.f.
Programs to compute X^Q functions.
- File "rxq300.dat.txt". 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.txt". The first of two programs used to generate rxq300.dat.
- File "axist3.res.txt". The output from axist3.f.
- File "axist4.f.txt". The second of these programs. Its results are in a form suitable for the subroutines RXQ,RQX.
- File "rqx.f.txt". A program that reads rxq300.dat and prints a table of values for RQX.
- File "rqx.res.txt". The results from running rqx.f.
- File "rxq.f.txt". A program that reads rxq300.dat and prints a table of values for RXQ.
- File "rxq.res.txt". The results from running rxq.f.
- File "rxqtst.f.txt". A program to test the coefficients against the asymptotic behaviour predicted. It uses the data rxq300.dat.
- File "rxqtst.res.txt". The results from running rxqtst.f.