[Users] boundary conditions with Llama

[Miguel Zilhão]

hi all,

i've been implementing a scalar wave evolution code to work with the Llama multipatch 
infrastructure. i think i have it working except for the boundary conditions. as a first test, i 
tried to call the NewRad infrastructure, in the same way that i usually do for my non-multipatch codes:

   ierr = NewRad_Apply(cctkGH, phi, rhs_phi, phi_0, one, n_phi)
   ierr = NewRad_Apply(cctkGH, Kphi, rhs_Kphi, Kphi_0, one, n_Kphi)

but this does not seem to work as i get instabilities very early on... am i missing something, or is 
NewRad not designed to work with Llama? if so, what is the "standard" way of imposing BCs for 
multipatch evolutions?

thanks,
Miguel