REMORA_InitAnalyticHMix_Upwelling_ML.H

Go to the documentation of this file.

    for ( MFIter mfi((mf_visc2_p), TilingIfNotGPU()); mfi.isValid(); ++mfi )
    {
        Array4<Real> const& visc2_p = (mf_visc2_p).array(mfi);
        Array4<Real> const& visc2_r = (mf_visc2_r).array(mfi);
        Array4<Real> const& diff2   = mf_diff2.array(mfi);
        Box bx = mfi.tilebox();
        bx.grow(IntVect(NGROW,NGROW,0));
        Gpu::streamSynchronize();
 
        int ncomp = mf_diff2.nComp();
 
        ParallelFor(makeSlab(bx,2,0), [=] AMREX_GPU_DEVICE (int i, int j, int )
        {
            visc2_p(i,j,0) = zero; //Real(5.0);
            visc2_r(i,j,0) = zero; // Real(5.0);
 
            for (int n = 0; n < ncomp; n++) {
                diff2(i,j,0,n) = zero;
            }
        });
    }