lib/field/vfield.cc
/**
 ********************************************************************************************************************************************
 * \brief   Operator is used to calculate time step #dt_out using CFL Condition
 *           
 *          When the parameters specify that time-step must be adaptively
 *          calculated using the Courant Number given in the parameters file,
 *          this function will provide the \f$ dt \f$ using the maximum values
 *          of the components of the vfield.
 *
 * \param   dt_out is a reference to the real value into which the calculated value of time-step is written
 *********************************************************************************************************************************************
 */
\brief   Operator is used to calculate time step #dt_out using CFL Condition

When the parameters specify that time-step must be adaptively
calculated using the Courant Number given in the parameters file,
this function will provide the \f$ dt \f$ using the maximum values
of the components of the vfield.

\param   dt_out is a reference to the real value into which the calculated value of time-step is written
\ brief Operator is used to calculate time step #dt_out using CFL Condition When the parameters specify that time - step must be adaptively calculated using the Courant Number given in the parameters file this function will provide the \ f$ dt \ f$ using the maximum values of the components of the vfield . \ param dt_out is a reference to the real value into which the calculated value of time - step is written
void vfield::computeTStp(real &dt_out) {
    real Umax, Vmax, Wmax;
    real delx, dely, delz; 
    delx = gridData.dXi;
    Umax = Vx.fieldMax();
#ifdef PLANAR
    dely = 0.0;
    Vmax = 1.0;
#else
    dely = gridData.dEt;
    Vmax = Vy.fieldMax();
#endif
    delz = gridData.dZt;
    Wmax = Vz.fieldMax();
    dt_out = real(gridData.inputParams.courantNumber*(delx/Umax + dely/Vmax + delz/Wmax));
}
void vfield :: computeTStp ( real & dt_out ) { real Umax , Vmax , Wmax ; real delx , dely , delz ; delx = gridData . dXi ; Umax = Vx . fieldMax ( ) ; #ifdef PLANAR dely = 0.0 ; Vmax = 1.0 ; #else dely = gridData . dEt ; Vmax = Vy . fieldMax ( ) ; #endif delz = gridData . dZt ; Wmax = Vz . fieldMax ( ) ; dt_out = real ( gridData . inputParams . courantNumber * ( delx / Umax + dely / Vmax + delz / Wmax ) ) ; }
\brief   Operator is used to calculate time step #dt_out using CFL Condition
\ brief Operator is used to calculate time step #dt_out using CFL Condition