<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div><div class="h5"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><br></div><div>  ! end</div><div>  call VecScatterEnd(  a_ts%from_petsc(0),a_XX,xxvec<wbr>(0),INSERT_VALUES,SCATTER_<wbr>FORWARD,ierr)</div><div>  a_n1 = n1</div></div></blockquote><div><br></div></div></div><div>Why is the line above in here?</div></div></div></div></blockquote><div><br></div><div>a_n1 is the return value of the subroutine.  I just did this to keep the types clear, real kind=8 vs PetscScalar. </div></div></div></div>