This is making a binary presence absence raster for anywhere Gelisols were mapped in STATSGO (even if just 1% of the comp_r). # Here we are building to generate a difference raster between the NSP ...