############################ 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.6-0 (2020-12-14) 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

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.SO4_ob, d.SO4_mod ,d.POCode,s.state from  ARDEQ_JUL_TRIAL  as d, site_metadata as s  WHERE d.SO4_ob is not NULL and d.network='IMPROVE' and s.stat_id=d.stat_id and d.ob_dates BETWEEN 20160701 and 20160731 and d.ob_datee BETWEEN 20160701 and 20160731 and ob_hour between 0 and 23  ORDER BY stat_id,ob_dates,ob_hour"
Error in round(stats.df$MEAN_OBS, 1) : 
  non-numeric argument to mathematical function
In addition: There were 11 warnings (use warnings() to see them)
Execution halted