diff --git a/README.md b/README.md
index 44ebc8307d896b325ac289f9704a4f57e6fa3e57..c1aff676bb6574af644a92fc264549f412fdfabf 100644
--- a/README.md
+++ b/README.md
@@ -23,7 +23,8 @@ Python requires all of the NCM databases:
 
 * MineralPhysicsDatabase.nc
 
-These codes are part of an effort to produce a three dimensional national crustal model for use with seismic hazard studies. (<https://doi.org/10.5066/P9T96Q67>)
+These codes are part of an effort to produce a three dimensional national crustal mode
+for use with seismic hazard studies. (<https://doi.org/10.5066/P9T96Q67>)
 
 ## Python
 
@@ -36,10 +37,14 @@ See the [README](matlab/README.md) file in the [matlab](matlab) directory for ru
 
 ## Help
 
-Any questions regarding the code can be directed to olboyd@usgs.gov.
+Any questions regarding the code can be directed to <olboyd@usgs.gov>.
 
 ## References
 
 *Boyd, O.S., 2020, Calibration of the U.S. Geological Survey
 National Crustal Model: U.S. Geological Survey Open-File
 Report 2020–1052, 23 p., <https://doi.org/10.3133/ofr20201052>*
+
+**Suggested citation:** Boyd, O. S., Clayton, B. S., 2019, GeoPhys—Set of utilities is used
+to extract geophysical information from the National Crustal Model: U.S. Geological Survey
+software release, <https://doi.org/10.5066/P9NQ5LNU>, IP-116829, BAO Approval 02/24/20.
diff --git a/geophys/README.md b/geophys/README.md
index 2ec8a2412213ffacd9de6e3b6aba57daeb173fce..186e5ee5619d6c2b29e5d05d71ae12530fd9f66e 100644
--- a/geophys/README.md
+++ b/geophys/README.md
@@ -1,8 +1,16 @@
 # GeoPhys Python routines ![Shear-wave velocity model in the LA area, looking north. Color scale is 100 to 1000 m/s Vs30 on the top and 100 to 5000 m/s shear-wave velocity on the sides. Vertical exageration is a factor of 10.](../images/LACrossSectionVelocitySmall.png)
 
-This set of utilities is used to extract geophysical profiles from the NCM (Boyd, 2020). Running the code requires the NCM Calibration dataset (<https://doi.org/10.5066/P9GO3CP8>), NCM Geologic Framework databases (<https://doi.org/10.5066/P9SBQENM>), NCM thermal model, NCM_TemperatureGrids.nc (<https://doi.org/10.5066/P935DT1G>), and the NCM Petrologic and Mineral Physics database (<https://doi.org/10.5066/P9HN170G>). These codes and datasets are part of an effort to produce a three dimensional national crustal model for use with seismic hazard studies. (<https://doi.org/10.5066/P9T96Q67>)
-
-Profiles can be returned from roughly -126 degrees to -100 degrees longitude and 22 to 53 degrees latitude. If the requested location is outside of the bounds of the NCM, a message to that effect will be returned.
+This set of utilities is used to extract geophysical profiles from the NCM (Boyd, 2020).
+Running the code requires the NCM Calibration dataset (<https://doi.org/10.5066/P9GO3CP8>),
+NCM Geologic Framework databases (<https://doi.org/10.5066/P9SBQENM>), NCM thermal model,
+NCM_TemperatureGrids.nc (<https://doi.org/10.5066/P935DT1G>), and the NCM Petrologic and
+Mineral Physics database (<https://doi.org/10.5066/P9HN170G>). These codes and datasets
+are part of an effort to produce a three dimensional national crustal model for use with
+seismic hazard studies. (<https://doi.org/10.5066/P9T96Q67>)
+
+Profiles can be returned from roughly -126 degrees to -100 degrees longitude and 22 to 53 degrees
+latitude. If the requested location is outside of the bounds of the NCM, a message to that effect
+will be returned.
 
 ## Prerequisites
 
@@ -103,7 +111,7 @@ Please see or query the code for help.
 poetry run help
 ```
 
-Any questions regarding the code and data can be directed to olboyd@usgs.gov.
+Any questions regarding the code and data can be directed to <olboyd@usgs.gov>.
 
 ## Development
 
diff --git a/matlab/README.md b/matlab/README.md
index 52fcd5fccd25fac51beee31cd022ef677c09e756..ad0bded50abb8058390f0167324c4e7a0f59dd7f 100644
--- a/matlab/README.md
+++ b/matlab/README.md
@@ -1,12 +1,24 @@
 # GeoPhys MATLAB routines ![Shear-wave velocity model in the LA area, looking north. Color scale is 100 to 1000 m/s Vs30 on the top and 100 to 5000 m/s shear-wave velocity on the sides. Vertical exageration is a factor of 10.](images/LACrossSectionVelocitySmall.png)
 
-This program is used to extract geophysical profiles from the NCM (Boyd, 2020). Running the code requires the NCM Calibration dataset (<https://doi.org/10.5066/P9GO3CP8>), NCM Geologic Framework databases (<https://doi.org/10.5066/P9SBQENM>), NCM thermal model, NCM_TemperatureGrids.nc (<https://doi.org/10.5066/P935DT1G>), and the NCM Petrologic and Mineral Physics database (<https://doi.org/10.5066/P9HN170G>). These codes and datasets are part of an effort to produce a three dimensional national crustal model for use with seismic hazard studies. (<https://doi.org/10.5066/P9T96Q67>)
-
-Profiles can be returned from roughly -126 degrees to -100 degrees longitude and 22 to 53 degrees latitude. If the requested location is outside of the bounds of the NCM, a message to that effect will be returned.
+This program is used to extract geophysical profiles from the NCM (Boyd, 2020).
+Running the code requires the NCM Calibration dataset (<https://doi.org/10.5066/P9GO3CP8>),
+NCM Geologic Framework databases (<https://doi.org/10.5066/P9SBQENM>), NCM thermal model,
+NCM_TemperatureGrids.nc (<https://doi.org/10.5066/P935DT1G>), and the NCM Petrologic and
+Mineral Physics database (<https://doi.org/10.5066/P9HN170G>). These codes and datasets
+are part of an effort to produce a three dimensional national crustal model for use with
+seismic hazard studies. (<https://doi.org/10.5066/P9T96Q67>)
+
+Profiles can be returned from roughly -126 degrees to -100 degrees longitude and 22 to 53
+degrees latitude. If the requested location is outside of the bounds of the NCM, a message
+to that effect will be returned.
 
 ## Installation
 
-Within matlab, add to your path the location of the GeoPhys matlab scripts. The LoadNCMData.m script assumes that the databases (mentioned above) are installed in NCMBase/database where NCMBase is ../../(Location of LoadNCMData.m). For example, GeoPhys.m might be located in Path/NCM/GeoPhys/matlab. The database directory would then be in Path/NCM/database. LoadNCMData.m can be edited to allow for a different location of the database.
+Within matlab, add to your path the location of the GeoPhys matlab scripts. The LoadNCMData.m
+script assumes that the databases (mentioned above) are installed in NCMBase/database where
+NCMBase is ../../(Location of LoadNCMData.m). For example, GeoPhys.m might be located in
+Path/NCM/GeoPhys/matlab. The database directory would then be in Path/NCM/database.
+LoadNCMData.m can be edited to allow for a different location of the database.
 
 ## Example
 
@@ -16,9 +28,9 @@ See example.m for an example.
 
 Please see or query the code for help.
 
-**matlab>> help GeoPhys**
+### Matlab Help
 
-Any questions regarding the code and data can be directed to olboyd@usgs.gov.
+Any questions regarding the code and data can be directed to <olboyd@usgs.gov>.
 
 ## References