Projection and New Discharge Name

[ChuliangXiao]

• GetProj function in geospatial.R

  map_proj <- ncdf4::ncatt_get(coordNC, varid=0, attname="MAP_PROJ")$value
  cen_lat <- ncdf4::ncatt_get(coordNC, varid=0, attname="MOAD_CEN_LAT")$value
  cen_lon <- ncdf4::ncatt_get(coordNC, varid=0, attname="STAND_LON")$value
  truelat1 <- ncdf4::ncatt_get(coordNC, varid=0, attname="TRUELAT1")$value
  truelat2 <- ncdf4::ncatt_get(coordNC, varid=0, attname="TRUELAT2")$value

  CEN_LAT and CEN_LON are specific to the nested domain. If available, MOAD_CEN_LAT and STAND_LON are should be used.

• TransUsgsProdStat function in usgs_ingest.R

  prodStatLookup <- c( X_00060_00011 ='Discharge (cfs)', ##instantaneous is 00011 but not worth saying IMO
  X_00060_00011_cd ='Discharge code',
  X_00065_00011 ='Stage (ft)',
  X_00065_00011_cd ='Stage code' )
  code2Name <- any(names %in% names(prodStatLookup))

  'Discharge (cfs)' in new USGS data is denoted by X_00060_00000 instead of X_00060_00011.