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)