Tiedosto:East sibiria at lgm 2.svg

Wikipediasta
Siirry navigaatioon Siirry hakuun

Alkuperäinen tiedosto (SVG-tiedosto; oletustarkkuus 1 053 × 744 kuvapistettä; tiedostokoko 930 KiB)

Yhteenveto

Kuvaus
English: East Siberia at Last Glacial Maximum (LGM). Almost all sites are above 9.5 ° C Tjuly isotherm. There was no sure habitation above 58 °N latitude. This means that tundra north from treeline was unhabitated.
Päiväys
Lähde Oma teos
Tekijä Merikanto
SVG kehittely
InfoField
 
Tämän SVG-tiedoston lähdekoodi on validoitu.
 
Tämä vektorigrafiikkatiedosto luotiin käyttäen apuna ohjelmaa unknown tool

LGM human occupation sites data is from Kuzmin&Keates 2016

Kuzmin, Yaroslav & G. Keates, Susan. (2016). Siberia and neighboring regions in the Last Glacial Maximum: did people occupy northern Eurasia at that time?. Archaeological and Anthropological Sciences.

doi:10.1007/s12520-016-0342-z

Image reproduction scripts in Ubuntu etc.

Init script

    1. get LGM climate data geotiff
  1. http://www.worldclim.org/paleo-climate1

mkdir data cd data

wget http://biogeo.ucdavis.edu/data/climate/cmip5/lgm/cclgmtn_2-5m.zip wget http://biogeo.ucdavis.edu/data/climate/cmip5/lgm/cclgmtx_2-5m.zip wget http://biogeo.ucdavis.edu/data/climate/cmip5/lgm/cclgmpr_2-5m.zip wget http://biogeo.ucdavis.edu/data/climate/cmip5/lgm/cclgmbi_2-5m.zip

unzip cclgmtn_2-5m.zip unzip cclgmtx_2-5m.zip unzip cclgmpr_2-5m.zip unzip cclgmbi_2-5m.zip

rm *.zip

cd ..

mkdir plot mkdir teksti

    1. init mapping utilzzz

sudo apt install gdal-bin python-gdal nco cdo ctioga2 libnetcdf libnetcdf-dev -y

    1. sumenna.bash

sudo apt install saga

  1. gcc -o ncblit1 ncblit1.c -lm -lnetcdf
  2. gcc -o ncprobe1 ncprobe1.c -lm -lnetcdf

Data downloading script, in ccsm dir

sorckanta1="cclgmtn" sorckanta2="cclgmtx" sorckanta3="cclgmpr" sorckanta4="cclgmtk" psorckanta1="cclgmbi"

cd data

for i in {1..12..1}

 do
  
    sorcnimi1=$sorckanta1$i.tif
    sorcnimi2=$sorckanta2$i.tif
    sorcnimi3=$sorckanta3$i.tif
    sorcnimi4=$sorckanta4$i.tif
    nesorcnimi1=$sorckanta1$i.nc
    nesorcnimi2=$sorckanta2$i.nc
    nesorcnimi3=$sorckanta3$i.nc
    nesorcnimi4=$sorckanta4$i.nc
   gdal_calc.py -A $sorcnimi1 -B $sorcnimi2 --outfile=$sorcnimi4 --calc="(A+B)/20"
  gdal_translate -ot Float32 -of netCDF $sorcnimi1 $nesorcnimi1
  gdal_translate -ot Float32 -of netCDF $sorcnimi2 $nesorcnimi2
  gdal_translate -ot Float32 -of netCDF $sorcnimi3 $nesorcnimi3
  gdal_translate -ot Float32 -of netCDF $sorcnimi4 $nesorcnimi4
done

for i in {1..15..1}

 do
    psorcnimi1=$psorckanta1$i.tif
   
    pesorcnimi1=$psorckanta1$i.nc
    
    gdal_calc.py -A $psorcnimi1 --outfile=temp.tif --calc="A"
  
    gdal_translate -ot Float32 -of netCDF temp.tif $pesorcnimi1
  
done

cd ..

Gaussian interpolation script, w/saga GIS

ingrid1="../ccsm/data/cclgmtk7.nc" ingrid2="../ccsm/data/cclgmtk1.nc" ingrid3="../ccsm/data/cclgmbi1.nc"

outgridbase1="tk7_gauss1" outgridbase2="tk1_gauss1" outgridbase3="ta_gauss1"

saga_cmd grid_filter 1 -INPUT $ingrid3 -RESULT $outgridbase3 -SIGMA 6 -RADIUS 6 gdal_translate -of GMT $outgridbase3".sdat" temp.nc

gdal_calc.py -A temp.nc --outfile="temp2.tif" --calc="(A/10)"

gdal_translate -of netCDF temp2.tif $outgridbase3".nc"

exit

saga_cmd grid_filter 1 -INPUT $ingrid1 -RESULT $outgridbase1 -SIGMA 6 -RADIUS 6 gdal_translate -of GMT $outgridbase1".sdat" temp.nc

gdal_calc.py -A temp.nc --outfile="temp2.tif" --calc="(A)"

gdal_translate -of netCDF temp2.tif $outgridbase1".nc"

saga_cmd grid_filter 1 -INPUT $ingrid2 -RESULT $outgridbase2 -SIGMA 6 -RADIUS 6 gdal_translate -of GMT $outgridbase2".sdat" temp.nc

gdal_calc.py -A temp.nc --outfile="temp2.tif" --calc="(A)"

gdal_translate -of netCDF temp2.tif $outgridbase2".nc" CertainLGM sites loydot.txt, charchoal

091.43   52.97  24 0 14 BL   Ui 1
087.95  55.90   24 0 14 BL   Shestakovo 
091.10  55.96   24 0 14 BL   Novoselovo
092.52  55.33   24 0 14 BL   Derbina 4 
093.23  56.05   24 0 14 BL   Gosudarev

103.50 52.90 24 0 14 BL Malta 103.42 52.95 24 0 14 BL Buret 100.33 58.33 24 0 14 BL Ust’-Kova 108.33 57.83 24 0 14 BL Alekseevsk 1 129.58 51.83 24 0 14 BL Ust’-Ulma 1 127.66 50.33 24 0 14 BL Khodulikha 142.40 46.77 24 0 14 BL Ogonki 5

Unvertain sites, bones 083.52 57.75 BL Mogochino 084.60 51.30 BL Anui 2, layer 3 090.95 55.01 BL Sabanikha 091.42 55.13 BL Kashtanka 1 092.05 55.20 BL Shlenka 108.50 50.18 BL Studenoe 2 112.50 57.50 BL Tesa 117.47 59.93 BL Khayrgas 133.75 63.17 BL Ikhine 2

ann takayrat.txt -15

jan tk1kayrat.txt -40 jul tk7kayrat.txt 9.5

ta.cpt

-30 100/200/100 10.0 100/200/100 ; arctic 10.0 255/127/50 20 255/127/50 ; temporal, subtropic. 20.0 255/100/100 30 255/100/100 ; subtropic. tropic

B blue F white N 127.5

Image gen

  1. siberia

lona=80 lonb=140 lata=40 latb=80

projektio=-JQ120/10.0i

  1. paletti="zones.cpt"

paletti="ta.cpt"

tk7="../ccsm/data/cclgmtk7.nc" tk1="../ccsm/data/cclgmtk1.nc" tka="../ccsm/data/cclgmtbi12.nc"

tk7s1="tk7_gauss1.nc" tk1s1="tk1_gauss1.nc" tas1="ta_gauss1.nc"

ta="tk7kayrat.txt" tb="tk1kayrat.txt" tc="takayrat.txt"

gmt makecpt -Crainbow -T-20/20/5 > t.cpt

  1. gmt makecpt -Crainbow -T9/10/5 > tk.cpt


gmt grdimage $tk7 -R$lona/$lonb/$lata/$latb $projektio -C$paletti -K > out.ps

gmt pscoast -K -O $projektio -R$lona/$lonb/$lata/$latb -B10g5/5g5 -S180/180/255 -W1p,200/200/255 \ -I1/1p,200/200/255 -I2/1p,200/200/255 -I3/0.5p,200/200/255 -I4/0.25p,200/200/255 \

-P -N1/1.5p,127/0/0,.- -Dc >>out.ps

gmt grdcontour $tk7s1 -R$lona/$lonb/$lata/$latb $projektio -W0.75,red@25 -C$ta -K -O >> out.ps

  1. gmt grdcontour $tk1s1 -R$lona/$lonb/$lata/$latb $projektio -W2,blue@25 -C$tb -K -O >> out.ps
  2. gmt grdcontour $tas1 -R$lona/$lonb/$lata/$latb $projektio -W3,darkgreen@25 -C$tc -K -O >> out.ps

cat loydot.txt | gmt psxy -R$lona/$lonb/$lata/$latb $projektio -O -K -G100/0/0 -Sc0.75 >> out.ps

cat loydote.txt | gmt psxy -R$lona/$lonb/$lata/$latb $projektio -O -K -G0/0/0 -Sc0.5 >> out.ps

gmt psxy -R$lona/$lonb/$lata/$latb $projektio -O -K -G255/0/0 -St0.75i << EOF >> out.ps 103.5 52.9 EOF

gmt pstext -R -J -O -K -F+f60,Helvetica-Bold,white=thin >> out.ps << END

100 65 Siberia

END

gmt pstext -R -J -O -K -F+f40,Times-Italic,blue >> out.ps << END

109 56 Baikal

END

gmt pstext -R -J -O -K -F+f75,Helvetica-Bold >> out.ps << END

120 75 LGM sites

END

gmt pslegend -R -J -O -K -DjTL+w2.5i+o0.4i -F+pthick+ithinner+gwhite --FONT_ANNOT_PRIMARY=17p,Helvetica-Bold<< EOF >> out.ps S 0.1i t0.75 0.08i red thinnest 0.3i Baikal site S 0.1i c0.5 0.08i 100/0/0 thinnest 0.3i Bone sites S 0.1i c0.5 0.08i black thinnest 0.3i Charcoal sites S 0.1i - 0.08i red thick,red 0.3i Tjuly 9.5 oC

  1. S 0.1i f0.1i+l+t 0.25i blue 0.25p 0.3i This is a fault

EOF

  1. gmt psscale -Dx14c/1.25c+w10c/0.5c+h -B+tcontinuous -Ct.cpt -By+loC -O -Ac >> out.ps

ps2pdf out.ps out.pdf pdf2svg out.pdf out.svg

  1. svgo out.svg

Lisenssi

Minä, tämän teoksen tekijänoikeuksien haltija, julkaisen täten tämän teoksen seuraavalla lisenssillä:
w:fi:Creative Commons
nimeäminen jaa samoin
Tämä tiedosto on lisensoitu Creative Commons Nimeä-JaaSamoin 4.0 Kansainvälinen -lisenssillä.
Voit:
  • jakaa – kopioida, levittää ja esittää teosta
  • remiksata – valmistaa muutettuja teoksia
Seuraavilla ehdoilla:
  • nimeäminen – Sinun on mainittava lähde asianmukaisesti, tarjottava linkki lisenssiin sekä merkittävä, mikäli olet tehnyt muutoksia. Voit tehdä yllä olevan millä tahansa kohtuullisella tavalla, mutta et siten, että annat ymmärtää lisenssinantajan suosittelevan sinua tai teoksen käyttöäsi.
  • jaa samoin – Jos muutat tai perustat tähän työhön, voit jakaa tuloksena syntyvää työtä vain tällä tai tämän kaltaisella lisenssillä.

Kuvatekstit

Lisää yhden rivin pituinen kuvaus tästä tiedostosta

Kohteet, joita tässä tiedostossa esitetään

esittää

Jotkut arvot ilman kohdetta Wikidata

7. joulukuu 2018

Tiedoston historia

Päiväystä napsauttamalla näet, millainen tiedosto oli kyseisellä hetkellä.

PäiväysPienoiskuvaKokoKäyttäjäKommentti
nykyinen7. joulukuuta 2018 kello 09.14Pienoiskuva 7. joulukuuta 2018 kello 09.14 tallennetusta versiosta1 053 × 744 (930 KiB)MerikantoUser created page with UploadWizard

Seuraava sivu käyttää tätä tiedostoa:

Metatieto