Data CRS ellipsoid added

Author: peterson-tim-j

R/file.summary.R

@@ -35,7 +35,9 @@ file.summary <- function(ncfile) {
                           from = as.Date(character()),
                           to = as.Date(character()),
                           timestep = character(),
-                          units = character()  )
+                          units = character(),
+                          ellipsoid.crs = character()
+                          )
 
   # Get number of variables in all groups
   for (igrp in grps$grps) {
@@ -58,12 +60,15 @@ file.summary <- function(ncfile) {
         grp.string = RNetCDF::grp.inq.nc(igrp)$name
         var.string = paste(grp.string, '/', var.tmp$name, sep='' )
 
+        ellipsoid.crs = RNetCDF::att.get.nc(igrp, 'NC_GLOBAL', 'CRS')
+
         summary.df.new =  data.frame( group = grp.string,
                                       var.string = var.string,
                                       from = date.from,
                                       to = date.to,
                                       timestep = timestep,
-                                      units = units)
+                                      units = units,
+                                      ellipsoid.crs = ellipsoid.crs)
 
         row.names(summary.df.new) = var.tmp$name
 

R/get.variableSource.R

@@ -20,41 +20,47 @@ get.variableSource <- function() {
                         timestep = rep('', nvars),
                         data.URL = rep('', nvars),
                         ncdf.name = rep('', nvars),
+                        ellipsoid.crs = rep('', nvars),
                         infillGaps = rep(F, nvars),
                         row.names = vars)
   var.data['precip',] <- c('Total daily precipitation',
                            'mm/day',
                            'daily',
                            'https://www.bom.gov.au/web03/ncc/www/awap/rainfall/totals/daily/grid/0.05/history/nat/',
                            'precip',
+                           '+proj=longlat +ellps=GRS80',
                            F)
 
   var.data['tmin',] <- c(  'Min daily temperature',
                            'deg_C',
                            'daily',
                            'https://www.bom.gov.au/web03/ncc/www/awap/temperature/minave/daily/grid/0.05/history/nat/',
                            'tmin',
+                           '+proj=longlat +ellps=GRS80',
                            F)
 
   var.data['tmax',] <- c(  'Max daily temperature',
                            'deg_C',
                            'daily',
                            'https://www.bom.gov.au/web03/ncc/www/awap/temperature/maxave/daily/grid/0.05/history/nat/',
                            'tmax',
+                           '+proj=longlat +ellps=GRS80',
                            F)
 
   var.data['vprp',] <- c(  '3pm daily vapour pressure',
                            'hpa',
                            'daily',
                            'https://www.bom.gov.au/web03/ncc/www/awap/vprp/vprph15/daily/grid/0.05/history/nat/',
                            'vprp',
+                           '+proj=longlat +ellps=GRS80',
                            F)
 
   var.data['solarrad',] <- c(  'Total daily solar radiation',
                            'MJ/m^2',
                            'daily',
                            'https://www.bom.gov.au/web03/ncc/www/awap/solar/solarave/daily/grid/0.05/history/nat/',
                            'solarrad',
+                           '+proj=longlat +ellps=GRS80',
                            T)
 
   return(var.data)