File:Pueblo area precipitation 1156 ad 1 1 1 1.png
Original file (1,600 × 1,600 pixels, file size: 1.42 MB, MIME type: image/png)
Captions
Summary
[edit]DescriptionPueblo area precipitation 1156 ad 1 1 1 1.png |
English: Pueblo area precipitation 1156 AD |
Date | |
Source | Own work |
Author | Merikanto |
Source of data is NOAA paleoclimate study "SW USA 2000 Year Growing Degree Days and Precipitation Reconstructions"
SW USA 2000 Year Growing Degree Days and Precipitation Reconstructions
Bocinsky, R.K.; Rush, J.; Kintigh, K.W.; Kohler, T.A.
Citation Information:
R. Kyle Bocinsky, Johnathan Rush, Keith W. Kintigh and Timothy A. Kohler. 2016. Exploration and exploitation in the macrohistory of the pre-Hispanic Pueblo Southwest. Science Advances, 2(4), e1501532. doi: 10.1126/sciadv.1501532
NOAA Study Page:
https://www.ncei.noaa.gov/access/paleo-search/study/19783
JSON Metadata:
https://www.ncei.noaa.gov/access/paleo-search/study/search.json?xmlId=17591
DIF Metadata:
http://www1.ncdc.noaa.gov/pub/data/metadata/published/paleo/dif/xml/noaa-recon-19783.xml
ISO Metadata:
http://www1.ncdc.noaa.gov/pub/data/metadata/published/paleo/iso/xml/noaa-recon-19783.xml
DOI:
https://doi.org/10.25921/8ctk-8v26
"R" code
- 3
-
- southwest usa 2000 years data temperature, precipitation curve
- "R"
-
- 26.01.2024 v 0000.0003
-
library(raster)
library(rgdal)
library(gdalUtils)
library(ncdf4)
library(terra)
library(pals)
create_ppt_map<-function(ext1)
{
files1<- list.files(path='./indata1',pattern=glob2rx("*PPT.nc4") , full.names = TRUE )
#print("Mosa ...")
print(files1)
#quit("yes")
len1<-length(files1)
print(len1)
r1<-raster(files1[1])
for (n in 2:len1)
{
print(files1[n])
}
template <- raster(ext1)
projection(template) <- '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs'
writeRaster(template, file="./indata1/PPT.tif", format="GTiff", overwrite=TRUE)
mosaic_rasters(gdalfile=files1,dst_dataset="./indata1/PPT.tif",of="GTiff")
gdalinfo("./indata1/PPT.tif")
#r1<-raster("./indata13/srtm.tif")
#plot(r1)
}
download_teering_reconstruction_slice_1<-function(palon1, palat1)
{
indir1<-"./indata1/"
dir.create(indir1)
slat1=as.character(as.integer(palat1))
slon1=as.character(as.integer(abs(palon1))+1)
nametrunc1=paste0(slon1, "W", slat1, "N")
iname1<-paste0(nametrunc1, "_GDD.nc4")
iname2<-paste0(nametrunc1, "_PPT.nc4")
inamedir1=paste0(indir1,iname1)
inamedir2=paste0(indir1,iname2)
print(slat1)
print(slon1)
print(iname1)
print(iname2)
urlbase1<-"https://www.ncei.noaa.gov/pub/data/paleo/treering/reconstructions/northamerica/usa/bocinsky2016/"
url1<-paste0(urlbase1, iname1)
url2<-paste0(urlbase1, iname2)
if(download_data==1)
{
try(download.file(url = url1,destfile = inamedir1), silent=FALSE)
try(download.file(url = url2,destfile = inamedir2), silent=FALSE)
}
}
-
download_data=0
make_mosaic=0
- ext1<-extent(c(-110,-105, 35, 39))
ext1<-extent(c(-111,-104, 34, 39))
- mesa verde
- sitename1<-"Mesa Verde"
- plat1<- 37.183889
- plon1<- -108.488611
- chaco
- sitename1="Chaco Canyon"
- plat1=36.058333
- plon1 =-107.958889
- casa grande
- sitename1="Casa Grande"
- plat1 = 32.997005
- plon1 = -111.532069
if(download_data==1)
{
palat1=33 ## 40
for (iy in 1:7)
{
palon1=-112 ## -100
for (ix in 1:8)
{
download_teering_reconstruction_slice_1(palon1, palat1)
palon1=palon1+1
}
palat1=palat1+1.0
}
}
if(make_mosaic==1)
{
create_ppt_map(ext1)
}
iname1<-"./indata1/PPT.tif"
- r1<-raster(iname1)
- r1<-rast(iname1)
rs1<-stack(iname1)
- str(rs1)
- quit("yes")
year1<-1156
- year1<-1050
- year1<-1
r0<-rs1year1
print(year1)
r1<-crop(r0, ext1)
- r1[r1>1000]=1000
r1[r1>1000]=1000
colo1="#001f00"
- plot ....
png(
"pueblo_area_precipitation_1156_ad_1_1_1_1.png",
width = 1600, height = 1600,
units = "px",
pointsize = 36, bg = "white"
#, res = NA,
#restoreConsole = TRUE
)
plot (r1, col=(brewer.blues(256)),
main="Pueblo area year precipitation in mm at 1156 AD",
xlab="lon", ylab="lat", cex=1
)
contour(r1, add=T, lwd=0.5)
points( c(-107.958889),c( 36.058333) , col = colo1, cex = 1.5, lwd=5)
text ( c(-107.958889),c( 36.058333+0.15) , "Chaco Canyon" , cex=2,col=colo1 )
points( c(-108.488611),( 37.183889), col = colo1, cex = 1.5, lwd=5)
text ( c(-108.488611),( 37.183889+0.15), "Mesa Verde" , cex=2, col=colo1 )
dev.off()
- str(r1)
quit("yes")
Licensing
[edit]- You are free:
- to share – to copy, distribute and transmit the work
- to remix – to adapt the work
- Under the following conditions:
- attribution – You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
- share alike – If you remix, transform, or build upon the material, you must distribute your contributions under the same or compatible license as the original.
File history
Click on a date/time to view the file as it appeared at that time.
Date/Time | Thumbnail | Dimensions | User | Comment | |
---|---|---|---|---|---|
current | 10:29, 26 January 2024 | 1,600 × 1,600 (1.42 MB) | Merikanto (talk | contribs) | Uploaded own work with UploadWizard |
You cannot overwrite this file.
File usage on Commons
There are no pages that use this file.