## NOTE resolution is fixed at SSE for the following
if (!exists("year.assessment")) {
year.assessment=lubridate::year(Sys.Date()) - 1 # choose one
year.assessment=lubridate::year(Sys.Date())
}
# ----------------------------------------------------------
# glue biological data sets together from various surveys and lookup environmental data where possible
p = aegis::aegis_parameters( DS="survey", yrs=1970:year.assessment )
survey.db( DS="set.init.redo", p=p )
survey.db( DS="cat.init.redo", p=p )
survey.db( DS="det.init.redo", p=p )
# the following requires the preformed 10.lookuptable.r to complete the data lookups of missing data
survey.db( DS="lengthweight.redo", p=p ) # # TODO:: parallelize me ... update the local tables (not necessary)
survey.db( DS="set.base.redo", p=p ) # adds temperature required for metabolism lookup in "det.redo"
survey.db( DS="det.redo", p=p ) # mass/length imputation and sanity checking
survey.db( DS="cat.redo", p=p ) # sanity checking and fixing mass estimates from det etc ...
survey.db( DS="set.redo", p=p ) # sanity checking and year filtering to 1999 - present
figure.bio.map.survey.locations(p=p) # see mpa/src/_Rfunctions/figure.trawl.density for more control
# ----------------------------------------------------------
# landings are like surveys -- but we do not want to interpolate this
# .. rather aggregate into meaningful areas, model patterns and map
# or modify as supplemntary data for distributional models
if (0) {
pl = aegis::aegis_parameters( DS="landings", yrs=1970:year.assessment ) # these are default years
landings.db( DS="rawdata", p=pl )
for ( vn in p$varstomodel) {
print(vn)
aegis_db ( DS="complete.redo", p=p )
aegis_db_map( p=p )
}
}
### end
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.