I just thought I would update my original question with the solution my colleague and I discovered in case anyone else has the same problem.
# Name: FinalKrigingCode.py
# Description: Uses an existing Geostatistical Layer (xml file) as a template
# to interpolate another variable, then exports the new Geostatistical
# Layer as a raster layer. Uses the ListFields function to iterate through
# all fields within the shapefile.
# Requirements: Geostatistical Analyst Extension
# Import system modules
import arcpy, os
# Set environment settings
arcpy.env.workspace = "C:\\Users\\courtney\\Desktop\\MundleStableIsotope"
arcpy.env.extent = "C:\\Users\\courtney\\Desktop\\MundleStableIsotope\\DetroitRiverBoundary.shp"
arcpy.env.mask = "C:\\Users\\courtney\\Desktop\\MundleStableIsotope\\DetroitRiverBoundary.shp"
Sediment1999 = "C:\\Users\\courtney\\Desktop\\MundleStableIsotope\\Sediment1999.shp"
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
#Execute FieldList
fieldList = arcpy.ListFields(Sediment1999, field_type="Double")
fieldNames = [f.name for f in fieldList ]
for name in fieldNames:
# Set local variables - CreateGeostatisticalLayer
inLayer = "C:\\Users\\courtney\\Desktop\\MundleStableIsotope\\Kriging45.xml"
inData = "C:\\Users\\courtney\\Desktop\\MundleStableIsotope\\Sediment1999.shp F1="+name
outLayer = name
# Execute CreateGeostatisticalLayer
arcpy.GACreateGeostatisticalLayer_ga(inLayer, inData, outLayer)
# Set local variables - GALayerToGrid
inLayer = name
outGrid = (os.path.join("C:\\Users\\courtney\\Desktop\\MundleStableIsotope\\Anisotropic\\Sediment\\1999", name+".tif"))
cellSize = 30
cellptsHor = 1
cellptsVer = 1
# Execute GALayerToGrid
arcpy.GALayerToGrid_ga(inLayer, outGrid, cellSize, cellptsHor, cellptsVer)