Only allow READ access to Constant coefficients in par_loop

Although in the underlying data model Constants are actually not, raise
an error if trying to write to one in a par_loop, since that's probably
not what was intended.
Fix assignment in direct parloop test

Should be = not +=.
Add some tests of par_loops accessing Constant coefficients