File:Ccsm trace21k tjuly 14559bp europe.svg
Original file (SVG file, nominally 1,650 × 1,275 pixels, file size: 1.72 MB)
Captions
Summary
[edit]DescriptionCcsm trace21k tjuly 14559bp europe.svg |
English: July temperature C at 14559 BP Late Glacial last interstadial in Europe. It is based on CCSM3 Trace 21k simulation and downscaling with current temparature map. |
|||
Date | ||||
Source | Own work | |||
Author | Merikanto | |||
SVG development InfoField |
|
Source of temperature data is NCAR Trace21k
http://www.cgd.ucar.edu/ccr/TraCE/
"TraCE-21ka: Simulation of Transient Climate Evolution over the last 21,000 years"
New image version: Obtained paltier 1ce 6 G map and cut it west europe size with "gdalwarp te". Used trace 21 k slices ca 1850 CE and 14995 CE. these called a and b. These was scalad to size of 10 min West europe map with "cdo ramapbil." glacial age c and current d. worldclim current climate was scaled to reso 10 min and europe size w/cdo remapbil to raster e. calculated difference of f=e - d with qgis raster calculator. or current high reso - current coarse reso. added g = c+f. then have downscaled raster.
Old image version.
Data is manually downloaded from ESG site extracted with bash and R script.
Old slice extraction script, Ubuntu.,Mint ......
variable=TS
taimi=3726
outi=ts_14585_raw.nc
filee=trace.10.14900-14351BP.cam2.h0.TS.0710101-0765012.nc
latvar=lat
lonvar=lon
datadir=/home/hirohito/ccsmdata/trace/
outdatadir=./data
cdo -selname,$variable,$latvar,$lonvar -seltimestep,$taimi $datadir$filee $outdatadir"/"$outi
Preprocessing script
want no data files in rundir
rm *.nc
- take snapshot from Trace21k TS file
- bash rainsnap.bash
- rain_infile=./data/trace_21k_7.nc
rain_infile=./data/ts_14585_raw.nc
orog_infile=./data/ice5g_v1.2_21.0k_10min.nc
- from etopo1 cropped and 2.5 min scaled orogeny
- orog_infile=/home/hirohito/grid/euro25b.nc
rain_inref=./data/curr_tmean7_10min.nc
ncks -O -C -vlong,lat,orog $orog_infile temp.nc
- ncks -O -C -vlon,lat,orog $orog_infile temp.nc
ncks -O -C -x -v bounds_lat,bounds_long temp.nc temp2.nc
ncap2 -O -s '{where( orog < 0 ) orog = 1; }' temp2.nc orog_in_10.nc
- ncks -O -C -v lon,lat,TS -x lat_bnds,long_bnds $rain_infile temp.nc
ncks -O -C -vlon,lat,TS $rain_infile temp.nc
ncks -O -C -x -v lat_bnds,long_bnds temp.nc temp2.nc
ncap2 -O -s '{TS=TS-273.13; }' temp2.nc rains.nc
cp rains.nc rain_coarse.nc
cdo remapdis,orog_in_10.nc,16 rains.nc rain_in_10.nc
ncks -O -C -vlon,lat,Band1 $rain_inref temp.nc
ncks -O -C -x -v lat_bnds,long_bnds temp.nc temp2.nc
ncrename -O -vBand1,TS temp2.nc temp3.nc
cdo remapdis,orog_in_10.nc,16 temp3.nc rain_ref_10.nc
- cdo latlonbox, -176.25 189.9
cp rain_coarse.nc c1.nc
cp orog_in_10.nc p1.nc
cp rain_ref_10.nc p2.nc
- cdo remapnn,lon=-68_lat=-18 ifile ofile
cdo sellonlatbox,-176.25,180.0,-90,90 c1.nc c2.nc
cdo sellonlatbox,-176.25,180.0,-90,90 p1.nc p1a.nc
cdo sellonlatbox,-176.25,180.0,-90,90 p2.nc p2a.nc
- ro dissever analysis
Rscript distek1.r
Dissever R 3.2 script
distek1.r
- basic script downscaling rainfall
- or temperature with R dissever geospatial package
- R version 3.4.4 (2018-03-15)
- you must install libs and data
- coarse_raster is coarse temperature or precipitation file, in netcdf,
- usually output from simulation
- p1 is orography/topography file, raster in netcdf format
- p2 training file for rainfall ot remparature
- fine rsolution is usually current temperature
- resolution of p1, p2 files must be same
- generic format for input file is netcdf w/vars lat,lon,z(lat, lon)
library(ncdf4)
library(dissever)
library(raster)
library(viridis)
library(dplyr)
library(lattice)
library(RColorBrewer)
coarse_raster <- raster("c2.nc")
p1<- raster("p1a.nc")
p2<- raster("p2a.nc")
raster_stack <- stack(p1,p2)
min_iter <- 5 # Minimum number of iterations
max_iter <- 10 # Maximum number of iterations
p_train <- 0.025 # Subsampling of the initial data
- pazka!
oma_juttu <- dissever(
coarse = coarse_raster,
fine = raster_stack,
method = "lm",
p = p_train,
min_iter = min_iter,
max_iter = max_iter,
verbose=1
)
- plot(oma_juttu)
writeRaster(oma_juttu$map, filename='out.nc', format="CDF", overwrite=TRUE)
Licensing
[edit]
I, the copyright holder of this work, hereby publish it under the following license:
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.
https://creativecommons.org/licenses/by-sa/4.0CC BY-SA 4.0 Creative Commons Attribution-Share Alike 4.0 truetrue
File history
Click on a date/time to view the file as it appeared at that time.
Date/Time | Thumbnail | Dimensions | User | Comment | |
---|---|---|---|---|---|
current | 15:27, 13 September 2019 | 1,650 × 1,275 (1.72 MB) | Merikanto (talk | contribs) | 10 yrs average 14559-14549 bp | |
14:46, 3 September 2019 | 1,650 × 1,275 (1.19 MB) | Merikanto (talk | contribs) | Tpo correct wrong downscaling | ||
13:30, 4 August 2019 | 1,650 × 1,275 (866 KB) | Merikanto (talk | contribs) | Diffenent time, time correction | ||
15:21, 15 December 2018 | 1,056 × 816 (758 KB) | Merikanto (talk | contribs) | User created page with UploadWizard |
You cannot overwrite this file.
File usage on Commons
The following page uses this file:
- File:Ccsm trace21k tjuly 146585bp europe.svg (file redirect)
Metadata
This file contains additional information such as Exif metadata which may have been added by the digital camera, scanner, or software program used to create or digitize it. If the file has been modified from its original state, some details such as the timestamp may not fully reflect those of the original file. The timestamp is only as accurate as the clock in the camera, and it may be completely wrong.
Width | 1650px |
---|---|
Height | 1275px |