I've set up a .dc sweep with Vgs stepped to attempt to reproduce the datasheet curves Carl posted, using Wimberleytech's J211 model. It seems to be somewhat better than the default NJF one, but it still doesn't look much like the curves in the datasheet. It has a sharper knee, and much lower output resistance in the saturation region.