{"id":605,"date":"2014-04-21T15:48:42","date_gmt":"2014-04-21T20:48:42","guid":{"rendered":"http:\/\/homepages.uc.edu\/~yaozo\/wordpress\/?p=605"},"modified":"2014-04-21T15:48:42","modified_gmt":"2014-04-21T20:48:42","slug":"monthly-total-precipitation-mm-for-august-1997-in-the-rocky-mountain-region-and-some-gridded-4km-elevation-data-sets-m","status":"publish","type":"post","link":"https:\/\/zhuoyao.net\/index.php\/2014\/04\/21\/monthly-total-precipitation-mm-for-august-1997-in-the-rocky-mountain-region-and-some-gridded-4km-elevation-data-sets-m\/","title":{"rendered":"Monthly total precipitation (mm) for August 1997 in the Rocky Mountain Region and some gridded 4km elevation data sets (m)."},"content":{"rendered":"\n\n\n
RMprecip<\/td>\nR Documentation<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Monthly total precipitation (mm) for August 1997 in the Rocky Mountain Region and some gridded 4km elevation data sets (m).<\/h2>\n

Description<\/h3>\n

RMprecip<\/code>\u00a0is a useful spatial data set of moderate size consisting of 806 locations. See www.image.ucar.edu\/Data for the source of these data.\u00a0PRISMelevation<\/code>\u00a0and\u00a0RMelevation<\/code>\u00a0are gridded elevations for the continental US and Rocky Mountain region at 4km resolution. Note that the gridded elevations from the PRISM data product are different than the exact station elevations. (See example below.)<\/p>\n

Format<\/h3>\n

The data set\u00a0RMprecip<\/code>\u00a0is a list containing the following components:<\/p>\n

\n
x<\/dt>\n
Longitude-latitude position of monitoring stations. Rows names are station id codes consistent with the US Cooperative observer network. The ranges for these coordinates are [-111, -99] for longitude and [35,45] for latitude.<\/p>\n<\/dd>\n
elev<\/dt>\n
Station elevation in meters.<\/p>\n<\/dd>\n
y<\/dt>\n
Monthly total precipitation in millimeters. for August, 1997<\/p>\n<\/dd>\n<\/dl>\n

The data sets\u00a0PRISMelevation<\/code>\u00a0and\u00a0RMelevation<\/code>\u00a0are lists in the usual R grid format for images and contouring<\/p>\n

They have the following components:<\/p>\n

\n
x<\/dt>\n
Longitude grid at approximately 4km resolution<\/p>\n<\/dd>\n
y<\/dt>\n
Latitude grid at approximately 4km resolution<\/p>\n<\/dd>\n
z<\/dt>\n
Average elevation for grid cell in meters<\/p>\n<\/dd>\n<\/dl>\n

These elevations and the companion grid formed the basis for the 103-Year High-Resolution Precipitation Climate Data Set for the Conterminous United States\u00a0ftp:\/\/ftp.ncdc.noaa.gov\/pub\/data\/prism100<\/a>\u00a0archived at the National Climate Data Center. This work was primarily authored by Chris Daly\u00a0www.prism.oregonstate.edu<\/a>\u00a0and his PRISM group but had some contribution from the Geophysical Statistics Project at NCAR. and is an interpolation of the observational data to a 4km grid that takes into account topography such as elevation and aspect.<\/p>\n

Details<\/h3>\n

The binary file\u00a0RData.USmonthlyMet.bin<\/code>\u00a0can be downwloaded from\u00a0http:\/\/www.image.ucar.edu\/Data\/US.monthly.met<\/a>\u00a0and also includes information on its source.<\/p>\n

# explicit source code to create the RMprecip data\ndir <- \"\" # include path to data file \nload(paste(dir, \"RData.USmonthlyMet.bin\", sep=\"\/\")\n#year.id<-  1963- 1895\nyear.id<- 103\n#pptAUG63<- USppt[ year.id,8,]\nloc<- cbind(USpinfo$lon, USpinfo$lat)\nxr<- c(-111, -99)\nyr<- c( 35, 45)\nstation.subset<-  (loc[,1]>= xr[1]) & (loc[,1] <= xr[2]) & (loc[,2]>= yr[1]) & (loc[,2]<= yr[2])\nydata<-  USppt[ year.id,8,station.subset]\nydata <- ydata*10 #  cm -> mm conversion\nxdata<- loc[station.subset,]\ndimnames(xdata)<- list( USpinfo$station.id[station.subset], c( \"lon\", \"lat\"))\nxdata<- data.frame( xdata)\ngood<- !is.na(ydata)\nydata<- ydata[good]\nxdata<- xdata[good,]\n     \ntest.for.zero.flag<- 1\ntest.for.zero( unlist(RMprecip$x), unlist(xdata), tag=\"locations\")\ntest.for.zero( ydata, RMprecip$y, \"values\")\n<\/pre>\n
Examples<\/h3>\n
# this data set was created  the \n# historical data  taken from \n# Observed monthly precipitation, min and max temperatures for the coterminous US \n# 1895-1997\n# NCAR_pinfill \n# see the Geophysical Statistics Project datasets page for the supporting functions \n# and details. \n\n# plot \nquilt.plot(RMprecip$x, RMprecip$y)\nUS( add=TRUE, col=2, lty=2)\n\n# comparison of station elevations with PRISM gridded values\n\ndata(RMelevation)\n\ninterp.surface( RMelevation, RMprecip$x)-> test.elev\n\nplot( RMprecip$elev, test.elev, xlab=\"Station elevation\", \nylab=\"Interpolation from PRISM grid\")\nabline( 0,1,col=\"blue\")\n\n# some differences  with high elevations probably due to complex\n# topography!\n\n#\n# view of Rockies looking from theSoutheast\n\nsave.par<- par(no.readonly=TRUE)\n\npar( mar=c(0,0,0,0))\n\n# fancy use of persp with shading and lighting.\npersp( RMelevation, theta=75, phi= 15, \n          box=FALSE, axes=FALSE, xlab=\"\", ylab=\"\", \n         border=NA,\n         shade=.95, lphi= 10, ltheta=80,\n         col= \"wheat4\", \n         scale=FALSE, expand=.00025)\n\n# reset graphics parameters and a more conventional image plot.\npar( save.par)\nimage.plot(RMelevation, col=topo.colors(256))\nUS( add=TRUE, col=\"grey\", lwd=2)\ntitle(\"PRISM elevations (m)\")\n<\/pre>\n
Results<\/h3>\n
R version 3.0.2 (2013-09-25) -- \"Frisbee Sailing\"\nCopyright (C) 2013 The R Foundation for Statistical Computing\nPlatform: x86_64-unknown-linux-gnu (64-bit)\n\nR is free software and comes with ABSOLUTELY NO WARRANTY.\nYou are welcome to redistribute it under certain conditions.\nType 'license()' or 'licence()' for distribution details.\n\nR is a collaborative project with many contributors.\nType 'contributors()' for more information and\n'citation()' on how to cite R or R packages in publications.\n\nType 'demo()' for some demos, 'help()' for on-line help, or\n'help.start()' for an HTML browser interface to help.\nType 'q()' to quit R.\n\n> library(fields)\nLoading required package: spam\nSpam version 0.29-3 (2013-04-23) is loaded.\nType 'help( Spam)' or 'demo( spam)' for a short introduction \nand overview of this package.\nHelp for individual functions is also obtained by adding the\nsuffix '.spam' to the function name, e.g. 'help( chol.spam)'.\n\nAttaching package: 'spam'\n\nThe following objects are masked from 'package:base':\n\n    backsolve, forwardsolve\n\nLoading required package: maps\n> png(filename=\"\/home\/oogasawa\/snapshot\/RGM3\/R_CC\/result\/fields\/RMprecip.Rd_%03d_medium.png\", width=480, height=480)\n> ### Name: RMprecip\n> ### Title: Monthly total precipitation (mm) for August 1997 in the Rocky\n> ###   Mountain Region and some gridded 4km elevation data sets (m).\n> ### Aliases: RMprecip RMelevation PRISMelevation\n> ### Keywords: datasets\n> \n> ### ** Examples\n> \n> # this data set was created  the \n> # historical data  taken from \n> # Observed monthly precipitation, min and max temperatures for the coterminous US \n> # 1895-1997\n> # NCAR_pinfill \n> # see the Geophysical Statistics Project datasets page for the supporting functions \n> # and details. \n> \n> # plot \n> quilt.plot(RMprecip$x, RMprecip$y)\n> US( add=TRUE, col=2, lty=2)\n> \n> # comparison of station elevations with PRISM gridded values\n> \n> data(RMelevation)\n> \n> interp.surface( RMelevation, RMprecip$x)-> test.elev\n> \n> plot( RMprecip$elev, test.elev, xlab=\"Station elevation\", \n+ ylab=\"Interpolation from PRISM grid\")\n> abline( 0,1,col=\"blue\")\n> \n> # some differences  with high elevations probably due to complex\n> # topography!\n> \n> #\n> # view of Rockies looking from theSoutheast\n> \n> save.par<- par(no.readonly=TRUE)\n> \n> par( mar=c(0,0,0,0))\n> \n> # fancy use of persp with shading and lighting.\n> persp( RMelevation, theta=75, phi= 15, \n+           box=FALSE, axes=FALSE, xlab=\"\", ylab=\"\", \n+          border=NA,\n+          shade=.95, lphi= 10, ltheta=80,\n+          col= \"wheat4\", \n+          scale=FALSE, expand=.00025)\n> \n> # reset graphics parameters and a more conventional image plot.\n> par( save.par)\n> image.plot(RMelevation, col=topo.colors(256))\n> US( add=TRUE, col=\"grey\", lwd=2)\n> title(\"PRISM elevations (m)\")\n> \n> \n> \n> \n> \n> dev.off()\nnull device \n          1 \n>\n<\/pre>\n
 <\/p>\n\n\n\n
\"\"<\/a>\"\"<\/a>\"\"<\/a>\"\"<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n","protected":false},"excerpt":{"rendered":"
RMprecip R Documentation Monthly total precipitation (mm) for August 1997 in the Rocky Mountain Region and some gridded 4km elevation data sets (m). Description RMprecip\u00a0is… <\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[20],"tags":[],"class_list":["post-605","post","type-post","status-publish","format-standard","hentry","category-r"],"_links":{"self":[{"href":"https:\/\/zhuoyao.net\/index.php\/wp-json\/wp\/v2\/posts\/605","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/zhuoyao.net\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/zhuoyao.net\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/zhuoyao.net\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/zhuoyao.net\/index.php\/wp-json\/wp\/v2\/comments?post=605"}],"version-history":[{"count":0,"href":"https:\/\/zhuoyao.net\/index.php\/wp-json\/wp\/v2\/posts\/605\/revisions"}],"wp:attachment":[{"href":"https:\/\/zhuoyao.net\/index.php\/wp-json\/wp\/v2\/media?parent=605"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/zhuoyao.net\/index.php\/wp-json\/wp\/v2\/categories?post=605"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/zhuoyao.net\/index.php\/wp-json\/wp\/v2\/tags?post=605"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}