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 SCATTER PLOT. > ################################################################ > > amet_base <- Sys.getenv("AMETBASE") > dbase <- Sys.getenv("AMET_DATABASE") > out_dir <- Sys.getenv("AMET_OUT") > ametRinput <- Sys.getenv("AMETRINPUT") > source(ametRinput) > > ################################### > ### Does not need to be changed ### > ################################### > 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 > ################################### > > ###################################### > ####### Configuration Options ######## > ###################################### > run_name1 <- Sys.getenv("AMET_PROJECT") # AMET project name > run_name2 <- Sys.getenv("AMET_PROJECT2") # Additional run to include on plot > > ### These options are set in run_info file but setting here will override that setting ### > #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 or month=2 or month=3 or month=4 or month=5 or month=6 or month=7 or month=8 or month=9 or month=10 or month=11 or month=12") > #batch_names <- c("All") > #hourly_ozone_analysis <- 'y' # 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,"scatter_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_Scatterplot.R",sep="") > run_script_command2 <- paste(amet_base,"/R_analysis_code/AQ_Scatterplot_single.R",sep="") > run_script_command3 <- paste(amet_base,"/R_analysis_code/AQ_Scatterplot_density.R",sep="") > run_script_command4 <- paste(amet_base,"/R_analysis_code/AQ_Scatterplot_bins.R",sep="") > run_script_command5 <- paste(amet_base,"/R_analysis_code/AQ_Scatterplot_percentiles.R",sep="") > run_script_command6 <- paste(amet_base,"/R_analysis_code/AQ_Scatterplot_skill.R",sep="") > run_script_command7 <- paste(amet_base,"/R_analysis_code/AQ_Scatterplot_mtom.R",sep="") > run_script_command8 <- paste(amet_base,"/R_analysis_code/AQ_Scatterplot_soil.R",sep="") > run_script_command9 <- paste(amet_base,"/R_analysis_code/AQ_Scatterplot_Multisim_plotly.R",sep="") > run_script_command10 <- paste(amet_base,"/R_analysis_code/AQ_Scatterplot_bins_plotly.R",sep="") > run_script_command11 <- paste(amet_base,"/R_analysis_code/AQ_Scatterplot_density_ggplot.R",sep="") > > ####################################################################################### > ### This portion of the code will create monthly stat plots for the various species ### > ####################################################################################### > for (m in 1:length(batch_query)) { + # mkdir_command <- paste("mkdir -p ",out_dir,"/",batch_names[m],sep="") + # cat(mkdir_command) + # system(mkdir_command) + } > if (hourly_ozone_analysis == 'y') { + max_limit <- O3_max_limit + averaging <- ozone_averaging + for (m in 1:length(batch_query)) { + species_list <- c("O3") + 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") + network_label <- c("AQS_Hourly") + pid <- network_label + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (scatter_plot == 'y') { + try(source(run_script_command1)) + try(source(run_script_command9)) + } + if (scatter_single_plot == 'y') { try(source(run_script_command2)) } + if (scatter_density_plot == 'y') { + try(source(run_script_command3)) + try(source(run_script_command11)) + } + if (scatter_bins_plot == 'y') { + try(source(run_script_command4)) + try(source(run_script_command10)) + } + if (scatter_percentiles_plot == 'y') { try(source(run_script_command5)) } + if (scatter_skill_plot == 'y') { try(source(run_script_command6)) } + if ((scatter_mtom_plot == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command7)) } + } + } + } ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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 [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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Loading required package: ggplot2 [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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 > if (daily_ozone_analysis == 'y') { + max_limit <- O3_max_limit + averaging <- ozone_averaging + species_list <- c("O3_1hrmax","O3_8hrmax") + for (m in 1:length(batch_query)) { + 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_Daily_O3") + network_label <- c("AQS_Daily") + pid <- network_label + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (scatter_plot == 'y') { + try(source(run_script_command1)) + try(source(run_script_command9)) + } + if (scatter_single_plot == 'y') { try(source(run_script_command2)) } + if (scatter_density_plot == 'y') { + try(source(run_script_command3)) + try(source(run_script_command11)) + } + if (scatter_bins_plot == 'y') { + try(source(run_script_command4)) + try(source(run_script_command10)) + } + if (scatter_percentiles_plot == 'y') { try(source(run_script_command5)) } + if (scatter_skill_plot == 'y') { try(source(run_script_command6)) } + if ((scatter_mtom_plot == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command7)) } + } + } + } ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated 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 >= and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ) ORDER BY stat_id,ob_dates,ob_hour" Error in .local(conn, statement, ...) : could not run statement: You have an error in your SQL syntax; check the manual that corresponds to your MySQL server version for the right syntax to use near 'and d.ob_datee <= and ( d.ob_dates between and and d.ob_datee between and ' at line 1 ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened > > if (aerosol_analysis == 'y') { + max_limit <- PM_max_limit + averaging <- aerosol_averaging + for (m in 1:length(batch_query)) { + species_list <- c("SO4","NO3","EC","OC","TC","PM_TOT") + 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") + network_label <- c("IMPROVE") + pid <- network_label + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (scatter_plot == 'y') { + try(source(run_script_command1)) + try(source(run_script_command9)) + } + if (scatter_single_plot == 'y') { try(source(run_script_command2)) } + if (scatter_density_plot == 'y') { + try(source(run_script_command3)) + try(source(run_script_command11)) + } + if (scatter_bins_plot == 'y') { + try(source(run_script_command4)) + try(source(run_script_command10)) + } + if (scatter_percentiles_plot == 'y') { try(source(run_script_command5)) } + if (scatter_skill_plot == 'y') { try(source(run_script_command6)) } + if ((scatter_mtom_plot == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command7)) } + } + } + for (m in 1:length(batch_query)) { + species_list <- c("SO4","NO3","NH4","TC","PM_TOT") + 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("CSN") + network_label <- c("CSN") + pid <- network_label + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (scatter_plot == 'y') { + try(source(run_script_command1)) + try(source(run_script_command9)) + } + if (scatter_single_plot == 'y') { try(source(run_script_command2)) } + if (scatter_density_plot == 'y') { + try(source(run_script_command3)) + try(source(run_script_command11)) + } + if (scatter_bins_plot == 'y') { + try(source(run_script_command4)) + try(source(run_script_command10)) + } + if (scatter_percentiles_plot == 'y') { try(source(run_script_command5)) } + if (scatter_skill_plot == 'y') { try(source(run_script_command6)) } + if ((scatter_mtom_plot == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command7)) } + } + } + for (m in 1:length(batch_query)) { + species_list <- c("SO4","NO3","TNO3","NH4","SO2") + 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("CASTNET") + network_label <- c("CASTNET") + pid <- network_label + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (scatter_plot == 'y') { + try(source(run_script_command1)) + try(source(run_script_command9)) + } + if (scatter_single_plot == 'y') { try(source(run_script_command2)) } + if (scatter_density_plot == 'y') { + try(source(run_script_command3)) + try(source(run_script_command11)) + } + if (scatter_bins_plot == 'y') { + try(source(run_script_command4)) + try(source(run_script_command10)) + } + if (scatter_percentiles_plot == 'y') { try(source(run_script_command5)) } + if (scatter_skill_plot == 'y') { try(source(run_script_command6)) } + if ((scatter_mtom_plot == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command7)) } + } + } + for (m in 1:length(batch_query)) { + species_list <- c("PM_TOT") + 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) + system(mkdir_command) + dates <- batch_names[m] + network_names <- c("AQS_Daily") + network_label <- c("AQS_Daily") + pid <- network_label + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (scatter_plot == 'y') { + try(source(run_script_command1)) + try(source(run_script_command9)) + } + if (scatter_single_plot == 'y') { try(source(run_script_command2)) } + if (scatter_density_plot == 'y') { + try(source(run_script_command3)) + try(source(run_script_command11)) + } + if (scatter_bins_plot == 'y') { + try(source(run_script_command4)) + try(source(run_script_command10)) + } + if (scatter_percentiles_plot == 'y') { try(source(run_script_command5)) } + if (scatter_skill_plot == 'y') { try(source(run_script_command6)) } + if ((scatter_mtom_plot == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command7)) } + } + } + for (m in 1:length(batch_query)) { + species_list <- c("PM_TOT") + 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") + network_label <- c("AQS_Hourly") + pid <- network_label + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (scatter_plot == 'y') { + try(source(run_script_command1)) + try(source(run_script_command9)) + } + if (scatter_single_plot == 'y') { try(source(run_script_command2)) } + if (scatter_density_plot == 'y') { + try(source(run_script_command3)) + try(source(run_script_command11)) + } + if (scatter_bins_plot == 'y') { + try(source(run_script_command4)) + try(source(run_script_command10)) + } + if (scatter_percentiles_plot == 'y') { try(source(run_script_command5)) } + if (scatter_skill_plot == 'y') { try(source(run_script_command6)) } + if ((scatter_mtom_plot == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command7)) } + } + } + } ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL TO OBS SCATTERPLOT ############################## ### AMET CODE: R_Scatterplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species for multiple networks on a single ### plot. Additionally, summary statistics are also included on the plot. The script will ### also allow a second run to plotted on top of the first run. Output format is png, pdf ### or both. ### ### Last Updated 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_Scatterplot_Multisim_plotly.R': No such file or directory ########################### SINGLE NETWORK SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_single.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single network that includes more statistics than the multiple network scatter plot. ### Although limited to a single network, up to two simulations may be included. Output ### format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ######################### ### AMET CODE: R_Scatterplot_density.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### scatterplot. This script will plot a single species from up to three networks on a ### single plot. Additionally, summary statistics are also included on the plot. The ### script will also allow a second run to plotted on top of the first run. Output for ### is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ################################################################################## Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ####################### MODEL TO OBS DENSITY SCATTERPLOT ########################## ### AMET CODE: R_Scatterplot_density_ggplot.r ### ### This script is part of the AMET-AQ system. This script creates a single model-to-obs ### density scatterplot, where higher density areas of the plot are shaded differently ### from lower density areas. This script will plot a single species for a single network ### using the R ggplot package. ### ### Last Updated by Wyat Appel: June, 2019 #################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened #################### MODEL TO OBS BINNED SCATTERPLOT ###################### ### AMET CODE: R_Scatterplot_bins.r ### ### This script is part of the AMET-AQ system. This script creates a binned ### bias/RMSE plot. This script will create a binned bias/RMSE plot for a single ### species from a single network but for multiple simulations. Output format ### in png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ########################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened Error in library(plotly) : there is no package called ‘plotly’ ########################## MODEL PERCENTILE SCATTER PLOT ############################# ### AMET CODE: R_Scatterplot_percentile.R ### ### This script is part of the AMET-AQ system. This script creates a scatter plot for ### a single species, network and simulation that is color coded by the observation ### percentile, split by 5, 25, 50, 75 and 95. Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ###################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened ############################ MODEL SKILL SCATTER PLOT ############################## ### AMET CODE: R_Scatterplot_skill.R ### ### This script is part of the AMET-AQ system. This script creates a unique forecast ### skill scatter plot. The scipt was orinally designed to work specifically with 8-hr ### max ozone, as all the lines and stats were based on the criteria for 8-hr ozone. ### However, users can now specifiy the criteria line, and the forecast statistics are ### calculated based on the user specified criteria. So, it may be used with any species. ### Output format is png, pdf or both. ### ### Last Updated by Wyat Appel: June, 2019 ##################################################################################### Error in .local(drv, ...) : Cannot allocate a new connection: 16 connections already opened > if (dep_analysis == 'y') { + max_limit <- dep_max_limit + averaging <- deposition_averaging + species_list <- c("SO4_dep","NO3_dep","NH4_dep","Precip") + 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("NADP") + network_label <- c("NADP") + pid <- network_label + query <- paste(query_string,"and (",batch_query[m],")",sep=" ") + system(mkdir_command) + if (scatter_plot == 'y') { + try(source(run_script_command1)) + try(source(run_script_command9)) + } + if (scatter_single_plot == 'y') { try(source(run_script_command2)) } + if (scatter_density_plot == 'y') { + try(source(run_script_command3)) + try(source(run_script_command11)) + } + if (scatter_bins_plot == 'y') { + try(source(run_script_command4)) + try(source(run_script_command10)) + } + if (scatter_percentiles_plot == 'y') { try(source(run_script_command5)) } + if (scatter_skill_plot == 'y') { try(source(run_script_command6)) } + if ((scatter_mtom_plot == 'y') && (exists("run_name2")) && (nchar(run_name2) > 0)) { try(source(run_script_command7)) } + } + } + } Error: object 'dep_max_limit' not found Execution halted