R version 3.2.3 (2015-12-10) -- "Wooden Christmas-Tree" Copyright (C) 2015 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. [Previously saved workspace restored] > ################################################################ > ### THIS FILE CONTAINS CODE TO DRAW A CUSTOMIZED BOXPLOT DISPLAY. It > ### draws side-by-side boxplots for the various groups, without median > ### lines or whiskers. This means that to see the boxplot, we have o > ### give it a colored background. We also draw points at the medians > ### of the boxplots, and connect these with lines. > ### NOTE: The user should make sure the data is sorted by group before > ### using this code. > ################################################################ > > batch <- 'y' > > amet_base <- Sys.getenv("AMETBASE") > dbase <- Sys.getenv("AMET_DATABASE") > out_dir <- Sys.getenv("AMET_OUT") > ametRinput <- Sys.getenv("AMETRINPUT") > source(ametRinput) > > #states <- c("All") > #states <- c("AL","AK","AZ","AR","CA","CO","CT","DE","DC","FL","GA","HI","ID","IL","IN","IA","KS","KY","LA","ME","MD","MA","MI","MN","MS","MO","MT","NE","NV","NH","NJ","NM","NY","NC","ND","OH","OK","OR","PA","RI","SC","SD","TN","TX","UT","VT","VA","WA","WV","WI","WY") > #states <- c("AZ","CA","CO","ID","IA","KS","MT","NE","NV","NM","ND","OK","OR","SD","TX","UT","WA","WY") # Western States > #states <- c("CA") > > ###################################### > ####### Configuration Options ######## > ###################################### > run_name1 <- Sys.getenv("AMET_PROJECT") # AMET project name > run_name2 <- Sys.getenv("AMET_PROJECT2") # Additional run to include on plot > print(run_name2) [1] "" > #start_date <- "20060101" # Set the start date of the analysis > #end_date <- "20061231" # Set the end date of the analysis > #batch_query <- c("month=1","month=2","month=3","month=4","month=5","month=6","month=7","month=8","month=9","month=10","month=11","month=12", > # "(month=12 or month=1 or month=2)", > # "(month=3 or month=4 or month=5)", > # "(month=6 or month=7 or month=8)", > # "(month=9 or month=10 or month=11)") > #batch_names <- c("January","February","March","April","May","June","July","August","September","October","November","December", > # "Winter","Spring","Summer","Fall") > > #hourly_ozone_analysis <- 'n' # Flag to include hourly ozone analysis > #daily_ozone_analysis <- 'y' # Flag to include daily ozone analysis > #aerosol_analysis <- 'y' # Flag to include aerosol analysis > #dep_analysis <- 'y' # Flag to include analysis of deposition performance > #gas_analysis <- 'y' # Flag to include gas analysis > > #ozone_averaging <- 'n' # Flag to average ozone data; options are n (none), d (daily), m (month), s (season), y (all) > #aerosol_averaging <- 'n' # Flag to average aerosol data; options are n (none), d (daily), m (month), s (season), y (all) > #deposition_averaging <- 'n' # Flag to sum deposition data; options are n (none), d (daily), m (month), s (season), y (all) > #gas_averaging <- 'n' # Flag to average gas data; options are n (none), d (daily), m (month), s (season), y (all) > ####################################### > > ##################### > ### Other Options ### > ##################### > > ### Custom Title ### > custom_title <- "" > ############################## > > ### Main Database Query String. ### > query_string<-paste(" and s.stat_id=d.stat_id and d.ob_dates >=",start_date,"and d.ob_datee <=",end_date,additional_query,sep=" ") > > ### Set and create output directory ### > #out_dir <- paste(out_dir,"spatial_plots",sep="/") > mkdir_main_command <- paste("mkdir -p",out_dir,sep=" ") > system(mkdir_main_command) # This will create a subdirectory with the name of the project > ####################################### > > run_script_command1 <- paste(amet_base,"/R_analysis_code/AQ_Stats_Plots.R",sep="") > run_script_command2 <- paste(amet_base,"/R_analysis_code/AQ_Plot_Spatial.R",sep="") > run_script_command3 <- paste(amet_base,"/R_analysis_code/AQ_Plot_Spatial_Diff.R",sep="") > run_script_command4 <- paste(amet_base,"/R_analysis_code/AQ_Plot_Spatial_MtoM.R",sep="") > run_script_command5 <- paste(amet_base,"/R_analysis_code/AQ_Plot_Spatial_Ratio.R",sep="") > run_script_command6 <- paste(amet_base,"/R_analysis_code/AQ_Plot_Spatial_interactive.R",sep="") > run_script_command7 <- paste(amet_base,"/R_analysis_code/AQ_Plot_Spatial_Diff_interactive.R",sep="") > run_script_command8 <- paste(amet_base,"/R_analysis_code/AQ_Stats_Plots_interactive.R",sep="") > > ########################################################################################## > ### This portion of the code will create monthly spatial plots for the various species ### > ########################################################################################## > #for (m in 1:length(batch_query)) { > # mkdir_command <- paste("mkdir ",out_dir,"/",batch_names[m],sep="") > # system(mkdir_command) > #} > if (hourly_ozone_analysis == 'y') { + averaging <- ozone_averaging + for (m in 1:length(batch_query)) { + species <- "O3" + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + dates <- batch_names[m] + network_names <- c("AQS_Hourly","CASTNET_Hourly","NAPS","EMEP_Hourly") + network_label <- c("AQS","CASTNET","NAPS","EMEP") + pid <- "multi_network" + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + print(query) + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + } + } [1] " and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 )" ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Loading required package: RMySQL Loading required package: DBI Loading required package: fields Loading required package: spam Loading required package: dotCall64 Loading required package: grid Spam version 2.5-1 (2019-12-12) is loaded. Type 'help( Spam)' or 'demo( spam)' for a short introduction and overview of this package. Help for individual functions is also obtained by adding the suffix '.spam' to the function name, e.g. 'help( chol.spam)'. Attaching package: ‘spam’ The following objects are masked from ‘package:base’: backsolve, forwardsolve Loading required package: maps See https://github.com/NCAR/Fields for an extensive vignette, other supplements and source code Loading required package: mapdata [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.O3_ob, d.O3_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.O3_ob is not NULL and d.network='AQS_Hourly' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.O3_ob' in 'field list' In addition: Warning message: In FUN(X[[i]], ...) : failed to assign RegisteredNativeSymbol for toeplitz to toeplitz since toeplitz is already defined in the ‘spam’ namespace Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.O3_ob, d.O3_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.O3_ob is not NULL and d.network='AQS_Hourly' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.O3_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory > > if (daily_ozone_analysis == 'y') { + averaging <- ozone_averaging + for (m in 1:length(batch_query)) { + species <- "O3_1hrmax" + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + dates <- batch_names[m] + network_names <- c("AQS_Daily_O3","CASTNET_Daily","NAPS_Daily_O3","EMEP_Daily_O3") + network_label <- c("AQS","CASTNET","NAPS","EMEP") + pid <- "multi_network" + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + # print(query) + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + } + for (m in 1:length(batch_query)) { + species <- "O3_8hrmax" + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + dates <- batch_names[m] + network_names <- c("AQS_Daily_O3","CASTNET_Daily","NAPS_Daily_O3","EMEP_Daily_O3") + network_label <- c("AQS","CASTNET","NAPS","EMEP") + pid <- "multi_network" + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + } + } ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.O3_1hrmax_ob, d.O3_1hrmax_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.O3_1hrmax_ob is not NULL and d.network='AQS_Daily_O3' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.O3_1hrmax_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.O3_1hrmax_ob, d.O3_1hrmax_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.O3_1hrmax_ob is not NULL and d.network='AQS_Daily_O3' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.O3_1hrmax_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.O3_8hrmax_ob, d.O3_8hrmax_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.O3_8hrmax_ob is not NULL and d.network='AQS_Daily_O3' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.O3_8hrmax_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.O3_8hrmax_ob, d.O3_8hrmax_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.O3_8hrmax_ob is not NULL and d.network='AQS_Daily_O3' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.O3_8hrmax_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory > > if (aerosol_analysis == 'y') { + averaging <- aerosol_averaging + for (m in 1:length(batch_query)) { + species_list <- c("SO4","NO3","NH4") + for (i in 1:length(species_list)) { + species <- species_list[i] + dates <- batch_names[m] + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + network_names <- c("IMPROVE","CSN","CASTNET") + network_label <- c("IMPROVE","CSN","CASTNET") + pid <- "multi_network" + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + print(query) + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + } + } + for (m in 1:length(batch_query)) { + species_list <- c("HNO3","TNO3") + for (i in 1:length(species_list)) { + species <- species_list[i] + dates <- batch_names[m] + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + network_names <- c("CASTNET") + network_label <- c("CASTNET") + pid <- network_names + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + } + } + for (m in 1:length(batch_query)) { + species <- "PM_TOT" + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + dates <- batch_names[m] + network_names <- c("IMPROVE","CSN","AQS_Daily","AQS_Hourly","NAPS") + network_label <- c("IMPROVE","CSN","AQS_Daily","AQS_Hourly","NAPS") + pid <- "multi_network" + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + } + for (m in 1:length(batch_query)) { + species_list <- c("EC","OC","TC") + for (i in 1:length(species_list)) { + species <- species_list[i] + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + dates <- batch_names[m] + network_names <- c("IMPROVE","CSN") + network_label <- c("IMPROVE","CSN") + pid <- "multi_network" + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + print(query) + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + if (ratio_plots == 'y') { try(source(run_script_command5)) } + } + } + } [1] " and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 )" ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.SO4_ob, d.SO4_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.SO4_ob is not NULL and d.network='IMPROVE' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.SO4_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.SO4_ob, d.SO4_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.SO4_ob is not NULL and d.network='IMPROVE' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.SO4_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory [1] " and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 )" ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.NO3_ob, d.NO3_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.NO3_ob is not NULL and d.network='IMPROVE' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.NO3_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.NO3_ob, d.NO3_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.NO3_ob is not NULL and d.network='IMPROVE' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.NO3_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory [1] " and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 )" ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.NH4_ob, d.NH4_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.NH4_ob is not NULL and d.network='IMPROVE' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.NH4_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.NH4_ob, d.NH4_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.NH4_ob is not NULL and d.network='IMPROVE' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.NH4_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.HNO3_ob, d.HNO3_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.HNO3_ob is not NULL and d.network='CASTNET' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.HNO3_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.HNO3_ob, d.HNO3_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.HNO3_ob is not NULL and d.network='CASTNET' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.HNO3_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.TNO3_ob, d.TNO3_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.TNO3_ob is not NULL and d.network='CASTNET' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.TNO3_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## [1] "SELECT d.network,d.stat_id,s.lat,s.lon,d.ob_dates,d.ob_datee,d.ob_hour,d.month , d.TNO3_ob, d.TNO3_mod ,d.POCode,s.state from PXM_LandUse_AQ as d, site_metadata as s WHERE d.TNO3_ob is not NULL and d.network='CASTNET' and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: Unknown column 'd.TNO3_ob' in 'field list' Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory [1] " and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 )" ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory [1] " and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 )" ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory [1] " and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 )" ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory > > if (dep_analysis == 'y') { + averaging <- deposition_averaging + for (m in 1:length(batch_query)) { + species_list <- c("SO4_dep","NO3_dep","NH4_dep","Precip") + for (i in 1:length(species_list)) { + species <- species_list[i] + dates <- batch_names[m] + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + network_names <- c("NADP") + network_label <- c("NADP") + pid <- network_label + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + print(query) + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + } + } + } [1] " and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 )" ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory [1] " and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 )" ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory [1] " and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 )" ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory [1] " and s.stat_id=d.stat_id and d.ob_dates >= 20100310 and d.ob_datee <= 20100410 and ( d.ob_dates between 20100310 and 20100410 and d.ob_datee between 20100310 and 20100410 )" ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory > > if (gas_analysis == 'y') { + averaging <- gas_averaging + for (m in 1:length(batch_query)) { + species_list <- c("SO2","NO","NO2","NOY","NOX","CO") + for (i in 1:length(species_list)) { + species <- species_list[i] + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + dates <- batch_names[m] + network_names <- c("AQS_Hourly","NAPS") + network_label <- c("AQS","NAPS") + pid <- "multi_network" + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + } + } + for (m in 1:length(batch_query)) { + species <- "SO2" + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + dates <- batch_names[m] + network_names <- c("CASTNET") + network_label <- c("CASTNET") + pid <- network_names + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + } + } ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory > > if (AE6_analysis == 'y') { + species_list <- c("Na","Cl","Fe","Al","Si","Ti","Ca","Mg","K","Mn","soil","NaCl","other","ncom","other_rem") + for (m in 1:length(batch_query)) { + for (i in 1:length(species_list)) { + species <- species_list[i] + dates <- batch_names[m] + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + network_names <- c("IMPROVE","CSN") + network_label <- c("IMPROVE","CSN") + pid <- "multi_network" + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + if (ratio_plots == 'y') { try(source(run_script_command5)) } + } + } + } ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory > if (AOD_analysis == 'y') { + species_list <- c("AOD_500") + for (m in 1:length(batch_query)) { + for (i in 1:length(species_list)) { + species <- species_list[i] + dates <- batch_names[m] + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + network_names <- c("AERONET") + network_label <- c("AERONET") + pid <- paste(run_name1,dates,species,sep="_") + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + } + } + } ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory > if (PAMS_analysis == 'y') { + species_list <- c("Isoprene","Ethane","Ethylene","Toluene") + for (m in 1:length(batch_query)) { + for (i in 1:length(species_list)) { + species <- species_list[i] + dates <- batch_names[m] + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + network_names <- c("AQS_Hourly") + network_label <- c("AQS_Hourly") + pid <- paste(run_name1,dates,species,sep="_") + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + } + } + species_list <- c("Isoprene","Ethane","Ethylene","Toluene","Acetaldehyde","Formaldehyde") + for (m in 1:length(batch_query)) { + for (i in 1:length(species_list)) { + species <- species_list[i] + dates <- batch_names[m] + figdir <- paste(out_dir,species,sep="/") + if (batch_names[m] != "None") { + figdir <- paste(out_dir,batch_names[m],species,sep="/") + } + mkdir_command <- paste("mkdir -p",figdir) + network_names <- c("AQS_Daily") + network_label <- c("AQS_Daily") + pid <- paste(run_name1,dates,species,sep="_") + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (stat_plots == 'y') { + try(source(run_script_command1)) + try(source(run_script_command8)) + } + if (spatial_plots == 'y') { + try(source(run_script_command2)) + try(source(run_script_command6)) + } + if ((diff_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { + try(source(run_script_command3)) + try(source(run_script_command7)) + } + if ((mtom_plots == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command4)) } + } + } + } ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory ############################## INTERACTIVE STATS PLOTS ############################## ### AMET CODE: AQ_Stats_Plots.R ### ### This code is part of the AMET-AQ system. The Stats and Plots code takes a MYSQL ### database query for a single species from one or more networks and a single simulation ### and calculates summary statistics for each site and the entire domain by network. ### Spatial plots are also produced for several select statistics, specifically MB, ME, FB, ### FE, NMB, NME, RMSE and Correlation. Images can be output as png, pdf or both. ### ### Last modified by Wyat Appel, June 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Stats_Plots_interactive.R': No such file or directory ###################################### SPATIAL PLOT ###################################### ### AMET CODE: AQ_Plot_Spatial.R ### ### This code is part of the AMET-AQ system. The Plot Spatial code takes a MYSQL database ### query for a single species from one or more networks and plots the observation value, ### model value, and difference between the model and ob for each site for each corresponding ### network. Mutiple values for a site are averaged to a single value for plotting purposes. ### The map area plotted is dynamically generated from the input data. ### ### Last modified by Wyat Appel: June, 2019 ########################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in file(filename, "r", encoding = encoding) : cannot open the connection In addition: Warning message: In file(filename, "r", encoding = encoding) : cannot open file '/home/pkasper/AMET_v14/R_analysis_code/AQ_Plot_Spatial_interactive.R': No such file or directory > > > proc.time() user system elapsed 3.472 0.264 3.987