diff --git a/src/EulerEq/rhs.jl b/src/EulerEq/rhs.jl
index 3d0681a57bee09cea23e99f7e72f38251149df81..0c315f62de297b3e093c61059df2333cf9174506 100644
--- a/src/EulerEq/rhs.jl
+++ b/src/EulerEq/rhs.jl
@@ -200,7 +200,7 @@ function step_fv_muscl!(P::Project, mesh::Mesh1d{FVElement}, state)
         dvâ‚Š[k], dvâ‚‹[k] = vâ‚Š-vâ‚€, vâ‚€-vâ‚‹
       end
     else
-      TODO()
+      # TODO Atm we assume that evolution does not reach bdry
     end
 
     # limit slope and compute left (L) and (R) interface values
@@ -295,7 +295,7 @@ function step_fv_muscl!(P::Project, mesh::Mesh1d{FVElement}, state)
       nflx_E_R[k]   = LLF(nx*fâ‚Š, nx*fâ‚‹, Eâ‚Š, Eâ‚‹, max_v)
     end
   else
-    TODO()
+    # TODO Atm we assume that evolution does not reach bdry
   end
 
   # evaluate rhs
@@ -305,6 +305,12 @@ function step_fv_muscl!(P::Project, mesh::Mesh1d{FVElement}, state)
     rhs_q[k]   = -(nflx_q_R[k]   + nflx_q_L[k]  )/dl
     rhs_E[k]   = -(nflx_E_R[k]   + nflx_E_L[k]  )/dl
   end
+  if !periodic
+    # TODO Atm we assume that evolution does not reach bdry
+    rhs_rho[1] = rhs_rho[end] = 0
+    rhs_q[1]   = rhs_q[end]   = 0
+    rhs_E[1]   = rhs_E[end]   = 0
+  end
 
   return
 end