diff --git a/doc/source/basic_usage.rst b/doc/source/basic_usage.rst
index 24c3d79..350cab4 100644
--- a/doc/source/basic_usage.rst
+++ b/doc/source/basic_usage.rst
@@ -1027,12 +1027,13 @@ The surface levels to interpolate also need to be specified.
     from __future__ import print_function
 
     from netCDF4 import Dataset
-    from wrf import getvar, interpline, CoordPair  
+    from wrf import getvar, vinterp 
     
     ncfile = Dataset("wrfout_d01_2016-10-07_00_00_00")  
     
-    # Interpolate tk to theta-e levels                
-    interp_levels = [200, 300, 500, 1000]
+    tk = getvar(ncfile, "tk")
+    # Interpolate tk to theta-e levels                 
+    interp_levels = [200, 300, 500, 1000]
     
     interp_field = vinterp(ncfile, 
                    field=tk,