diff --git a/docs/requirements.txt b/docs/requirements.txt
index 2f34e05a..b614e9ec 100644
--- a/docs/requirements.txt
+++ b/docs/requirements.txt
@@ -3,6 +3,5 @@ sphinx_rtd_theme
 docutils==0.16
 numpy
 xarray
-matplotlib
 netCDF4
 matplotlib
diff --git a/ush/UFS_plot_domains.py b/ush/UFS_plot_domains.py
index a6e629bd..b4162145 100755
--- a/ush/UFS_plot_domains.py
+++ b/ush/UFS_plot_domains.py
@@ -1,5 +1,7 @@
 #!/usr/bin/env python
-
+"""
+Plot domains?
+"""
 import matplotlib.pyplot as plt
 from mpl_toolkits.basemap import Basemap
 from matplotlib.path import Path
@@ -108,15 +110,20 @@
 # Define a function to get the lambert points in the gnomonic space
 
 def get_lambert_points(gnomonic_map, lambert_map,pps):
-    print("Hello from a function")
-    
-    # This function takes the lambert domain we have defined, lambert_map, as well as 
-    # pps (the number of points to interpolate and draw for each side of the lambert "rectangle"), 
-    # and returns an array of two lists: one a list of tuples of the 4*ppf + 4 vertices mapping the approximate shape 
-    # of the lambert domain on the gnomonic map, the other a list of "draw" instructions to be used by
-    # the PathPatch function
+    """
+    This function takes the lambert domain we have defined, lambert_map, as well as 
+    pps (the number of points to interpolate and draw for each side of the lambert "rectangle"), 
+    and returns an array of two lists: one a list of tuples of the 4*ppf + 4 vertices mapping the approximate shape 
+    of the lambert domain on the gnomonic map, the other a list of "draw" instructions to be used by
+    the PathPatch function.
     
-    # pps is recommended 10 or less due to time of calculation
+    pps is recommended 10 or less due to time of calculation.
+
+    gnomonic_map: ???
+    lambert_map: the lambert domain we have defined
+    pps: the number of points to interpolate and draw for each side of the lambert "rectangle"
+    """
+    print("Hello from a function")
     
     # Start array with bottom left point, "MOVETO" instruction
     vertices = [gnomonic_map(*lambert_map(lambert_map.xmin, lambert_map.ymin, inverse= True))]