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 funciton 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.
Programs to compute Y^B functions.
Programs to compute X^C functions.
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.
Programs to compute Y^H functions.
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.