Skip to content

Commit

Permalink
Merge branch 'main' of https://github.com/Gulf-IEA/Caribbean-ESR into…
Browse files Browse the repository at this point in the history
… main

merging for unknown reason
  • Loading branch information
MandyKarnauskas-NOAA committed Feb 6, 2024
2 parents f444aec + 1f4a343 commit f664524
Show file tree
Hide file tree
Showing 40 changed files with 3,987 additions and 0 deletions.
28 changes: 28 additions & 0 deletions .Rhistory
Original file line number Diff line number Diff line change
@@ -0,0 +1,28 @@
rm(list = ls())
dat <- read.table("../indicator_data/surface_omega_series.txt", skip = 0, header = F)
head(dat)
tail(dat)
dat$dates <- paste0(month.abb[dat$V2], dat$V1)
datdata <- dat$dates
inddata <- data.frame(dat$V3)
labs <- c("Ocean acidification", "Surface aragonite saturation", "")
indnames <- data.frame(matrix(labs, nrow = 3, byrow = F))
s <- list(labels = indnames, indicators = inddata, datelist = datdata)
class(s) <-"indicatordata"
plotIndicatorTimeSeries(s)
library(plotTimeSeries)
library(spam)
rm(list = ls())
dat <- read.table("../indicator_data/surface_omega_series.txt", skip = 0, header = F)
head(dat)
tail(dat)
dat$dates <- paste0(month.abb[dat$V2], dat$V1)
datdata <- dat$dates
inddata <- data.frame(dat$V3)
labs <- c("Ocean acidification", "Surface aragonite saturation", "")
indnames <- data.frame(matrix(labs, nrow = 3, byrow = F))
s <- list(labels = indnames, indicators = inddata, datelist = datdata)
class(s) <-"indicatordata"
plotIndicatorTimeSeries(s)
inddata <- s
save(inddata, file = "OA_test.RData")
3 changes: 3 additions & 0 deletions .gitignore
Original file line number Diff line number Diff line change
@@ -1 +1,4 @@

.Rproj.user

/.quarto/
16 changes: 16 additions & 0 deletions .quarto/_freeze/OA/execute-results/html.json
Original file line number Diff line number Diff line change
@@ -0,0 +1,16 @@
{
"hash": "65457366e3373b6766a45354a36143ce",
"result": {
"markdown": "## Ocean acidification indicator\n\nSent by F. Gomez on 02/19/2023\n\nNotes from Fabian:\n\nderived surface omega series from the MOM-Topaz hindcast 0.10deg resolution model \\*txt files have this structure: column1: year, column2=month, and column3=carbon system variable The mean magnitude of the simulated omega trends, around 9e-3 year-1, is consistent with observed trends in the Subtropical North Atlantic\n\nQuestion: What is causing the acceleration after 2008?\\\nAnswer: made a Taylor decomposition to figure out what was driving that acceleration in the simulated ??Ar anomaly. Interannual ??Ar changes were mostly driven by the balance between dic and alkalinity. There was a positive trend from the 1980s until mid 2000s, which contributed to moderate the ??Ar decline. But that trend vanished in the last 15 years or so, which probably has to do with this accelerated decline in ??Ar in the last decade.\n\n\n::: {.cell}\n\n```{.r .cell-code}\nlibrary(plotTimeSeries)\nlibrary(spam)\n```\n:::\n\n::: {.cell}\n\n```{.r .cell-code}\nrm(list = ls())\n\ndat <- read.table(\"indicator_data/surface_omega_series.txt\", skip = 0, header = F)\nhead(dat)\ntail(dat)\ndat$dates <- paste0(month.abb[dat$V2], dat$V1)\n```\n:::\n\n\nFormat indicator object\n\n\n::: {.cell}\n\n```{.r .cell-code}\ndatdata <- dat$dates \ninddata <- data.frame(dat$V3) \nlabs <- c(\"Ocean acidification\", \"Surface aragonite saturation\", \"\") \nindnames <- data.frame(matrix(labs, nrow = 3, byrow = F)) \ns <- list(labels = indnames, indicators = inddata, datelist = datdata)\nclass(s) <-\"indicatordata\"\n```\n:::\n\n\nSave and plot\n\n\n::: {.cell}\n\n```{.r .cell-code}\nplotIndicatorTimeSeries(s)\n```\n\n::: {.cell-output-display}\n![](OA_files/figure-html/unnamed-chunk-4-1.png){width=672}\n:::\n\n```{.r .cell-code}\ninddata <- s \nsave(inddata, file = \"OA_test.RData\")\n```\n:::\n",
"supporting": [
"OA_files\\figure-html"
],
"filters": [
"rmarkdown/pagebreak.lua"
],
"includes": {},
"engineDependencies": {},
"preserve": {},
"postProcess": true
}
}
20 changes: 20 additions & 0 deletions .quarto/_freeze/OA/execute-results/tex.json
Original file line number Diff line number Diff line change
@@ -0,0 +1,20 @@
{
"hash": "65457366e3373b6766a45354a36143ce",
"result": {
"markdown": "## Ocean acidification indicator\n\nSent by F. Gomez on 02/19/2023\n\nNotes from Fabian:\n\nderived surface omega series from the MOM-Topaz hindcast 0.10deg resolution model \\*txt files have this structure: column1: year, column2=month, and column3=carbon system variable The mean magnitude of the simulated omega trends, around 9e-3 year-1, is consistent with observed trends in the Subtropical North Atlantic\n\nQuestion: What is causing the acceleration after 2008?\\\nAnswer: made a Taylor decomposition to figure out what was driving that acceleration in the simulated ??Ar anomaly. Interannual ??Ar changes were mostly driven by the balance between dic and alkalinity. There was a positive trend from the 1980s until mid 2000s, which contributed to moderate the ??Ar decline. But that trend vanished in the last 15 years or so, which probably has to do with this accelerated decline in ??Ar in the last decade.\n\n\n\n::: {.cell}\n\n```{.r .cell-code}\nlibrary(plotTimeSeries)\nlibrary(spam)\n```\n:::\n\n::: {.cell}\n\n```{.r .cell-code}\nrm(list = ls())\n\ndat <- read.table(\"indicator_data/surface_omega_series.txt\", skip = 0, header = F)\nhead(dat)\ntail(dat)\ndat$dates <- paste0(month.abb[dat$V2], dat$V1)\n```\n:::\n\n\n\nFormat indicator object\n\n\n\n::: {.cell}\n\n```{.r .cell-code}\ndatdata <- dat$dates \ninddata <- data.frame(dat$V3) \nlabs <- c(\"Ocean acidification\", \"Surface aragonite saturation\", \"\") \nindnames <- data.frame(matrix(labs, nrow = 3, byrow = F)) \ns <- list(labels = indnames, indicators = inddata, datelist = datdata)\nclass(s) <-\"indicatordata\"\n```\n:::\n\n\n\nSave and plot\n\n\n\n::: {.cell}\n\n```{.r .cell-code}\nplotIndicatorTimeSeries(s)\n```\n\n::: {.cell-output-display}\n![](OA_files/figure-pdf/unnamed-chunk-4-1.pdf){fig-pos='H'}\n:::\n\n```{.r .cell-code}\ninddata <- s \nsave(inddata, file = \"OA_test.RData\")\n```\n:::\n",
"supporting": [
"OA_files\\figure-pdf"
],
"filters": [
"rmarkdown/pagebreak.lua"
],
"includes": {},
"engineDependencies": {
"knitr": [
"{\"type\":\"list\",\"attributes\":{},\"value\":[]}"
]
},
"preserve": null,
"postProcess": false
}
}
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file not shown.
14 changes: 14 additions & 0 deletions .quarto/_freeze/index/execute-results/html.json
Original file line number Diff line number Diff line change
@@ -0,0 +1,14 @@
{
"hash": "e689c5f9263022dd5ee70a83347bdcd8",
"result": {
"markdown": "## ESR Lab Reports\n\nESR lab reports are quarto books that compile the R scripts for each indicator as well as any other descriptions of how the ESR will be put together. \n\n## Caribbean ESR\n\nThis lab report is for the Caribbean ESR that was compiled from 2023-2024. This is the first ESR for the Caribbean region. Here we can add more text about why this ESR is being compiled.\n\n\n::: {.cell warnings='false'}\n\n```{.r .cell-code}\nlibrary(plotTimeSeries)\nlibrary(spam)\n```\n:::\n",
"supporting": [],
"filters": [
"rmarkdown/pagebreak.lua"
],
"includes": {},
"engineDependencies": {},
"preserve": {},
"postProcess": true
}
}
7 changes: 7 additions & 0 deletions .quarto/_freeze/site_libs/clipboard/clipboard.min.js

Some generated files are not rendered by default. Learn more about how customized files appear on GitHub.

1 change: 1 addition & 0 deletions .quarto/cites/index.json
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
{"Lab-report\\OA.qmd":[],"index.qmd":[],"OA.qmd":[]}
Loading

0 comments on commit f664524

Please sign in to comment.