make_rtheta

Create a simple rtheta import file that can be used with Sirocco.

Synopsis:

Create an import file of an rtheta (polar) grid) with either angles in theta or angles set so that the size of cells at a single radius is the same

Command line usage:

usage:

make_rtheta.py file.pf

where file.pf is a parameter file that gives the information needed to create the grid.

If this file does not exist then one will be asked for the necessary information, and the parameter file will be created (for future use).

Description:

The routine produces a spherically expanding wind on an r-theta grid, that can be imported into Sirocco.

The model is intended for testing various aspects of Sirocco, rather than being a model that is physically realistic.

Specifically, although the model has a velocity that increases linearly from an inner to an outer radius.

The file that is created has a uniform density

The rtheta grid can be one that has equal angle bins, or one where the angular bins are adjusted so that the volume of the cells at a particular radius are essentiall constant.

Warning: The routine does not read or create a Sirocco .pf file

To make use of the file produced here, one needs to create an Sirocco pf file that uses the import option and has a coordinate sytem type of polar

Primary routines:

doit

Notes

This routine was developed for testing of the convergence properties of rad/hydro calculations

The routine makes use of a very simplified python implementation of rdpar.

It would be straightfoward to adapt this routine to produce a real stellar wind model with parameters that are consistent with what is used in Sirocco.

History:

230601 ksl

Coding begun

Functions

doit(pf)	Run the program
gen_model(x)	Generate a model given a set of parameters
get_inputs(pf)	Get the inputs