File:Mis3 france warm kriged tjuly 1.svg

From Wikimedia Commons, the free media repository
Jump to navigation Jump to search

Original file(SVG file, nominally 1,650 × 1,275 pixels, file size: 1.02 MB)

Captions

Captions

July temperature on MIS3 warm stage in france.

Summary

[edit]
Description
English: Temperature of July in France during warm instage on MIS3.
Date
Source Own work
Author Merikanto
SVG development
InfoField
 
The SVG code is valid.
 
This vector image was created with Adobe Illustrator by Merikanto.

Tjuly data is from Davies& Gollop 2003

https://www.researchgate.net/publication/236109423_The_Human_Presence_in_Europe_during_the_Last_Glacial_Period_II_Climate_Tolerance_and_Climate_Preferences_of_Mid-and_Late_Glacial_Hominids

Davies, William & Gollop, Piers. (2003). The Human Presence in Europe during the Last Glacial Period II: Climate Tolerance and Climate Preferences of Mid-and Late Glacial Hominids.

Original Tjuly data is from Pollard et al OIS3 project.

https://www.researchgate.net/figure/Simulated-patterns-of-late-Quaternary-mean-temperatures-in-Europe-for-summer-JJA_fig2_236109423

Barron, E., & Pollard, D. (2002). High-resolution climate simulations of oxygen isotope stage 3 in Europe. Quaternary Research, 58(3), 296-309. https://doi.org/10.1006/qres.2002.2374

Some points data is intuitively manually selected from plottings.

Kriging with etopo1 dataset and r. Visualization with NASA Panoply. Conversion to SVG on net.

"ID","lon","lat","TS","Elev" "Calais",1.85,50.95,10.0,0 "Marse",5.0,41.0,25.0,0.0 "Marse2",4.0,43.0,23.0,0.0 "P5",0.0,50.0,10.0,0.0 "Up1",15.0,47.5,12,1000.0 "Up2",15.0,46.5,8,2000.0 "Mare1",-10.0,46.5,10,00.0 "Aqvit1",0.5,44.5,17,200 "LesEyzies",1.018,44.936,17.0,75 "Loire",0,47,17.0,250

library(raster) require(rgdal) library(sp) library(gstat) library(automap) library(RColorBrewer) library(viridis)

elevname0<-"./eurodata/area_glacial.nc"

outname1="base_tk7.nc"

elev1<-raster(elevname0)

names(elev1)<-"Elev" print(names(elev1)) print(elev1)

pointable1 <- read.csv("./mis3warm4.txt")

projstr1<-"+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"

print (projstr1)

points1 <- SpatialPointsDataFrame(coords = pointable1[,c("lon","lat")],

                                    data        = pointable1,
                                    proj4string = CRS(projstr1))
                                    

names(points1)<-c("ID","lon","lat","TS","Elev")

knitr::kable(head(points1, n=10))

print("SP")

points1$extracted <- extract(x=elev1, y=points1)

newgrid <- as(elev1, "SpatialGridDataFrame") names(newgrid) <- "Elev"

print("Variogram ...")

tsvar1 <- variogram(TS~Elev, points1) print(".") sph.model <- vgm(psill=15, model="Exp", range=500, nugget=20) print(".") sph.fit <- fit.variogram(object = tsvar1, model = sph.model)

print("Kriging ...")

rk1 <- krige(TS~Elev, points1, newgrid, model = sph.fit)

rout0 <- as(rk1[1], "RasterLayer")

plot(rout0)

plot(rout0, col=viridis(100))

crs(rout0) <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0" writeRaster(rout0, filename=outname1, overwrite=TRUE, format="CDF", varname="T_July", varunit="Celsius", longname="Temperature of July", xname="lon", yname="lat")

Licensing

[edit]
I, the copyright holder of this work, hereby publish it under the following license:
w:en:Creative Commons
attribution share alike
This file is licensed under the Creative Commons Attribution-Share Alike 4.0 International license.
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/TimeThumbnailDimensionsUserComment
current08:52, 30 September 2019Thumbnail for version as of 08:52, 30 September 20191,650 × 1,275 (1.02 MB)Merikanto (talk | contribs)User created page with UploadWizard

There are no pages that use this file.

File usage on other wikis

The following other wikis use this file:

Metadata