This package contains a fortran77 code which calculates a Roche model
with the irradiation. It is described in :
Budaj J., 2011, AJ 141, 59
First you need to unzip and untar the archive:
gunzip roche.tar.gz
tar xvf roche.tar

Go to the newly created directory 'public'.
cd public

The directory contains the following files:
 readme                    this file
 roche                     executable version-may not work on your comp
 roche7p.f                 source code
 roche.gnu                 gnuplot plotting file
 roche.in                  input file
 roche.inexample1          example of input file
 roche.inexample2          example of input file
 roche.inexample3          example of input file
 roche.out                 output file
 roche.outexample1         example of output file
 roche.outexample2         example of output file 
 roche.outexample3         example of output file 
 roche.run                 batch file to run examples

compile the source code:
g77 roche7p.f -o roche
or
gfortran roche7p.f -o roche

check and then run the example batch file with three tests/examples
(exoplanet, secondary, contact system):
./roche.run

You may need to adjust the dimensation parameters ndim1,ndim2
according to the memory available.

Input:
Explanation of the input parameters is at the end of the input files.

Output:  
  Rsub         Rback        Rpole        Rside
radius at the substellar point, antistellar point, rotation pole, side point
Rmean- radius of a sphere with the same volume as the Roche surface
Rcross =sqrt(Rside*Rpole) i.e. radius of the circle with the same 
crossection as the Roche surface during the transit)
ff-Roche lobe filling parameter

roche.out contains the following columns in cgs units
x, y, Roche surface, normalized gravity, cosine of the irradiation angle, 
temperature

Jan Budaj (budaj@ta3.sk)
Tatranska Lomnica
Feb 24, 2011