Jump to content

File:Beringia 20000bp duration of snow cover days 1.png

Page contents not supported in other languages.
This is a file from the Wikimedia Commons
From Wikipedia, the free encyclopedia

Original file (3,021 × 1,661 pixels, file size: 437 KB, MIME type: image/png)

Summary

Description
English: Duration of snow cover in days, East Beringia, 20000 years ago. Chelsa Trace 21ka variable bio/scd 200.
Date
Source Own work
Author Merikanto
Camera location65° 00′ 00″ N, 150° 00′ 00″ W  Heading=1° Kartographer map based on OpenStreetMap.View this and other nearby images on: OpenStreetMapinfo

https://chelsa-climate.org/chelsa-trace21k/

https://envicloud.wsl.ch/#/?prefix=chelsa%2Fchelsa_V1%2Fchelsa_trace%2F

directory orog

Tex dem file

https://os.zhdk.cloud.switch.ch/envicloud/chelsa/chelsa_V1/chelsa_trace/orog/CHELSA_TraCE21k_dem_-80_V1.0.tif

https://os.zhdk.cloud.switch.ch/envicloud/chelsa/chelsa_V1/chelsa_trace/orog/CHELSA_TraCE21k_gle_-80_V1.0.tif

https://os.zhdk.cloud.switch.ch/envicloud/chelsa/chelsa_V1/chelsa_trace/orog/CHELSA_TraCE21k_glz_-80_V1.0.tif

Trace21ka selected data ownloading "R" script is under this link

File:Beringia_8000bp_1.png

rem rem chelsa trace21ka orography and glaciers of beringia rem rem CHELSA_TraCE21k_dem_-80_V1.0.tif rem rem uses gmt 6.x, gdal, R rem rem v 0000.0009 rem 11.10.2021 rem

set taimi=200


set kaption=-B+glightblue+t"Beringia "%taimi%"00 years ago \n Duration of snow cover in days" rem set kaption=-B+glightblue+t"Beringia "%taimi%"00 years ago \n SWE"

rem set kaption=-B+glightblue+t"North Siberia "%taimi%"00 years ago \n Duration of snow cover, days" rem set kaption=-B+glightblue+t"North Siberia "%taimi%"00 years ago \n Snow water content SWE"


rem whole beringia rem set rect1=-R-220.0/-120.0/40.0/80.0 rem eastern beringia set rect1=-R-180.0/-130.0/50.0/75.0 rem whole beringia rem set rect2=-R140/260/40/80

rem east asia rem set rect1=-R80.0/180.0/40.0/80.0


set size1=-JQ24 set annot1=-B10m10m

set namebase1=CHELSA_TraCE21k_ set namebase2=_V1.0.tif

set params1=%size1% %rect1%

set gridskale=3600

set inputdem=indata\%namebase1%dem_-%taimi%%namebase2% set inputgle=indata\%namebase1%gle_-%taimi%%namebase2% set inputglz=indata\%namebase1%glz_-%taimi%%namebase2% set inputtmin7=indata\%namebase1%tasmin_7_-%taimi%%namebase2% set inputtmax7=indata\%namebase1%tasmax_7_-%taimi%%namebase2% set inputtmin1=indata\%namebase1%tasmin_1_-%taimi%%namebase2% set inputtmax1=indata\%namebase1%tasmax_1_-%taimi%%namebase2% set inputbio1=indata\%namebase1%bio01_-%taimi%%namebase2% set inputbio12=indata\%namebase1%bio12_-%taimi%%namebase2% set inputswe=indata\%namebase1%swe_-%taimi%%namebase2% set inputscd=indata\%namebase1%scd_-%taimi%%namebase2%

rem upscale grid, because in 180 deg is gap

rem if you comment this, you do not (pre)process rem goto jump1

rem skale grids only rem goto jump2

rscript dl1.r %taimi%

del dem.nc del gle.nc del glz.nc del sign.nc del gsign.nc

del t7min.nc del t7max.nc del t7.nc del t1min.nc del t1max.nc del t1.nc

del bio1.nc del bio12.nc del swe.nc del scd.nc


jump2

gdalwarp -of netcdf -ts %gridskale% 0 %inputdem% dem.nc gdalwarp -of netcdf -ts %gridskale% 0 %inputgle% gle.nc gdalwarp -of netcdf -ts %gridskale% 0 %inputglz% glz.nc gdalwarp -of netcdf -ts %gridskale% 0 %inputtmin7% t7min.nc gdalwarp -of netcdf -ts %gridskale% 0 %inputtmax7% t7max.nc gdalwarp -of netcdf -ts %gridskale% 0 %inputtmin1% t1min.nc gdalwarp -of netcdf -ts %gridskale% 0 %inputtmax1% t1max.nc

gdalwarp -of netcdf -ts %gridskale% 0 %inputbio1% bio1.nc gdalwarp -of netcdf -ts %gridskale% 0 %inputbio12% bio12.nc gdalwarp -of netcdf -ts %gridskale% 0 %inputswe% swe.nc gdalwarp -of netcdf -ts %gridskale% 0 %inputscd% scd.nc

grdmath dem.nc SIGN = sign.nc grdmath glz.nc SIGN = gsign.ncgrd grdmath t7min.nc t7max.nc ADD 20 DIV 273.15 SUB = t7.nc grdmath t1min.nc t1max.nc ADD 20 DIV 273.15 SUB = t1.nc

jump1

set dem1=./dem.nc set gle1=./gle.nc set glz1=./glz.nc


gmt makecpt -Chot > my_sea.cpt gmt makecpt -Cetopo1 > my_topo.cpt gmt makecpt -Cgebco > my_glacier.cpt gmt makecpt -Clightblue,darkgreen,darkgreen -T-%gridskale%/%gridskale%/%gridskale% > colors1.cpt gmt makecpt -Crainbow -T-30/30/10 -Z > temp.cpt

gmt begin out png

psbasemap %params1% %kaption% -Ba10

rem ncimage t7.nc %params1% -Ctemp.cpt

grdimage %dem1% %params1% -Q -Ccolors1.cpt grdimage %glz1% %params1% -Q -Cmy_glacier.cpt

gmt grdcontour sign.nc %params1% -W0.15p,black -L-1/1 -C1 gmt grdcontour gsign.nc %params1% -W0.0075p,220/220/255 -L0/2 -C1

rem gmt grdcontour t7.nc %params1% -W0.0075p,blue -L-25/25 -C5 -A5 rem gmt grdcontour t7.nc %params1% -W0.0075p,red -L10/10 -C1 -A1 rem gmt grdcontour t1.nc %params1% -W0.0075p,blue -L-3/-3 -C1 -A1 rem gmt grdcontour t1.nc %params1% -W0.0075p,blue -L-38/-38 -C1 -A1 rem gmt grdcontour bio1.nc %params1% -W0.0075p,lightblue -L-10/-10 -C1 -A1 rem gmt grdcontour bio1.nc %params1% -W0.0075p,lightblue -L-5/-5 -C1 -A1

rem snow water content rem gmt grdcontour swe.nc %params1% -W0.01p,lightblue -L0/200 -C10 -A- rem gmt grdcontour swe.nc %params1% -W0.03p,blue -L0/200 -C20 -A20


rem snow water content gmt grdcontour scd.nc %params1% -W0.01p,lightblue -L0/300 -C20 -A- gmt grdcontour scd.nc %params1% -W0.03p,blue -L0/300 -C20 -A40

rem grdimage %dem1% %params1% -Q -Ccolors1.cpt rem grdimage %glz1% %params1% -Q -Cmy_glacier.cpt

rem gmt grdcontour sign.nc %params1% -W0.15p,black -L-1/1 -C1 rem gmt grdcontour gsign.nc %params1% -W0.0075p,220/220/255 -L0/2 -C1

gmt end


Downloader



    1. download some selected chelsa trace21ka files from chelsa repository
    2. R 4.03 script
    3. 11.10.2021
    4. 0000.0006
    1. set timeout greater than 60, so we can download big files

options(timeout=2400)


download_trace21ka<-function(taimi, variable1, variable2) { address1="https://os.zhdk.cloud.switch.ch/envicloud/chelsa" address2="/chelsa_V1/chelsa_trace/" namebase1="CHELSA_TraCE21k_" namebase2="_V1.0.tif" url1=paste0(address1,address2,variable2,"/") dolly1=paste0(url1,namebase1,variable1,"_-",taimi,namebase2) desti1=paste0("indata/",namebase1,variable1,"_-",taimi,namebase2) #dolly1="https://os.zhdk.cloud.switch.ch/envicloud/chelsa/chelsa_V1/chelsa_trace/orog/CHELSA_TraCE21k_dem_-100_V1.0.tif" #desti1="out.tif" download.file(dolly1,desti1,method = "libcurl", mode = "wb" ) }


download_selected_trace29ka<-function(taimi) { variable2="orog" variable1="dem" download_trace21ka(taimi, variable1, variable2) variable1="glz" download_trace21ka(taimi, variable1, variable2) variable1="gle" download_trace21ka(taimi, variable1, variable2)

variable1="tasmax_7" variable2="tasmax" download_trace21ka(taimi, variable1, variable2)

variable1="tasmin_7" variable2="tasmin" download_trace21ka(taimi, variable1, variable2)

variable1="tasmax_1" variable2="tasmax" download_trace21ka(taimi, variable1, variable2)

variable1="tasmin_1" variable2="tasmin" download_trace21ka(taimi, variable1, variable2)

variable1="pr_7" variable2="pr" download_trace21ka(taimi, variable1, variable2)

variable1="pr_1" variable2="pr" download_trace21ka(taimi, variable1, variable2)

variable1="bio01" variable2="bio" download_trace21ka(taimi, variable1, variable2)

variable1="bio12" variable2="bio" download_trace21ka(taimi, variable1, variable2)

variable1="swe" variable2="bio" download_trace21ka(taimi, variable1, variable2)


variable1="scd" variable2="bio" download_trace21ka(taimi, variable1, variable2)

print("Trace21ka stuff downloaded.")

}



          1. unlink("./indata", recursive=TRUE)

dir.create("./indata")

args = commandArgs(trailingOnly=TRUE)

  1. print (args[1])

taimi=args[1]

print(taimi)


download_selected_trace29ka(taimi)

print(".")

Licensing

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.

Captions

Duration of snow cover in days, East Beringia, 20000 years ago

Items portrayed in this file

depicts

image/png

24078f0b77167c16ea5a70021c6b72c34426ad85

447,618 byte

1,661 pixel

3,021 pixel

File history

Click on a date/time to view the file as it appeared at that time.

Date/TimeThumbnailDimensionsUserComment
current05:28, 11 October 2021Thumbnail for version as of 05:28, 11 October 20213,021 × 1,661 (437 KB)MerikantoUpdate
14:03, 10 October 2021Thumbnail for version as of 14:03, 10 October 20213,021 × 1,661 (624 KB)MerikantoUpdate
13:39, 10 October 2021Thumbnail for version as of 13:39, 10 October 20213,027 × 2,142 (665 KB)MerikantoUploaded own work with UploadWizard

The following page uses this file:

Global file usage

The following other wikis use this file:

Metadata