TALK=T;RUN(1,1)
L(218) TEXT(1D Y; test of u/r = constant TITLE DISPLAY This case differs form case 218 in the boundary conditions for u1 which corrspond to those of solid-body rotation, i.e. u1/r is given the same value at both inner and outer boundaries. Identical values should be found atvintermediate locations also. ENDDIS IURVAL=-1 ! this setting causes the 'val' in the coval to be i ! interpreted as the value of u1/r, not yhat of u1 IURPRN=-1 ! this ensures that it is values of u1/r which are printed IURINI=-1 ! this insures that fiinit(u1) is treated as fiinit(u1/r)