| station.obj.dm {clim.pact} | R Documentation |
Create a station object for use as predictand in empirical
downscaling on monthly data. Also see station.obj.
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)
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. |
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. |
R.E. Benestad
## 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)