objDS               package:clim.pact               R Documentation

_O_b_j_e_c_t_i_v_e _d_o_w_n_s_c_a_l_i_n_g _o_f _m_o_n_t_h_l_y _m_e_a_n_s

_D_e_s_c_r_i_p_t_i_o_n:

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

_U_s_a_g_e:

     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",
               plot.res=FALSE,plot.rate=FALSE,xtr.args="",
               swsm="step",predm="predict",lsave=TRUE,rmac=TRUE,
               silent=FALSE)

_A_r_g_u_m_e_n_t_s:

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.

     .

_V_a_l_u_e:

     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

_A_u_t_h_o_r(_s):

     R.E. Benestad

_E_x_a_m_p_l_e_s:

     ## 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)

