From edb827a9e618ce5f178d9e475adc3a131dec0bb5 Mon Sep 17 00:00:00 2001 From: jennifer Date: Wed, 16 Aug 2023 12:19:35 -0600 Subject: [PATCH 1/4] Add SNODAS folder with R tutorials --- .../Download_Process_masked_SNODAS_data.Rmd | 214 +++++ .../Download_Process_masked_SNODAS_data.ipynb | 464 +++++++++++ ...ownload_Process_masked_SNODAS_data.nb.html | 775 ++++++++++++++++++ notebooks/SNODAS/img/nsidc_logo.png | Bin 0 -> 17603 bytes 4 files changed, 1453 insertions(+) create mode 100644 notebooks/SNODAS/Download_Process_masked_SNODAS_data.Rmd create mode 100644 notebooks/SNODAS/Download_Process_masked_SNODAS_data.ipynb create mode 100644 notebooks/SNODAS/Download_Process_masked_SNODAS_data.nb.html create mode 100644 notebooks/SNODAS/img/nsidc_logo.png diff --git a/notebooks/SNODAS/Download_Process_masked_SNODAS_data.Rmd b/notebooks/SNODAS/Download_Process_masked_SNODAS_data.Rmd new file mode 100644 index 0000000..3f4b725 --- /dev/null +++ b/notebooks/SNODAS/Download_Process_masked_SNODAS_data.Rmd @@ -0,0 +1,214 @@ +--- +title: "How to download and work with SNDOAS data" +output: html_notebook +--- + +## 1. Tutorial Overview +This notebook demonstrates how to download and unzip masked SNODAS files that are from after October 2013 and select one variable and compile it into a RasterStack for quick visualization and finally save as a .grd file. + +### Credits +This tutorial is based on an R script written in 2021 by Cara Thompson at New Mexico State University and has been converted to a tutorial with some minor additions by Jennifer Roebuck of NSIDC. + +For questions regarding the tutorial or to report problems, please create a new issue in the [NSIDC-Data-Tutorials repo](https://github.com/nsidc/NSIDC-Data-Tutorials/issues) + +### Learning Objectives + +By the end of this demonstration you will be able to: +1. Download and unzip masked SNODAS files +2. Convert the files to .bil format and put them in a RasterStack +3. Quickly visualize the rasters of the data +4. Save the RasterStack as .grd file + +### Prerequisites +1. The following libraries should already be installed: R.utils, rvest, dplyr, raster, sp, maptools, and data.table. + + +### Time requirement +Allow approximately 20 minutes to complete this tutorial. + +## 2. Tutorial steps + +### Load necessary libraries +```{r} +library(R.utils) +library(rvest) +library(dplyr) +library(raster) +library(sp) +library(maptools) +library(data.table) +``` + +### Download SNODAS data +We will download all the SNODAS masked files for June 2023. If you wish to download masked files from a different time period you can modify the path in the url variable below. + +```{r} +# set the URL for the SNODAS directory we want to download files from +url <- "https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun" + +# read html content from url +page <- read_html(url) + +# Get a list of the files listed in the directory at this url +files <- page %>% html_nodes("a") %>% html_attr("href") + +for(i in 2:length(files)){ + # generate the url for each of the files + u <- paste(url,files[i], sep="/") + # download each of the files + download.file(u,files[i], mode = "wb") +} +``` + +We will save the working directory path to a variable. +```{r} +WD <- getwd() +``` + +### Read in the SNODAS files +Next, we read in the .tar files as a list containing .tar +```{r} +files_1 <- list.files(path = WD, pattern = "*.tar", full.names = T, recursive = T) +``` + +Decompress the .tar files into the working directory +```{r} +asd <- sapply(files_1, untar) +``` + +### Select the variable we wish to use +There are eight parameters from SNODAS, each has an associated ID code which is used the filename. A full description of the ID codes can be found in Table 1 of the [SNODAS user guide](https://nsidc.org/sites/default/files/g02158-v001-userguide_2_1.pdf). We will choose one of the parameters to compile all the files related to this parameter. + +We can uncomment one of the variables below that corresponds to the parameter we are interested in. In this tutorial we will use the snow depth parameter as an example. + +```{r} +id_code = "us_ssmv11036tS" # Snow Depth +#id_code = "us_ssmv11034tS" # SWE +#id_code = "us_ssmv11044bS" # Snow melt runoff at base of snow pack +#id_code = "us_ssmv11050IL00" # Sublimation from the snow pack +#id_code = "us_ssmv11039IL00" # sublimation of blowing snow +#id_code = "us_ssmv11038wS" # snow pack average temperature +#id_code = "us_ssmv01025SIL01" # Solid precipitation +#id_code = "us_smsv01025SIL00" # Liquid precipitation +``` + +Now we will compile all the files that contain this id_code in the filename i.e. a .dat.gz file and .txt.gz file for each day of data that we downloaded. + +```{r} +files_2 <- list.files(path = WD, pattern = id_code, full.names = T, recursive = T) +files_2 +``` + +Decompress all the .gz files. +```{r} +dcompSD <- sapply(files_2, gunzip, remove=T) +``` + +We will create a new directory for these uncompressed files and move them over. +```{r} +dir.create("Snow Depth") +file.copy(dcompSD, "Snow Depth") +``` + +### Convert the .dat files to .bil +Next we need to convert the .dat files to .bil so they can be read as rasters. + +```{r} +files <- list.files(path=paste(WD,"/Snow Depth",sep=""), pattern="*.dat", full.names = T) +newfiles <- gsub(".dat$", ".bil", files) +file.rename(files, newfiles) +newfiles +``` + +### Create a text file header for each file +Now we will create a corresponding text header file for each .bil file using the associated .txt file that was downloaded. + +**Please note** that the parameters used in the m_post_2013 variable below are specific to SNODAS masked files for October 2013 onwards. If you wish to process masked files prior to October 2013 then you would need to change ulxmap and ulymap to the following: +ulxmap -124.72958333333334 +ulymap 52.870416666666664 + +If you wish to process unmasked files prior to October 2013 then you would need to change nrow, ncols, ulxmap and ulymap to the following: +nrows 4096 +ncols 8192 +ulxmap -130.51291666666665 +ulymap 58.228750000000005 + +If you wish to process unmasked files from October 2013 onwards then you would need to change the nrows, ncols, ulxmap and ulymap to the following: +nrows 4096 +ncols 8192 +ulxmap -130.51250000000002 +ulymap 58.229166666666664 +```{r} + +mylist<- list.files(path=paste(WD,"/Snow Depth",sep=""), pattern = "*.bil", full.names = T) +m_post_2013<-"SNODAS Snow Depth Description = +nrows 3351 +ncols 6935 +nbands 1 +nbits 16 +pixeltype signedint +byteorder M +layout dat +ulxmap -124.72916666666667 +ulymap 52.87083333333334 +xdim 0.00833333333 +ydim 0.00833333333 +pixeltype SIGNEDINT +nodata -9999" + + +mylist.hdr<-paste(mylist,".hdr",sep="") + +for(i in mylist.hdr) { + write(m_post_2013,i) +} + +``` + +Change the file type of the header text files from .bil.hdr to .hdr since the .bil file is already separate. + +```{r} +files4 <- list.files(path=paste(WD,"/Snow Depth",sep=""), pattern="*.bil.hdr", full.names = T) +newfiles <- gsub(".bil.hdr$", ".hdr", files4) +file.rename(files4, newfiles) +newfiles +``` + +### Create a RasterStack of all the rasters +Create a RasterStack of all the snow depth raster files. +```{r} +list.ras <- list.files(path=paste(WD,"/Snow Depth",sep=""), pattern = ".bil", full.names = T) +list.ras +ras_stack <- stack(list.ras) +ras_stack +``` + +### Plot RasterStack +Now we define a projection for the RasterStack and plot each of the rasters (must be in the format of what it is in, not what you want it to be). Please note that SNODAS data are point estimates in latitude/longitude coordinates with the horizontal datum WGS 84 thus the decision to project the data and the choice of projection is up to the user. In this tutorial we will use Geographic longitude/latitude. + +```{r} +crs.geo <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84") +proj4string(ras_stack) <- crs.geo +ras_stack +plot(ras_stack) +``` + +### Save RasterStack +Save the RasterStack to a new directory. We will write in the "raster"/native format to preserve layer names (i.e. dates). There is no attribute table because there is only one attribute per file (snow depth in mm). + +```{r} +#create a new directory Rasters +dir.create(paste(WD,"/Snow Depth","/Rasters",sep="")) +fp <- paste(WD,"/Snow Depth","/Rasters",sep="") +#save the RasterStack as a .grd file +writeRaster(ras_stack,filename=file.path(fp,"ras_stack.grd"), format="raster", progress="text", options=c("INTERLEAVE=BAND","COMPRESS=LZW")) + +``` + +## 3. Learning outcomes recap + +We have learned how to: +1. Download and unzip masked SNODAS files +2. Convert the files to .bil format and put them in a raster stack +3. Quickly visualize the rasters of the data +4. Save the raster stack as .grd file \ No newline at end of file diff --git a/notebooks/SNODAS/Download_Process_masked_SNODAS_data.ipynb b/notebooks/SNODAS/Download_Process_masked_SNODAS_data.ipynb new file mode 100644 index 0000000..8ff03aa --- /dev/null +++ b/notebooks/SNODAS/Download_Process_masked_SNODAS_data.ipynb @@ -0,0 +1,464 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "id": "c0b71b3c", + "metadata": {}, + "source": [ + "
\n", + "\n", + "\n", + "# How to download and work with SNODAS data\n", + "\n", + "
" + ] + }, + { + "cell_type": "markdown", + "id": "c55a92ad", + "metadata": {}, + "source": [ + "## 1. Tutorial Overview \n", + "This notebook demonstrates how to download and unzip masked SNODAS files that are from after October 2013 and select one variable and compile it into a RasterStack for quick visualization and finally save as a .grd file.\n", + "\n", + "### Credits\n", + "This tutorial is based on an R script written in 2021 by Cara Thompson at New Mexico State University and has been converted to a tutorial by Jennifer Roebuck of NSIDC.\n", + "\n", + "For questions regarding the tutorial or to report problems, please create a new issue in the [NSIDC-Data-Tutorials repo](https://github.com/nsidc/NSIDC-Data-Tutorials/issues)\n", + "\n", + "### Learning Objectives\n", + "\n", + "By the end of this demonstration you will be able to:\n", + "1. Download and unzip masked SNODAS files\n", + "2. Convert the files to .bil format and put them in a RasterStack\n", + "3. Quickly visualize the rasters of the data\n", + "4. Save the RasterStack as .grd file\n", + "\n", + "### Prerequisites\n", + "1. The following libraries should already be installed: R.utils, rvest, dplyr, raster, sp, maptools, and data.table.\n", + "\n", + "\n", + "### Time requirement\n", + "Allow approximately 20 minutes to complete this tutorial." + ] + }, + { + "cell_type": "markdown", + "id": "f8733fb0", + "metadata": {}, + "source": [ + "## 2. Tutorial steps \n", + "\n", + "### Load necessary libraries" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "69083e0b", + "metadata": {}, + "outputs": [], + "source": [ + "library(R.utils)\n", + "library(rvest)\n", + "library(dplyr)\n", + "library(raster)\n", + "library(sp)\n", + "library(maptools)\n", + "library(data.table)" + ] + }, + { + "cell_type": "markdown", + "id": "034f3efa", + "metadata": {}, + "source": [ + "### Download SNODAS data\n", + "We will download all the SNODAS masked files for June 2023. If you wish to download masked files from a different time period you can modify the path in the url variable below. " + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "9c7dc73f", + "metadata": {}, + "outputs": [], + "source": [ + "# set the URL for the SNODAS directory we want to download files from\n", + "url <- \"https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun\" \n", + "\n", + "# read html content from url\n", + "page <- read_html(url)\n", + "\n", + "# Get a list of the files listed in the directory at this url\n", + "files <- page %>% html_nodes(\"a\") %>% html_attr(\"href\")\n", + "\n", + "for(i in 2:length(files)){\n", + " # generate the url for each of the files \n", + " u <- paste(url,files[i], sep=\"/\")\n", + " # download each of the files\n", + " download.file(u,files[i], mode = \"wb\")\n", + "}" + ] + }, + { + "cell_type": "markdown", + "id": "4bfdbe57", + "metadata": {}, + "source": [ + "We will save the working directory path to a variable." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "64af70ee", + "metadata": {}, + "outputs": [], + "source": [ + "WD <- getwd()" + ] + }, + { + "cell_type": "markdown", + "id": "b951d2bf", + "metadata": {}, + "source": [ + "### Read in the SNODAS files\n", + "Next, we read in the .tar files as a list containing .tar" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "9a8ba949", + "metadata": {}, + "outputs": [], + "source": [ + "files_1 <- list.files(path = WD, pattern = \"*.tar\", full.names = T, recursive = T)" + ] + }, + { + "cell_type": "markdown", + "id": "20e4077c", + "metadata": {}, + "source": [ + "Decompress the .tar files into the working directory" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "f34d53fa", + "metadata": {}, + "outputs": [], + "source": [ + "asd <- sapply(files_1, untar)" + ] + }, + { + "cell_type": "markdown", + "id": "f91f1493", + "metadata": {}, + "source": [ + "### Select the variable we wish to use\n", + "There are eight parameters from SNODAS, each has an associated ID code which is used the filename. A full description of the ID codes can be found in Table 1 of the [SNODAS user guide](https://nsidc.org/sites/default/files/g02158-v001-userguide_2_1.pdf). We will choose one of the parameters to compile all the files related to this parameter.\n", + "\n", + "We can uncomment one of the variables below that corresponds to the parameter we are interested in. In this tutorial we will use the snow depth parameter as an example. " + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "5b7c1b18", + "metadata": {}, + "outputs": [], + "source": [ + "id_code = \"us_ssmv11036tS\" # Snow Depth\n", + "#id_code = \"us_ssmv11034tS\" # SWE\n", + "#id_code = \"us_ssmv11044bS\" # Snow melt runoff at base of snow pack\n", + "#id_code = \"us_ssmv11050IL00\" # Sublimation from the snow pack\n", + "#id_code = \"us_ssmv11039IL00\" # sublimation of blowing snow\n", + "#id_code = \"us_ssmv11038wS\" # snow pack average temperature\n", + "#id_code = \"us_ssmv01025SIL01\" # Solid precipitation \n", + "#id_code = \"us_smsv01025SIL00\" # Liquid precipitation" + ] + }, + { + "cell_type": "markdown", + "id": "52cfb398", + "metadata": {}, + "source": [ + "Now we will compile all the files that contain this id_code in the filename i.e. a .dat.gz file and .txt.gz file for each day of data that we downloaded." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "328a5934", + "metadata": {}, + "outputs": [], + "source": [ + "files_2 <- list.files(path = WD, pattern = id_code, full.names = T, recursive = T)\n", + "files_2" + ] + }, + { + "cell_type": "markdown", + "id": "15e5313e", + "metadata": {}, + "source": [ + "Decompress all the .gz files." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "50fb8483", + "metadata": {}, + "outputs": [], + "source": [ + "dcompSD <- sapply(files_2, gunzip, remove=T)" + ] + }, + { + "cell_type": "markdown", + "id": "b0747e19", + "metadata": {}, + "source": [ + "We will create a new directory for these uncompressed files and move them over." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "a7504997", + "metadata": {}, + "outputs": [], + "source": [ + "dir.create(\"Snow Depth\")\n", + "file.copy(dcompSD, \"Snow Depth\")" + ] + }, + { + "cell_type": "markdown", + "id": "7c11e739", + "metadata": {}, + "source": [ + "### Convert the .dat files to .bil\n", + "Next we need to convert the .dat files to .bil so they can be read as rasters." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "914a2108", + "metadata": {}, + "outputs": [], + "source": [ + "files <- list.files(path=paste(WD,\"/Snow Depth\",sep=\"\"), pattern=\"*.dat\", full.names = T)\n", + "newfiles <- gsub(\".dat$\", \".bil\", files)\n", + "file.rename(files, newfiles)\n", + "newfiles" + ] + }, + { + "cell_type": "markdown", + "id": "2d4b6101", + "metadata": {}, + "source": [ + "### Create a text file header for each file\n", + "Now we will create a corresponding text header file for each .bil file using the associated .txt file that was downloaded. \n", + "\n", + "**Please note** that the parameters used in the m_post_2013 variable below are specific to SNODAS masked files for October 2013 onwards. If you wish to process masked files prior to October 2013 then you would need to change ulxmap and ulymap to the following:\n", + "\n", + "`ulxmap -124.72958333333334`\n", + "\n", + "`ulymap 52.870416666666664`\n", + "\n", + "If you wish to process unmasked files prior to October 2013 then you would need to change nrow, ncols, ulxmap and ulymap to the following:\n", + "\n", + "`nrows 4096`\n", + "\n", + "`ncols 8192`\n", + "\n", + "`ulxmap -130.51291666666665`\n", + "\n", + "`ulymap 58.228750000000005`\n", + "\n", + "\n", + "If you wish to process unmasked files from October 2013 onwards then you would need to change the nrows, ncols, ulxmap and ulymap to the following:\n", + "\n", + "`nrows 4096`\n", + "\n", + "`ncols 8192`\n", + "\n", + "`ulxmap -130.51250000000002`\n", + "\n", + "`ulymap 58.229166666666664`\n" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "1ecf55dd", + "metadata": {}, + "outputs": [], + "source": [ + "mylist<- list.files(path=paste(WD,\"/Snow Depth\",sep=\"\"), pattern = \"*.bil\", full.names = T)\n", + "m_post_2013<-\"SNODAS Snow Depth Description = \n", + "nrows 3351\n", + "ncols 6935\n", + "nbands 1\n", + "nbits 16\n", + "pixeltype signedint\n", + "byteorder M\n", + "layout dat\n", + "ulxmap -124.72916666666667\n", + "ulymap 52.87083333333334\n", + "xdim 0.00833333333\n", + "ydim 0.00833333333\n", + "pixeltype SIGNEDINT\n", + "nodata -9999\"\n", + "\n", + "\n", + "mylist.hdr<-paste(mylist,\".hdr\",sep=\"\")\n", + "\n", + "for(i in mylist.hdr) {\n", + " write(m_post_2013,i)\n", + "}" + ] + }, + { + "cell_type": "markdown", + "id": "1b0f96ca", + "metadata": {}, + "source": [ + "Change the file type of the header text files from .bil.hdr to .hdr since the .bil file is already separate." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "d9c5996e", + "metadata": {}, + "outputs": [], + "source": [ + "files4 <- list.files(path=paste(WD,\"/Snow Depth\",sep=\"\"), pattern=\"*.bil.hdr\", full.names = T)\n", + "newfiles <- gsub(\".bil.hdr$\", \".hdr\", files4)\n", + "file.rename(files4, newfiles)\n", + "newfiles" + ] + }, + { + "cell_type": "markdown", + "id": "458f61c5", + "metadata": {}, + "source": [ + "### Create a RasterStack of all the rasters\n", + "Create a RasterStack of all the snow depth raster files. " + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "91fec11f", + "metadata": {}, + "outputs": [], + "source": [ + "list.ras <- list.files(path=paste(WD,\"/Snow Depth\",sep=\"\"), pattern = \".bil\", full.names = T)\n", + "list.ras\n", + "ras_stack <- stack(list.ras)\n", + "ras_stack" + ] + }, + { + "cell_type": "markdown", + "id": "d7cdc584", + "metadata": {}, + "source": [ + "### Plot RasterStack\n", + "Now we define a projection for the RasterStack and plot each of the rasters (must be in the format of what it is in, not what you want it to be). Please note that SNODAS data are point estimates in latitude/longitude coordinates with the horizontal datum WGS 84 thus the decision to project the data and the choice of projection is up to the user. In this tutorial we will use Geographic longitude/latitude." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "9a6fed2e", + "metadata": {}, + "outputs": [], + "source": [ + "crs.geo <- CRS(\"+proj=longlat +ellps=WGS84 +datum=WGS84\")\n", + "proj4string(ras_stack) <- crs.geo\n", + "ras_stack\n", + "plot(ras_stack)" + ] + }, + { + "cell_type": "markdown", + "id": "e5584358", + "metadata": {}, + "source": [ + "### Save RasterStack\n", + "Save the RasterStack to a new directory. We will write in the \"raster\"/native format to preserve layer names (i.e. dates). There is no attribute table because there is only one attribute per file (snow depth in mm)." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "6fec3702", + "metadata": {}, + "outputs": [], + "source": [ + "#create a new directory Rasters\n", + "dir.create(paste(WD,\"/Snow Depth\",\"/Rasters\",sep=\"\"))\n", + "fp <- paste(WD,\"/Snow Depth\",\"/Rasters\",sep=\"\")\n", + "#save the RasterStack as a .grd file\n", + "writeRaster(ras_stack,filename=file.path(fp,\"ras_stack.grd\"), format=\"raster\", progress=\"text\", options=c(\"INTERLEAVE=BAND\",\"COMPRESS=LZW\"))" + ] + }, + { + "cell_type": "markdown", + "id": "83a04c93", + "metadata": {}, + "source": [ + "## 3. Learning outcomes recap\n", + "\n", + "We have learned how to:\n", + "1. Download and unzip masked SNODAS files\n", + "2. Convert the files to .bil format and put them in a raster stack\n", + "3. Quickly visualize the rasters of the data\n", + "4. Save the raster stack as .grd file" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "646d3b9d", + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3 (ipykernel)", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.9.15" + } + }, + "nbformat": 4, + "nbformat_minor": 5 +} diff --git a/notebooks/SNODAS/Download_Process_masked_SNODAS_data.nb.html b/notebooks/SNODAS/Download_Process_masked_SNODAS_data.nb.html new file mode 100644 index 0000000..ec8a266 --- /dev/null +++ b/notebooks/SNODAS/Download_Process_masked_SNODAS_data.nb.html @@ -0,0 +1,775 @@ + + + + + + + + + + + + + +How to download and work with SNDOAS data + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +
+ + + + + + + + +
+

1. Tutorial Overview

+

This notebook demonstrates how to download and unzip masked SNODAS files that are from after October 2013 and select one variable and compile it into a RasterStack for quick visualization and finally save as a .grd file.

+
+

Credits

+

This tutorial is based on an R script written in 2021 by Cara Thompson at New Mexico State University and has been converted to a tutorial with some minor additions by Jennifer Roebuck of NSIDC.

+

For questions regarding the tutorial or to report problems, please create a new issue in the NSIDC-Data-Tutorials repo

+
+
+

Learning Objectives

+

By the end of this demonstration you will be able to: 1. Download and unzip masked SNODAS files 2. Convert the files to .bil format and put them in a RasterStack 3. Quickly visualize the rasters of the data 4. Save the RasterStack as .grd file

+
+
+

Prerequisites

+
    +
  1. The following libraries should already be installed: R.utils, rvest, dplyr, raster, sp, maptools, and data.table.
    2. +
+
+
+

Time requirement

+

Allow approximately 20 minutes to complete this tutorial.

+
+
+
+

2. Tutorial steps

+
+

Load necessary libraries

+ + + +
library(R.utils)
+library(rvest)
+ + +
Loading required package: xml2
+ + +
library(dplyr)
+library(raster)
+library(sp)
+library(maptools)
+library(data.table)
+ + + +
+
+

Download SNODAS data

+

We will download all the SNODAS masked files for June 2023. If you wish to download masked files from a different time period you can modify the path in the url variable below.

+ + + +
# set the URL for the SNODAS directory we want to download files from
+url <- "https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun" 
+
+# read html content from url
+page <- read_html(url)
+
+# Get a list of the files listed in the directory at this url
+files <- page %>% html_nodes("a") %>% html_attr("href")
+
+for(i in 2:length(files)){
+  # generate the url for each of the files 
+  u <- paste(url,files[i], sep="/")
+  # download each of the files
+  download.file(u,files[i], mode = "wb")
+}
+ + +
trying URL 'https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun/SNODAS_20230601.tar'
+Content type 'application/octet-stream' length 5724160 bytes (5.5 MB)
+==================================================
+downloaded 5.5 MB
+
+trying URL 'https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun/SNODAS_20230602.tar'
+Content type 'application/octet-stream' length 6543360 bytes (6.2 MB)
+==================================================
+downloaded 6.2 MB
+
+trying URL 'https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun/SNODAS_20230603.tar'
+Content type 'application/octet-stream' length 6737920 bytes (6.4 MB)
+==================================================
+downloaded 6.4 MB
+
+trying URL 'https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun/SNODAS_20230604.tar'
+Content type 'application/octet-stream' length 6021120 bytes (5.7 MB)
+==================================================
+downloaded 5.7 MB
+
+trying URL 'https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun/SNODAS_20230605.tar'
+Content type 'application/octet-stream' length 5621760 bytes (5.4 MB)
+==================================================
+downloaded 5.4 MB
+
+trying URL 'https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun/SNODAS_20230606.tar'
+Content type 'application/octet-stream' length 6307840 bytes (6.0 MB)
+==================================================
+downloaded 6.0 MB
+
+trying URL 'https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun/SNODAS_20230607.tar'
+Content type 'application/octet-stream' length 5038080 bytes (4.8 MB)
+==================================================
+downloaded 4.8 MB
+
+trying URL 'https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun/SNODAS_20230608.tar'
+Content type 'application/octet-stream' length 5345280 bytes (5.1 MB)
+==================================================
+downloaded 5.1 MB
+
+trying URL 'https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun/SNODAS_20230609.tar'
+Content type 'application/octet-stream' length 5171200 bytes (4.9 MB)
+==================================================
+downloaded 4.9 MB
+
+trying URL 'https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun/SNODAS_20230610.tar'
+Content type 'application/octet-stream' length 5427200 bytes (5.2 MB)
+==================================================
+downloaded 5.2 MB
+
+trying URL 'https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun/SNODAS_20230611.tar'
+Content type 'application/octet-stream' length 5806080 bytes (5.5 MB)
+==================================================
+downloaded 5.5 MB
+
+trying URL 'https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun/SNODAS_20230612.tar'
+Content type 'application/octet-stream' length 5724160 bytes (5.5 MB)
+==================================================
+downloaded 5.5 MB
+
+trying URL 'https://noaadata.apps.nsidc.org/NOAA/G02158/masked/2023/06_Jun/SNODAS_20230613.tar'
+Content type 'application/octet-stream' length 5478400 bytes (5.2 MB)
+==================================================
+downloaded 5.2 MB
+ + + +

We will save the working directory path to a variable.

+ + + +
WD <- getwd()
+ + + +
+
+

Read in the SNODAS files

+

Next, we read in the .tar files as a list containing .tar

+ + + +
files_1 <- list.files(path = WD, pattern = "*.tar", full.names = T, recursive = T)
+ + + +

Decompress the .tar files into the working directory

+ + + +
asd <- sapply(files_1, untar)
+ + + +
+
+

Select the variable we wish to use

+

There are eight parameters from SNODAS, each has an associated ID code which is used the filename. A full description of the ID codes can be found in Table 1 of the SNODAS user guide. We will choose one of the parameters to compile all the files related to this parameter.

+

We can uncomment one of the variables below that corresponds to the parameter we are interested in. In this tutorial we will use the snow depth parameter as an example.

+ + + +
id_code = "us_ssmv11036tS" # Snow Depth
+#id_code = "us_ssmv11034tS" # SWE
+#id_code = "us_ssmv11044bS" # Snow melt runoff at base of snow pack
+#id_code = "us_ssmv11050IL00" # Sublimation from the snow pack
+#id_code = "us_ssmv11039IL00" # sublimation of blowing snow
+#id_code = "us_ssmv11038wS" # snow pack average temperature
+#id_code = "us_ssmv01025SIL01" # Solid precipitation 
+#id_code = "us_smsv01025SIL00" # Liquid precipitation
+ + + +

Now we will compile all the files that contain this id_code in the filename i.e. a .dat.gz file and .txt.gz file for each day of data that we downloaded.

+ + + +
files_2 <- list.files(path = WD, pattern = id_code, full.names = T, recursive = T)
+files_2
+ + +
 [1] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060105HP001.dat.gz"
+ [2] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060105HP001.txt.gz"
+ [3] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060205HP001.dat.gz"
+ [4] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060205HP001.txt.gz"
+ [5] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060305HP001.dat.gz"
+ [6] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060305HP001.txt.gz"
+ [7] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060405HP001.dat.gz"
+ [8] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060405HP001.txt.gz"
+ [9] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060505HP001.dat.gz"
+[10] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060505HP001.txt.gz"
+[11] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060605HP001.dat.gz"
+[12] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060605HP001.txt.gz"
+[13] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060705HP001.dat.gz"
+[14] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060705HP001.txt.gz"
+[15] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060805HP001.dat.gz"
+[16] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060805HP001.txt.gz"
+[17] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060905HP001.dat.gz"
+[18] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023060905HP001.txt.gz"
+[19] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023061005HP001.dat.gz"
+[20] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023061005HP001.txt.gz"
+[21] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023061105HP001.dat.gz"
+[22] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023061105HP001.txt.gz"
+[23] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023061205HP001.dat.gz"
+[24] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023061205HP001.txt.gz"
+[25] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023061305HP001.dat.gz"
+[26] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/us_ssmv11036tS__T0001TTNATS2023061305HP001.txt.gz"
+ + + +

Decompress all the .gz files.

+ + + +
dcompSD <- sapply(files_2, gunzip, remove=T)
+ + + +

We will create a new directory for these uncompressed files and move them over.

+ + + +
dir.create("Snow Depth")
+file.copy(dcompSD, "Snow Depth")
+ + +
 [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
+[16] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
+ + + +
+
+

Convert the .dat files to .bil

+

Next we need to convert the .dat files to .bil so they can be read as rasters.

+ + + +
files <- list.files(path=paste(WD,"/Snow Depth",sep=""), pattern="*.dat", full.names = T)
+newfiles <- gsub(".dat$", ".bil", files)
+file.rename(files, newfiles)
+ + +
 [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
+ + +
newfiles
+ + +
 [1] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060105HP001.bil"
+ [2] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060205HP001.bil"
+ [3] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060305HP001.bil"
+ [4] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060405HP001.bil"
+ [5] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060505HP001.bil"
+ [6] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060605HP001.bil"
+ [7] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060705HP001.bil"
+ [8] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060805HP001.bil"
+ [9] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060905HP001.bil"
+[10] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023061005HP001.bil"
+[11] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023061105HP001.bil"
+[12] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023061205HP001.bil"
+[13] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023061305HP001.bil"
+ + + +
+
+

Create a text file header for each file

+

Now we will create a corresponding text header file for each .bil file using the associated .txt file that was downloaded.

+

Please note that the parameters used in the m_post_2013 variable below are specific to SNODAS masked files for October 2013 onwards. If you wish to process masked files prior to October 2013 then you would need to change ulxmap and ulymap to the following: ulxmap -124.72958333333334 ulymap 52.870416666666664

+

If you wish to process unmasked files prior to October 2013 then you would need to change nrow, ncols, ulxmap and ulymap to the following: nrows 4096 ncols 8192 ulxmap -130.51291666666665 ulymap 58.228750000000005

+

If you wish to process unmasked files from October 2013 onwards then you would need to change the nrows, ncols, ulxmap and ulymap to the following: nrows 4096 ncols 8192 ulxmap -130.51250000000002 ulymap 58.229166666666664

+ + + +
#setwd("/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth")
+mylist<- list.files(path=paste(WD,"/Snow Depth",sep=""), pattern = "*.bil", full.names = T)
+x<-"SNODAS Snow Depth Description = 
+nrows 3351
+ncols 6935
+nbands 1
+nbits 16
+pixeltype signedint
+byteorder M
+layout dat
+ulxmap -124.72958333333334
+ulymap 52.870416666666664
+xdim 0.00833333333
+ydim 0.00833333333
+pixeltype SIGNEDINT
+nodata -9999"
+
+mylist.hdr<-paste(mylist,".hdr",sep="")
+
+for(i in mylist.hdr) {
+  write(x,i)
+}
+
+ + + +

Change the file type of the header text files from .bil.hdr to .hdr since the .bil file is already separate.

+ + + +
#setwd("/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth")
+files4 <- list.files(path=paste(WD,"/Snow Depth",sep=""), pattern="*.bil.hdr", full.names = T)
+newfiles <- gsub(".bil.hdr$", ".hdr", files4)
+file.rename(files4, newfiles)
+ + +
 [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
+ + +
newfiles
+ + +
 [1] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060105HP001.hdr"
+ [2] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060205HP001.hdr"
+ [3] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060305HP001.hdr"
+ [4] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060405HP001.hdr"
+ [5] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060505HP001.hdr"
+ [6] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060605HP001.hdr"
+ [7] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060705HP001.hdr"
+ [8] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060805HP001.hdr"
+ [9] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060905HP001.hdr"
+[10] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023061005HP001.hdr"
+[11] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023061105HP001.hdr"
+[12] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023061205HP001.hdr"
+[13] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023061305HP001.hdr"
+ + + +
+
+

Create a RasterStack of all the rasters

+

Create a RasterStack of all the snow depth raster files.

+ + + +
#create new directory to save the RasterStack in
+#dir.create("Rasters")
+#setwd("/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/Rasters")
+
+list.ras <- list.files(path=paste(WD,"/Snow Depth",sep=""), pattern = ".bil", full.names = T)
+list.ras
+ + +
 [1] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060105HP001.bil"
+ [2] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060205HP001.bil"
+ [3] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060305HP001.bil"
+ [4] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060405HP001.bil"
+ [5] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060505HP001.bil"
+ [6] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060605HP001.bil"
+ [7] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060705HP001.bil"
+ [8] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060805HP001.bil"
+ [9] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023060905HP001.bil"
+[10] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023061005HP001.bil"
+[11] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023061105HP001.bil"
+[12] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023061205HP001.bil"
+[13] "/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/us_ssmv11036tS__T0001TTNATS2023061305HP001.bil"
+ + +
ras_stack <- stack(list.ras)
+ras_stack
+ + +
class      : RasterStack 
+dimensions : 3351, 6935, 23239185, 13  (nrow, ncol, ncell, nlayers)
+resolution : 0.008333333, 0.008333333  (x, y)
+extent     : -124.7338, -66.94208, 24.94958, 52.87458  (xmin, xmax, ymin, ymax)
+crs        : NA 
+names      : us_ssmv11//60105HP001, us_ssmv11//60205HP001, us_ssmv11//60305HP001, us_ssmv11//60405HP001, us_ssmv11//60505HP001, us_ssmv11//60605HP001, us_ssmv11//60705HP001, us_ssmv11//60805HP001, us_ssmv11//60905HP001, us_ssmv11//61005HP001, us_ssmv11//61105HP001, us_ssmv11//61205HP001, us_ssmv11//61305HP001 
+min values :                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768 
+max values :                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767
+ + + +
+
+

Plot RasterStack

+

Now we define a projection for the RasterStack and plot each of the rasters (must be in the format of what it is in, not what you want it to be). Please note that SNODAS data are point estimates in latitude/longitude coordinates with the horizontal datum WGS 84 thus the decision to project the data and the choice of projection is up to the user. In this tutorial we will use Geographic longitude/latitude.

+ + + +
crs.geo <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")
+proj4string(ras_stack) <- crs.geo
+ras_stack
+ + +
class      : RasterStack 
+dimensions : 3351, 6935, 23239185, 13  (nrow, ncol, ncell, nlayers)
+resolution : 0.008333333, 0.008333333  (x, y)
+extent     : -124.7338, -66.94208, 24.94958, 52.87458  (xmin, xmax, ymin, ymax)
+crs        : +proj=longlat +datum=WGS84 +no_defs 
+names      : us_ssmv11//60105HP001, us_ssmv11//60205HP001, us_ssmv11//60305HP001, us_ssmv11//60405HP001, us_ssmv11//60505HP001, us_ssmv11//60605HP001, us_ssmv11//60705HP001, us_ssmv11//60805HP001, us_ssmv11//60905HP001, us_ssmv11//61005HP001, us_ssmv11//61105HP001, us_ssmv11//61205HP001, us_ssmv11//61305HP001 
+min values :                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768,                -32768 
+max values :                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767,                 32767
+ + +
plot(ras_stack)
+ + +

+ + + +
+
+

Save RasterStack

+

Save the RasterStack to a new directory. We will write in the “raster”/native format to preserve layer names (i.e. dates). There is no attribute table because there is only one attribute per file (snow depth in mm).

+ + + +
#setwd("/Users/jero7025/Documents/NOAA/Datasets/G02158/R_code/Snow Depth/Rasters")
+dir.create(paste(WD,"/Snow Depth","/Rasters",sep=""))
+fp <- paste(WD,"/Snow Depth","/Rasters",sep="")
+writeRaster(ras_stack,filename=file.path(fp,"ras_stack.grd"), format="raster", progress="text", options=c("INTERLEAVE=BAND","COMPRESS=LZW"))
+ + +

+  |                                                                            
+  |                                                                      |   0%
+  |                                                                            
+  |===                                                                   |   4%
+  |                                                                            
+  |======                                                                |   8%
+  |                                                                            
+  |========                                                              |  12%
+  |                                                                            
+  |===========                                                           |  16%
+  |                                                                            
+  |==============                                                        |  20%
+  |                                                                            
+  |=================                                                     |  24%
+  |                                                                            
+  |====================                                                  |  28%
+  |                                                                            
+  |======================                                                |  32%
+  |                                                                            
+  |=========================                                             |  36%
+  |                                                                            
+  |============================                                          |  40%
+  |                                                                            
+  |===============================                                       |  44%
+  |                                                                            
+  |==================================                                    |  48%
+  |                                                                            
+  |====================================                                  |  52%
+  |                                                                            
+  |=======================================                               |  56%
+  |                                                                            
+  |==========================================                            |  60%
+  |                                                                            
+  |=============================================                         |  64%
+  |                                                                            
+  |================================================                      |  68%
+  |                                                                            
+  |==================================================                    |  72%
+  |                                                                            
+  |=====================================================                 |  76%
+  |                                                                            
+  |========================================================              |  80%
+  |                                                                            
+  |===========================================================           |  84%
+  |                                                                            
+  |==============================================================        |  88%
+  |                                                                            
+  |================================================================      |  92%
+  |                                                                            
+  |===================================================================   |  96%
+  |                                                                            
+  |======================================================================| 100%
+ + + +
+
+
+

3. Learning outcomes recap

+

We have learned how to: 1. Download and unzip masked SNODAS files 2. Convert the files to .bil format and put them in a raster stack 3. Quickly visualize the rasters of the data 4. Save the raster stack as .grd file

+ +
+ +
LS0tCnRpdGxlOiAiSG93IHRvIGRvd25sb2FkIGFuZCB3b3JrIHdpdGggU05ET0FTIGRhdGEiCm91dHB1dDogaHRtbF9ub3RlYm9vawotLS0KCiMjIDEuIFR1dG9yaWFsIE92ZXJ2aWV3IApUaGlzIG5vdGVib29rIGRlbW9uc3RyYXRlcyBob3cgdG8gZG93bmxvYWQgYW5kIHVuemlwIG1hc2tlZCBTTk9EQVMgZmlsZXMgdGhhdCBhcmUgZnJvbSBhZnRlciBPY3RvYmVyIDIwMTMgYW5kIHNlbGVjdCBvbmUgdmFyaWFibGUgYW5kIGNvbXBpbGUgaXQgaW50byBhIFJhc3RlclN0YWNrIGZvciBxdWljayB2aXN1YWxpemF0aW9uIGFuZCBmaW5hbGx5IHNhdmUgYXMgYSAuZ3JkIGZpbGUuCgojIyMgQ3JlZGl0cwpUaGlzIHR1dG9yaWFsIGlzIGJhc2VkIG9uIGFuIFIgc2NyaXB0IHdyaXR0ZW4gaW4gMjAyMSBieSBDYXJhIFRob21wc29uIGF0IE5ldyBNZXhpY28gU3RhdGUgVW5pdmVyc2l0eSBhbmQgaGFzIGJlZW4gY29udmVydGVkIHRvIGEgdHV0b3JpYWwgd2l0aCBzb21lIG1pbm9yIGFkZGl0aW9ucyBieSBKZW5uaWZlciBSb2VidWNrIG9mIE5TSURDLgoKRm9yIHF1ZXN0aW9ucyByZWdhcmRpbmcgdGhlIHR1dG9yaWFsIG9yIHRvIHJlcG9ydCBwcm9ibGVtcywgcGxlYXNlIGNyZWF0ZSBhIG5ldyBpc3N1ZSBpbiB0aGUgW05TSURDLURhdGEtVHV0b3JpYWxzIHJlcG9dKGh0dHBzOi8vZ2l0aHViLmNvbS9uc2lkYy9OU0lEQy1EYXRhLVR1dG9yaWFscy9pc3N1ZXMpCgojIyMgTGVhcm5pbmcgT2JqZWN0aXZlcwoKQnkgdGhlIGVuZCBvZiB0aGlzIGRlbW9uc3RyYXRpb24geW91IHdpbGwgYmUgYWJsZSB0bzoKMS4gRG93bmxvYWQgYW5kIHVuemlwIG1hc2tlZCBTTk9EQVMgZmlsZXMKMi4gQ29udmVydCB0aGUgZmlsZXMgdG8gLmJpbCBmb3JtYXQgYW5kIHB1dCB0aGVtIGluIGEgUmFzdGVyU3RhY2sKMy4gUXVpY2tseSB2aXN1YWxpemUgdGhlIHJhc3RlcnMgb2YgdGhlIGRhdGEKNC4gU2F2ZSB0aGUgUmFzdGVyU3RhY2sgYXMgLmdyZCBmaWxlCgojIyMgUHJlcmVxdWlzaXRlcwoxLiBUaGUgZm9sbG93aW5nIGxpYnJhcmllcyBzaG91bGQgYWxyZWFkeSBiZSBpbnN0YWxsZWQ6IFIudXRpbHMsIHJ2ZXN0LCBkcGx5ciwgcmFzdGVyLCBzcCwgbWFwdG9vbHMsIGFuZCBkYXRhLnRhYmxlLgoKCiMjIyBUaW1lIHJlcXVpcmVtZW50CkFsbG93IGFwcHJveGltYXRlbHkgMjAgbWludXRlcyB0byBjb21wbGV0ZSB0aGlzIHR1dG9yaWFsLgoKIyMgMi4gVHV0b3JpYWwgc3RlcHMgCgojIyMgTG9hZCBuZWNlc3NhcnkgbGlicmFyaWVzCmBgYHtyfQpsaWJyYXJ5KFIudXRpbHMpCmxpYnJhcnkocnZlc3QpCmxpYnJhcnkoZHBseXIpCmxpYnJhcnkocmFzdGVyKQpsaWJyYXJ5KHNwKQpsaWJyYXJ5KG1hcHRvb2xzKQpsaWJyYXJ5KGRhdGEudGFibGUpCmBgYAoKIyMjIERvd25sb2FkIFNOT0RBUyBkYXRhCldlIHdpbGwgZG93bmxvYWQgYWxsIHRoZSBTTk9EQVMgbWFza2VkIGZpbGVzIGZvciBKdW5lIDIwMjMuIElmIHlvdSB3aXNoIHRvIGRvd25sb2FkIG1hc2tlZCBmaWxlcyBmcm9tIGEgZGlmZmVyZW50IHRpbWUgcGVyaW9kIHlvdSBjYW4gbW9kaWZ5IHRoZSBwYXRoIGluIHRoZSB1cmwgdmFyaWFibGUgYmVsb3cuIAoKYGBge3J9CiMgc2V0IHRoZSBVUkwgZm9yIHRoZSBTTk9EQVMgZGlyZWN0b3J5IHdlIHdhbnQgdG8gZG93bmxvYWQgZmlsZXMgZnJvbQp1cmwgPC0gImh0dHBzOi8vbm9hYWRhdGEuYXBwcy5uc2lkYy5vcmcvTk9BQS9HMDIxNTgvbWFza2VkLzIwMjMvMDZfSnVuIiAKCiMgcmVhZCBodG1sIGNvbnRlbnQgZnJvbSB1cmwKcGFnZSA8LSByZWFkX2h0bWwodXJsKQoKIyBHZXQgYSBsaXN0IG9mIHRoZSBmaWxlcyBsaXN0ZWQgaW4gdGhlIGRpcmVjdG9yeSBhdCB0aGlzIHVybApmaWxlcyA8LSBwYWdlICU+JSBodG1sX25vZGVzKCJhIikgJT4lIGh0bWxfYXR0cigiaHJlZiIpCgpmb3IoaSBpbiAyOmxlbmd0aChmaWxlcykpewogICMgZ2VuZXJhdGUgdGhlIHVybCBmb3IgZWFjaCBvZiB0aGUgZmlsZXMgCiAgdSA8LSBwYXN0ZSh1cmwsZmlsZXNbaV0sIHNlcD0iLyIpCiAgIyBkb3dubG9hZCBlYWNoIG9mIHRoZSBmaWxlcwogIGRvd25sb2FkLmZpbGUodSxmaWxlc1tpXSwgbW9kZSA9ICJ3YiIpCn0KYGBgCgpXZSB3aWxsIHNhdmUgdGhlIHdvcmtpbmcgZGlyZWN0b3J5IHBhdGggdG8gYSB2YXJpYWJsZS4KYGBge3J9CldEIDwtIGdldHdkKCkKYGBgCgojIyMgUmVhZCBpbiB0aGUgU05PREFTIGZpbGVzCk5leHQsIHdlIHJlYWQgaW4gdGhlIC50YXIgZmlsZXMgYXMgYSBsaXN0IGNvbnRhaW5pbmcgLnRhcgpgYGB7cn0KZmlsZXNfMSA8LSBsaXN0LmZpbGVzKHBhdGggPSBXRCwgcGF0dGVybiA9ICIqLnRhciIsIGZ1bGwubmFtZXMgPSBULCByZWN1cnNpdmUgPSBUKQpgYGAKCkRlY29tcHJlc3MgdGhlIC50YXIgZmlsZXMgaW50byB0aGUgd29ya2luZyBkaXJlY3RvcnkKYGBge3J9CmFzZCA8LSBzYXBwbHkoZmlsZXNfMSwgdW50YXIpCmBgYAoKIyMjIFNlbGVjdCB0aGUgdmFyaWFibGUgd2Ugd2lzaCB0byB1c2UKVGhlcmUgYXJlIGVpZ2h0IHBhcmFtZXRlcnMgZnJvbSBTTk9EQVMsIGVhY2ggaGFzIGFuIGFzc29jaWF0ZWQgSUQgY29kZSB3aGljaCBpcyB1c2VkIHRoZSBmaWxlbmFtZS4gQSBmdWxsIGRlc2NyaXB0aW9uIG9mIHRoZSBJRCBjb2RlcyBjYW4gYmUgZm91bmQgaW4gVGFibGUgMSBvZiB0aGUgW1NOT0RBUyB1c2VyIGd1aWRlXShodHRwczovL25zaWRjLm9yZy9zaXRlcy9kZWZhdWx0L2ZpbGVzL2cwMjE1OC12MDAxLXVzZXJndWlkZV8yXzEucGRmKS4gV2Ugd2lsbCBjaG9vc2Ugb25lIG9mIHRoZSBwYXJhbWV0ZXJzIHRvIGNvbXBpbGUgYWxsIHRoZSBmaWxlcyByZWxhdGVkIHRvIHRoaXMgcGFyYW1ldGVyLgoKV2UgY2FuIHVuY29tbWVudCBvbmUgb2YgdGhlIHZhcmlhYmxlcyBiZWxvdyB0aGF0IGNvcnJlc3BvbmRzIHRvIHRoZSBwYXJhbWV0ZXIgd2UgYXJlIGludGVyZXN0ZWQgaW4uIEluIHRoaXMgdHV0b3JpYWwgd2Ugd2lsbCB1c2UgdGhlIHNub3cgZGVwdGggcGFyYW1ldGVyIGFzIGFuIGV4YW1wbGUuIAoKYGBge3J9CmlkX2NvZGUgPSAidXNfc3NtdjExMDM2dFMiICMgU25vdyBEZXB0aAojaWRfY29kZSA9ICJ1c19zc212MTEwMzR0UyIgIyBTV0UKI2lkX2NvZGUgPSAidXNfc3NtdjExMDQ0YlMiICMgU25vdyBtZWx0IHJ1bm9mZiBhdCBiYXNlIG9mIHNub3cgcGFjawojaWRfY29kZSA9ICJ1c19zc212MTEwNTBJTDAwIiAjIFN1YmxpbWF0aW9uIGZyb20gdGhlIHNub3cgcGFjawojaWRfY29kZSA9ICJ1c19zc212MTEwMzlJTDAwIiAjIHN1YmxpbWF0aW9uIG9mIGJsb3dpbmcgc25vdwojaWRfY29kZSA9ICJ1c19zc212MTEwMzh3UyIgIyBzbm93IHBhY2sgYXZlcmFnZSB0ZW1wZXJhdHVyZQojaWRfY29kZSA9ICJ1c19zc212MDEwMjVTSUwwMSIgIyBTb2xpZCBwcmVjaXBpdGF0aW9uIAojaWRfY29kZSA9ICJ1c19zbXN2MDEwMjVTSUwwMCIgIyBMaXF1aWQgcHJlY2lwaXRhdGlvbgpgYGAKCk5vdyB3ZSB3aWxsIGNvbXBpbGUgYWxsIHRoZSBmaWxlcyB0aGF0IGNvbnRhaW4gdGhpcyBpZF9jb2RlIGluIHRoZSBmaWxlbmFtZSBpLmUuIGEgLmRhdC5neiBmaWxlIGFuZCAudHh0Lmd6IGZpbGUgZm9yIGVhY2ggZGF5IG9mIGRhdGEgdGhhdCB3ZSBkb3dubG9hZGVkLgoKYGBge3J9CmZpbGVzXzIgPC0gbGlzdC5maWxlcyhwYXRoID0gV0QsIHBhdHRlcm4gPSBpZF9jb2RlLCBmdWxsLm5hbWVzID0gVCwgcmVjdXJzaXZlID0gVCkKZmlsZXNfMgpgYGAKCkRlY29tcHJlc3MgYWxsIHRoZSAuZ3ogZmlsZXMuCmBgYHtyfQpkY29tcFNEIDwtIHNhcHBseShmaWxlc18yLCBndW56aXAsIHJlbW92ZT1UKQpgYGAKCldlIHdpbGwgY3JlYXRlIGEgbmV3IGRpcmVjdG9yeSBmb3IgdGhlc2UgdW5jb21wcmVzc2VkIGZpbGVzIGFuZCBtb3ZlIHRoZW0gb3Zlci4KYGBge3J9CmRpci5jcmVhdGUoIlNub3cgRGVwdGgiKQpmaWxlLmNvcHkoZGNvbXBTRCwgIlNub3cgRGVwdGgiKQpgYGAKCiMjIyBDb252ZXJ0IHRoZSAuZGF0IGZpbGVzIHRvIC5iaWwKTmV4dCB3ZSBuZWVkIHRvIGNvbnZlcnQgdGhlIC5kYXQgZmlsZXMgdG8gLmJpbCBzbyB0aGV5IGNhbiBiZSByZWFkIGFzIHJhc3RlcnMuCgpgYGB7cn0KZmlsZXMgPC0gbGlzdC5maWxlcyhwYXRoPXBhc3RlKFdELCIvU25vdyBEZXB0aCIsc2VwPSIiKSwgcGF0dGVybj0iKi5kYXQiLCBmdWxsLm5hbWVzID0gVCkKbmV3ZmlsZXMgPC0gZ3N1YigiLmRhdCQiLCAiLmJpbCIsIGZpbGVzKQpmaWxlLnJlbmFtZShmaWxlcywgbmV3ZmlsZXMpCm5ld2ZpbGVzCmBgYAoKIyMjIENyZWF0ZSBhIHRleHQgZmlsZSBoZWFkZXIgZm9yIGVhY2ggZmlsZQpOb3cgd2Ugd2lsbCBjcmVhdGUgYSBjb3JyZXNwb25kaW5nIHRleHQgaGVhZGVyIGZpbGUgZm9yIGVhY2ggLmJpbCBmaWxlIHVzaW5nIHRoZSBhc3NvY2lhdGVkIC50eHQgZmlsZSB0aGF0IHdhcyBkb3dubG9hZGVkLiAKCioqUGxlYXNlIG5vdGUqKiB0aGF0IHRoZSBwYXJhbWV0ZXJzIHVzZWQgaW4gdGhlIG1fcG9zdF8yMDEzIHZhcmlhYmxlIGJlbG93IGFyZSBzcGVjaWZpYyB0byBTTk9EQVMgbWFza2VkIGZpbGVzIGZvciBPY3RvYmVyIDIwMTMgb253YXJkcy4gSWYgeW91IHdpc2ggdG8gcHJvY2VzcyBtYXNrZWQgZmlsZXMgcHJpb3IgdG8gT2N0b2JlciAyMDEzIHRoZW4geW91IHdvdWxkIG5lZWQgdG8gY2hhbmdlIHVseG1hcCBhbmQgdWx5bWFwIHRvIHRoZSBmb2xsb3dpbmc6CnVseG1hcCAtMTI0LjcyOTU4MzMzMzMzMzM0CnVseW1hcCA1Mi44NzA0MTY2NjY2NjY2NjQKCklmIHlvdSB3aXNoIHRvIHByb2Nlc3MgdW5tYXNrZWQgZmlsZXMgcHJpb3IgdG8gT2N0b2JlciAyMDEzIHRoZW4geW91IHdvdWxkIG5lZWQgdG8gY2hhbmdlIG5yb3csIG5jb2xzLCB1bHhtYXAgYW5kIHVseW1hcCB0byB0aGUgZm9sbG93aW5nOgpucm93cyA0MDk2Cm5jb2xzIDgxOTIKdWx4bWFwIC0xMzAuNTEyOTE2NjY2NjY2NjUKdWx5bWFwIDU4LjIyODc1MDAwMDAwMDAwNQoKSWYgeW91IHdpc2ggdG8gcHJvY2VzcyB1bm1hc2tlZCBmaWxlcyBmcm9tIE9jdG9iZXIgMjAxMyBvbndhcmRzIHRoZW4geW91IHdvdWxkIG5lZWQgdG8gY2hhbmdlIHRoZSBucm93cywgbmNvbHMsIHVseG1hcCBhbmQgdWx5bWFwIHRvIHRoZSBmb2xsb3dpbmc6Cm5yb3dzIDQwOTYKbmNvbHMgODE5Mgp1bHhtYXAgLTEzMC41MTI1MDAwMDAwMDAwMgp1bHltYXAgNTguMjI5MTY2NjY2NjY2NjY0CmBgYHtyfQoKbXlsaXN0PC0gbGlzdC5maWxlcyhwYXRoPXBhc3RlKFdELCIvU25vdyBEZXB0aCIsc2VwPSIiKSwgcGF0dGVybiA9ICIqLmJpbCIsIGZ1bGwubmFtZXMgPSBUKQptX3Bvc3RfMjAxMzwtIlNOT0RBUyBTbm93IERlcHRoIERlc2NyaXB0aW9uID0gCm5yb3dzIDMzNTEKbmNvbHMgNjkzNQpuYmFuZHMgMQpuYml0cyAxNgpwaXhlbHR5cGUgc2lnbmVkaW50CmJ5dGVvcmRlciBNCmxheW91dCBkYXQKdWx4bWFwIC0xMjQuNzI5MTY2NjY2NjY2NjcKdWx5bWFwIDUyLjg3MDgzMzMzMzMzMzM0CnhkaW0gMC4wMDgzMzMzMzMzMwp5ZGltIDAuMDA4MzMzMzMzMzMKcGl4ZWx0eXBlIFNJR05FRElOVApub2RhdGEgLTk5OTkiCgoKbXlsaXN0LmhkcjwtcGFzdGUobXlsaXN0LCIuaGRyIixzZXA9IiIpCgpmb3IoaSBpbiBteWxpc3QuaGRyKSB7CiAgd3JpdGUobV9wb3N0XzIwMTMsaSkKfQoKYGBgCgpDaGFuZ2UgdGhlIGZpbGUgdHlwZSBvZiB0aGUgaGVhZGVyIHRleHQgZmlsZXMgZnJvbSAuYmlsLmhkciB0byAuaGRyIHNpbmNlIHRoZSAuYmlsIGZpbGUgaXMgYWxyZWFkeSBzZXBhcmF0ZS4KCmBgYHtyfQpmaWxlczQgPC0gbGlzdC5maWxlcyhwYXRoPXBhc3RlKFdELCIvU25vdyBEZXB0aCIsc2VwPSIiKSwgcGF0dGVybj0iKi5iaWwuaGRyIiwgZnVsbC5uYW1lcyA9IFQpCm5ld2ZpbGVzIDwtIGdzdWIoIi5iaWwuaGRyJCIsICIuaGRyIiwgZmlsZXM0KQpmaWxlLnJlbmFtZShmaWxlczQsIG5ld2ZpbGVzKQpuZXdmaWxlcwpgYGAKCiMjIyBDcmVhdGUgYSBSYXN0ZXJTdGFjayBvZiBhbGwgdGhlIHJhc3RlcnMKQ3JlYXRlIGEgUmFzdGVyU3RhY2sgb2YgYWxsIHRoZSBzbm93IGRlcHRoIHJhc3RlciBmaWxlcy4gCmBgYHtyfQpsaXN0LnJhcyA8LSBsaXN0LmZpbGVzKHBhdGg9cGFzdGUoV0QsIi9Tbm93IERlcHRoIixzZXA9IiIpLCBwYXR0ZXJuID0gIi5iaWwiLCBmdWxsLm5hbWVzID0gVCkKbGlzdC5yYXMKcmFzX3N0YWNrIDwtIHN0YWNrKGxpc3QucmFzKQpyYXNfc3RhY2sKYGBgCgojIyMgUGxvdCBSYXN0ZXJTdGFjawpOb3cgd2UgZGVmaW5lIGEgcHJvamVjdGlvbiBmb3IgdGhlIFJhc3RlclN0YWNrIGFuZCBwbG90IGVhY2ggb2YgdGhlIHJhc3RlcnMgKG11c3QgYmUgaW4gdGhlIGZvcm1hdCBvZiB3aGF0IGl0IGlzIGluLCBub3Qgd2hhdCB5b3Ugd2FudCBpdCB0byBiZSkuIFBsZWFzZSBub3RlIHRoYXQgU05PREFTIGRhdGEgYXJlIHBvaW50IGVzdGltYXRlcyBpbiBsYXRpdHVkZS9sb25naXR1ZGUgY29vcmRpbmF0ZXMgd2l0aCB0aGUgaG9yaXpvbnRhbCBkYXR1bSBXR1MgODQgdGh1cyB0aGUgZGVjaXNpb24gdG8gcHJvamVjdCB0aGUgZGF0YSBhbmQgdGhlIGNob2ljZSBvZiBwcm9qZWN0aW9uIGlzIHVwIHRvIHRoZSB1c2VyLiBJbiB0aGlzIHR1dG9yaWFsIHdlIHdpbGwgdXNlIEdlb2dyYXBoaWMgbG9uZ2l0dWRlL2xhdGl0dWRlLgoKYGBge3J9CmNycy5nZW8gPC0gQ1JTKCIrcHJvaj1sb25nbGF0ICtlbGxwcz1XR1M4NCArZGF0dW09V0dTODQiKQpwcm9qNHN0cmluZyhyYXNfc3RhY2spIDwtIGNycy5nZW8KcmFzX3N0YWNrCnBsb3QocmFzX3N0YWNrKQpgYGAKCiMjIyBTYXZlIFJhc3RlclN0YWNrClNhdmUgdGhlIFJhc3RlclN0YWNrIHRvIGEgbmV3IGRpcmVjdG9yeS4gV2Ugd2lsbCB3cml0ZSBpbiB0aGUgInJhc3RlciIvbmF0aXZlIGZvcm1hdCB0byBwcmVzZXJ2ZSBsYXllciBuYW1lcyAoaS5lLiBkYXRlcykuIFRoZXJlIGlzIG5vIGF0dHJpYnV0ZSB0YWJsZSBiZWNhdXNlIHRoZXJlIGlzIG9ubHkgb25lIGF0dHJpYnV0ZSBwZXIgZmlsZSAoc25vdyBkZXB0aCBpbiBtbSkuCgpgYGB7cn0KI2NyZWF0ZSBhIG5ldyBkaXJlY3RvcnkgUmFzdGVycwpkaXIuY3JlYXRlKHBhc3RlKFdELCIvU25vdyBEZXB0aCIsIi9SYXN0ZXJzIixzZXA9IiIpKQpmcCA8LSBwYXN0ZShXRCwiL1Nub3cgRGVwdGgiLCIvUmFzdGVycyIsc2VwPSIiKQojc2F2ZSB0aGUgUmFzdGVyU3RhY2sgYXMgYSAuZ3JkIGZpbGUKd3JpdGVSYXN0ZXIocmFzX3N0YWNrLGZpbGVuYW1lPWZpbGUucGF0aChmcCwicmFzX3N0YWNrLmdyZCIpLCBmb3JtYXQ9InJhc3RlciIsIHByb2dyZXNzPSJ0ZXh0Iiwgb3B0aW9ucz1jKCJJTlRFUkxFQVZFPUJBTkQiLCJDT01QUkVTUz1MWlciKSkKCmBgYAoKIyMgMy4gTGVhcm5pbmcgb3V0Y29tZXMgcmVjYXAKCldlIGhhdmUgbGVhcm5lZCBob3cgdG86CjEuIERvd25sb2FkIGFuZCB1bnppcCBtYXNrZWQgU05PREFTIGZpbGVzCjIuIENvbnZlcnQgdGhlIGZpbGVzIHRvIC5iaWwgZm9ybWF0IGFuZCBwdXQgdGhlbSBpbiBhIHJhc3RlciBzdGFjawozLiBRdWlja2x5IHZpc3VhbGl6ZSB0aGUgcmFzdGVycyBvZiB0aGUgZGF0YQo0LiBTYXZlIHRoZSByYXN0ZXIgc3RhY2sgYXMgLmdyZCBmaWxl
+ + + +
+ + + + + + + + + + + + + + + + diff --git a/notebooks/SNODAS/img/nsidc_logo.png b/notebooks/SNODAS/img/nsidc_logo.png new file mode 100644 index 0000000000000000000000000000000000000000..35a06e02d3e9d47b3a392d05a550364d8cfd8aaa GIT binary patch literal 17603 zcmXtAWmH?;5=Be!;O-XO-J!T!f#UA&?ry~$iWe>JE$&)U+})kx^78$7S?lITR+4i@ z_w1S67*%B%R3t(qC@3gYIoZ$ZkhAxH2jCOrp5L5j4swEXlaSK{KyLm3^JvI(L}yuD zHz+99o&OH#MJbO{$cqH-QabJ$PL}RoCaxAxUS3`-HjZ{~W+u)SEKaUgxfeo&P*CJh za-YRDz4OlUy|PHO2(E&x9DDrQ9x}|q$j9Nc6v48HLomg%T0k{bRk)h(oDPdYBk@SqZ{p2ffzCuC~o={!KE@wwd*HIMtrZr;R8IFRO#UZFz0^oL3*)9|IrA z`iQR}5E)K+4Wd zk`LCt4;(;-IuDE$^aylRYQ$_EwCJg4Ym=iD(dw%25ed?)AeY2qaN))_RFN!Gg>vCY z#NvKw}mwl+vQl%fNbfvn^+JTWe#Wwx9v*dwiqJWUhR};?(V>L%&Jewg*2H=-wxK zNJ}Caf1nRQN9GUZC-AboHkS(@o_wpEgI!=Ol27eeu{L*%IHvpcovb z<3uJSc|Qa)mF1!PBR>kM0wSXW&chsp9;T%W0_-HwBqD%{mDR=Y! zu@PK#teJ`HGCLV_q&~MG*qgS(kincKtz|=@x7w-O)|CIRYsKb*P)t6i7%Uk1veQ$Ot)jaeoht zFm<3X!o1BXGdWTYia9Za4x|qEen>5Bxb6!g=~Ew<46zGIg=Qy4xm~ieYrh@CvA3PZ z$w!ze7}Omh9xXj_<+(W7ikQ#!hG?p!wc&LOR;R~~)!|o%>=*35F$MIkCokSmY*}b> ztPeGVucQ2h0@3rUD@P=xCJJ_VodxQ3befs2eKAotwKL2m&1{u~-SBIa;g2JMO5-Dd z4{kUGdb<$4U^gh=D44aYS>&~sB;J_89sWx)&t7vD;eWC`++H`MvGZrgX9(XLmt-^? z`19Wa72s&Mnk5gfFBXLFP(A_iO%0#`k9Q%QK0f~zYTgMTChWR#t zmzMJDZQu^2JY!={0{yVi8K9l#!jzy@7PC47z4_|IVE*9(SqaWob#JYh8XM+rUB#E1 z-nh$~W11!%0Xxk6w*hpN-B)`xBM2)V2BG$_n2(_K!P>BHp31@f}byrimn zhKgEp6n(O3`F?&bD-Oi`P7KWJu>SWyBK$wjurN>EM~P2oB7hwyuXHFMYLHcC!t@~E z{>Yhzn9(#U`mb+&h&)n7zl+14G*0+-e?&|cs7@4I1RCl#r7tdCnkwim{%}yNqPKFU z9(87e!OyD#sD}6M*SpQl!9I?x63y}HYniW=m8gEoX|fQSr01k3i{#cNJEPO!3n)I647mS-rBcBm#POU zkUV@=q+pqZMC3nxPq3d2rD`f2-j5h=(>;jFa801d1+VUJ9sJDe=X7#1AL=CGpD|Dt zl$psLhM4;I=Gm!DL!%kwODD>5%jwJ!lZvY5Z;%ifD=WznYJqH3-&KP5b3)D}1X*d` z`UlhQRD$4j6w0(p^a&U(cO=!RM!~g~MSKwzhVNx*sM#I;>g!k0YI*-O9NhY*5w}>ulQ!0TUz379euw&v~e8C z&@3U+u(!JuU-W3qUr4oMyo#c{L#=3ky6VvIhWG+lQx9!6EM!mnA;5A3y%0aNEDg`J zrY-SsV|jF$cP(jC9}haV`y&1UHW?)PxN}{j(-y|do8WkF9tJi!ETw6Bnrj^D?x~IF z>2A&GO9qlwgze(Zl-6}e3md#aoX^yn+C19Q&EQj`LCW6@qWh~*dQ`4?t$NpWG(-KC z?7$$Pf&BE>)+{YOy`{Bep3AzVv@aIQssFCOs;mBlR;)IE$d~~I)da1m|+ElK7(uj@KKmnSh z2$E!W%#fDAC+VeC-od(QksaIjvI(0umzz0$|Ldcn?gkUOciA*Ap7U}iKI>?kWSuy< zx-moz`UYGx25m&FvHg0N*&5o!Jl^zV)$oN4h%)LJ zc>_ICEU+9XFldlzc-Z5PO|VY-t%zf##Rf&1<{89Y7i zvEoMDzS&`A6SkySGL)%=r*Vw+`F7sdD7aHASSgx5^&wnE;cb{I*zt-apW5$;jCqe+lIGEQz;iecSxCiDCnZi|A_yv&+IaF!-4a@o6nQanL zdwr-Y_aVOeENulzS#S^P!rOK~u#GkRd)3m&pWBKC*LBbSKkZ{v@<<(s=DD)UOC$7X z%8F9u(aJJ)8L!$Ev!~m08faG9R;<-t%}pjM2vy}J$G&sYt)ZqriL@3G%OlthmFvug ze%8s7KVfW5km1oWn)Pt7v?srbvGm#?ko6jA2R4((Q<7sRTrWg*AxfDM6!V3JDdjHF z!g;1boyD5IGfpb>IclWrY?PNqQyRcLhM719o>`p5CWMu96H@XkxlAAo6ajdo`p)a= zE*;j^`&_i&bD)EOTXDTkB6H!TeKH5e^!qfu>jQmkh3^+m-cHz@Afh!i!A-$k9Y_B= zi&$i#G;Hq_Qo5-1jJ^J>DvX5cL9GO5CS^WcD>k6Cd<-zsthBTATano4)>n#@Xax!e z1AU@P>JFXKY8`skFSLgzauMr0ZtRziY$*O51hDLU4U=+0iRF)_bgUmth1j?Ve z&Roj+3iJsIR!Tejt1K$*d_9ce?QbPpwCv??pXm2(dQ)Arn)0uHb$LUun$IRP)cs>ubf>e2u*AKd&c}ZK8h^V)1vRgE3M5 zq3remMYa7SW+}~k`cV&%#?_Etr$zT=JKQ>>)ZUP}lKC_Pu;4zXt)ef#xCL+s7yt#c z(S$5o&;>+zEuHfz?RW^1d^RIoz4vVDsLd@B%NjUT-Y(N)9XR!r--4kSge8NQc}alC zlrp|P>t)VJuQjE6{F{u2zv%1h+(D(`1D+an+zaa(v$-G=^t4(CF}^83fQOtWv06Fm<_m;CLlAmgcIcPsRP(s=~EES%z9Fg)#b?| zhY#P+!5;?~RvX5)4*YI(^E>s+5bHWSI@<5h(vxn{V~q7DL+mhC{js`Fjvr&Hk(Tz9 zN*k(SskZW!frFpF!Cl^lTvt{Mp%PGvrY>6wl%@ES`fqHnr>gFD*3ORaqn&2$rfm_} zDeK)V10Mt$ZHdvme5z(P=a3RH6`$6W3CH8I!8!wKZrvaWVD#tZ*0#*`b=b%^`-5*~emBo= z;?uvk12@`l1wpZoWob7`6oXMO95XKr#X$(zYY(rh_8aScV&vUPLmEJB4|^C~6`Svq zhvJz6ML)(jLMS%Y7jjTq5|Cm%uVL^9(^)*VMnh&V2(9=P;=Qh|8XcmsU0; zHM9s$a|EFPpH1!8S0@&J{TlP?WW;%=gUK^Cf?B|D%Xo2~*N_KHjSE23t zbePc*vZ+9yw3JDLw<}$rc7xo@e*R0+|KWm8*RvP9AkeEY8ZmQVx_JfZL>@UAAdiH! z1pbxJ-F0863dInIY+c1-#fDEyP$~OwZHqNqEcIHyTsYM-o)=@Q_it$ehtGlf*AfF5 zKZM9idV1oPnIh{v)Tw;A{VIs#4OlAoJG3>G{mMPhf0MKEI43LE-jhT#mN(cp!RNO} zPjs>zjxpA-`%SVHnoDL50)C3S%EIA0Lvkams$^!SCSF|n^Jk*G)v+{Hih*Ax69eRt_1*Gokc2nkHLQ*87@~zddFkgPfjPFwR9f{pGKU4z zr4j0>ALbWd2$!=`rO5HkJbz_US@9Mf3#T?^Slew!oYKlM84b*1cNh;PYm&$ZU z)5dyGU%tpH?!h2YhZj&!edIA^x&~5j$Q?S^mb z$lm(clsaY8ceZRg|b8DI9V5;?0cyoCCI~2O8oP)E` zh3{A^|JD$><3Cf{ujgiS9=I-L3l`--0MNmP`2F3U^vvfRp2FZ0Q!<4L?#6~+QA+5q z_(>B*Ft-uF+*@UR)4*Y0==M*}|JAsuqY;W1xthBDlbsB&MF1xU2U-qgJ#QIp(bVek z$Jxoxik{=;Y}LP0Op|g$#In4+n!j2qib?nro}Sb^%iBwcGcy^J10~oy7PfRWh|OGo zqJad|RhD|`)pALs#`li~FWry7JsW86PyszS*f7gS3Nfetw2z+{_(5mssP#AUE1|$q z+S{MYOMkDyFZSc>ZK`f#t85J46nznZ)W$U=W=*~PW>eQch?&alp*-RH)ef7Rf>Ijj zj@tP(j*#IhndwW*oT8aY0$waQm~+<$kvN|?c%d)Y zO|>QHL3|QLBT_{ZL#m)Wlhc*Yw1<{{4Rf3VjhrBjyu9N^O$z}T1#AjF_pekRE6usb z+!kpsNK4&qbeB2|NWX|@Rfojb`lB4yOCkxOAy&71y?9E0*U4>M4m3SHy5$t`$8+9z z+KZ)jjR=M%tcIWji6qfc<@r^-(QJSqxwHmx3XQK620ibU&isWE7Dsf3Q`^|j5w9f zo+iSiMgU2kPKXl1SYiEDYF?p&$OjQ`E4E#io(`W$xE@dw1`-M?#3{10w8yr=biml* zr^R7|K}`p);;GanMh zo7g)-*(*zzBVS~I{yn-|{P%FbwswkGB!KdFO)A_@)iA_E9|tgSdFR!RJarSDF25Jz zGu|aA*9c$>jnW*xX<4FJQPcOZL3e8>(?plJOeQ>h0zCmk#QFKsvpOibsX({ zLr2M?Yl)6XmZ65YhYPgc{rw!C>^XZC!@Zi?k%gO|%0^=ZDPCO8$r6}T%x=?p`gU-| zwvC>ct0lx>9Tu+tPibMt=&7lz1wmIQf~p5|HXJJ^L3*XL`uG7&nGOx?2X3Ce?#@6F z=}ko__C#t`c0a+Q!mC zN{D`N!8ViyJgI`M6adY-#`22+4J2$dd}hHf08zjpC_5T!_@Re|IDZ;>(r59Fgl#wQ zWI#f`;v>OzS)bg@KCO@0d9hU7jX1W7Xmw)}_h#hTdoN>2ZD2h^}CS;Lk3) z`Zf_sA4Hx}N`+z#BSSrHJYe=L(fh8}vebyw6OId<_IzS)H_8AiijB>A5nJH8qlOMg zS)vCs$)jX4in0b=+Vk2O$qo1reUAM|X5O8UQT_)D%iy6)aYf9DW-~}5%}#x%dD*%K zRQ|gLpp^;{`u7$mfL5j^-V4Wbb3<)zR#l_1X>HjK7j|rYLrX5bmIc2B*u+m;l+~Kg;qrfkmZDE;p!>oc@jY!yyQ3~ z!C`$DnHY-mgaIIB+Kl>`ZK$W5Ju`8#?l;18*|^H*>?&M6>6PlET`tk`y)xsRB;ZGm zBZGFNDR)bG9t1^s=b_1Wvxt2!qH>SU9XsPG+VT{@N0bo*?Bq~vV1(YkZ^xHDNGW;P zpC3)Rtff5(DOS^4&DZ_%nfjiSB7vlj=D9D|0WTS(yQHZwqZ-@_Eh>@}S6AUw)O(Hs z!z~dlrc4m$n*mG1xadHJ(*d<@Ekq>9ToM=db3>OV>(FXQ*%SA@~|CT_{AxgQ{tnZn&A2%(@R9e+>)6KVBF;Cqk1Sk%}G;A|c9E+hyv zK>7J1WT_gW;5)X5Rq&Wl&j%Fqx{A)kncZSh;5gFB>5rU-1==$+NX;!}x^yl-yt1|?I2wLajDnuaLPPYm zD{E0F9lR*UH5|FXlEF1!R(PEJq|7&dUl(Y(pta>9aFgigHADXEogA5V;c7 z*B53InzFs+Hl95Pa_By5;-3?JrLmZ}$~6`S3Upn}uhnLuK*-jL%LhkJSMTc~MlYp% z>cX?+fP_r9&{*Cvpsq4oW-jeOle^N%q*T%fgZ}GXEQf5lyHr|;szZ5>oA~?M;S^hT zXJowJGLJiT-Q0@kWZFq%o=AkL%R-xHLwhA!R|FU~9%}U9p7dlQ7Bs!`JrfPGNQl^c z%+~M2Q*E+s-ps|M8P&%td&40dnDr%zAmP5}QYF*Pwx@u%m(^-}!0a(3e)5^yQo&N0 z@H7*8I*?SJ6C4z^#GNsap*z~Nu0EGwrP4g1j+5=&l$?N%2h)q6Y~XFJ3dz0 zqBMrQ7Vu9v%iF4H>jtp)YwOUe@RQ?}tf!YH&>O^2+Lh()-=47SyaN>6BA~UbFp;#wYfzKOfS|YK&NjRl`bG;R8y+cO z;t+VjEc>AGZ;ahCzxMHt0qsC&KnU=^cuA*Z_rcaO6p@ifsKvd|;?5Cd>hUSXz*Fy2 zU+W#}w&zMES8KD51Xtb7%tMN)ukQWr*+M3CcTWqyQ$)-UCApBe2hUb4EiKBb*og+G zWl>ZJRIg*%&ChrEY`T!QDT-OXN@m^}V=e_aP^vi?xG%^fc{^4$`uh2B3-C2(ILBGI zs#94jIq4ocNckqO>U5BMYHv$1)0_{t>5iuD9v!Un#dJBzAwNK|tuc&UJp-a~_O=Fr z?5o52RoCV-G!48;`Y|WzH$jj7-YkF!%vQeAf?pgl50L`xA|zrqv{2!k+Fy{AR?hi; z>mbiEFp5|;ISdmfF_TvL6lSu55SguO)VDTl+uHhYAp8!ELJ3JM#9DH0OX*x_JcI5M z@u!{ChVV~#YUtuN12Uq;=0E& zY^uDzy)t7X8A=fB`6~&&$eK7;6~$=|-|*_V1;_rhUL+45B}mfKockmfty&} ze!ieyWyL3}mCPOW4_*ShotXc?ym3P(?-%RZ<10CsHm>_+kBVc4THMF8x>a?&i?uL2}i5z(yy(acXOw0BBu;wij^!S zX0tO8aG8eZ&5VE;dHDqZ_IGqXPsz%Qgg%zM{-I>4HLVS; zSZYKMa|=XCrzE%d&vFtEnrMTVuguq2UAuI9+?;l1j*n+tf!$>^Z(DeMlbuM^5G*=_ zug6ZkE-7BcG$)K86F)z zBSeD4GmvvxY1z!KdBAuHAJyBNdj~to_w_Z4pPT>f}{zp zRt6+Sq{yqvh2+wwTcC{Wa%rYhmMN0r3gALgE$FUOcUDJA(!%B$@h2dDZAXyQFTF@L z98#o+apjm``K>G4HkQUxVs)fEQT(za%?vAFLCUcmhmgrJGpFTNp#O-P7Qq#3_ROY? z-{O`FU#?8LxUb<#*dZ9-o7l781*dCji7HRl-kJYwvT$ScFi}U$ry*&4PQuhB9sR-m zKf3=*-D-i*{!fU;AVJu)N?IQby52k7a)+CNvOXs4EeB2Jsty!vKu0RtW~#gX!%I5c z(Lhsqf7=xm9{I0d#{MniH$}N-dQ($Z;dt?X2}WmiZsoG28iANC}o^4J0P)cdh%OD7-||+%Jf(QM!?q8Ab=)e2?NRX2uPn007dQNniZ$TIS{8kSQLIBPM); z>^NC}uB7Pxp@uv zS-W1aw!~9g6w5dxl<_m)l@HP{ssr66hJ#$Y5UX8GT_k1D-a>+cV`aCxiGe$ zzO2;tBTZBL)}SQBX4>=V-Q}nq-UbwFD|^Wj{}3sk;AqU{D#|#s_ik)Zgl=0;OM}cp zbN+b9J$x@*xiJj~zCJf{>RhE4T4k|Fk>a7{O6cWO-dz3@Y`1&5f749ReBl{S4jgR4jtbG zV0FhUtByfpY=mWAAh~vef+BaCYBu5oPJDnaC0@LFYmr#Ak!T0yHmTV`IANXe3wX@S zRsq>XkQhP&^zXWFu~#})jE&{LmB3@bK}dkdQn0z3n-!WXc>8|DsH~_tx-~=-Q;e$) zGk!SQj?gfZYPYmQxFkyjN)P*>+N6Cs;lH^^1go zP((;`{M3c8wD>VJRKk8Qkvf6rm+`eiw2D>CheQ=*F(Wo3r|$UC34jnK z<>rgxu%H#G8~LAhX%OxVMcQSa4{>_=stX8supslnZR=P2CiyovcB8ch6^6S^l_0fy^ZRT?RgY{gOnI_+6A3 zzVt~SdyLq=`UIcV^wdW;po6~Xb$pu^|6)aA$+e85&MbzQfgE%J*eZ6#v>(n#G;~mR>AEOs!jMA)5)Gq;)~35 z;V_OBah#Tk(&8~Oit3rit_Oc5%grBQY_Uc2;c6$RwV^+mo5>h=KbTn$E*N`cM>Kfd z_OnNFuFc}mEXAK``D2LF$vO1XVqS3n(eZSZN!l$~`R#O;GBIu$73>e6KOg|4cf1wJ zo~^&z3T(H_C$El`R+)PKJMe`)uUv9fKUJh_NM>-mCPS10ISYg z%yN%uJNmLcA$Ab42qY|H*|`bWxXV_|kbPg;mB0wK7AjapOZg|bnjJ9xZ|!+>fGD!Ing~c@uLm=$+*liMn<#2T{DR$7Fvd#n|vk@k*%)<93@@QN+zXRb{E7Z27RI zwf^NjwVzhvbl-=#!1VUpU_JJA(3J^u8{*;@S1*c6fCw!>_~qU0y-V}oC^+N{tJ|8{ z9=cDT*1?RNgX(5l^ht(ZYG2K=r?nlQIjJne_Ye*XCR4JpIB6?uW@SE{$@xce z>cG6jz31+#drI0o99F#8F|GrJg^ zu`Dri#h8(N_6M7dCKNONfX{+ZEM; z{Vm1xXH-!#IaMDUxtO4J!Sq>HSDz(QfkduEFv_ zcS2h0kX1RzELBzAe)G>$flfsXhZW3Z@i%^eVXL$BaPLH1>=6k*EijhDEY##V4abL5 zZfVE%{Niw2kW2w>uL3WB_y_qmueN^-Os05xXE{3;COKw;hT$HkgC;CLj*Z~;dyB2x zaK&{BtOi`wUy!m!pzV|c&)BCXKrsh1x5xe9$}4fa2_b`o^WV_p*R$&3rNEcJ@`kdG zH0+4tCkN21!#FWRe20;uG?|G$8id~WtxoX*61X)A9V6Szt?11GJl}m?1WJtx?*j(#tZW>x&(>P_mEEX2{rDiIs z8Dyw*++RgqG2i=P=+vgRX>-i*ayxzZag1*7mo^BGDD-jv>9^nFZU!8N#ypDL z&G0jc{u5-v4!5r%&Z_Y_Y4h+_hPAkKc0!g~o(*%%+HE-8wxDwxbJy;`@~tz2@86eg z$QT?}{v<$&2GxKfoS|&(sD2e!3;PKlEnEzVGFe9*D_?Z}dx{h%?M^*oNdY`+&WI*1J_MG{Oa9mLZx1@!vh!k{R%=?(t|P9CKwp~eAIUV3vQNWE zw?7%V+^ugcirvdq7&01q`f-Qu=)KiZV;-tWv~3M>U2!p`V=^z1YlfWi+iIu$i1l)| zmZDg?z48cexadwTdw-c_mPC8gy0V;Z#p z$rz0j_vSTlxKo8ZAzHbfyXbK7v86>AWG&90!l=Z>oS&29ym#icGO*IKoRU_jrdwHl zzg5dT?zc^IEbhJHzyI6)*coRfoe8zUYIlCx)^C5k%JXtJxH)3HbelI&!H2OSKw7R?wD5R*>0cqXPC2Sy~vMfbvrHfLKl7)%sl7EJCSiJ`}G}>NBR5>v@(C-z3->~>0jNxOqH?E zM^@Pqb6L6WwBRi7bNEu02|?a{5Bfc&?bOZc^bEp(>PG7dlPw@Nmr}cGH-H;)Ay(&s zUC>HnYbehvf3_{>)>PN{qDW3`BBRF(5mwAW{Le(BH5NLS3X@oZXddiu5=o3z4ognZ7z%HR?2szx|;t9Y3JHLsABZdMg#i{4M8nV92eUgF;b!=&Rl_QrGDAI z-9R~%B98i>5UDT)97`tF&((w~o55WMcJ9QChO|8|1OXB>(`)vdR&Z=rzA`0ywUdZd zYY_-3!jKg}^$jaVdPIo82oN%JYC@ew>rz`9uDbL-I+L}aO2h8?WR+*0N>LmE98e!` zv5yh*~;_LkKlp%_ne#0P5 zin8Vlt+RWI&1!>JuE3H~i8KFtx$J+{LN!2e;jeq9=j6_ZDq} zG2!t6iacEd-T9~C^=RWnnFGhPxH0AHgk{}Q24+FkAh1bteE7ZtW38;R(ieyQLNYBa zHOyoA5;);G))K{gHFW8&!+dY+onfP!fDsC^P_j>VKl44|cq(8^mAQxqFFYy2P5Ci+ zqOiCv$))3l{_HY%MsYHJ%Jr~7$23bRtL#pQB2R&*S+e(R)^RMgz3b$lW5%mN>?e;6 zpaK|x1S}5i+T+wOyS0Ry5*_LFY3X$u2A4?wOOGJvBKrYv2utp%lRSB|6 z1Xo=`lq84+uisID$I`8RwwQe)=8qYBX&MNvf}LQjZ;w(H0B?|Rn7yZ8x)*ym>3hm3 zpS7SfHB`ZB_9tG>|KP{a!u|+7MdJh^fJjul6qw;c=dpjBS!X9rDP~HFO(e>9IQSXf zD=Bn9Z6>Z2I=3thToH5~b+*<(*Yte?wKGvuD|{fNc^peVsRt>fep&n^5-=S&?>H!4 zmg9&Cm?+Hn4gw5iF!gvZn7?~0&Y#&|T8^EK#fMXRuPY}$p3Hb34y@Hu_B+bZrq^8O zJ1}7rsmmg7TD{-If+HOHwVx)Ca#6!mh*u{oQ5(0rjl7FeV0L_U22Px-|4NH|AqOl3b83~s7s7@?x%#MT8Wb63$ zss{7c?F#0*&Zh}vr`wJ4{4oxB34Mt#NYjsedh)_EC4bu0hV1#=rlsi_wV1h$Vz)Yx zRpbrvN(|qIPq2#0! zfj1JlR(-!Ki37r*<*g4wQfA`YbHZJYp3cWjzn;GLn=r9}7oiMu-dO1hQl$242;*;W zN}M4`PibmBnO(po`QPOd95g28oLwHAN%!|uw%Qy*JASuiN;cP@?_8hWrDf`~fgg;) zn+6_^gm}5Ml_e9U{<93Snkbdkbu70?;&mk!`v%^57#w{(Q161b47di{^x^8H^~|SGz-+S7n1qA4%~lFl!X)U7#(43DF=R>nytXfM^|yAjyu|C<#XE(vh(9z`)O2p*x~18_ zgY`<3(fAmZsG7jzxvRFMTyn+YCA`pY{l?TJi8^elX1G&ON4=V%q8y=9P&esN4~79O zJU z{>^Sm6oo4snL1cMXDXNw~{8>MS%xJs$3{bX&Q*l!kkjQwI`bwH5F6Bug_ z4^y);Ee#{EOB>Dpvzb!XthSXPzI183pgQF+Fh~Z#UdL*eq+$S9nl=MV_jZ*uIJzQQ>`CK0@D;p~~XlQ!Qu-~g2N)W#m_wo*1Gwf4LA;?lGt#lr_k;a@>o<{HVj zf%@|)dP2nw-W)jAw-=rV{_j<<0ZNx{(xT^)h7GL7URi_oS`@>t%agtyB;i|)tcCsg zo6j)4M7?UeZ!pKd^C+L0$$i)1Cv0qOxrW?$s6S_-HPm1K97wNbr(+p10qMsy5+LPS zf?U23_0RqI)lrg^Ud?-(ze=`q^wV&L@zMceDs;&&j}w=(bITK{KYt^X-d{pSL#Nj` z>xxgZ;-3DN@m`<2kP+HQnbpRXFcFk6%G(gLQ0!yL`mu)Szzv#3O3>2$`^UZJI%Wl} zxt05Jvm!&FC>a#kx{{b z*oS?Kt%}~22ryyC&g%@m{2y@K#*J`BUkPirrPk8@MjAWA|$33^VyNxlWDy*VWa}&3?bj zI$Ey;;Le7-$VTp8ASfSfxc>AuA3ep>06jxgl(f^fQJLPxVVu>0%r5{$iFn892W40C(EE(F% zY!sQ@RMGH0BieUdFB2RV0X%)Eb$bahhN9np>V12ACvTeR%qUpkIjy>xX9+iXpQ09tiUmL$9|9=sHb4yxMt9c`+X8W7?2v9G(yxlf(Qas~-&H?Uh z1RP~)rN4YqA4dprTR68r{O%Sg=5KzMWKG?7!8MeuxZIOtslR)Y{FJgj9U%Hnh?n=9 zY4>}E!xYJnoZ&qB1b&*ANNK(Ygr4og%%B)}PFiTTBlMedz4yV?pWOZTYx zli;Sg=ACd3nx@iRSsWGDK$az8_->cYBrKwp0-!{m&C?E_@yGNI@0R-30K5j#LkKh62i~Z}R`G z0>UXZ)36&$d{}!P!Qe5DGAVfdj|5gWr{U5x%*@AweBd|P*6A!3u%QgE2W+Ys_k1%H zmn-_MESX9{z*kSYhmX5c#63rNmL00zK!`&Dz!f1Gx5|83xphLD&f197Zd3^8$oKyYHA}^UU&`oFy8h zojY|jb}hvu@DIXFCu&v1^kmT6)<3-JI?rc>t|pP73O*ie(XvIuHedIBrhh-9$xrRx zC+P^q)%vaQzvkIhLR`$*JAA~*KaN*)ssBd7$$0q2X-V>!&uXn?0UqS} z+jGIy@pG#=tH$(abNrTca zH~h&swj0wr8|o}y${(ccl*K~*;Db#+%_<2F3fcbI^>m%m7M+=1U?J^(T3vg2RDc^$ z3v^vHB9O~7;$Km2_7B{N_fo_=(|J?YOigyS0%>X+x#IE?nerp~%K@(c9wl55`PTqG zeN$(Fen-^oetfr7{ zykzj4{pR~R2h5WMB4#p_-o^cO8q(_}tNZwbSKton1v*y&cThzvAD1vt$jkEJyHU98 zms&Da4DOMAzviSSV@QS(o0;AU$gzFP3dwpMa9%ps_u016`L~T`=;tK4@-la((0`uo z8N?;)B#-eOne(V{9l^#R9hmla;pWC*5uXl7kWD}bfaC(CAki%>akk(zSIR4{`uY;5 zEY{aYxzkG$$OQ&K|6Lb=bi;plUGs`Ue0^7lj016unP&Onp)w^cKUd7Yv0KqEs&u%( zR>fS?%zCs0DE&!%DXYs_g!SWje9#7lD0l{vZ{9Pu0{iJAA#t?~`4Qf^BQ3?RARu)j zD3JBKD9}sD^iNo)1EMYqAoT-1v=BnS)iwV&keEmonq@$ zBKnW$+$+X5mXtG0X_^&m4}!nm*Urbud9Ox?T}yD7RQUtQETQ$r=@K?{L}1(PuD|4+ z|Kf&6P>HkrrY)BhLLW&71Drd6JLqIZZZ9iwDx{Mh{|im-U|uY)a?2zC;^+4~ zL8Pxc4XFu^P8K4Q!z{RNG(YjEU?EmNR)P>Hz9PPN@frN|!c!872AhDhgb*KERPc$w ze+#B1T;uzq0;4&e^o_%llqr-m%!3mI3o(kNG`8l?C zHtpu;4?NHFue^oy4`w+L(wqkAh#nV=#~r#KWlcpd412Ee<=b zq7u$nzKmbmk_Vvg25tfVB!t-aZ{B~*@Lxx80@4uuW8gqAPO|&Q+Z%b{`8RoW z<-2I1fRVdo3F&UZ^z}MnHzzpW^giXRpH>=s*55Weir2@&c{BOuStoGTu?ON#9XP%T zW1anrxv%Sg0{&+ZoB;D1PRoGvfX@OxFurN@+&de2_@8gH^8NJ$hlUN#&Wmu%+DRff z5yThyls60_usiS+6><2|`CRn*a|hv0OOQhI>p zz`4MQ@PFDnW0i)XC=7pBqjnIfwJ`+SqNtSWR>f7QZ|M{G5V{HCP+TfzvKWjsMIugZ z65H!=BE+C#7j2*)-05aS5YPV6$q|+~48&aF3(L7R}vSHhs1YoC0S`sR_JFf6r74oJMPH11HI8 z0k!n~Fp7Eb16)7A^L4S!S47YrveLj&N(IcAkr& Date: Wed, 16 Aug 2023 13:10:12 -0600 Subject: [PATCH 2/4] Add README file --- notebooks/SNODAS/README.md | 13 +++++++++++++ 1 file changed, 13 insertions(+) create mode 100644 notebooks/SNODAS/README.md diff --git a/notebooks/SNODAS/README.md b/notebooks/SNODAS/README.md new file mode 100644 index 0000000..5996590 --- /dev/null +++ b/notebooks/SNODAS/README.md @@ -0,0 +1,13 @@ +## Download and Process SNODAS data + +## Summary + +In this tutorial we demonstrate how to download SNODAS data from the NOAA@NSIDC HTTPS File System and select one variable and compile it into a RasterStack for quick visualization and save as a .grd file. + +### Set up + +### Key Learning Objectives +1. Download and unzip masked SNODAS data files +2. Convert the files to .bil format and put them in a RasterStack +3. Quickly visualize rasters of the data +4. Save the RasterStack as a .grd file From dcfcaa1176ae6e076375eef1be69773890ed9f1f Mon Sep 17 00:00:00 2001 From: jennifer Date: Wed, 16 Aug 2023 13:12:13 -0600 Subject: [PATCH 3/4] Update README --- notebooks/SNODAS/README.md | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/notebooks/SNODAS/README.md b/notebooks/SNODAS/README.md index 5996590..025070b 100644 --- a/notebooks/SNODAS/README.md +++ b/notebooks/SNODAS/README.md @@ -2,7 +2,7 @@ ## Summary -In this tutorial we demonstrate how to download SNODAS data from the NOAA@NSIDC HTTPS File System and select one variable and compile it into a RasterStack for quick visualization and save as a .grd file. +In this tutorial we demonstrate how to download SNODAS data from the NOAA@NSIDC HTTPS File System and select one variable and compile it into a RasterStack for quick visualization and save as a .grd file. Please note that the tutorial demonstrate processing masked files from after October 2013. If you wish to process any masked files prior to October 2013 or any unmasked files you will need to modify the tutorial slightly, details on how to do this are included in the tutorial. ### Set up From 37a5d97c9ec2b4a2d3e7a2aff77f5b5cdb29ae46 Mon Sep 17 00:00:00 2001 From: jennifer Date: Wed, 16 Aug 2023 13:18:42 -0600 Subject: [PATCH 4/4] Add credits to README file --- notebooks/SNODAS/README.md | 2 ++ 1 file changed, 2 insertions(+) diff --git a/notebooks/SNODAS/README.md b/notebooks/SNODAS/README.md index 025070b..2171bde 100644 --- a/notebooks/SNODAS/README.md +++ b/notebooks/SNODAS/README.md @@ -2,6 +2,8 @@ ## Summary +**Credits** This tutorial is based on code provided by Cara Thompson at New Mexico State University. It has been modifed slightly and converted into a tutorial by NSIDC. + In this tutorial we demonstrate how to download SNODAS data from the NOAA@NSIDC HTTPS File System and select one variable and compile it into a RasterStack for quick visualization and save as a .grd file. Please note that the tutorial demonstrate processing masked files from after October 2013. If you wish to process any masked files prior to October 2013 or any unmasked files you will need to modify the tutorial slightly, details on how to do this are included in the tutorial. ### Set up