NLDAS NOAH Processing in R

# Libraries
  # Tidyverse
  library(package = "tidyverse")
  library(package = "lubridate")
  # NetCDF/OPeNDAP support
  require("ncdf4")          # netCDF file access
  require("ncdf4.helpers")  # netCDF support files

# directories 

  root_loc = "http://kyrill.ias.sdsmt.edu:8080/thredds/dodsC/NLDAS/NOAH_OUTPUT/"

  yyyy = 2016
  mm   =   07
  dd   =   18
  hh   =   18

# directories 
  root_loc = "http://kyrill.ias.sdsmt.edu:8080/thredds/dodsC/NLDAS/NOAH_OUTPUT/"
  yyyy = 2016
  mm   =   07
  dd   =   18
  hh   =   18
  
  

   # Pulling NLDAS Coordinate Metadata
   yyyy = 2016
   mm   =   07
   dd   =   18
   hh   =   18
   file_loc = str_c(root_loc,
                    str_pad(string = yyyy, 
                            width  =    4, 
                            pad     = "0"),
                    "/",
                    str_pad(string =   mm, 
                            width  =    2, 
                            pad     = "0"),
                    "/",
                    str_pad(string =   dd, 
                            width  =    2, 
                            pad     = "0"),
                    "/NLDAS_NOAH0125_H.",
                    str_pad(string = yyyy, 
                            width  =    4, 
                            pad     = "0"),
                    str_pad(string =   mm, 
                            width  =    2, 
                            pad     = "0"),
                    str_pad(string =   dd, 
                            width  =    2, 
                            pad     = "0"),
                    str_pad(string =   hh, 
                            width  =    2, 
                            pad     = "0"),
                    ".grb",
                    sep = "")
   
   
   nc.file = nc_open(filename = file_loc)  # name of the file
   
   # Extracting Latitude and Longitude Arrays
  
   lon =  ncvar_get(nc    =  nc.file, # netCDF file handle
                    varid =    "lon") # netCDF variable name
  
   lat =  ncvar_get(nc    =  nc.file, # netCDF file handle
                    varid =    "lat") # netCDF variable name   

      file_loc = str_c(root_loc,
                       yyyy,
                       "/",
                       mm,
                       "/",
                       dd,
                       "/NLDAS_NOAH0125_H.",
                       yyyy,
                       mm,
                       dd,
                       hh,
                       ".grb",
                       sep = "")

LS0tCnRpdGxlOiAiTkxEQVMgTk9BSCBQcm9jZXNzaW5nIGluIFIiCm91dHB1dDogaHRtbF9ub3RlYm9vawotLS0KCgoKYGBge3J9CgojIExpYnJhcmllcwoKICAjIFRpZHl2ZXJzZQoKICBsaWJyYXJ5KHBhY2thZ2UgPSAidGlkeXZlcnNlIikKICBsaWJyYXJ5KHBhY2thZ2UgPSAibHVicmlkYXRlIikKCiAgIyBOZXRDREYvT1BlTkRBUCBzdXBwb3J0CgogIHJlcXVpcmUoIm5jZGY0IikgICAgICAgICAgIyBuZXRDREYgZmlsZSBhY2Nlc3MKICByZXF1aXJlKCJuY2RmNC5oZWxwZXJzIikgICMgbmV0Q0RGIHN1cHBvcnQgZmlsZXMKCgoKYGBgCgoKCmBgYHtyfQoKIyBkaXJlY3RvcmllcyAKCiAgcm9vdF9sb2MgPSAiaHR0cDovL2t5cmlsbC5pYXMuc2RzbXQuZWR1OjgwODAvdGhyZWRkcy9kb2RzQy9OTERBUy9OT0FIX09VVFBVVC8iCgogIHl5eXkgPSAyMDE2CiAgbW0gICA9ICAgMDcKICBkZCAgID0gICAxOAogIGhoICAgPSAgIDE4CgoKYGBgCgoKCgoKYGBge3J9CgojIGRpcmVjdG9yaWVzIAoKICByb290X2xvYyA9ICJodHRwOi8va3lyaWxsLmlhcy5zZHNtdC5lZHU6ODA4MC90aHJlZGRzL2RvZHNDL05MREFTL05PQUhfT1VUUFVULyIKCiAgeXl5eSA9IDIwMTYKICBtbSAgID0gICAwNwogIGRkICAgPSAgIDE4CiAgaGggICA9ICAgMTgKICAKICAKCmBgYAoKCgoKYGBge3J9CgogICAjIFB1bGxpbmcgTkxEQVMgQ29vcmRpbmF0ZSBNZXRhZGF0YQoKICAgeXl5eSA9IDIwMTYKICAgbW0gICA9ICAgMDcKICAgZGQgICA9ICAgMTgKICAgaGggICA9ICAgMTgKCgogICBmaWxlX2xvYyA9IHN0cl9jKHJvb3RfbG9jLAogICAgICAgICAgICAgICAgICAgIHN0cl9wYWQoc3RyaW5nID0geXl5eSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB3aWR0aCAgPSAgICA0LCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHBhZCAgICAgPSAiMCIpLAogICAgICAgICAgICAgICAgICAgICIvIiwKICAgICAgICAgICAgICAgICAgICBzdHJfcGFkKHN0cmluZyA9ICAgbW0sIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgd2lkdGggID0gICAgMiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICBwYWQgICAgID0gIjAiKSwKICAgICAgICAgICAgICAgICAgICAiLyIsCiAgICAgICAgICAgICAgICAgICAgc3RyX3BhZChzdHJpbmcgPSAgIGRkLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHdpZHRoICA9ICAgIDIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgcGFkICAgICA9ICIwIiksCiAgICAgICAgICAgICAgICAgICAgIi9OTERBU19OT0FIMDEyNV9ILiIsCiAgICAgICAgICAgICAgICAgICAgc3RyX3BhZChzdHJpbmcgPSB5eXl5LCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHdpZHRoICA9ICAgIDQsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgcGFkICAgICA9ICIwIiksCiAgICAgICAgICAgICAgICAgICAgc3RyX3BhZChzdHJpbmcgPSAgIG1tLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHdpZHRoICA9ICAgIDIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgcGFkICAgICA9ICIwIiksCiAgICAgICAgICAgICAgICAgICAgc3RyX3BhZChzdHJpbmcgPSAgIGRkLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHdpZHRoICA9ICAgIDIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgcGFkICAgICA9ICIwIiksCiAgICAgICAgICAgICAgICAgICAgc3RyX3BhZChzdHJpbmcgPSAgIGhoLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHdpZHRoICA9ICAgIDIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgcGFkICAgICA9ICIwIiksCiAgICAgICAgICAgICAgICAgICAgIi5ncmIiLAogICAgICAgICAgICAgICAgICAgIHNlcCA9ICIiKQogICAKICAgCiAgIG5jLmZpbGUgPSBuY19vcGVuKGZpbGVuYW1lID0gZmlsZV9sb2MpICAjIG5hbWUgb2YgdGhlIGZpbGUKICAgCiAgICMgRXh0cmFjdGluZyBMYXRpdHVkZSBhbmQgTG9uZ2l0dWRlIEFycmF5cwogIAogICBsb24gPSAgbmN2YXJfZ2V0KG5jICAgID0gIG5jLmZpbGUsICMgbmV0Q0RGIGZpbGUgaGFuZGxlCiAgICAgICAgICAgICAgICAgICAgdmFyaWQgPSAgICAibG9uIikgIyBuZXRDREYgdmFyaWFibGUgbmFtZQogIAogICBsYXQgPSAgbmN2YXJfZ2V0KG5jICAgID0gIG5jLmZpbGUsICMgbmV0Q0RGIGZpbGUgaGFuZGxlCiAgICAgICAgICAgICAgICAgICAgdmFyaWQgPSAgICAibGF0IikgIyBuZXRDREYgdmFyaWFibGUgbmFtZSAgIAoKCmBgYAoKCgoKYGBge3J9CmYKICAgICAgZmlsZV9sb2MgPSBzdHJfYyhyb290X2xvYywKICAgICAgICAgICAgICAgICAgICAgICB5eXl5LAogICAgICAgICAgICAgICAgICAgICAgICIvIiwKICAgICAgICAgICAgICAgICAgICAgICBtbSwKICAgICAgICAgICAgICAgICAgICAgICAiLyIsCiAgICAgICAgICAgICAgICAgICAgICAgZGQsCiAgICAgICAgICAgICAgICAgICAgICAgIi9OTERBU19OT0FIMDEyNV9ILiIsCiAgICAgICAgICAgICAgICAgICAgICAgeXl5eSwKICAgICAgICAgICAgICAgICAgICAgICBtbSwKICAgICAgICAgICAgICAgICAgICAgICBkZCwKICAgICAgICAgICAgICAgICAgICAgICBoaCwKICAgICAgICAgICAgICAgICAgICAgICAiLmdyYiIsCiAgICAgICAgICAgICAgICAgICAgICAgc2VwID0gIiIpCgoKCgpgYGAKCg==