From 53562369d69016c04bac19148d7da8d80e80cc49 Mon Sep 17 00:00:00 2001 From: Charles Perretti Date: Wed, 24 Jan 2018 10:12:47 -0500 Subject: [PATCH] Fixed extrapolation area plot for NE shelf --- R/Plot_data_and_knots.R | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) 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)) ){