File:Beringia 20000bp duration of snow cover days 1.png
Original file (3,021 × 1,661 pixels, file size: 437 KB, MIME type: image/png)
|
This is a file from the Wikimedia Commons. The description on its description page there is shown below.
|
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 location |   | View this and other nearby images on: OpenStreetMap |  |
|---|
https://chelsa-climate.org/chelsa-trace21k/
https://envicloud.wsl.ch/#/?prefix=chelsa%2Fchelsa_V1%2Fchelsa_trace%2F
directory orog
Tex dem file
Trace21ka selected data ownloading "R" script is under this link
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
- download some selected chelsa trace21ka files from chelsa repository
- R 4.03 script
- 11.10.2021
- 0000.0006
- 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.")
}
- unlink("./indata", recursive=TRUE)
dir.create("./indata")
args = commandArgs(trailingOnly=TRUE)
- print (args[1])
taimi=args[1]
print(taimi)
download_selected_trace29ka(taimi)
print(".")
Licensing
- 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 | Dimensions | User | Comment | |
|---|---|---|---|---|
| current | 21:28, 10 October 2021 | 3,021 × 1,661 (437 KB) | Merikanto | Update |
File usage
The following page uses this file:
