         //Compute parallel dissipative energy for N/////////////////////////////
         if( p.nu_parallel != 0)
         {
-                if (p.pollim==1) dzN_.set_boundaries( p.bc, 0, 0);                
+                if (p.pollim==1) dzN_.set_boundaries( p.bc, 0, 0);             
+                #ifndef MPI_VERSION
                 if (p.pardiss==0)
                 {
                     dzN_.forward(y[i], omega); 

                    dzDIR_.backwardT(omega,lambda);
                    dg::blas1::axpby( 0.5*p.nu_parallel, lambda, 1., yp[i+2]); 
                }
                if (p.pardiss==1)
                {
                    if (p.pollim==1) dzN_.set_boundaries( p.bc, 0, 0);

                    dzDIR_.backwardT(omega,lambda);
                    dg::blas1::axpby( 0.5*p.nu_parallel, lambda,  1., lambda,lambda); 
                }
                #endif

                if (p.pardiss==1)
                {
/*                    dzDIR_.dzz(y[i+2],omega);                                          //dz^2 U 
