diff --git a/R/Plot_data_and_knots.R b/R/Plot_data_and_knots.R
index 839fd19..1dfb988 100644
--- a/R/Plot_data_and_knots.R
+++ b/R/Plot_data_and_knots.R
@@ -25,10 +25,10 @@ Plot_data_and_knots = function( Extrapolation_List, Spatial_List, Data_Geostat,
   # Plot data and grid
   png( file=paste0(PlotDir,Plot1_name), width=6, height=6, res=200, units="in")
     par( mfrow=c(2,2), mar=c(3,3,2,0), mgp=c(1.75,0.25,0) )
-    plot( Extrapolation_List$Data_Extrap[which(Extrapolation_List$Area_km2_x>0),c('Lon','Lat')], cex=0.01, main="Extrapolation (Lat-Lon)" )
+    plot( Extrapolation_List$Data_Extrap[which(Extrapolation_List$a_el > 0),c('Lon','Lat')], cex=0.01, main="Extrapolation (Lat-Lon)" )
     map( "world", add=TRUE )
     if( !any(is.na(Extrapolation_List$Data_Extrap[,c('E_km','N_km')])) ){
-      plot( Extrapolation_List$Data_Extrap[which(Extrapolation_List$Area_km2_x>0),c('E_km','N_km')], cex=0.01, main="Extrapolation (North-East)" )
+      plot( Extrapolation_List$Data_Extrap[which(Extrapolation_List$a_el > 0),c('E_km','N_km')], cex=0.01, main="Extrapolation (North-East)" )
     }
     plot( Spatial_List$loc_x, col="red", pch=20, main="Knots (North-East)")
     if( all(c('E_km','N_km')%in%names(Data_Geostat)) ){