This code supports the manuscript "Gaussian process forecasts Pseudogymnoascus destructans will cover coterminous United States by 2030." The code is used to fit Gaussian process models to publicly accessible monitoring data on the spread of white-nose syndrome in North America. These models are used to make predictions on a fine spatial grid, giving a forecast (and hindcast) of the spread of white-nose syndrome at any location. Also contained in the code is a retrospective cross validation experiment, producing parameter estimates and model scoring over time. The code also relies on the GRTS grid for model prediction, which is publicly accessible at https://doi.org/10.5066/p9o75ydv. Shapefiles such as administrative boundaries can be used to add to plots are not required for the analysis.