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)