src/solvers/hydro/hydro.cc
/**
 ********************************************************************************************************************************************
 * \brief   Function to initialize the forcing terms for velocity
 *
 *          The forcing terms for the velocity field are initialized here.
 *          Out of the different forcings available in the force class,
 *          the appropriate forcing is chosen according to the parameters set by the user.
 ********************************************************************************************************************************************
 */
\brief   Function to initialize the forcing terms for velocity
The forcing terms for the velocity field are initialized here.
Out of the different forcings available in the force class,
the appropriate forcing is chosen according to the parameters set by the user.
\ brief Function to initialize the forcing terms for velocity The forcing terms for the velocity field are initialized here . Out of the different forcings available in the force class the appropriate forcing is chosen according to the parameters set by the user .
void hydro::initVForcing() {
    switch (inputParams.forceType) {
        case 0:
            if (mpiData.rank == 0) std::cout << "Running hydrodynamics simulation with zero velocity forcing" << std::endl << std::endl;
            vForcing = new zeroForcing(mesh, V);
            break;
        case 1:
            if (mpiData.rank == 0) std::cout << "Running hydrodynamics simulation with random velocity forcing" << std::endl << std::endl;
            vForcing = new randomForcing(mesh, V);
            break;
        case 2:
            if (mpiData.rank == 0) std::cout << "Running hydrodynamics simulation with rotation" << std::endl << std::endl;
            vForcing = new coriolisForce(mesh, V);
            break;
        case 5:
            if (mpiData.rank == 0) std::cout << "Running hydrodynamics simulation with constant pressure gradient along X-direction" << std::endl << std::endl;
            vForcing = new constantPGrad(mesh, V);
            break;
        default:
            if (mpiData.rank == 0) std::cout << "WARNING: Chosen velocity forcing is incompatible with hydrodynamics runs. Defaulting to zero forcing" << std::endl << std::endl;
            vForcing = new zeroForcing(mesh, V);
    }
}
void hydro :: initVForcing ( ) { switch ( inputParams . forceType ) { case 0 : if ( mpiData . rank == 0 ) std :: cout << " " << std :: endl << std :: endl ; vForcing = new zeroForcing ( mesh , V ) ; break ; case 1 : if ( mpiData . rank == 0 ) std :: cout << " " << std :: endl << std :: endl ; vForcing = new randomForcing ( mesh , V ) ; break ; case 2 : if ( mpiData . rank == 0 ) std :: cout << " " << std :: endl << std :: endl ; vForcing = new coriolisForce ( mesh , V ) ; break ; case 5 : if ( mpiData . rank == 0 ) std :: cout << " " << std :: endl << std :: endl ; vForcing = new constantPGrad ( mesh , V ) ; break ; default : if ( mpiData . rank == 0 ) std :: cout << " " << std :: endl << std :: endl ; vForcing = new zeroForcing ( mesh , V ) ; } }
\brief   Function to initialize the forcing terms for velocity
          The forcing terms for the velocity field are initialized here.
\ brief Function to initialize the forcing terms for velocity The forcing terms for the velocity field are initialized here .