make_rtheta

Create a simple rtheta model file.

Synopsis:

Create a model which can be imported into Python of an rtheta grid, with either either equal angles in theta or angles set so that the size of cells at a single radius is the same

Command line usage (if any):

usage:

make_rtheta.py file.pf

Description:

The routine reads inputs from a parameter file, or from the screen, and produces a spherically expanding wind on an r-theta grid, that can be imported into Python.

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

Specifically, although the model has a velocity that increases linear from an inner to an outer radius, this is a unifrom density model

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.

The routine does not read a Python .pf file so to make use of it with a Python run, one needs to make sure the inputs to the Python run are consistent with what is generated here.

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 Python.

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