def METHOD_NAME(output_file, epsg_code,
                           latitude, longitude, elevation, resistivity_data,
                           **kwargs):
    """ Resistivity_data in (elevation, latitude, longitude) grid. """
    with create_dataset(output_file) as dataset:
        dataset.description = 'Resistivity Model'
        z_label = kwargs.get('z_label', 'elevation')
        # dimensions
        dataset.createDimension('latitude', latitude.shape[0])
        dataset.createDimension('longitude', longitude.shape[0])
        dataset.createDimension(z_label, elevation.shape[0])
        # variables
        x = dataset.createVariable('latitude', 'f4', ('latitude',))
        y = dataset.createVariable('longitude', 'f4', ('longitude',))
        z = dataset.createVariable(z_label, 'f4', (z_label,))
        for var in [x, y]:
            var.units = 'degree'
        z.units = 'm'
        resistivity = dataset.createVariable('resistivity', 'f4', (z_label, 'latitude', 'longitude'))
        resistivity.grid_mapping = 'crs'
        resistivity.long_name = 'resistivity'
        resistivity.units = "ohm-m"
        # populate variables
        x[:] = latitude
        y[:] = longitude
        z[:] = elevation
        resistivity[:, :, :] = resistivity_data
        # attach crs info
        crs_var = dataset.createVariable('crs', 'i4', ())
        set_grid_mapping_attrs_geographic(crs_var, epsg_code)