Hi Paul,<br><br>Thanks for pointing me to the taylor example.  I'm now successfully running simulations on an annulus, but the part of my code in userchk that computed volume averages and inner products is no longer working properly.<br>

<br>e.g., (glsc3(vx,bm1,vx,n) returns NaN.  Is this to be expected?<br><br>Best,<br><br>David<br>