I'm still messing about with the resistor model. I'm sure it isn't a good general purpose solution. It may not be any good at all. I dunno, but think it might work for the circuit I'm working on. It is low voltage, low current DC.
.SUBCKT 1N4728A 1 2
* Node 1 -> A
* Node 2 -> K
* Curve fit to a single, used sample. YMMV
* sampled 0.5mA - 100mA. Very likely bogus outside that range.
* y=a*I**b
***************************************
VS 1 3 0 ;shunt
RZ 3 2 R ={ if( I(VS) > 0, 0.885*abs(I(VS))**-0.972 , 4.48*abs(I(VS))**-0.878) }
.ENDS
Tom.