| 87 | } |
| 88 | |
| 89 | void |
| 90 | CNS::initData () |
| 91 | { |
| 92 | BL_PROFILE("CNS::initData()"); |
| 93 | |
| 94 | const Real* dx = geom.CellSize(); |
| 95 | const Real* prob_lo = geom.ProbLo(); |
| 96 | MultiFab& S_new = get_new_data(State_Type); |
| 97 | Real cur_time = state[State_Type].curTime(); |
| 98 | |
| 99 | #ifdef AMREX_USE_OMP |
| 100 | #pragma omp parallel |
| 101 | #endif |
| 102 | for (MFIter mfi(S_new); mfi.isValid(); ++mfi) |
| 103 | { |
| 104 | const Box& box = mfi.validbox(); |
| 105 | cns_initdata(&level, &cur_time, |
| 106 | BL_TO_FORTRAN_BOX(box), |
| 107 | BL_TO_FORTRAN_ANYD(S_new[mfi]), |
| 108 | dx, prob_lo); |
| 109 | } |
| 110 | |
| 111 | MultiFab& C_new = get_new_data(Cost_Type); |
| 112 | C_new.setVal(1.0); |
| 113 | } |
| 114 | |
| 115 | void |
| 116 | CNS::computeInitialDt (int finest_level, |