def METHOD_NAME(self):
    fileName = str(self.datadir.join('example_jan.nc'))
    simulationStartTime = '0001-01-01'
    variableMap = {
        'avgSurfaceTemperature':
            ['time_avg_avgValueWithinOceanRegion_avgSurfaceTemperature',
             'other_string',
             'yet_another_string'],
        'daysSinceStartOfSim':
            ['time_avg_daysSinceStartOfSim',
             'xtime',
             'something_else'],
        'avgLayerTemperature':
            ['time_avg_avgValueWithinOceanLayerRegion_avgLayerTemperature',
             'test1',
             'test2'],
        'Time': [['xtime_start', 'xtime_end'],
                 'time_avg_daysSinceStartOfSim']}
    variableList = ['avgSurfaceTemperature', 'avgLayerTemperature',
                    'refBottomDepth', 'daysSinceStartOfSim']
    config = self.setup_config()
    for calendar in ['gregorian', 'noleap']:
        # preprocess_mpas will use variableMap to map the variable names
        # from their values in the file to the desired values in
        # variableList
        ds = open_multifile_dataset(
            fileNames=fileName,
            calendar=calendar,
            config=config,
            simulationStartTime=simulationStartTime,
            timeVariableName='Time',
            variableList=variableList,
            variableMap=variableMap)
        # make sure the remapping happened as expected
        dsVarList = list(ds.data_vars.keys()) + list(ds.coords.keys())
        assert(numpy.all([var in dsVarList for var in variableList]))