objDS {clim.pact} | R Documentation |
'Objective' downscaling based on DS
. The
function selects region according to a correlation analysis, setting
the borders where the correlation with the station series become
zero. A fit to a truncated Fourier expansion is used to describe
profiles of correlation coefficients in zonal and meridional
directions from where the station is located.
Version 2.1-5: smallest spatial domain for objDS is set to +- 10deg N and E of the station location in order to avoid problems with dimensions if one only contains one data point (the domain becomes a line instead of a rectangle).
objDS(field.obs,field.gcm,station,plot=TRUE,positive=NULL, mon=NULL,direc="output/",cal.id=NULL, ldetrnd=TRUE,i.eofs=seq(1,8,by=1),ex.tag="", method="lm",leps=FALSE,param="t2m",failure.action=NULL, plot.res=FALSE,plot.rate=FALSE,xtr.args="", swsm="step",predm="predict",lsave=FALSE,rmac=TRUE, silent=FALSE,qualitycontrol=TRUE,LINPACK=TRUE,wOBS=0.25)
field.obs |
The gridded observation predictor retrieve.nc . |
field.gcm |
The climate simulation predictor retrieve.nc . |
station |
A climate.station object (station.obj or
station.obj.dm ).
[e.g. from getnacd , getnordklim or station.obj ]. |
plot |
'TRUE' produces figures. |
positive |
'TRUE': only consider the region where correlations are positive (important for temperature predictors). |
mon |
month or season to downscale,loops though the 12 calendar months if NULL. |
direc |
name of directory inwhich the output is dumped (e.g. figures, tables). |
cal.id |
ID tag used for calibration. By default use the
first field (catFields ) for calibration. |
ldetrnd |
F for no detrending; T for removing linear trends before model calibration. |
i.eofs |
select which EOFs to include in the setp-wise screening. |
ex.tag |
Extra labelling tag for file names for experiments. |
method |
Sets the method to use for regression. Method is set
to "lm" by default, but "anm" allows the incorporation of an
analog model, see anm . "anm.weight" weights the
principal components according to the eigenvalues, whereas "anm"
uses unweighted series. |
leps |
'TRUE' produces EPS figures (files). |
param |
Name of parameter (for plot labels). |
plot.res |
'TRUE' shows statistics for residuals. |
plot.rate |
'TRUE' shows analysis of rate-of-change. |
xtr.args |
Extra/additional arguments in the formula. |
swsm |
Step-wise screening method, default=='step'; 'none' skips stepwise sceeening. |
predm |
Prediction method, default is "predict" |
lsave |
TRUE -> saves the result on disc |
rmac |
TRUE -> subtracts (removes) the annual cycle in station data. |
silent |
TRUE -> no output to screen. |
qualitycontrol |
TRUE: perform a qualitry control consisting of comparing the smoothness of the monthly trend estimates throughout the year and repeat the downscaling with a smaller domain if adjacent trend estimates vary significantly (diff greater than 3*variance of rates). |
LINPACK |
'TRUE': svd; 'FALSE':La.svd |
wOBS |
Used for weighting down GCM results in the common EOF
analysis. The GCM results are re-scaled after the
analysis.wOBS=NULL skips this scaling/re-scaling. |
failure.action |
If stepwise screening in
DS fails to select any variables, then call the
function specified by a character string, unless set to NULL. |
An 'objDS' object - a list of objects:
station | a 'station' object (see getnacd ) |
Jan | a 'ds' (see DS ) object |
Feb | a 'ds' object |
Mar | a 'ds' object |
Apr | a 'ds' object |
May | a 'ds' object |
Jun | a 'ds' object |
Jul | a 'ds' object |
Aug | a 'ds' object |
Sep | a 'ds' object |
Oct | a 'ds' object |
Nov | a 'ds' object |
Dec | a 'ds' object |
R.E. Benestad
## Not run: library(clim.pact) oslo<-getnordklim("Oslo-Blindern") slp.obs <- retrieve.nc("ncep_slp.nc") # Get gridded observations/analysis from NCEP slp.gcm <- retrieve.nc("EH4OPYC_B2_slp.nc") # Get results from climate models ds <- objDS(field.obs=slp.obs,field.gcm=slp.gcm,station=oslo) t2m <- retrieve.nc("~/data/analysis/DNMI_t2m.nc") t2m.gcm <- retrieve.nc("~/data/mpi/mpi-gsdio_t2m.nc", x.rng=range(t2m$lon),y.rng=range(t2m$lat)) ds <- objDS(field.obs=t2m,field.gcm=t2m.gcm,station=oslo) ## End(Not run)