Simulation functions set up (though probably need tweaking to get quadratic ice effect to better match up with observed data range), separable and space-time proc conv models set up.  About to switch 'CPIF' model from a funny Poisson distribution to a multinomial sort of framework.