Create a gist now

Instantly share code, notes, and snippets.

What would you like to do?
# Devide the reanalysis data
sea_temp <- filter(all_jan1_0.5, variable == "BRAN/temp")
air_temp <- filter(all_jan1_0.5, variable == "ERA/temp")
currents <- filter(all_jan1_0.5, variable == "BRAN/u" | variable == "BRAN/v") %>%
select(-date, -index) %>%
spread(key = variable, value = value) %>%
rename(u = "BRAN/u", v = "BRAN/v")
winds <- filter(all_jan1_0.5, variable == "ERA/u" | variable == "ERA/v") %>%
select(-date, -index) %>%
spread(key = variable, value = value) %>%
rename(u = "ERA/u", v = "ERA/v")
# Reduce wind/ current vectors
lon_sub <- seq(10, 40, by = 1)
lat_sub <- seq(-40, -15, by = 1)
# currents <- currents[(currents$lon %in% lon_sub & currents$lat %in% lat_sub),]
winds <- winds[(winds$lon %in% lon_sub & winds$lat %in% lat_sub),]
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment