PyObject * pyCitcom_ic_initVelocity(PyObject *self, PyObject *args)
{
    PyObject *obj;
    struct All_variables* E;

    if (!PyArg_ParseTuple(args, "O:initVelocity", &obj))
        return NULL;

    E = (struct All_variables*)(PyCObject_AsVoidPtr(obj));

    initial_velocity(E);

    Py_INCREF(Py_None);
    return Py_None;
}
/*
Initialize values which are not problem dependent.
NOTE: viscosity may be a function of all previous input fields (temperature,
pressure, velocity, chemistry) and so is always to be done last.
*/
void common_initial_fields(
    struct All_variables *E
)
{
    void initial_qterm();
    void read_viscosity_option();
    void initial_velocity();

    report(E,"Initialize velocity");
    initial_velocity(E);
    report(E,"Initialize pressure field");
    initial_qterm(E);

    
    return;
}