File:Ccsm trace21k tjuly 14559bp europe.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.72 MB)

Captions

Captions

Add a one-line explanation of what this file represents

Summary

[edit]
Description
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
 
The SVG code is valid.
 
This map was created with Adobe Illustrator by Merikanto.
This map is not a real SVG file and is just a wrapper containing one or more raster graphics without vector coding.
Using SVG as just a wrapper is undesirable.
Note: This template should be used if the SVG file contains only raster graphics. If the SVG file contains both vector and raster graphics, then use {{BadSVG}}.

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

  1. take snapshot from Trace21k TS file
  1. bash rainsnap.bash
  1. rain_infile=./data/trace_21k_7.nc

rain_infile=./data/ts_14585_raw.nc

orog_infile=./data/ice5g_v1.2_21.0k_10min.nc

  1. from etopo1 cropped and 2.5 min scaled orogeny
  2. orog_infile=/home/hirohito/grid/euro25b.nc

rain_inref=./data/curr_tmean7_10min.nc

ncks -O -C -vlong,lat,orog $orog_infile temp.nc

  1. 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

  1. 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

  1. 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

  1. 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

    1. ro dissever analysis

Rscript distek1.r

Dissever R 3.2 script

distek1.r

    1. basic script downscaling rainfall
    2. or temperature with R dissever geospatial package
    3. R version 3.4.4 (2018-03-15)
  2. you must install libs and data
    1. coarse_raster is coarse temperature or precipitation file, in netcdf,
    2. usually output from simulation
    3. p1 is orography/topography file, raster in netcdf format
    4. p2 training file for rainfall ot remparature
    5. fine rsolution is usually current temperature
    6. resolution of p1, p2 files must be same
    7. 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

    1. 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

)

  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:
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
current15:27, 13 September 2019Thumbnail for version as of 15:27, 13 September 20191,650 × 1,275 (1.72 MB)Merikanto (talk | contribs)10 yrs average 14559-14549 bp
14:46, 3 September 2019Thumbnail for version as of 14:46, 3 September 20191,650 × 1,275 (1.19 MB)Merikanto (talk | contribs)Tpo correct wrong downscaling
13:30, 4 August 2019Thumbnail for version as of 13:30, 4 August 20191,650 × 1,275 (866 KB)Merikanto (talk | contribs)Diffenent time, time correction
15:21, 15 December 2018Thumbnail for version as of 15:21, 15 December 20181,056 × 816 (758 KB)Merikanto (talk | contribs)User created page with UploadWizard

You cannot overwrite this file.

The following page uses this file:

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.

Retrieved from "https://commons.wikimedia.org/w/index.php?title=File:Ccsm_trace21k_tjuly_14559bp_europe.svg&oldid=852038983"