station.obj.dm {clim.pact}R Documentation

Make daily climate station series object.

Description

Create a station object for use as predictand in empirical downscaling on monthly data. Also see station.obj.

Usage

station.obj.dm(t2m,precip,dd,mm,yy,
                       obs.name=NULL,unit=NULL,ele=NULL,
                       station=NULL,lat=NULL,lon=NULL,alt=NULL,
                       location="unspecified",wmo.no=NULL,
                       start=NULL,yy0=NULL,country=NULL,ref=NULL)

Arguments

t2m a vector holding daily mean temperature.
precip a vector holding daily precipitation.
dd a vector holding day of month.
mm a vector holding the month.
yy a vector holding the year.
obs.name the name of observation: eg c("Daily mean temperature","Daily precipitation").
unit the unite of observation: eg c("deg C","mm/day").
ele element code: eg c("tam","rr").
station local (national) station number.
lat latitude.
lon longitude.
alt altitude.
location name of location.
wmo.no WMO number of station.
start start of measurements.
yy0 first year of record.
country name of country.
ref reference to the data.

Value

a "daily.station.record"-class object.
t2m a vector holding daily mean temperature.
precip a vector holding daily precipitation.
day a vector holding day of month.
month a vector holding the month.
year a vector holding the year.
obs.name the name of observation:
eg c("Daily mean temperature","Daily precipitation").
unit the unite of observation: eg c("deg C","mm/day").
ele element code: eg c("tam","rr").
station local (national) station number.
lat latitude.
lon longitude.
alt altitude.
location name of location.
wmo.no WMO number of station.
start start of measurements.
yy0 first year of record.
country name of country.
ref reference to the data.

Author(s)

R.E. Benestad

Examples

## Not run: 
blindern.raw <-read.table("~/data/stations/blindern_rr_day.dat",header=TRUE)
blindern.raw$rr[blindern.raw$rr < 0] <- NA
yy <- floor(blindern.raw$yyyymmdd/10000)
mm <- floor(blindern.raw$yyyymmdd/100) - 10000*yy
dd <- blindern.raw$yyyymmdd - 100*mm  - 10000*yy
blindern <- station.obj.dm(t2m=rep(NA,length(blindern.raw$rr)),
                           precip=blindern.raw$rr,
                           dd=dd,mm=mm,yy=yy,
                           obs.name=c("T(2m)","recip"),
                           unit=c("deg C","mm/day"),ele=NULL,
                           station=18700,lat=59.95,lon=10.71,alt=94,
                           location="Oslo-Blindern",wmo.no=NULL,
                           start=NULL,yy0=1937,country="Norway",
                           ref="www.met.no")
## End(Not run)

[Package clim.pact version 2.2-15 Index]