diff --git a/utility_scripts/make_mpas_to_lat_lon_mapping.py b/utility_scripts/make_mpas_to_lat_lon_mapping.py
new file mode 100755
index 000000000..00d508356
--- /dev/null
+++ b/utility_scripts/make_mpas_to_lat_lon_mapping.py
@@ -0,0 +1,39 @@
+#!/usr/bin/env python
+
+'''
+Creates a mapping file that can be used with ncremap (NCO) to remap MPAS files
+to a latitude/longitude grid.
+
+Usage: Copy this script into the main MPAS-Analysis directory (up one level).
+Modify the grid name, the path to the MPAS grid file and the output grid
+resolution.
+'''
+
+from mpas_analysis.shared.interpolation import Remapper
+from mpas_analysis.shared.grid import MpasMeshDescriptor
+from mpas_analysis.shared.climatology import get_comparison_descriptor
+from mpas_analysis.configuration import MpasAnalysisConfigParser
+
+
+# replace with the MPAS mesh name
+inGridName = 'oQU240'
+
+# replace with the path to the desired mesh or restart file
+inGridFileName = '/media/xylar/extra_data/analysis/edison/G-QU240-master-intel/run/mpaso.rst.0001-01-06_00000.nc'
+
+config = MpasAnalysisConfigParser()
+config.read('mpas_analysis/config.default')
+# replace 1.0 with the desired resolution of the output mesh
+config.set('climatology', 'comparisonLatResolution', '1.0')
+config.set('climatology', 'comparisonLonResolution', '1.0')
+
+inDescriptor = MpasMeshDescriptor(inGridFileName, inGridName)
+
+outDescriptor = get_comparison_descriptor(config, 'latlon')
+outGridName = outDescriptor.meshName
+
+mappingFileName = 'map_{}_to_{}.nc'.format(inGridName, outGridName)
+
+remapper = Remapper(inDescriptor, outDescriptor, mappingFileName)
+
+remapper.build_mapping_file(method='bilinear')