diff --git a/DESCRIPTION b/DESCRIPTION
index ad606306956ae9ba51162b987f8f0a8e9a7a8db3..6119e508b663e0d00b45e2eb801b4f6dc12dea65 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -31,6 +31,8 @@ Imports:
     stats,
     plyr,
     ncdf4,
+    rgdal,
+    sp,
     multiApply (>= 2.0.0)
 Suggests:
     easyVerification,
diff --git a/NAMESPACE b/NAMESPACE
index 627bf51ce445554a9aa98ee491eb9163808fc2dc..34b6b1a527b4220428d912ef9be98ef638251cc3 100644
--- a/NAMESPACE
+++ b/NAMESPACE
@@ -65,6 +65,8 @@ importFrom(grDevices,rgb)
 importFrom(grDevices,svg)
 importFrom(grDevices,tiff)
 importFrom(plyr,take)
+importFrom(rgdal,readOGR)
+importFrom(sp,spTransform)
 importFrom(stats,acf)
 importFrom(stats,confint)
 importFrom(stats,cor)
diff --git a/R/PlotEquiMap.R b/R/PlotEquiMap.R
index 37847dcc40166713a346b3859c38a75f85699fcc..cd9858dbc706ad47a198d52d6a099a7e81b813c1 100644
--- a/R/PlotEquiMap.R
+++ b/R/PlotEquiMap.R
@@ -148,6 +148,12 @@
 #'  the corresponding device.
 #'@param res Resolution of the device (file or window) to plot in. See 
 #'  ?Devices and the creator function of the corresponding device.
+#'@param input_shapefile A character string describing the pathway to your
+#'  shapefile of the region of interest.
+#'@param lwd_shapefile Line width of the shapefile boundary lines. 
+#'  Defaults to 2.
+#'@param col_shapefile Colour of the shapefile boundary lines. 
+#'  Defaults to 'black'.
 #'@param \dots Arguments to be passed to the method. Only accepts the following 
 #'  graphical parameters:\cr
 #'  adj ann ask bg bty cex.sub cin col.axis col.lab col.main col.sub cra crt 
@@ -206,6 +212,8 @@
 #'@import graphics GEOmap geomapdata maps
 #'@importFrom grDevices dev.cur dev.new dev.off gray
 #'@importFrom stats cor
+#'@importFrom rgdal readOGR
+#'@importFrom sp spTransform
 #'@export
 PlotEquiMap <- function(var, lon, lat, varu = NULL, varv = NULL,
                         toptitle = NULL, sizetit = NULL, units = NULL, 
@@ -234,7 +242,8 @@ PlotEquiMap <- function(var, lon, lat, varu = NULL, varv = NULL,
                         margin_scale = rep(1, 4), title_scale = 1, 
                         numbfig = NULL, fileout = NULL, 
                         width = 8, height = 5, size_units = 'in', 
-                        res = 100, ...) {
+                        res = 100, input_shapefile = NULL, lwd_shapefile = 2, 
+                        col_shapefile = 'black', ...) {
   # Process the user graphical parameters that may be passed in the call
   ## Graphical parameters to exclude
   excludedArgs <- c("cex", "cex.axis", "cex.lab", "cex.main", "col", "din", "fig", "fin", "lab", "las", "lty", "lwd", "mai", "mar", "mgp", "new", "oma", "ps", "tck")
@@ -603,6 +612,13 @@ PlotEquiMap <- function(var, lon, lat, varu = NULL, varv = NULL,
     }
   }
 
+  # Check input_shapefile
+  if (!is.null(input_shapefile)) {
+    if (!is.character(input_shapefile)) {
+      input_shapefile <- NULL
+      warning("The 'input_shapefile' must be character string.")
+    }
+  }
   #library(GEOmap)
   #library(geomapdata)
   #library(maps)
@@ -797,6 +813,27 @@ PlotEquiMap <- function(var, lon, lat, varu = NULL, varv = NULL,
       counter <- counter + 1
     }
   }
+
+  #
+  #  shapefile
+  # ~~~~~~~~~~
+  #
+  if (!is.null(input_shapefile)) {
+  if (!is.numeric(lwd_shapefile)) {
+    lwd_shapefile = 2
+    warning(" Parameter 'lwd_shapefile' must be numeric, it has been set to 2.")
+  }
+  # plot the boundary of the region of interest from shapefile
+  # Read the shapefile 
+  shp1 <- readOGR(dsn = input_shapefile, layer = 'RDD_ETRS89')
+  warning("The input_shapefile was read with the layer named RDD_ETRS89.")
+  s1 <- spTransform(shp1, CRS("+proj=longlat"))
+  warning("The input_shapefile was transformed with CRS: +proj=longlat.")
+  roi <- SpatialPolygons2map(s1)
+  # Plot the region of interest
+  map(roi, interior = FALSE, add = TRUE, lwd = lwd_shapefile,  col = col_shapefile)
+  }
+
   #
   #  PlotWind
   # ~~~~~~~~~~
diff --git a/data/shapefile_sample/RDD_ETRS89.dbf b/data/shapefile_sample/RDD_ETRS89.dbf
new file mode 100644
index 0000000000000000000000000000000000000000..685e26bf8024facb258d3204eb66820f21393e59
Binary files /dev/null and b/data/shapefile_sample/RDD_ETRS89.dbf differ
diff --git a/data/shapefile_sample/RDD_ETRS89.prj b/data/shapefile_sample/RDD_ETRS89.prj
new file mode 100644
index 0000000000000000000000000000000000000000..f2bf688c9ca2ebddc29739980ff80223cf5ff443
--- /dev/null
+++ b/data/shapefile_sample/RDD_ETRS89.prj
@@ -0,0 +1 @@
+PROJCS["Datum_73_Hayford_Gauss_IGeoE",GEOGCS["GCS_Datum_73",DATUM["D_Datum_73",SPHEROID["International_1924",6378388.0,297.0]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",200180.598],PARAMETER["False_Northing",299913.01],PARAMETER["Central_Meridian",-8.131906111111112],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",39.66666666666666],UNIT["Meter",1.0]]
\ No newline at end of file
diff --git a/data/shapefile_sample/RDD_ETRS89.sbn b/data/shapefile_sample/RDD_ETRS89.sbn
new file mode 100644
index 0000000000000000000000000000000000000000..191fcd71c057ca7d19043b035029d56803a7dfaf
Binary files /dev/null and b/data/shapefile_sample/RDD_ETRS89.sbn differ
diff --git a/data/shapefile_sample/RDD_ETRS89.sbx b/data/shapefile_sample/RDD_ETRS89.sbx
new file mode 100644
index 0000000000000000000000000000000000000000..7aa46006b1499ff60cb275c65e99880bf43991f1
Binary files /dev/null and b/data/shapefile_sample/RDD_ETRS89.sbx differ
diff --git a/data/shapefile_sample/RDD_ETRS89.shp b/data/shapefile_sample/RDD_ETRS89.shp
new file mode 100644
index 0000000000000000000000000000000000000000..f280ab7de7191ad8fec5941f1474ef2e95150830
Binary files /dev/null and b/data/shapefile_sample/RDD_ETRS89.shp differ
diff --git a/data/shapefile_sample/RDD_ETRS89.shp.xml b/data/shapefile_sample/RDD_ETRS89.shp.xml
new file mode 100644
index 0000000000000000000000000000000000000000..89aba32a7f7a9151b74b13afee69cd9d29196a69
--- /dev/null
+++ b/data/shapefile_sample/RDD_ETRS89.shp.xml
@@ -0,0 +1,2 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<metadata xml:lang="pt"><Esri><CreaDate>20141103</CreaDate><CreaTime>09591000</CreaTime><SyncOnce>FALSE</SyncOnce><SyncDate>20110204</SyncDate><SyncTime>11112200</SyncTime><ModDate>20110204</ModDate><ModTime>11112200</ModTime><DataProperties><lineage><Process ToolSource="C:\Programas\ArcGIS\ArcToolbox\Toolboxes\Data Management Tools.tbx\CopyFeatures" Date="20081104" Time="141222">CopyFeatures \\oceanus\e$\Ortos\SIG_IVDP.gdb\FreguesiasRDD_polygon "Database Connections\sdeAdmin@OCEANUS.sde\sde.SDE.FreguesiasRDD_polygon" # 0 0 0</Process><Process ToolSource="C:\Program Files\ArcGIS\ArcToolbox\Toolboxes\Data Management Tools.tbx\CopyFeatures" Date="20110201" Time="124337">CopyFeatures "Database Connections\ACanha@oceanus.sde\sde.SDE.FreguesiasRDD" "C:\Documents and Settings\Antonio Canha\Desktop\ESRI\sde_SDE_FreguesiasRDD.shp" # 0 0 0</Process><Process ToolSource="C:\Program Files\ArcGIS\Desktop10.0\ArcToolbox\Toolboxes\Data Management Tools.tbx\CopyFeatures" Date="20110204" Time="111127">CopyFeatures "C:\Documents and Settings\cmatias\Desktop\ESRI\sde_SDE_FreguesiasRDD.shp" "C:\Documents and Settings\cmatias\Desktop\IVDP - Instituto dos Vinhos do Douro e do Porto\Dados.gdb\sde_SDE_FreguesiasRDD" # 0 0 0</Process><Process ToolSource="c:\program files (x86)\arcgis\desktop10.0\ArcToolbox\Toolboxes\Data Management Tools.tbx\CalculateField" Date="20130531" Time="100156">CalculateField "Database Connections\Connection to geodouro.sde\IVDP_SDE.SDE.FreguesiasRDD" Dicofre Left( [Dicofre], 6) VB #</Process></lineage><itemProps><itemName Sync="TRUE">sde_SDE_FreguesiasRDD</itemName><imsContentType Sync="TRUE">002</imsContentType><itemLocation><linkage Sync="TRUE">Server=PROMETEU; Service=5151; Database=IVDP_SDE; User=sde; Version=sde.DEFAULT</linkage><protocol Sync="TRUE">ArcSDE Connection</protocol></itemLocation></itemProps><coordRef><type Sync="TRUE">Projected</type><geogcsn Sync="TRUE">GCS_Datum_73</geogcsn><projcsn Sync="TRUE">Datum_73_Hayford_Gauss_IGeoE</projcsn><peXml Sync="TRUE">&lt;ProjectedCoordinateSystem xsi:type='typens:ProjectedCoordinateSystem' xmlns:xsi='http://www.w3.org/2001/XMLSchema-instance' xmlns:xs='http://www.w3.org/2001/XMLSchema' xmlns:typens='http://www.esri.com/schemas/ArcGIS/10.0'&gt;&lt;WKT&gt;PROJCS[&amp;quot;Datum_73_Hayford_Gauss_IGeoE&amp;quot;,GEOGCS[&amp;quot;GCS_Datum_73&amp;quot;,DATUM[&amp;quot;D_Datum_73&amp;quot;,SPHEROID[&amp;quot;International_1924&amp;quot;,6378388.0,297.0]],PRIMEM[&amp;quot;Greenwich&amp;quot;,0.0],UNIT[&amp;quot;Degree&amp;quot;,0.0174532925199433]],PROJECTION[&amp;quot;Transverse_Mercator&amp;quot;],PARAMETER[&amp;quot;False_Easting&amp;quot;,200180.598],PARAMETER[&amp;quot;False_Northing&amp;quot;,299913.01],PARAMETER[&amp;quot;Central_Meridian&amp;quot;,-8.131906111111112],PARAMETER[&amp;quot;Scale_Factor&amp;quot;,1.0],PARAMETER[&amp;quot;Latitude_Of_Origin&amp;quot;,39.66666666666666],UNIT[&amp;quot;Meter&amp;quot;,1.0],AUTHORITY[&amp;quot;ESRI&amp;quot;,102160]]&lt;/WKT&gt;&lt;XOrigin&gt;-5423200&lt;/XOrigin&gt;&lt;YOrigin&gt;-14095100&lt;/YOrigin&gt;&lt;XYScale&gt;10000&lt;/XYScale&gt;&lt;ZOrigin&gt;-100000&lt;/ZOrigin&gt;&lt;ZScale&gt;10000&lt;/ZScale&gt;&lt;MOrigin&gt;-100000&lt;/MOrigin&gt;&lt;MScale&gt;10000&lt;/MScale&gt;&lt;XYTolerance&gt;0.001&lt;/XYTolerance&gt;&lt;ZTolerance&gt;0.001&lt;/ZTolerance&gt;&lt;MTolerance&gt;0.001&lt;/MTolerance&gt;&lt;HighPrecision&gt;true&lt;/HighPrecision&gt;&lt;WKID&gt;102160&lt;/WKID&gt;&lt;/ProjectedCoordinateSystem&gt;</peXml></coordRef></DataProperties><ArcGISFormat>1.0</ArcGISFormat></Esri><idinfo><native Sync="TRUE">Microsoft Windows XP Version 5.1 (Build 2600) Service Pack 3; ESRI ArcCatalog 9.3.1.3000</native><descript><langdata Sync="TRUE">pt</langdata><abstract>REQUIRED: A brief narrative summary of the data set.</abstract><purpose>REQUIRED: A summary of the intentions with which the data set was developed.</purpose></descript><citation><citeinfo><origin>REQUIRED: The name of an organization or individual that developed the data set.</origin><pubdate>REQUIRED: The date when the data set is published or otherwise made available for release.</pubdate><title Sync="TRUE">sde_SDE_FreguesiasRDD</title><ftname Sync="TRUE">sde_SDE_FreguesiasRDD</ftname><geoform Sync="TRUE">vector digital data</geoform><onlink Sync="TRUE">\\RSCAD03\C$\Documents and Settings\Antonio Canha\Desktop\ESRI\sde_SDE_FreguesiasRDD.shp</onlink></citeinfo></citation><timeperd><current>REQUIRED: The basis on which the time period of content information is determined.</current><timeinfo><sngdate><caldate>REQUIRED: The year (and optionally month, or month and day) for which the data set corresponds to the ground.</caldate></sngdate></timeinfo></timeperd><status><progress>REQUIRED: The state of the data set.</progress><update>REQUIRED: The frequency with which changes and additions are made to the data set after the initial data set is completed.</update></status><spdom><bounding><westbc Sync="TRUE">REQUIRED: Western-most coordinate of the limit of coverage expressed in longitude.</westbc><eastbc Sync="TRUE">REQUIRED: Eastern-most coordinate of the limit of coverage expressed in longitude.</eastbc><northbc Sync="TRUE">REQUIRED: Northern-most coordinate of the limit of coverage expressed in latitude.</northbc><southbc Sync="TRUE">REQUIRED: Southern-most coordinate of the limit of coverage expressed in latitude.</southbc></bounding></spdom><keywords><theme><themekt>REQUIRED: Reference to a formally registered thesaurus or a similar authoritative source of theme keywords.</themekt><themekey>REQUIRED: Common-use word or phrase used to describe the subject of the data set.</themekey></theme></keywords><accconst>REQUIRED: Restrictions and legal prerequisites for accessing the data set.</accconst><useconst>REQUIRED: Restrictions and legal prerequisites for using the data set after access is granted.</useconst><natvform Sync="TRUE">Shapefile</natvform></idinfo><dataIdInfo><envirDesc Sync="TRUE">Microsoft Windows XP Version 5.1 (Build 2600) Service Pack 3; ESRI ArcGIS 10.0.0.2414</envirDesc><dataLang><languageCode Sync="TRUE" value="pt" country="PT"/></dataLang><idCitation><resTitle Sync="TRUE">sde_SDE_FreguesiasRDD</resTitle><presForm><PresFormCd Sync="TRUE" value="005"/></presForm></idCitation><spatRpType><SpatRepTypCd Sync="TRUE" value="001"/></spatRpType><descKeys><thesaName uuidref="723f6998-058e-11dc-8314-0800200c9a66"/><keyword Sync="TRUE">002</keyword></descKeys></dataIdInfo><metainfo><langmeta Sync="TRUE">pt</langmeta><metstdn Sync="TRUE">FGDC Content Standards for Digital Geospatial Metadata</metstdn><metstdv Sync="TRUE">FGDC-STD-001-1998</metstdv><mettc Sync="TRUE">local time</mettc><metc><cntinfo><cntorgp><cntper>REQUIRED: The person responsible for the metadata information.</cntper><cntorg>REQUIRED: The organization responsible for the metadata information.</cntorg></cntorgp><cntaddr><addrtype>REQUIRED: The mailing and/or physical address for the organization or individual.</addrtype><city>REQUIRED: The city of the address.</city><state>REQUIRED: The state or province of the address.</state><postal>REQUIRED: The ZIP or other postal code of the address.</postal></cntaddr><cntvoice>REQUIRED: The telephone number by which individuals can speak to the organization or individual.</cntvoice></cntinfo></metc><metd Sync="TRUE">20110201</metd></metainfo><mdLang><languageCode Sync="TRUE" value="pt"/></mdLang><mdStanName Sync="TRUE">ISO 19115 Geographic Information - Metadata</mdStanName><mdStanVer Sync="TRUE">DIS_ESRI1.0</mdStanVer><mdChar><CharSetCd Sync="TRUE" value="004"/></mdChar><mdHrLv><ScopeCd Sync="TRUE" value="005"/></mdHrLv><mdHrLvName Sync="TRUE">dataset</mdHrLvName><distinfo><resdesc Sync="TRUE">Downloadable Data</resdesc><stdorder><digform><digtinfo><transize Sync="TRUE">0,000</transize><dssize Sync="TRUE">0,000</dssize></digtinfo></digform></stdorder></distinfo><distInfo><distributor><distorTran><onLineSrc><orDesc Sync="TRUE">002</orDesc><linkage Sync="TRUE">file://\\RSCAD03\C$\Documents and Settings\Antonio Canha\Desktop\ESRI\sde_SDE_FreguesiasRDD.shp</linkage><protocol Sync="TRUE">Local Area Network</protocol></onLineSrc><transSize Sync="TRUE">0,000</transSize></distorTran><distorFormat><formatName Sync="TRUE">Shapefile</formatName></distorFormat></distributor><distFormat><formatName Sync="TRUE">File Geodatabase Feature Class</formatName></distFormat></distInfo><spdoinfo><direct Sync="TRUE">Vector</direct><ptvctinf><esriterm Name="sde_SDE_FreguesiasRDD"><efeatyp Sync="TRUE">Simple</efeatyp><efeageom Sync="TRUE">4</efeageom><esritopo Sync="TRUE">FALSE</esritopo><efeacnt Sync="TRUE">0</efeacnt><spindex Sync="TRUE">TRUE</spindex><linrefer Sync="TRUE">TRUE</linrefer></esriterm></ptvctinf></spdoinfo><spref><horizsys><cordsysn><geogcsn Sync="TRUE">GCS_Datum_73</geogcsn><projcsn Sync="TRUE">Datum_73_Hayford_Gauss_IGeoE</projcsn></cordsysn><planar><planci><plance Sync="TRUE">coordinate pair</plance><plandu Sync="TRUE">meters</plandu><coordrep><absres Sync="TRUE">0.000000</absres><ordres Sync="TRUE">0.000000</ordres></coordrep></planci></planar><geodetic><horizdn Sync="TRUE">D_Datum_73</horizdn><ellips Sync="TRUE">International_1924</ellips><semiaxis Sync="TRUE">6378388.000000</semiaxis><denflat Sync="TRUE">297.000000</denflat></geodetic></horizsys><vertdef><altsys><altenc Sync="TRUE">Explicit elevation coordinate included with horizontal coordinates</altenc><altres Sync="TRUE">0.000100</altres></altsys></vertdef></spref><refSysInfo><RefSystem><refSysID><identCode Sync="TRUE" code="102160">Datum_73_Hayford_Gauss_IGeoE</identCode><idCodeSpace Sync="TRUE">ESRI</idCodeSpace><idVersion Sync="TRUE">10.0.0</idVersion></refSysID></RefSystem></refSysInfo><spatRepInfo><VectSpatRep><topLvl><TopoLevCd Sync="TRUE" value="001"/></topLvl><geometObjs Name="sde_SDE_FreguesiasRDD"><geoObjTyp><GeoObjTypCd Sync="TRUE" value="002"/><geoObjCnt Sync="TRUE">0</geoObjCnt></geoObjTyp><geoObjCnt Sync="TRUE">0</geoObjCnt></geometObjs></VectSpatRep></spatRepInfo><eainfo><detailed Name="sde_SDE_FreguesiasRDD"><enttyp><enttypl Sync="TRUE">sde_SDE_FreguesiasRDD</enttypl><enttypt Sync="TRUE">Feature Class</enttypt><enttypc Sync="TRUE">0</enttypc></enttyp><attr><attrlabl Sync="TRUE">OBJECTID</attrlabl><attalias Sync="TRUE">OBJECTID</attalias><attrtype Sync="TRUE">OID</attrtype><attwidth Sync="TRUE">4</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale><attrdef Sync="TRUE">Internal feature number.</attrdef><attrdefs Sync="TRUE">ESRI</attrdefs><attrdomv><udom Sync="TRUE">Sequential unique whole numbers that are automatically generated.</udom></attrdomv></attr><attr><attrlabl Sync="TRUE">ID</attrlabl><attalias Sync="TRUE">ID</attalias><attrtype Sync="TRUE">Double</attrtype><attwidth Sync="TRUE">8</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale></attr><attr><attrlabl Sync="TRUE">SUB_REGIÃO</attrlabl><attalias Sync="TRUE">SUB_REGIÃO</attalias><attrtype Sync="TRUE">String</attrtype><attwidth Sync="TRUE">50</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale></attr><attr><attrlabl Sync="TRUE">ÁREA</attrlabl><attalias Sync="TRUE">ÁREA</attalias><attrtype Sync="TRUE">Double</attrtype><attwidth Sync="TRUE">8</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale></attr><attr><attrlabl Sync="TRUE">PERIMETRO</attrlabl><attalias Sync="TRUE">PERIMETRO</attalias><attrtype Sync="TRUE">Double</attrtype><attwidth Sync="TRUE">8</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale></attr><attr><attrlabl Sync="TRUE">SECÇÃO</attrlabl><attalias Sync="TRUE">SECÇÃO</attalias><attrtype Sync="TRUE">String</attrtype><attwidth Sync="TRUE">50</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale></attr><attr><attrlabl Sync="TRUE">SECTOR</attrlabl><attalias Sync="TRUE">SECTOR</attalias><attrtype Sync="TRUE">String</attrtype><attwidth Sync="TRUE">50</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale></attr><attr><attrlabl Sync="TRUE">Freguesia</attrlabl><attalias Sync="TRUE">Freguesia</attalias><attrtype Sync="TRUE">String</attrtype><attwidth Sync="TRUE">200</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale></attr><attr><attrlabl Sync="TRUE">Concelho</attrlabl><attalias Sync="TRUE">Concelho</attalias><attrtype Sync="TRUE">String</attrtype><attwidth Sync="TRUE">200</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale></attr><attr><attrlabl Sync="TRUE">Distrito</attrlabl><attalias Sync="TRUE">Distrito</attalias><attrtype Sync="TRUE">String</attrtype><attwidth Sync="TRUE">200</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale></attr><attr><attrlabl Sync="TRUE">Dicofre</attrlabl><attalias Sync="TRUE">Dicofre</attalias><attrtype Sync="TRUE">String</attrtype><attwidth Sync="TRUE">50</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale></attr><attr><attrlabl Sync="TRUE">OBSERVAÇÕE</attrlabl><attalias Sync="TRUE">OBSERVAÇÕE</attalias><attrtype Sync="TRUE">String</attrtype><attwidth Sync="TRUE">200</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale></attr><attr><attrlabl Sync="TRUE">ID1</attrlabl><attalias Sync="TRUE">ID1</attalias><attrtype Sync="TRUE">Double</attrtype><attwidth Sync="TRUE">8</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale></attr><attr><attrlabl Sync="TRUE">Shape</attrlabl><attalias Sync="TRUE">Shape</attalias><attrtype Sync="TRUE">Geometry</attrtype><attwidth Sync="TRUE">0</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale><attrdef Sync="TRUE">Feature geometry.</attrdef><attrdefs Sync="TRUE">ESRI</attrdefs><attrdomv><udom Sync="TRUE">Coordinates defining the features.</udom></attrdomv></attr><attr><attrlabl Sync="TRUE">Shape_Length</attrlabl><attalias Sync="TRUE">Shape_Length</attalias><attrtype Sync="TRUE">Double</attrtype><attwidth Sync="TRUE">8</attwidth><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale><attrdef Sync="TRUE">Length of feature in internal units.</attrdef><attrdefs Sync="TRUE">ESRI</attrdefs><attrdomv><udom Sync="TRUE">Positive real numbers that are automatically generated.</udom></attrdomv></attr><attr><attrlabl Sync="TRUE">Shape_area</attrlabl><attalias Sync="TRUE">Shape_Area</attalias><attrtype Sync="TRUE">Double</attrtype><attwidth Sync="TRUE">8</attwidth><attrdef Sync="TRUE">Area of feature in internal units squared.</attrdef><attrdefs Sync="TRUE">ESRI</attrdefs><attrdomv><udom Sync="TRUE">Positive real numbers that are automatically generated.</udom></attrdomv><atprecis Sync="TRUE">0</atprecis><attscale Sync="TRUE">0</attscale></attr></detailed></eainfo><mdDateSt Sync="TRUE">20110204</mdDateSt></metadata>
diff --git a/data/shapefile_sample/RDD_ETRS89.shx b/data/shapefile_sample/RDD_ETRS89.shx
new file mode 100644
index 0000000000000000000000000000000000000000..c786c6fb851414c63d276cf094741793b9fa8a4e
Binary files /dev/null and b/data/shapefile_sample/RDD_ETRS89.shx differ
diff --git a/man/PlotEquiMap.Rd b/man/PlotEquiMap.Rd
index cf45ead4b3555ba1417ebcc18d7b01140550b11c..86dcc40529907c7f569dda9b8c0fb21450783ef4 100644
--- a/man/PlotEquiMap.Rd
+++ b/man/PlotEquiMap.Rd
@@ -22,7 +22,7 @@ PlotEquiMap(var, lon, lat, varu = NULL, varv = NULL, toptitle = NULL,
   bar_extra_margin = rep(0, 4), boxlim = NULL, boxcol = "purple2",
   boxlwd = 5, margin_scale = rep(1, 4), title_scale = 1, numbfig = NULL,
   fileout = NULL, width = 8, height = 5, size_units = "in", res = 100,
-  ...)
+  input_shapefile = NULL, lwd_shapefile = 2, col_shapefile = "black", ...)
 }
 \arguments{
 \item{var}{Array with the values at each cell of a grid on a regular 
@@ -212,6 +212,15 @@ the corresponding device.}
 \item{res}{Resolution of the device (file or window) to plot in. See 
 ?Devices and the creator function of the corresponding device.}
 
+\item{input_shapefile}{A character string describing the pathway to your
+shapefile of the region of interest.}
+
+\item{lwd_shapefile}{Line width of the shapefile boundary lines. 
+Defaults to 2.}
+
+\item{col_shapefile}{Colour of the shapefile boundary lines. 
+Defaults to 'black'.}
+
 \item{\dots}{Arguments to be passed to the method. Only accepts the following 
 graphical parameters:\cr
 adj ann ask bg bty cex.sub cin col.axis col.lab col.main col.sub cra crt