{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Coursework Support\n", "\n", "
\n", "\n", "### Imports" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "library(RNetCDF)\n", "library(tidync)\n", "library(ncmeta)\n", "library(tidyverse)\n", "library(lubridate)\n", "library(ggplot2)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "### Loading Data\n", "\n", "We'll begin by loading the NetCDF file using `tidync`. When we do so we retrieve both the values of the data-array it contains as well as the metadata that describes it. " ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "Data Source (1): era5.nc ...\n", "\n", "Grids (4) : \n", "\n", "[1] D0,D1,D2 : u10, v10, ssrd **ACTIVE GRID** ( 745488 values per variable)\n", "[2] D0 : longitude\n", "[3] D1 : latitude\n", "[4] D2 : time\n", "\n", "Dimensions 3 (all active): \n", " \n", " dim name length min max start count dmin dmax unlim coord_dim \n", " \n", "1 D0 longi~ 501 90 1.40e2 1 501 9.00e1 1.40e2 FALSE TRUE \n", "2 D1 latit~ 186 -12 6.50e0 1 186 -1.20e1 6.50e0 FALSE TRUE \n", "3 D2 time 8 1043144 1.05e6 1 8 1.04e6 1.05e6 FALSE TRUE " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "nc_filepath = '../data/Coursework-Support/era5.nc'\n", "\n", "nc_file = tidync(nc_filepath)\n", "\n", "nc_file" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "We can see that there are 3 dimemensions in the file, calling `hyper_dims` lets us inspect them in more detail.\n", "\n", "We can see that longitude, latitude and time are not 'unlimited', this means that the extent of the data is fixed. 'coord_dim' as True denotes that the dimension is being used is attached to the dataset, in this case they all are but think about examples where this maybe different (e.g. providing multiple versions of coordinates for different CRS')." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\t\n", "\t\n", "\t\n", "\n", "
namelengthstartcountidunlimcoord_dim
longitude501 1 501 0 FALSE TRUE
latitude 186 1 186 1 FALSE TRUE
time 8 1 8 2 FALSE TRUE
\n" ], "text/latex": [ "\\begin{tabular}{r|lllllll}\n", " name & length & start & count & id & unlim & coord\\_dim\\\\\n", "\\hline\n", "\t longitude & 501 & 1 & 501 & 0 & FALSE & TRUE \\\\\n", "\t latitude & 186 & 1 & 186 & 1 & FALSE & TRUE \\\\\n", "\t time & 8 & 1 & 8 & 2 & FALSE & TRUE \\\\\n", "\\end{tabular}\n" ], "text/markdown": [ "\n", "| name | length | start | count | id | unlim | coord_dim |\n", "|---|---|---|---|---|---|---|\n", "| longitude | 501 | 1 | 501 | 0 | FALSE | TRUE |\n", "| latitude | 186 | 1 | 186 | 1 | FALSE | TRUE |\n", "| time | 8 | 1 | 8 | 2 | FALSE | TRUE |\n", "\n" ], "text/plain": [ " name length start count id unlim coord_dim\n", "1 longitude 501 1 501 0 FALSE TRUE \n", "2 latitude 186 1 186 1 FALSE TRUE \n", "3 time 8 1 8 2 FALSE TRUE " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "hyper_dims(nc_file)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "We can see that there are 4 grids and that 3 of them are the dimensions we just explored, what's left over? - the variables. \n", "\n", "Here we can see that there are 3 variables, two for the u and v components of a wind vector, the other for surface solar radiation. From these we estimate solar and wind power production at different locations in Indonesia." ] }, { "cell_type": "code", "execution_count": 32, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\t\n", "\t\n", "\t\n", "\n", "
idnametypendimsnattsdim_coord
3 u10 NC_SHORT3 6 FALSE
4 v10 NC_SHORT3 6 FALSE
5 ssrd NC_SHORT3 7 FALSE
\n" ], "text/latex": [ "\\begin{tabular}{r|llllll}\n", " id & name & type & ndims & natts & dim\\_coord\\\\\n", "\\hline\n", "\t 3 & u10 & NC\\_SHORT & 3 & 6 & FALSE \\\\\n", "\t 4 & v10 & NC\\_SHORT & 3 & 6 & FALSE \\\\\n", "\t 5 & ssrd & NC\\_SHORT & 3 & 7 & FALSE \\\\\n", "\\end{tabular}\n" ], "text/markdown": [ "\n", "| id | name | type | ndims | natts | dim_coord |\n", "|---|---|---|---|---|---|\n", "| 3 | u10 | NC_SHORT | 3 | 6 | FALSE |\n", "| 4 | v10 | NC_SHORT | 3 | 6 | FALSE |\n", "| 5 | ssrd | NC_SHORT | 3 | 7 | FALSE |\n", "\n" ], "text/plain": [ " id name type ndims natts dim_coord\n", "1 3 u10 NC_SHORT 3 6 FALSE \n", "2 4 v10 NC_SHORT 3 6 FALSE \n", "3 5 ssrd NC_SHORT 3 7 FALSE " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "hyper_vars(nc_file)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "This is all very useful but what about if we wanted to store additional, more general metadata about the dataset?\n", "\n", "Luckily NetCDFs allow us to store further information as 'attributes' to both the top-level of the file as well as for individual variables. By specifying `'NC_GLOBAL'` as the variable we can list only the attributes for the top-level of the file, this includes information such as when the file was downloaded." ] }, { "cell_type": "code", "execution_count": 36, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\t\n", "\t\n", "\n", "
idnamevariablevalue
0 ConventionsNC_GLOBAL CF-1.6
1 history NC_GLOBAL 2020-12-10 14:47:04 GMT by grib_to_netcdf-2.16.0: /opt/ecmwf/eccodes/bin/grib_to_netcdf -S param -o /cache/data0/adaptor.mars.internal-1607611614.7236085-5982-17-ac6b25d3-2c5a-4a80-8c51-bece0fdfbc32.nc /cache/tmp/ac6b25d3-2c5a-4a80-8c51-bece0fdfbc32-adaptor.mars.internal-1607611614.7241852-5982-5-tmp.grib
\n" ], "text/latex": [ "\\begin{tabular}{r|llll}\n", " id & name & variable & value\\\\\n", "\\hline\n", "\t 0 & Conventions & NC\\_GLOBAL & CF-1.6 \\\\\n", "\t 1 & history & NC\\_GLOBAL & 2020-12-10 14:47:04 GMT by grib\\_to\\_netcdf-2.16.0: /opt/ecmwf/eccodes/bin/grib\\_to\\_netcdf -S param -o /cache/data0/adaptor.mars.internal-1607611614.7236085-5982-17-ac6b25d3-2c5a-4a80-8c51-bece0fdfbc32.nc /cache/tmp/ac6b25d3-2c5a-4a80-8c51-bece0fdfbc32-adaptor.mars.internal-1607611614.7241852-5982-5-tmp.grib\\\\\n", "\\end{tabular}\n" ], "text/markdown": [ "\n", "| id | name | variable | value |\n", "|---|---|---|---|\n", "| 0 | Conventions | NC_GLOBAL | CF-1.6 |\n", "| 1 | history | NC_GLOBAL | 2020-12-10 14:47:04 GMT by grib_to_netcdf-2.16.0: /opt/ecmwf/eccodes/bin/grib_to_netcdf -S param -o /cache/data0/adaptor.mars.internal-1607611614.7236085-5982-17-ac6b25d3-2c5a-4a80-8c51-bece0fdfbc32.nc /cache/tmp/ac6b25d3-2c5a-4a80-8c51-bece0fdfbc32-adaptor.mars.internal-1607611614.7241852-5982-5-tmp.grib |\n", "\n" ], "text/plain": [ " id name variable \n", "1 0 Conventions NC_GLOBAL\n", "2 1 history NC_GLOBAL\n", " value \n", "1 CF-1.6 \n", "2 2020-12-10 14:47:04 GMT by grib_to_netcdf-2.16.0: /opt/ecmwf/eccodes/bin/grib_to_netcdf -S param -o /cache/data0/adaptor.mars.internal-1607611614.7236085-5982-17-ac6b25d3-2c5a-4a80-8c51-bece0fdfbc32.nc /cache/tmp/ac6b25d3-2c5a-4a80-8c51-bece0fdfbc32-adaptor.mars.internal-1607611614.7241852-5982-5-tmp.grib" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ncmeta::nc_atts(nc_filepath, 'NC_GLOBAL')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "To access the attributes for a specific variable we just pass that variable's name to `nc_atts`. Doing this for `'u10'` provides us information such as the units and the value for missing values. \n", "\n", "The `scale_factor` and `add_offset` enable intelligent methods for compressing data, these are automatically applied to the data on loading and saving to redue the size of the stored data, you can find out more [here](https://www.unidata.ucar.edu/software/netcdf/workshops/2010/bestpractices/Packing.html)." ] }, { "cell_type": "code", "execution_count": 37, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\n", "
idnamevariablevalue
0 scale_factor u10 0.0003157431
1 add_offset u10 1.074366
2 _FillValue u10 -32767
3 missing_valueu10 -32767
4 units u10 m s**-1
5 long_name u10 10 metre U wind component
\n" ], "text/latex": [ "\\begin{tabular}{r|llll}\n", " id & name & variable & value\\\\\n", "\\hline\n", "\t 0 & scale\\_factor & u10 & 0.0003157431 \\\\\n", "\t 1 & add\\_offset & u10 & 1.074366 \\\\\n", "\t 2 & \\_FillValue & u10 & -32767 \\\\\n", "\t 3 & missing\\_value & u10 & -32767 \\\\\n", "\t 4 & units & u10 & m s**-1 \\\\\n", "\t 5 & long\\_name & u10 & 10 metre U wind component\\\\\n", "\\end{tabular}\n" ], "text/markdown": [ "\n", "| id | name | variable | value |\n", "|---|---|---|---|\n", "| 0 | scale_factor | u10 | 0.0003157431 |\n", "| 1 | add_offset | u10 | 1.074366 |\n", "| 2 | _FillValue | u10 | -32767 |\n", "| 3 | missing_value | u10 | -32767 |\n", "| 4 | units | u10 | m s**-1 |\n", "| 5 | long_name | u10 | 10 metre U wind component |\n", "\n" ], "text/plain": [ " id name variable value \n", "1 0 scale_factor u10 0.0003157431 \n", "2 1 add_offset u10 1.074366 \n", "3 2 _FillValue u10 -32767 \n", "4 3 missing_value u10 -32767 \n", "5 4 units u10 m s**-1 \n", "6 5 long_name u10 10 metre U wind component" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ncmeta::nc_atts(nc_filepath, 'u10')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "Let's dig a little deeper into the dimensions we have, what does 'time' look like? \n", "\n", "Before we can access the values we need to activate the correct grid, which currently is the one containing 'u10', 'v10', and 'ssrd'. We'll switch over to 'time' by activating it, then extract the values by calling `hyper_tibble`." ] }, { "cell_type": "code", "execution_count": 71, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
    \n", "\t
  1. 1043144
    2. \n", "\t
  2. 1043148
    3. \n", "\t
  3. 1043152
    4. \n", "\t
  4. 1043156
    5. \n", "\t
  5. 1047488
    6. \n", "\t
  6. 1047492
    7. \n", "\t
  7. 1047496
    8. \n", "\t
  8. 1047500
    9. \n", "
\n" ], "text/latex": [ "\\begin{enumerate*}\n", "\\item 1043144\n", "\\item 1043148\n", "\\item 1043152\n", "\\item 1043156\n", "\\item 1047488\n", "\\item 1047492\n", "\\item 1047496\n", "\\item 1047500\n", "\\end{enumerate*}\n" ], "text/markdown": [ "1. 1043144\n", "2. 1043148\n", "3. 1043152\n", "4. 1043156\n", "5. 1047488\n", "6. 1047492\n", "7. 1047496\n", "8. 1047500\n", "\n", "\n" ], "text/plain": [ "[1] 1043144 1043148 1043152 1043156 1047488 1047492 1047496 1047500" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "times <- (nc_file %>% activate(time) %>% hyper_tibble)$time\n", "\n", "times" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "Hmm these don't like like any dates we're used to working with, what's going on?\n", "\n", "Let's inspect the attributes of the time variable to find out more." ] }, { "cell_type": "code", "execution_count": 51, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\t\n", "\t\n", "\t\n", "\n", "
idnamevariablevalue
0 units time hours since 1900-01-01 00:00:00.0
1 long_nametime time
2 calendar time gregorian
\n" ], "text/latex": [ "\\begin{tabular}{r|llll}\n", " id & name & variable & value\\\\\n", "\\hline\n", "\t 0 & units & time & hours since 1900-01-01 00:00:00.0\\\\\n", "\t 1 & long\\_name & time & time \\\\\n", "\t 2 & calendar & time & gregorian\\\\\n", "\\end{tabular}\n" ], "text/markdown": [ "\n", "| id | name | variable | value |\n", "|---|---|---|---|\n", "| 0 | units | time | hours since 1900-01-01 00:00:00.0 |\n", "| 1 | long_name | time | time |\n", "| 2 | calendar | time | gregorian |\n", "\n" ], "text/plain": [ " id name variable value \n", "1 0 units time hours since 1900-01-01 00:00:00.0\n", "2 1 long_name time time \n", "3 2 calendar time gregorian " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ncmeta::nc_atts(nc_filepath, 'time')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "Lets create a function that maps from the time dimension to the relevant date string.\n", "\n", "We can now see that we have data for 4 different hours of the day, from two separate days in the year." ] }, { "cell_type": "code", "execution_count": 69, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[1] \"2019-01-01 08:00:00 UTC\" \"2019-01-01 12:00:00 UTC\"\n", "[3] \"2019-01-01 16:00:00 UTC\" \"2019-01-01 20:00:00 UTC\"\n", "[5] \"2019-07-01 08:00:00 UTC\" \"2019-07-01 12:00:00 UTC\"\n", "[7] \"2019-07-01 16:00:00 UTC\" \"2019-07-01 20:00:00 UTC\"" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "hours_to_datetime <- function(hour) {\n", " start_date <- ymd('1900-01-01')\n", " datetime <- start_date + dhours(hour)\n", " return(datetime)\n", "}\n", "\n", "datetime <- hours_to_datetime(times)\n", "\n", "datetime" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "Now we'll start looking at the actual data, to make things easier lets take a slice of it, filtering for only the first hour we have data for. To do this we can use the `hyper_filter` function which will filter over multiple dimensions if passed. We also have to make sure to reactivate the relevant data grid, in this case the solar radiation, using `activate`. Finally we convert the tinync object into a hyper-array using `hyper_array`." ] }, { "cell_type": "code", "execution_count": 142, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Class: tidync_data (list of tidync data arrays)\n", "Variables (1): 'ssrd'\n", "Dimension (2): longitude,latitude,time (501, 186)\n", "Source: C:/Users/Spatial/Downloads/ESDA-Spatial-master/ESDA-Spatial-master/tutorials/era5.nc" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "radiation_slice <- nc_file %>% \n", " activate(ssrd) %>% \n", " hyper_filter(time=time<1043145) %>% \n", " hyper_array\n", "\n", "radiation_slice" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "We're now in a position where we can start visualising the data!" ] }, { "cell_type": "code", "execution_count": 152, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAMAAADKOT/pAAAAVFBMVEUAAABNTU1oaGh8fHx9\nACWMjIyampqiBwanp6eysrK9vb3DIgDHx8fQ0NDZ2dnhPADh4eHp6entYgDw8PDyhAD1oQD3\nujz40HT75Jr/9Lf//8j///+x2akoAAAACXBIWXMAABJ0AAASdAHeZh94AAAgAElEQVR4nO3d\niVbiyhpAYXBo9dA2ijjdvP97XpLKUGPGvyoV2Hute66CSQrIZ5uBsCuIaHG7tQdAdA0BiUgg\nIBEJBCQigYBEJBCQiAQCEpFAQCISCEhEAgGJSCAgEQkEJCKBgEQkEJCIBAISkUBAIhIISEQC\nAYlIICARCQQkIoGARCQQkIgEAhKRQEAiEghIRAIBiUggIBEJBCQigYBEJBCQiAQCEpFAQCIS\nCEhEAgGJSCAgEQkEJCKBgEQkEJCIBAISkUBAIhIISEQCAYlIICARCQQkIoGARCQQkIgEAhKR\nQEAiEghIRAIBiUggIBEJBCQigYBEJBCQiAQCEpFAQCISCEhEAgGJSCAgEQkEJCKBgEQkEJCI\nBAISkUBAIhIISEQCAYlIICARCQQkIoGARCQQkIgEAhKRQEAiEghIRAIBiUggIBEJBCQigYBE\nJBCQiAQCEpFAQCISCEhEAgGJSCAgEQkEJCKBgEQkEJCIBJoP6d/L467s8fmf4HiINtlcSB93\nu6570SERba+5kJ53+79v1Vfvr/vds9yAiLbYXEj73Vv79dtuLzMYoq02F9JuF/qG6AbjXyQi\ngRZsI72+V1+xjUQ0f/f3vbbX7u5DckhE22vBcaTn6jjS/vGF40h087GbgEggIBEJlADSjmhj\nzVjLZahEXwRR2a/qR/Vd9eX/yQVLyQfSQt5E/q4b0uqLoFsJSEQCAYlIIM1RrehS6EdnLwVI\ndOUBiUggIBENpW/4lLk/oTlSNwQQLSwhpAkHsIBE47pFSP8BiaS7RUjF237sJU+AROP66QpQ\nahhFHkjSbaS3sW/nAxKN6zYhXf66exv+oWWLoJvqRiFltAi6ioC09iLoOhqUlMYRkGjbAWnl\nRdB1BKSVF0HX0UVOAFJ9vs+vwhZ7HECiTQeklRdB1xGQVl4EXUlKj36ekLr9cktzYtAMR5Mn\nABJtOyCtuwi6koC07iLoSgLSuougq+mrO3TUOlJ2yv9G382gAhJtPSCtuQi6moC05iLoatLf\n8todRmoQAYloVEBacxF0NQFpzUUQSQYkIoGARCQQkIgEAhKRQEAiEghIRAIBiWJ1PK49goQB\niWIFJPFJMlwERQ9I4pNkuAiK3vHS6bT2KBIFJIoVkMQnyXARFD0giU+S4SIoekASnyTDRVCK\njjezwwFIFDEgyU6S4SIoRUCSnSTDRZB8h8PBvfFGJAGJxAJS7EkyXATJB6TYk2S4CJLPC+lG\nzhQCEokFpNiTZLgIkg9IsSfJcBEkH5BiT5LhIogkAxKRQEAiEghIFLnDoTwN/No3lIBEkQOS\n3CQZLoJSBSS5STJcBC3pcOnp0sgfP9ZFHdO6AYlmBCQ7INGMgGQHJJoRkOyARNM71JAexk9S\nQaqm851GtP2ARNMDkhOQaHpAcgISTe3QOHp6eBhJqdlIOgJp0SQZLoJmByRPQKKpAckTkGha\nh8Nh+ibS8XjtlIBE0wKSNyDRtIDkDUg0LSB5AxJN6XDQII09Q8iAdKWnCgGJpgSkQECiKQEp\nEJBoTCUbBejYQJoyeann1HaNloBEYwLSQECiMQFpICDRcA8PlSNtF/bkWWiOlKQrowQkGg5I\ngwGJhgPSYECi4YA0GJBoOLWnYdmBoBLQWVV+KT7GlQMSDQekwYBEwwFpMCDRqBSkcv2fu3HT\nOjpf47EkINGogNQfkGhUQOoPSDSxBeu/6QhI0QNSxgHJF5BoYkDyBSSaGJB8AYnS1jJqrgR+\nHZdCARKlDUgLJslwEbRSQFowSYaLoBXTHNWaNh+QKH1AmjlJhougFQPSzEkyXARZNVcnSbOs\nhtGiE8pzCkikAtKigEQqIC0KSKQC0qKAdNuVa/CxaBg1Z/Acx2//f35+lc1YdKuouybKjLnk\nEpBuOyAJBaTbDkhCAenW0j5mz9g8qRR91pTqD9ZTh3gO7gpeoitKRcrR93d54/f3JFAHDdH5\nlG4DLU5AurWAFCUg3VpAihKQbqKHS/WX+sexdJDUBlIJo7RU3mftn66nuMCq1/9z/fMKUlP5\nXXnriDHVis5d7eWMBR95qlJCev+z278UxX93u/1zpEWQPyDFLiGkj/3u0n8v5X9391EWQYGA\nFLuEkJ53l3+Hnve7Px/FR/W1/CIoEJBilxDSvppwt/uo/m8fYxHkqfzUym4Pw4MH0mUjv4NR\nUrLur1ftQ3fbqZNUO/q5VEsqvkYeWjIW+qldF0XokacsIaTdrvtv83/CiyBPQErRCv8ilf/9\n4F+kZAEpRStsIz1/1F/LL4KsypX+cIHU3VJ+9l75KS3t3QpSWbtOq59srlHSbSKp/ztqGzf1\nRpK11O/vojRVfVF/VdW/7VRSas4V2p4l9tpddUBKFceRrjogpYozG66pes91XbNSHg41pHIn\neC+k6rSf+suTkbmc5sSeylJgLD8/5X+qL34qS2rDSd8E85flyULtr4DQph+QrikgxQpINxWQ\nYpUtJI4jxQhIsdoQpJ2exCKuvfKlPVcr/7loPjHlrG64vOzdeTfl/oaiVtRCaufiPyOnnpFa\nof1r9eDK/vv7W/xc/ldyqlRVlgp1ClG5k2LWo16pZsxd1g/wp91mA1LCgHS9ASlhQLrC9POs\nWy+f/pSkQ3meUHn+d+2ob+7tvJYO87dC5GSvid++H6rzbqAF+/rum9fSLErVE9SeKl+ck0L6\n9/JYbQE9Pv+LtYhbCEj+bgXSx522N4FThOYHJH+3Aul5t//7Vn31/rrnpNXZaedZf+lmvrTX\nWN2gKJW7wA8tnqenMZesk4D0v//9b+ksRmyJacVk1ELSazYli7SQ9ru39us33kYxOyAFuhVI\nxtEhDsjODkiBbgUS/yLJBKRAtwLpso30+l59xTbSvKq3oTYpSy4ju1JS3m9KmLePwEPsO8Xe\nBr3mUN2llLu/77W9dncfURZx3QFJ64YhFf+eq+NI+8cXjiPNCUhatwwpp0VssK9vm5EfknmL\nOpS09tj7K9f/sZfDq/MfX4pJyXOwrjstC0gbCkhaQMpkERsMSFpAymQR+eaejl3uYtWvH+fb\n/W0pKtpb1TVFVnkko2sGPvCWv3EzkhiQP+dULO2UeyBlF5CWzUhiQP6AtKmAtGxGEgPyB6RN\nBaRlM5IYkD8gbaT6YiXHY7Xel3y6F6pU8a1nOmok2avj5ZaVHsykyrW/WTHnz6V+HsRGZee8\nXwVIeQYkIMkGpAJI0wPSCovIr+pwhNZZT60ihiKDkVqH1Iu79gOZ3yJG9RyqLcRYWZKAlGNA\nApJ0QALS3DkAKe0ickqhaF4Sn6P2Y4+dE4E+zdd2w4oWVx41OC/czBrI3kjqXgQgZRCQRALS\nGovIKSCJBKQ1FpFTQBIJSGssIoe0k31aE97Uq+VjBCSt6tyq2J8Joz3dZ+1FANKKAUk2IK2w\niBwCkmxAWmERq2e91/XTpOFxZF/f5NMt6vbBVkrzcX+fxksDpPUCUqSAlHQRqwekSAEp6SJW\n7+f397eh5Nvu6dlACmIqzxcvqlXpmN/nGSfL+xmEojm/44C0XkCKFZBSLmL1gBQrIKVcxOoB\nKVZASrmI9futqvc3WDYsRxYf7a0UGqbqc86Nd2K0ux5i7YGI+0Eqc4vuqDA+RF4FpNUCUqSA\nlHIR6wekSAEp5SLWrfzAbwOSbsmQZDj6djL9Fcb7bJtZrP1gry+HEZBWCkibDki5BKRNB6Rc\n+u3qThTqXGiUKkblR7qUVdN+ez7dpZ6ovFtN1uOouivx472uPI6AtE5A2nJAyiYgbTkgZROQ\nthyQ1kyt/wrP/1Q6JXP3Qd8lDn+6+VQTFrkez8m4ZW+jBdKaASmjgLTdgJRRQNpw3+qM70aR\nLcmg1DsjbQPrMmmi0V9Xi6+QAqT1AlI+AWnDASmfgLTV6l3XpqOW0o9NqXdeaiYpRn21nWpJ\nM88Tdy8pfQJSmoCUVUDaakDKKiBtNSBlFZC2Xw+kztKIWaQZ7bWm3mgyE5LnY3fY2ZA8IOUQ\nkDYfkHIISJvvfz5JGqWfQUdqJgmG2l/EDw+P36mRNJ2S9f5jIK0UkHIISJsPSDkEpM33v15J\na49uQpuG1FwhZvqUfkdASh6QcghImw9IOQSkzQekHALS5gtAKi8Zua02DklRmjGp1xGQkgek\nHALS5gNSDgFp8+mI1PfbM1S1eUglpelTehkBKX1AyiEgbT4g5RCQtl8J6H/dKT4ZnOwzr01D\nKipLMzeRPI6AlD4g5RGQNh6Q8ghIGw9IeQSkLfU5cLHHDbcVSKGjRedZkgL7GoAUOSCtHpCu\nISCtHpC2X/URLfoNkY8YDV0RT7StQAo1B9LxGIDEde2iBqSMA9J2AlLGAWk7tYzqy5ls8STv\nzRf45Gn1adcT51VC8m0hnYAUNSBlEJC2H5AyCEhbr772cHsN4o1dk8FXddm9cT96OMQdyoQ+\n/Rt00yEFtpDKW4EULyBlEpC2HZAyCUjbDkiZBKRtVwFqP4Wv/GrjkLRPnimK/iNWx+OMq2rH\nzDvcyfsaOkifZZ0jIEUMSBkFpO0GpIwCUsYde1YX/bOVtbZMyfxQQd+aqZ4Odbp1+dxktJ20\n/Nyp49F21FAC0sKAZAWkKQGpCUhWQJoSkJr6tgN+TEnVeeBfzY7wbWb+ZvD9RAfpdDoqSRWl\np6Tj9Lf0HZZHBUlX1EgC0sKAZAWkKQGpCUhWQJoSkJqAZAWkKQFpRPpmea2olrT2yBY0wKj9\nzdLtbcjrWNKCjsqRzajb3wCkOAEJSEMBaURAAtJQQBrM2TjSJK09tgVpj8b/AzUc9THHV8To\nYG0gda9niwlIMQISkIYD0mBAAtJwQOrPJ0h79tce3qK6Vafvonancuf31UA6HDpHpiL91U0J\n6eN5f/nvy91ud/830iJyCEhAGtXctfx9v9sVH5f/lN1HWUQWAQlIo5q7lv/ZPX5c/vPn/WLq\nz+45xiKyCEhAGtXctXy3+6j/c/krb7ePsYgc0hDZRxw+Nw+pOWwSvGxIfRypOoS0dUiHqvoI\nUrWnoXkNjXOl1GudFNLlP/ud9o34InIISEAa1/w/7d6K4qX8T/kvUu9GEpAyDUgZQHrb7Z/f\nisf9RdLr3e41xiJWT50W1Bo6m2+n3Dqk+lof6nRn9+72nKDjsYGUzztkp2coao8hOY7qUu7+\nfq332JW9xFnE2gGp/QJIQy1Yy//+uSsVPb68R1vEugGp/QJIQ3FmQ7DqjG/F6GzWSVp7jIs6\n1Xl2bddrXf1d+ZXaxkg8QpG6raNTp6h8EV1F9RVhfoAkGZCAND4gBQMSkMYnspZf53EkIAFp\nfJEg7fQkFpG+ZldD6ejUpVHaOqRjl3XPwYRU37RdRo6i5iBSZpBWX0SMgASkCQEpFJCANCEg\neeo+wEVtIJ0sSFdxHKnQJQXvaO/ZHKRytAdz+8jZ0HWPI60A6d/LY7UF9Pj8L9Yi1gpIQJra\n3LX8407bm3Btb+wDEpCmNnctf97t/1anfhfvr/tre2Nfo6jZ9e2DdAVv7GvWs0MfpOZnVhng\nrA5m6mF4Pna5exWbWkoJIe3VOyiq3q7tjX1AAtLUFr2xz/uN2CLWC0hAmhr/InkCEpCmtmAb\n6VW9feL6tpF0R4qReRjiSnY16Nvi9j0upM1Q8ihqLgfpkdSlY0q5+/te22t39xFlEWsFJCBN\nbsFxpOfqONL+8eXajiMBCUiT48wGq/YTXPQNJIfRVVzXrnvbqOfCJp5tpBw+qq8r9K5eJ+PU\noCrzFo8nIC0PSCogTQtIVkBSAWlaQLKqIZUbSO1zbl0yQ9ttutowe68yPLaDe5J3W/0BzMdu\n++jpksBC4zXk6NNwpB/PsDecgLQ8IKmANC0gWQFJBaRpAckKSCogTQtIWiWi398aknZukHYS\nvsFoVUhlxk2hK3iH69Y16476UTeQlgxUvvrwXndDNT710TN+Rc3bmjtH1v02KyAtC0h1QJoc\nkLSAVAekyQFJq3J0kVS/M7Z7xv0bSGtKag57aDedp0rSVibzjmataiVlRanefNVuMU8G8jgx\ntoK6uwPsTicgLQtIKiBNHxWQtICkAtL0UQFJy4B07l6Q5sVwGK0H6ayf2tLc5Dlrpi9t5bFm\nfTYhzR7lxAF10/VMWD7r1a+5+vvL8DpILQX9a/v0oKN+Trt1nrgKSMsCUj1rIE0OSFpAqmcN\npMkBSQtI9ayBNDkgadVXfq7fG9usSd2LkZWj7rDWZazlbY6I/sxVzZqxBmnuEE/NIRtzqcXg\nIAd+wIIUoKD/CtQlVY/KZWTNBUjLAhKQgCQQkIAEpOXVFzXTzg9yr3mSCaTibDcZUne46HRu\nj8po61sNaeb41Fybzazuds9RKyu1/vf8QA3pXP9sLyTTUvOY+iBVAWlRQAISkAQCEpCAtDTt\noz2cz3LRt0eiQxo+z8f5w19JOo2T1PyQ+ejqZbuQDrP2frcznzQHbfXv+akG0uVhqJ/1UNK/\nORmpO59UQJIPSAWQgLQ8IBVAAtLygFQACUjLa6/37XyWy1l3ZEj6FLkEiVWzWp+td0W0TPTV\nQl/xqzNyxkG6/LA1A23hJqQ5D8FebV1IBvl2IcrAAKSvDtJJf6eHL8eS+umnNiCJB6Rm4UAC\n0oKA1CwcSECaWUtIvUaf2rqkX/PEn+xQ6tW6WZ21e9T6YDByGrMEE9AQpFlnfpvrbX2j+v9j\nUXspCmubqNu66V9m96tOewtzS6D5KcfSsUP3BKRIAQlIQBIISEAC0uS+7e81R5+f2mvUrXA9\njnogTb7OnJrKXJnLm47WRduWOPJDMu/VVtHB9Tq4BDULc922H4ingQUav+xcSIW2MFtS+4N+\nSPotQBoTkIAEJIGABCQgCQQkIAFJoB/re5ORsRbbx5DCkOyLM5ZzONs3jsqApLbC3RViviTz\nN4TFqGWmQZrxEKr5FKfQgagFkL7qJ717iTRI5sfOeB0NQ6puBNKYgAQkIAkEJCABaWk/Px5H\n34FDSGMcKUjWGtlNPjigk3a+THsAqRlFu444K0W3Fk2EVLgXBDGG03Vcdk0777J7NpHqXxq9\n02uMDEc+SGXOu2YtR8ZxJO1mIA0GJCABSSAgAQlIAnkhdZtIhqNuVetXVJ+xbe9EPtkrqVV7\nknezPp8NSN267KPUrUPm9tLwM9BMctbT7tccRfm4PmOTr+i2mAYnbAnZijQgzrI8kEKnCAFp\nSkACEpAEAhKQgCQQkIAEpOV5HOnvn/CcHdQnqaheXXddLLxrqN3Zqp6ptq60L/0YSGqawWeg\nnoc2SfVVc3dHctTKnazu+XQQtYzcqQYhme+jANLogAQkIAkEJCABSaBvB5Lzztizx1GAUlFC\n8pNppuoZix+nfpJSw6h7wXU8JqRhuKpmNie7Ztyao6w+zLx6TsovTubRI9Xldg8k8726R+35\n7N9MAtJQQAISkAQCEpCAtLjvb/v7r+AWknN5EK8jfUPHmPWIlToEqRtGCJK5ckyDVFbNyIbU\nPAxt5qOf2MTZikpI3m2koqVk0XvyWwLSuICkAhKQFgUkFZCAtCggqYAEJNEsSPp65EJyJFWz\nmLL+WqnLiAtDGvVWCi+kbg7ZQ6r2mBTVNScbR0FIhfZcafCeeikBaVpAAhKQBAISkIC0OOtq\ndtYmkgvJoVQscvTlZo7DOO6hQ+rG2H6rIRiz7HJOvZDUsuc8rEJd3zty7dhaSOG8kPopAWlK\nQAISkAQCEpCANLXyEifOLc6+77O9bjmqjFksYOSD9OlAal5bnZGxddN8JwipvabdvIcVOcfM\nECMNUg+jDlM1CZCCAUkPSECaGZD0gASkmQFJD0hAmpnjqLDWXw8j/6qmzWIJJO3j0T3HkU4G\npKAjY2TtN8ML7zuO1Cx71qMqi0WwWsvrp0J9X1SHj54Gd4qUz4jvHRQeSGoCIAUDkh6QgDQz\nIOkBCUgzu6yv7i3WFlK9QgU/EE9fS7W1buaIyqV6FBmfReeBFEQ+WlHVAKSZj0nNesnEobwf\nQF6v/9X9o+fTU/tTQAoGJDMgAWlWQDIDEpCE0h11jI49GyLtGnYyNkpmLt8VZEKyd9j2j8uG\nPpR+npELaXDq5v/9H8g3/kmwpwxO64V0sAEMNgJRGZAmBCQgAUkgIAEJSAIBCUhAEujThGQ6\nctfZ6oMZ60lnrbr24vsceSD1OdIOgI1eifVHZp9nVN7bOyMNUvm8lE/Bqf3uqL6bVjsS/931\n2h6g1D9nfffHCERlQJoQkIAEJIGABCQgjc8+fFSnVlprE+nopK9s5ou8DJL5btuzXbNIy5Er\nSPv5epqxIziZvzy0+U0QqT9LgS2m8fVdb/wpJKnfkXYS/ZMxn15GQPIFJH9AAtKkgOQPSECa\nUoXFe7u+GutaetJnUH7r3Dg6A5KPkX25Do8jZwq1ITeykzalO++RD6yDNPyzYy6lEl5qH6SA\npPKJbSE1Nz4MIioDkhOQAgEJSFMCUiAgAWlKQAoEJCCN7qzWV8895lo8nZE6EWY8pLPzPtoA\nI3Ndril166vf0Xni+l+PyHzoDiS1pjY/r96eWowD4andXTB90oeHdvV3IA3sbegbSs9PAMkM\nSL0jAlIoIJkBqXdEQAoFJKtz4HPFP521cJCR9Ed8W5BOVu7yA4YcfiOXr0/qLqND3F1k23NR\nudENbpUEe6jqhzR1roOTAMkKSMGA1BOQrIAUDEg9AcmsWlN8d3xa66AXTlPzjezYgo4CjoOI\nPApHPTXGspulFO0ube05qCYo7zS3miY2dX1/0PNAasY3i+dAQDIDUs9TA6RwQDIDUs9TA6Rw\nQDIDUs9TA6Rwq0DaDc1iLUhqrfPcoR280Vai9nCNyahLaFDdIGxIAUKDjuZAKs7GUbQaUZln\n4d3NhfgvlEAVjzAk7VUCUuyA1BeQekoIaWcWYxFLA1JfQOopIaR/+y1A8q9T7taB+fL4IUkN\n6lwPwTpFqU/RBEjV92MuPnJuF18vvb3Ht/jqlKjj6INUnh4eJvxwyUNhMhh5IVWUemY1Z7Ap\n/7T7eNzdv1dzyPhPOyAFA1JPabeR/u52fwsgTR0UkIa7LUjF+/3u8WN7kOxTvnUtPZCktpHK\nFf0cgmR7thj5HSlMFq3BkWhLr620acvsJCWs0lF4L1VivlLOC3N25zNn+yn5XruX3f4VSNMG\nBaTBbg9S8XY3sKdh+SJmB6TegBRujeNIf4A0bVBAGuwmIWWxCG/DkFwyAUiSg6pWdAeSudCQ\no1aLdijMqtuDMDCSHkiFuVDteUxzEKk5F8kDqaU0DtKFEZCWBqTegBRuLUjbOiALpOZpAFKg\nfCCNPu1hhc7aumt60dZdeUJl5QppOtIO5Lgf4mJL0qx0Ci1EzR3Dkronwr2rdaVZKg8kiT8j\nbtqz7oWkH0vqfY1mbh+V8afdmIBkPxHuXUBKMEmGi5gUkOwnwr0LSAkmyXARUzr7V11n9Y3g\nqNlCciHZG2ZHOw3S54iahzj0TJz8jowht5COyyA9NGcJ9Z8spOHwQzK2kwIvUnnz0/yTwpNC\n+vfyWG0BPT7/i7WIKAFJeyaA5C0hpI87bW/CfZRFRApI2jMBJG8JIT3v9n/fqq/eX/e75xiL\niBSQtGcCSN4SQtrv3tqv33b7GIuIlPvuiRAk+WWfCt/ugkmQ+vh8dR/UXiPp3+Fwubuc96iR\nL4FUyikVqf+f8n4KVc3HeLNsOY5j8MDWwl+CSd9qHvpGbBGRAlI3HCD541+kEQGpGw6Q/KXd\nRnqt3mm+sW0k8+99r6P2jGfztRBx5dlEasfTLtWL+2SfXPTlEdRBaoQOjGeso+aH5z4JD42f\nGYzKnvQroVSWZg5kXCl3f99re+3uPqIsIkZAMgKSt7THkZ6r40j7x5ctHUcCkhGQvHFmQ1No\nDdLW2mo11fal1uussVbrky77q1st3N197b+knYeSvolkmGm//P7+dih5L9hsPqjgo5J8N18H\naeYMDEhzZzI6IDUBCUgLAlITkIC0ICA1AQlICwJSfzajg3W+Sb3C6qt1N/FySO4ZPs3pSs7l\nId0dDhqkL7fvLlPSwJD6IQ3vrBhXu/aPNeAdVSJEZUDqD0h2QPIGpP6AZAckb0Dq8qwC2hqr\nOeqHJHiikOfkbeMEoZCkdlxhR35Il4bG1AdJrmblH+dAvSLq2t9RhxUOSF1AAtLsgNQFJCDN\nDkhdLiTbkfsusXqF7VZnY/pqnVuy2nn2ffd9pIsB6TQaUrsTfAyk0Mk2MXgFt3CMhdUjql8b\n+VGMCkhdQALS7IDUBSQgzQ5IXUAC0uyA1OQ7ANIPqVl9Q4eRBLZ9bUfDH3up72uo91T4GPkh\nfS2AFKeHEQurRhTns8rHB6QmIAFpQUBqAhKQFgSkvgKQzE0ke0XWpl/8wmrbOH5Ip5P3fKF2\nC6mDpB8vcizpwobG9JTkQFL1xthpMtZTVACpPyB5ApIvIPUFJE9A8gUkrXIF7L4z11EXkp+R\n7Bk0ClK3AWNBOuk5jvSd3x0XBdKQ9DXFUSpI9cJW3OyZFJC0gASkuQFJC0hAmhuQtIAEpLkB\nSUtz1K2kBpHu0ic+Rs0PyY3IOMNHh3QZ18mXDuxsQNJmolFyGxrSk+wjHFoakJYEpDogAWlJ\nQKoDEpCWtPJxJGMDxDTSZhFSEwkzOpln+JiQms7uTdo9FhnjdKOZkmY8yOOMC94pQTqjhH6n\nByRPQOoLSL6A5AlIfQHJF5A8GZD8ktzto0J+53cA0tm0Uud1dDbJfLaz6JE0NKpJD1JtUba/\nbSZmLChnRkDyBqS+gOQLSJ6A1BeQfAHJkxeS+7nH1g/IX/KtWuPLdVs7KeHckweSRubzU58+\nvJ00YmCH5ldG+5QU+jfVXfqbIM2tydHVz3YhvvEZISB5AtJAQHICkicgDQQkJyB5AtJAQHIC\nknvNE5NJ4Doj5kpbQxIeWnvqd32GtpfCp3Grzcj6OfdWF9Lw9U90Ms5vmTb9VHntiZsAQpsw\nd0dAAhKQJAISkIAk0M1D8lzO7uQi6VWkQZI8WfmsO9LflJsSUMwAABIOSURBVNfVrfv9G0++\nPLNp5zU0NJuRmQnpyfrZsSIuP9fNculzGT8gAQlIAgEJSEAS6OYheT58uF+Rg6iD1LwPrfoA\nhcWfAey74MKXb9Wv7xijx9iA8qGsv7ksvX9w1TaL35FryaY05tFXT+flSVVP7sKnMkFAAhKQ\nBAISkIAkEJCABCSBbhxSuT45N7qMPHTadbOof7JQ58WUM2ghjfpI7nAaI+vic1/6hbvbg0y9\neNzB27sujKrn5hRcjasHap8rZfzaKbloOxy0o0KTDgmFh5BXQAISkAQCEpCAJNBtQ6rXKOtG\nZ7UIIionrVa35i//+jCSAvSwGJIJxqPo5+dHl+Ra8igaB0mb+hhekb2/ZfrW+nYDaszjV4f4\nTtrzm29AAhKQBAISkIAk0E1DUuuUe5vlqFtX/Nsgl4nqF7o7P6jeQloG6SsAqVNUV99rIDAc\nWQ9kIqS+MbbzVM/d4GOqIY05k6qea/OLKm9JQAISkAQCEpCAJBCQgAQkgW4eku82B9JRX//0\n+tedpScJFebhIgvSz48jqcOkM3J/IwxCMhz2DLCb6+jHVB1xO4x5v0k1UvUFkOYFJBWQgLQo\nIKmABKRFJYGk1gD7Vq+jo5+RWtPcOS8G1NSg0aXYkH79lloHzuMwr3jX5Tg8DzEqK+c68VHV\nJ4YP/lw91BZS1pSAZAUkIM0JSFZAAtKcbhfS0b8KGJDad6NZGxY9mxAPD/Ve78Vv7CtsSNo5\nQi2kX5VpqR+Sl1E3a5Ph+G2f8TXv+Rv6uQY9kGYHJBWQgLQoIKmABKRFAUkFJCAtKtFxpB5H\nar3TLtlhfJSkWtPq98bqkz89PXQtH2C3W8GGVBRK0m+XZkkN0GEUcvRl9NlNvfwRhBq1u+HU\nWQbSrICkAhKQFgUkFZCAtKj1PrHP3kJqKRmQqi89Zzs/dZBErgKuHTHSIJX3NDf/GtU/WGoY\ngBR2pEFa/gCCjWKhQSr4VPM5AUkFJCAtCkgqIAFpUSt+GHMYUrUONpA8Uz6ZCY2nEeO5p4P0\nv7Ka0nfzHj8/JHfX95cLSf2eiNqYDR4TUt4ByQpIQJoTkKyABKQ5AckKSECaE5D01NpWr3WH\ngy7JguSuZ09xIBWjISlL5e6GEKRmX0Mfo/YglNT4l7QdR0AyAhKQZgYkPSABaWZA0jI3KExI\nxw5S4UB6chjJfbh5Ccl/x097/rcFqTlJaAQk11E1qdjol7UZRkAyAhKQ5gYkLSABaW5A6go6\n0iCdPwvnc/4cRKKQLvkh/XSQfr3bSNqnsjSMLEgeRllBql6TtYcwLiB1AQlIswNSF5CANDsg\ndQEJSLO7VUj2GSrNG2E9iqxzhMqfNg+zRIfk70eX9OuDZFz+5NQHybr4ST6QCgUpf09Aar4H\nEpAWBKTmeyABaUG3Csl8w6XG6ODpWL2Op/pK1OZcqiu0reGok2S8TbaDpF2cTkuHpNHpUDVb\nVykewehyv/QJkOqvgQSkZQGp+hpIQFrW7UJqurw+2vZRkJH/VOSQoxSQfrQ0R+0b+8wLK+uQ\nzg6d9rvPzywhHYA0KyCNCEg5BSQgAUkgIAEJSAKlhPTxZ7e7f61n0juXaJDcl6O1EoLUOPLP\nzQ8p1vDtXEo6JK+kIKRPE1KqhzAuIGl97Hdlj2omQBIISNmUENLz7r+Lpv/299VMgCQQkLIp\nIaS9mvB9f/e+DiTvi2GeGmRLqm8Mzi4gKcrwnX58aR/a11E6W5AUpe60oC93GynNQxgdkLTp\n6gk/7u+BJBGQMiohpLvdR/PVPZAEAlJGJYT03+5P/dX77j45pNArcQxDam/zTab9nIko5Ste\nf5LfdxBSa8mWZEHSRQWuyLx2eTtKuvv7udXzugOSREDKpqQHZN8em6/e/wBJICBl082c2QAk\nIMXsNiCpFd57lw2po9QYciE5+ySMhMcermGkS/p2IDk7HEZASvYQricgAQlIAq0FKe3OBiAB\nKXL5QNrpSSyiqzXhqbynWcVcVEfvAZUeRSmPdnx/qw/06yh1H33eA6l9l2wQ0meyh5Db4ar5\n3cSfdkACUuyABCQgCXQTkMKdzFxGLqReRUOOJkCbuI61ohxIraSzLqmD1E3UTDltyVSWFNK/\nl0f1lqTnf7EWMTEgAUmmlG/su9P2JtxHWcTkgAQkmZK+sW//96366v11v3uOsYjJAQlIMiV9\nY99b+/Xbbh9jERM7npwsRZ61ORWkyemORkP6/rYEAmleK7yxz/1GbBETAxKQpOJfJCABSaC0\n20iv79VXuWwjeSCddEWBzRSLTvnjYyElaDakWtLa499oKXd/32t77e4++n4SSPMD0iqlPY70\nXB1H2j++5HEcCUhAkuqGz2w4+hzp+SdzFakmQYr7AXRf3lOEPJC+rWA0OyABCUgCAQlIQBII\nSEACkkC3C+k4BCkwXcBRI2nswsUehycN0tnDKAzpG0hzAxKQgCQQkIAEJIFuGJLlxoQVnCyg\nqLYUF8jozI91cRz1QLpsJSU+RehK3iQLJCABSSAgAQlIAt0qJHuDSD8xqM9R4VF00k8ajz7y\ncemO7MeqINknCFVdpky9mQSkiAFpaUBKHJCABCSBgAQkIAl0o5B8ilpKfRO6OxismcQe+dhq\nSC4iDZLPUeKy+4zNuQEJSEASCEhAApJAtwkpoEix6JkuzKiekcTYBOZR9Rl0dCoPIxmQ2m8S\nX/1k4NfWhgISkIAkEJCABCSBbg2Sd/VvGzVxwFE1L4ERyq1avZtIn/Y17TRIq35i39OaC58f\nkIAEJIGABCQgCQQkIAFJoBuDdLDWf8tF77SaHS8jEUhqRotno5p6GKmBJLT4mwpIQAKSQEAC\nEpAEujVI9WnbMxxZm1fa1922yGJIUh6N+bmQyrfHlmzaj0LXIYkt/KYCEpCAJBCQgAQkgW4K\nkvaOvImIypytq3a9bCWNms/QQhbOwUzno0Z57nZ+//w4ktQpQucVz9xh97dgQJIKSIkCEpCA\nJBCQgAQkgW4IUvCyjuMdudvt9ZopBumynMVz0NPH2FZfHPLnx5GkPvty1XOEgCQYkKQCUqKA\nBCQgCXQzkPou2T1G0nEkpGyuftJkIWog1ZtIvwal+gPR1Xnj6zyUp0trLHdpQAISkAQCEpCA\nJNCtQDI+Z89naWgGjaOzc0q18DaSdF5Gl02kilGZJqlyVLRvCVzl98ImGQEJSEASCUhAApJA\nQAISkAS6DUgHT9McDUFq1zzhkQtU49GrjyL9/raS2l0NX80HwtSPZ+THtN98QAISkAQCEpCA\nJNBNQPI5MiyNmMd2IRXq8/sMRxWkXw+kz+7j0DN+PBkGJCABSSAgAQlIAt0CpJCjwwRH24Zk\nUar3fnsgfRqQ1h70lgISkIAkEJCABCSBgAQkIAl0/ZCCjKZCOgXSIK301oNxGbsaLEiefQ2r\nfrbL9gISkIAkEJCABCSBrh6SDKOxkOTGHaPqA1scSPUm0pe1iQSkSQEJSEASCEhAApJAtwSp\nkaN/PZrRNUAqN5Oq/7cg+Rhpe/Q5AXxEQBoVkKg/II0KSNQfkEYFJOrvhiAdvY2ayXEkJLlh\nx07f1RBy1EK61E25zes3Rg9IYwJSNyWQvAFpTEDqpgSSt2uHZDEqmg8nr78cO5sQIg3SeeOQ\nvIxO7UblwdhSApMdkEYFJCD1B6RRAQlI/d0GpAk7un31MdIhiY06ct/dpU8qR6Fd3zWkA5CG\nA9KIgASkoYA0IiABaSggjQhIQBrqJiDFZLRZSIpSDcm7o6He1eCDRHZAGgxIQBoOSIMBCUjD\nXTeky0ow5UygcCMcbekdpQ0k9WF92me5OIeQSkZrj3YTAWlUQKL+gDQqIFF/1w2pqM8IWjqT\nAULbUtTWnh8UPjnoGHXj6Jo2vIA0JiBFCUixA1KKgCQYkMYEpCgBKXYJFjGpMKANOyrfS9FA\nOtuMTkXNKObygRQ7ICUJSHIBaUxAihKQYie3CJmPWelXtFVHlSR1gpBnA6mI7qi4pjOPgDQm\nIEUJSLEDUpqAJBaQxnSF20dV7d5vc/voNOGYwcOc5XanlD9dy2XygDQmIAUDkgpIYwJSMCCp\ngDQmIAUDkuraIcl07tnRsHFI5XsodEPVMaRi/G+gh0vFtKuhHLtDVECKG5DSBCSxgDQmIAUD\nkioppH8vj7uyx+d/sRYRJ7+g7UMqigbS6WRBGp2CNL7ueuJFtYl0JZISQvq423XdR1lErIAU\nDkhVCSE97/Z/36qv3l/3u+cYi4gVkMIBqSohpP3urf36bbePsYhIXSmiKnWCkM5o+fsge2ou\nqXJs934DafJ0u9A3YouIFJDEAtKiSar4FynHgCRT2m2k1/fqq61tIwFJLCAtmkR1r+21u/uI\nsog4XS+jCtLZutb3aEmTTw86dnWQps4kz9IeR3qujiPtH1+2dRwJSP6A1MWZDSMCkj8gdQEp\nXG3Fa2jy4ZaMqx+OtpbPn1fvG14dR9dzhhCQegLS9IAUdRLPTDZwHAlI0wNS1Ek8M3HmstOT\nWIRQAUXX46iGJKCoUNczCd7phdRI2rgo/rQbCkhTAlLUSSYvIqOCkNYemFgWpEXzKnUE7/RD\nqu98Ck62jWas5fJwVlnGuBiJUzYD2fhIEryxb+PPUJyyGUk2A9n4SBK8sW/jz1CcshlJNgPZ\n+EgSvLFv489QnLIZSTYD2fhIEryNYuPPUJyyGUk2A9n4SBK8sW/jz1CcshlJNgPZ+Ej4F2md\nshlJNgPZ+EgSvLFv489QnLIZSTYD2fhIEryxb+PPUJyyGUk2A9n4SBK8sW/jz1CcshlJNgPZ\n+Eg4s2GdshlJNgPZ+EiAtE7ZjCSbgWx8JPmMnmjDAYlIICARCQQkIoGARCQQkIgEAhKRQEAi\nEghIRAIBiUggIBEJBCQigYBEJBCQiAQCEpFAQCISKBqk5/1u//zRd0OqnAX/d5fLSC79W+NX\nmTOQtz+73Z/3DEbysdpqclkvzJdi2khivYzq2ih3PTekylnwc3XDPv2L5XsKPvYrQHIG8prL\nU/K+VyNZw/SbeXnGiStspJfx327/Vrztd/+CN6TKWfDb7s9H+dvnz+ojKXtc4XPY3IHsLzd8\nPA5cVS3FSP5UY3hO/+IU5TD0l2LqChvpZXzevV7++3f3ErwhVc6CH9VDTr8C+56Cv7M+i0d6\nIH+r1fdj4DqfKUayW+vFufxmvTeWOnWFjTTix135j/Pb7jF4Q6pCC07/WnlG8m69eisN5I92\n2dx1R1L/oZuedHH5XWK8FFNX2Egvo/ObZbVfNYEFfwx9Fk2Skdzv3leA5Azkble87Ku/eNce\nyUv9p13yP1yKN//17Ee/PLcK6b/qX+6VR/Ky+7vGXzGeF0dd63P9kRT/lXsb9v8lH4k1DCAN\nj6TqfZ/8b0x3JNXfDXlAKnc2/En/74Dvd0tZ+n+QrGEAaXgkZR/75H/Y+f6iKvc35wGp3EZ6\nT39wwhnJf+WfdhfSq/yTlCOkvT0M54ZUeRd8v8LxLGckf6q/LleA5Dwlq/2Wc0Zytys31D7W\nON5oPf6pK2zUvXbv9l6797X22hkLfr+7X+OAnz2SJR9GLzqQ9Y4IOCNZb/e3vdSpK2ykEb9U\nv21fu0N8zg2pchf8mn6HnXckq0EKvDjv6Z8XZyTq34EVjmiVGS/E1BX29s5sWGF9CYykKocz\nG97LT7i6bJn8XX0kz7vy7Lbn9L9vy3I8s+Hyx25Ztcqq8Wk3pM0eyZ+V/h3wPCfmVysO5CWX\nF6c+w22dX3XNSzFrhY31MqqzeLVxaTekzR7JWn9QeZ4T86s1B/J6n8eLU59zvcJIChvSxBV2\nla06omsLSEQCAYlIICARCQQkIoGARCQQkIgEAhKRQEAiEghIRAIBiUggIBEJBCQigYBEJBCQ\niAQCEpFAQCISCEhEAgGJSCAgEQkEJCKBgEQkEJCIBAISkUBAIhIISEQCAYlIICARCQQkIoGA\nRCQQkIgEAhKRQEAiEghIRAIBiUggIBEJBCQigYBEJBCQiAQCEpFAQCISCEhEAgGJSCAgEQkE\nJCKBgEQkEJCIBAISkUBAIhIISEQCAYlIICARCQQkIoGARCQQkIgEAhKRQEAiEghIRAIBiUgg\nIBEJBCQigYBEJBCQiAQCEpFAQCISCEhEAgGJSCAgEQkEJCKBgEQkEJCIBAISkUBAIhIISEQC\nAYlIICARCQQkIoGARCQQkIgE+j9Y3SeteWNHJgAAAABJRU5ErkJggg==", "text/plain": [ "plot without title" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "image(radiation_slice[[1]], useRaster=TRUE)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "Unfortunately the plot isn't particularly helpful as its both oriented incorrectly and the x/y ratio is wrong.\n", "\n", "Digging a little deeper we can see that the latitude values descend (from a higher to lower value) as expected, it is the image function that reverses them." ] }, { "cell_type": "code", "execution_count": 162, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
    \n", "\t
  1. 6.5
    2. \n", "\t
  2. 6.40000009536743
    3. \n", "\t
  3. 6.30000019073486
    4. \n", "\t
  4. 6.19999980926514
    5. \n", "\t
  5. 6.09999990463257
    6. \n", "
\n" ], "text/latex": [ "\\begin{enumerate*}\n", "\\item 6.5\n", "\\item 6.40000009536743\n", "\\item 6.30000019073486\n", "\\item 6.19999980926514\n", "\\item 6.09999990463257\n", "\\end{enumerate*}\n" ], "text/markdown": [ "1. 6.5\n", "2. 6.40000009536743\n", "3. 6.30000019073486\n", "4. 6.19999980926514\n", "5. 6.09999990463257\n", "\n", "\n" ], "text/plain": [ "[1] 6.5 6.4 6.3 6.2 6.1" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "longitudes = (nc_file %>% activate(longitude) %>% hyper_tibble)$longitude\n", "latitudes = (nc_file %>% activate(latitude) %>% hyper_tibble)$latitude\n", "\n", "latitudes[1:5]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "We'll create a new variable where we flip the data matrix in the y direction.\n", "\n", "We'll also add a value for `asp` which is an abbreviation of aspect ratio, if we provide the lat/lon ratio we can achieve a more realistic visualisation." ] }, { "cell_type": "code", "execution_count": 175, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAMAAADKOT/pAAAAVFBMVEUAAABNTU1oaGh8fHx9\nACWMjIyampqiBwanp6eysrK9vb3DIgDHx8fQ0NDZ2dnhPADh4eHp6entYgDw8PDyhAD1oQD3\nujz40HT75Jr/9Lf//8j///+x2akoAAAACXBIWXMAABJ0AAASdAHeZh94AAAgAElEQVR4nO3d\niXqaWgNGYUyb1samZtAMh/u/zyPznpjkYwO63ucfUhNla1hhFJMUwGTJ0gMAbgEhAQKEBAgQ\nEiBASIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBIgAAhAQKEBAgQEiBASIAAIQEChAQIEBIg\nQEiAACEBAoQECBASIEBIgAAhAQKEBAgQEiBASIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBI\ngAAhAQKEBAgQEiBASIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBIgAAhAQKEBAgQEiBASIAA\nIQEChAQIEBIgQEiAACEBAoQECBASIEBIgAAhAQKEBAgQEiBASIAAIQEChAQIEBIgQEiAACEB\nAoQECBASIEBIgAAhAQKEBAgQEiBASIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBIgAAhAQKE\nBAgQEiBASIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBIgECEkBJgY66Yy/XhLDAJQImQAAFC\nAgQICRAgJECAkAABQgIECAkQICRAgJAAAUICBAgJECAkQICQAAFCAgQICRAgJECAkAABQgIE\nCAkQICRAgJAAAUICBAgJECAkQICQAAFCAgQICRAgJECAkAABQgIECAkQICRAgJAAAUICBAgJ\nECAkQICQAAFCAgTih/T0kCSPz7NOAogtYkhJfscfSW4/yySAhcQOaZ/s39P0dZ88zTEJYCGx\nQ9ol79nX78nDHJMAFhI7pCQx/iGfBLCQ2CH9qkLazTEJYCFRQ3r88/Sc/L18+b7v3ttASNiY\nqCEV8i9373NMAlhIzONILy9PT4+P+S6HfWdHhISt4cwGQICQAAFCAgSWConjSLgp6wkpMSkm\nAcTDqh0gQEiAACEBAoQECBASIEBIgMACJ60O2MNNSNiYiCE9ERJuVtSzv3c/5p4EsIyo20gv\nPRcPEkwCWETcnQ1PycvckwCWwF47QICQAAFCAgQICRAgJECAkAABQgIECAkQICRAgJAAAUIC\nBAgJECAkQICQAAFCAgQICRAgJECAkAABQgIECAkQICRAgJAAAUICBAgJECAkQICQAAFCAgQI\nCRAgJECAkAABQgIECAkQICRAgJAAAUICBAgJECAkQICQAAFCAgQICRAgJECAkAABQgIECAkQ\nICRAgJAAAUICBAgJECAkQICQAAFCAgQICRAgJECAkAABQgIECAkQICRAgJAAAUICBAgJECAk\nQICQAAFCAgQICRAgJECAkAABQgIECAkQICRAgJAAAUICBAgJECAkQICQAAFCAgQICRAgJECA\nkAABQgIECAkQICRAgJAAAUICBAgJECAkQICQAAFCAgQICRAgJECAkAABQgIECAkQICRAgJAA\nAUICBAgJECAkQICQAAFCAgQICRAgJECAkAABQgIEoob0789jknnc/5trEsAiIob0/pA0fswy\nCWAhEUPaJ7u/L/lXr8+7ZD/HJICFRAxpl7zUX78kuzkmASwkYkhJ0vYP2SSAhbBEAgTibiM9\nv+ZfsY2EWxNz9/cPY6/dw/sskwCWEfc40j4/jrR7/MNxJNwWzmwABAgJECAkQGCpkDiOhJuy\nnpASk2ISGqeAt8zSA9PJn9CxNuWhDpm2bzaTOB4Kvy/Kb9ZfbBSrdu2yZtL2lJYenkzxdI6K\nlNorSltCun5a60JI7QhpNEKa9y4rnMQIoZJOZWWb93E+l+uqTUsD7/pz9MTclH7fTkmE1K81\npFso6RLSuQqpzGnGqb35Id1ISoQ0QHtJW4/p8/Pj4+N8OlkpjXuInyMXTG91SCkhzY2QIiEk\nlajvRxq8h3t1IeUz2O2F9FWHdHJCGrqh9PNnHtKYGJqQ0mzlbvyoVyliSE8LhKRZ32/msdsL\n6VLS+Vw/GXOnw6BHGLs8Kpj77a65/wrFXLV72XVf8kQwCZc6pPBSSTKRJXx9fdWLJDel9Nhx\nbFWj6+jt1kTdRnrpfjufYhIOQupESDJxdzY8Ge82n2kSc7D2Dt/Q6t1XFtJl5e5sPZlyO6nY\nkJl1+jfU0c3vtZt++ljqhhQMSjLYuL6/v6uQzufgQqnYkplvBIQ0N0KKgJCUbjuk40FwRnMm\nnJKdk2LAkXxmK3VFSF/Ffru2ki4tLT3YbbjtkIrT+g/TUmqPyC5pM2exXhZCZUjmIimQUrNU\nuqVlxzzuIqRZS9pqSN9GSB9lSu4Bpbe2kG7l6I8OIfUjJELqdR8hlS1lMRlvFRhR16CQNlNS\nvuP7+7tJKVCS+SYlQup34yG5JXkGPUjxo90hZV/KRj07P6T2kuyKbuUkUzVCGoCQmnsSUtj9\nhNSUZP5j4MMMC2ntKV06Kb5oNpKMzaTgdtL8pzfciDsKyTV2obT5kIqOzpdknJDCJS093G25\n9ZC6UhpT0sCQ5nyb9lRlJm5IWUqfZUkfhHQtQiIkQhK4/ZDaSxoZ0rHMpj0kwblIczkbspLs\nkOqUzHW7da+nrs0dhNSa0ph9d5sO6WzLF0mfwZDMRVL+hJYe+mbcRUihkkbuBd9uSPno3JQ+\nA7sbmpU74wmxy24gQiIkQhK4j5DSrogGldQX0ttaQzp58pU76wTwMqRy3S6/T/l0Fng+Gz3e\neychGSkFKhoYUnjHnR3S6lIKvDO+WSR9lR19ftYhXVIqQ1rkmWz2OneEREiEJHA3IaWTOkrb\nzxISr9qJZ197Q85Yt8v3gH99mR2VKaXnRY8hbbOjOwopLVoqy3mbWpLxtTIk9RLNKr4ebbUH\n3O2oPJa06DEkQhKaaxJVSHUMY2qyV+6Mfxhvg5s4vvJhJj6K+3j+MrTY3VCU9PnphCSb+l0h\nJEIiJIF7C8nZZzCmpGqlsLmbM4MqNpHyB5r8KCV3T30dUrVu58t3NmC8+wrJXay4JQ24r0Uf\nUqrryD/kZYRU7Lgz8jFCWvZ0VbaRhAhJgJBiureQcu7sb5XRdcfOlPpLHDS26Q9RcD8+zEqp\nPnU1FFK8dTvdX43F3WVIzka400XXHdtLGrRMi8g4Yy60TDqHtpIWGCYhzWr+SThzllNF+/1C\nJ0iYKc4+8GFC16jzQ/L3N6Rp7L12t1ISIRESIQnca0hpaOe18a+2u3WeQ365bf6BD+FfzNuu\n6eRuJDUbS5cIow6VkOYUJaS2izAMT+ngZHQY/O6deYOrrwcUujB+d0izjuuW3W1Il5m5O6S2\nlFre2lT+e+i0Zc8ixAjJTKk3JE5ruB4hERIhCdxxSOGSzN0O4bsF3rPe/CvKwDs1l9VyU/JD\n+nA7oqRrEVIoo2YPnn+vQ6eeSc5ZWhPDVSE195xviDeMkAiJkATuOaS07XCS2ZJ3l0ghjTx5\nxgliaEhlSx8NQrrOXYTU8Yl95mFY7/1J2feCjzchJaHslJ4qn+a9rm5I/r6G5u3mvvwu0Z7B\n7Zxqdx8hpe1zeOioarMTofhgv8CdVlGS8Z488+oLfSGVHbWHxFW/r0BIhERIAvcRUlq0FLi5\ntaND50X2rZ/7bYlZUvlG8eZ6QF9eSIHto96Q1vhJFCs4sNDpbkJqWSx5p8xZC5f2kpoNpd8B\ncwzf9xVghxTePvJDav5V3mt12y5xtz6vQEiEREgCdxRSmv86/FsO/ok+TksdDxfsKNa7pY2A\nnA+EbVmtq0I6WyGdjX+uM6TIG5/XuK+QQg5pZ0hNSoGjSm0hxSjJWhB9B0IKvX+iPSRnibSu\nHdOEdKUIk7B+L4eWVTurpHw29B9muUWSW1FdkvPhLG5I9Y6IQEgfqwyp+LUsPYpOhJT/g5AI\naZp7Dcn5rRzsnQ6Bko7VqTXZj5+s+WyZkuzVuv9yZkhNRvXw3ZCMepz93/MPf6h8rbpjO3Ut\n7jWkgPa9d+WcWMyEafHhXcYd5w7p6yt8q7E4+i8Ykn3eU0dIzqGk9YS0gYJKhGRoK6kJKZv7\nikvGmff77cekG1R2vLXt9jKk//6zQ/psQjKfkR9SsKQVhbSyNcwuhGQgJEK6FiGZ8tktGFI5\nEzYHOt27zrZyVwTTdvt3KKRi1c47oT0QUiilskHV+KcIHHFYK0JyvDUb6OYCyQkpdM9ZQqpO\npQt9x9vTYB1H8kM62seR2koq9/dpxt9qyD648ncx80g0CMlBSIR0DUKy1MdbAiGdrFW7U+qv\nwl/y+VmSlFS+Fbx6f0T5Rtb8W9Xt35avJiR3n521365z7e58nn/lbtiZClZIHEe6QqRrNvg3\nWSEd7aNI5gxYlJT/oHX/OqSL6QNsQrIuq5DdnPohVb19lHu/u0LqKqkOacaUinOr+n7KCGn1\nR2TvNqRi1vJvdxdJ9nb6qemoOk5j372JSBJSeXUSP6TqTUhGRPYbzavxHb2UToGU3DNYW56d\nRPHnqb+j+ghy+nZc/6kNhOQgJEK6xt2GVF5G37vVWrcLzYCm8KymWKvL2SF9WB1Vb4stI/Iu\nfNIeUriks9mReZ5e34s4+tyDMqTen6sySgnpanG2kYKzyVuwpLdwST3XwZka1EfownNOSMGK\nzsYSM3CxvsAOB0f7ItdWP+jwJ3Usz/Tt/cFioPlXRUjDp7EAQgrcSEiENNa9h3RyD5iEQqpv\nC6zdZbuJi190tvu7epRsM2nyvjujmkBIRkrlt90MjDdRjA3JLLFz3a15zHTQZerKlbQhBwfK\nB81GQUhXinQcqZyj3Btb5r6GuUw61b/o7M7l/FHsbchbmjC8c1tIzW3mjnErIS8l/wkMD6k9\npeDL0xFetQtnWBV5l+UQ1p3RnYfkvSEiN6AkI6a0nnmyX3UeUlFP0dKkkKzrc9sHe6zPkShX\n6kLLy3BK/SEVr035IgTGlj3Rg/fyGEvx7NVo3vVY/Ku4w7hlS/sQ1oWQCImQBO48pLQzpDd3\n73FbVcdq9SffiG5Cmja0s1WSuS/hwxNarWuvKhCSdRjpnPada5fVEPgU3aqi+lhRXpH5bpRq\nJbjPsfyft56NtLW495De/A3kYCXdMRUhHcqQRKzPXrFmc2fm/zC2aoYLPE71WH0jc96pFc6o\nPnvX7Wjg0z8Y7w+b+ErGcO8hBT5Wu2N9pa2lOiSpk5GSH0xo54B7ZWLv5/xbAwu33oG573kM\nXjLGysjMbPgLYF0nbcILGQMhERIhCRBSgF3K0ZtZjsEfkF9Wt5jh85Ca1a7+DSC7pCaXUIZ+\nRkM+RfZQ78GuXpE0/48hNVfsrL9AI16A8sVOR+wxXwwhBQQ7seYTfwmVVr9u3e/7rVkk2RkF\nljvhK9g511RtlkrtGfWGNO5JmiGNfwXMCRHSFQgpQ0iENNGKQmrbCHA7an5INoyWkEK9uBXV\n33F6MbeXwhmJQ8pdFVG++9PaCbrqkgjJYPy+6wWSHVJw87mmXiQZx5K8kILHs6rR+iH5H3d+\nXUbFIibaLC2+1OaMCMlgrX+EQjL3RIVKOspDqnY3uCHZa5/mSI3vmCHZx6LKhwycNNE7pPzp\nC59iz8QIaQpCKhESIU2x8DZSpTOk1pKMB5g8FxgpnP2Qju6WnJPSydjIanrxOnLO4+sZUqyQ\nfqZjO1q0OUKqBN5L0xaSs0iy52PzASb/QQ0dUG0pxxhnVdLJ3lth8856LSLrG5H2esx9fv7s\nn1gxooP6UtEjEVKFkAhpAkKqhU9fbQ/p0LyNIJxReY70lDG55/j0hGTum3c2koaElP1s34ha\n59Y5Vvda32JsrT8XI5Je4HY8Qqq1XcDBnElbQvLeiaMa0cktqUqpoyJ7kTQkpLqiASHF2fld\n9TPw3frFL+T3kvsmCKlbuKTOkJo7T5/lWkPySzocnJKadbtgSsGQsjt0j6gzpOAFzq5hXB5w\n4HsjQ8OSXfJ2AELqRkgOQgojpErLHGDsXm62k6xVqPZtpOkhBa7UHQ7p4HVkX67cKOhsxdTs\nZqjO5OseUXdIU5+voUrg6hDM67DPXxMh9TIO1RwOgRm2+Wwyd/tIsjHhluSHFB6U/yEu5pl2\n7vLJWNh1D8f7a9Hz01f/MWlCuiqC37/NkObfDUFI/YxZ1y3J/Sx0+bQv87VTkr23ITySaoUz\n8Al9NiOk+iy+7vGk/oK3beRv178oWTzFtQHT65ZJ5V6huqJ8NSI/m7JtPBP3ohBSP0IyxkNI\nYYTU79RRkr+dop32m78LfHRI3RfectYY+16J4oXoG7O1Nnz9s69W63pSMl73323qF6ntd5Td\nPmEFkJAGqN/z46cU2FSRTvqtXCQFZnlnMMGOvE+eONncRx38Qvjfq256ezNCSqNcath41bsq\nqjtqHdGEkghpAELyXgj/e4QU4S4rnERQyyaCu2USLskKTDmoUAvuUMIZNSWd/IismoZklJ7t\nF8L8ljPV5sFiveWifNU7MjJeMON+9u7+Q/7BTdfEREim3pD8pVJLSMKc8k/M8EPyP+c2mFJb\nPdlD2P/qP5TaHpIxzWoAkuc+VP4iGGfctSyOhoR05VKJkEytO628+fd48ObgmTpKyw8Jagmp\nVWdJxpKjuKHteLTp5KxZNt8JTD7NV/am9DTuAFI+/9c76n62dVSfn9LxSFf93gjJREhdCKkD\nIZneqs+IcxlbGMYsbK9RzbZq14zB2F/QrGdODMnaoul0MlbtzHW70NTr26NtJRXZWOcytIX0\ne8Ju7jaEZCn/PPu8j8ELzDz2hpNukVQPwQ+pqyZvf8O0kE7NDgtrok5Ixh+QY3n12atfhrEz\nvHlunReS8WshpNkRUtdLQ0jtCMmRzS6Bm+2dx/1bKMVMJR2Zdcwn0EJHSn5M1t0GDiD07J3J\nW2u0h3TKebvNVs7Ye7aGVLrmQXvvQkiW4s9u8G2i1mzcuSho2I+Qb4APnW/9BaMTkt+EsUzw\nS3rz9nY3dxr+4lhP3e/Y3jI8pJNialIYf9/ekK4bS8dPEJLrnJ/KGf6Gc3LByJKO5ZWNrxuX\nc3JDIKamo6olZ5kVDGn4CN7skAILxP7HaAofMMEhBbVPtCWkriTq8xiNPwfmHsCOmAjJRUht\nCGlVIT09JMnj86yTmKbtEon5LPwRuuKpxZyx7F9yPc9dNayWTwxztlvsLWv/QzpDq3ZDR/Bm\nP2nj8cpbBzyG9SJdd3V9c0TFUw59q3XVrmcVrQ6pOFMifPJe4H4RQ0ryO/5IcvtZJjEv5+LB\nRi/H0FxbLIDK+1rfuXbyLaqZ2Q6ppaSyJ3OYQwdgLXysLrNn1P1A1TfLBXX+GrxV/zoeg9dw\nGjCc9sVg966G7m2kN2M5FwwpcPfYIe2T/Xuavu6TpzkmMa/Q+7KNksKzbDl/vJnrRdcOoLOk\nwLnpXcMa3bX1iVFWRtm0++9d/X/oxPGrF0xd65OBiOqOhu9sCIfkP0TskHbJe/b1e/IwxyTm\nRUiEtJqQksT4h3wSSv4eh+ClQvIZa8h8OmKua1NspLWGVK/cD1q7e6uHP2ziHRtcpyl/HdJp\nZ+S1C4VkNDA0pvaSrJhih/SrCmk3xySUPj8/nVvskFoOi4ZmNeMxTt4tg53dz7C0U2oPKVyS\nOfoBUz/0hTRhx8HUfQ5tDuYHRuc3pEUYvSdb5C+JdXHdsEP9SFFDevzz9Jz8vXz5vu/e20BI\nAYQ01s2GVMi/3L3PMQmlT68k53KKoU90bdaX6rnMfowJIZklu0WVIzE+orM1JXeH27COstmx\nN6S4b+YbyovmcOgLqXo+odVCN6TysWIeR3p5eXp6fMx3Oew7O1pFSCGfzkex2iWdfO6SY1pJ\nfkgfbkjNG27Cm0nNP3UhNY989Zl1ERKszwBsFlCtmleqoyT7G5zZMErwGr/tIXkZZaas3QXX\n64yi20JyPtnPyaB/wsUDdYa03kVSrogn+5idOiRvz1vFeKn6Oyr/SUijEBIhEdI1xmwm9a7Y\nFSkZP2A+9JC4Qg1NDWnQyxAMqXoWGwjJX087tIZkv1ThjPy8lgpp/ceRSp+fX84trbvuejMq\nSqrPk7MftrpX12A6MjpVG0nhkOzPnhkVUv047gKpGrgx58X7vPMBitck+yq0gElbDiW9GU/n\naOy9CVW0ypASk2ISGl9fbkj+2l1wkRSe51PjPYLWowaXUq5wqsaSwfu7a1xWwj8JI79T7ytw\nMM8kr59uM9CmrHUtk+oXNHR+cdmAf6+jo7uj5lZW7XoQEiERkoRXkrUX/FTPWlZKrVsz5UMU\nyRjV1HfvGsqpumvVXPHfJqR63W5QSOYKWpfqLmbExretkMq1JinzKphpM6/33s/9M2PGVGXg\nTasvpHBHhNQvFJLxcZHGr6mZ11o7alpy5sf6/r3jucz7TQAnazHVeUXwek4yMxp0DMm5kLhT\nu/FgR/1n27gXV7GeUHn1vA7WQttsqSrBvYN70Vp3ebSKkP79ecy3gB73/+aaxAwIiZBWFdL7\ng7E34ccsk5hJYDPp09lMOjn77vpDSv29dm89a3YtzJDKmcurqPUE9f6Hr3+0a83Oeo/uFc8h\ne5xq55r/rbaU6lm942H9s1DaQvL/9tQ/GA7JvC1iSPtk9/cl/+r1ebf+k1YN3umrn3ZJ9obS\n0JACrjtMa4eUp3EsLw/R19GwCXgZWSWV3zRCGl9SyzKtmK27O+qdoLdbyAjJvGdwMT4kpPzr\niCHtkpf665f1v42ik1FS2867a0K6SjlfVyUZ38kDct7mfk1KoZDqlKyIj9ceSjJCah7hUDyJ\n+g3qqbVsSo05vfOxy19SOCSrWiej5uBTy6qdIfo1G0L/kE0iGkIipMVCuqUlUubz09zpYITk\n7XKwz9aWD6Saq0/ummEVilWOOeMXOy56J1AWZ9+/LaTrDiU142v5UD07evv6/YfuqTYh1UMM\nfeKBnZFd9mFNIV22kZ5f86+2to0U9tkoUvLnVjejj7ZLfU0yYMPKaaBKYeBCyQ3SCKl5MPOP\n/RVPwixp+A4LcwHVoQ6pvIffTXabsSiy5N/szihqSNWVuHIPq39jXy9CSglpkZDSf/v8ONLu\n8c+WjiN1uDT0dZGlVF7p7midfxMnpIEC2ziDMjJYfymch21CunJ81b4/p8XiBe28Zz7vd/1A\nFVL5s35FNuM3aKzddePMhmm+voqSsl9VNRcYZ1evJ6S6n+qoSnbjyJDMv9XN9pgX0tUl1TOu\nPdG0d5TdP1BvI5U/FFgchSuqnldPSPndCGmar6/vi0tKVkj1GZ7ezoZFU/KPCY+8ipbx5JyH\ndi8/fu0gr76sV8f93JDSQxGTsxZn/5kwvPlXlfEzJKRpCKl4aEIaP2RCMpghnY39qx2bSUsN\nNXAY6zTyVArjz4T9DWNmPM5zBvgk+Yq39VyNEoyBGy15JXWtB14Q0jTZFtJ3tb/hVB9LOjbv\nO3AWSssNNXAwePwZSfUc5Ny+6o7S8kBf/e9ifOHzYc2QQh/1+eYoYyOk6aodd8bRc2+ZdF5B\nSIoHqf6OB75lhZTPqaH3+6xLx3LGCMm5eqG7+zP//RLSdISUI6SRCMlRhvT5UZ4rVJRkvOpm\nSUsPdqJ628L/Vj0fVlvjl9vW1VFgh0TXpo+zmWTcYp87mf9eCUmhPsGhPuvu6Jd0voGOmvfs\neN85eiGtbEOpVXBnXGDJEy7o4yP/3ROSAiER0vjpElKAmZKzTl3/Bm5l1S6f3ZzbzRWjjXXk\nltQcYAp2ZB3KKH7r2bo9IWl8Ni05568a75i9gZDqWtpuN761qZjScqe4ndJbKKNPS3GSGCEJ\nVfscgmeCVysDSw9ykkAtpUNgd97WQio1Kdl76cqOPoMdEZIQIRHSGITUpjqe9NFcltDZSt14\nSPUfh8DOBjuk7IttbSiZwi0RUjzG8aTzuc7J3Muz9AinKI+otHw4WRNS+cUx8Gl529HswTNT\n8reRSoQkRUj1V4TUh5A6lSl9NDmdbyWktPnb0HK6axHSWx3Sph3sllr225W/a0KSs1/guqhb\nOEXI+OPQ9WPZRtRsH1UeW7OxFDqQVP96CWkO/h+s+gje0kObwngy4R8o08lX/W4lpLQ6U9c4\nS8VKqFhME9IcCImQ+hFSv+qaKG5M3lXEN8X4+xD+gTKefF9E8HDTZjklGYpNRt4hO0XHjGKF\nVLXUPgdugvl8ws+jeEHyq7JWp91lN6/iXUkTt07LfQ5eSafyzEpCmqDjT66dUcP7rKUNsRaz\nwbW7JqTq+nGRh9hu+l4e8+w7Z3FESNMQkoOQxiCkYYoVoHpV6DM/nWTpQU1ihJT27Mhf3fZR\naLijP47KOCHc6YiQ5lNuSTSbFF9bD6lczg760RUtjc7nNHjgq/24chvvnRUskSIgpJUgpM2r\nDriUM192Lcklh3OfwrGcu850atXyDnRCmpl36HLmjj5jnoGk/8S0qIoCRt6pZZHEcaS5nc/b\nPrmuy1ZCCp6qfmVIaVtInNkwL0JaHiHdqP/+y/77X/3PJccyxVZCapFv2lxxP0Jah/9M+Q1b\n3QOx6ZCqfQTj79myTCKkyAhpFQhp6/5zfZeWHthomw8p/I75XuGSCCmytpDKi9EsPbwRth1S\neXmg8fe0r9FFSEvpzOhr4KkDxr6K5Ww6pKqjK04IDJdESJER0ioQ0tb1dNT5BlTzMeKM9maV\n63VXnaFOSAsq0/j+7gqpeTtt10N95w8SY9A3q94+uvKtHn5InCIUzWUxU3wEektH1tvSux7o\n27SlvRPrceX+ugZ77ZZDSOtBSNtVvEE2L+k7tF5nXSXlo2NDvvjx+p7Fu3DjPY+bMC0kvyJC\nWoq5TLEuhFBfLa38/RQXUUs/mn0Qn80nINSfAV3dK/t+screrMN7086+F/np3hZCWg9C2jBC\nWo/shCAnI/tS0tXvx77YtH9VIvuSn01G5m6lpZ/s7SGkNTFDsiuyPqcscN12O7z6svaX/5xO\nzmcF5lOaq6WVbptFuHwRIa3Hl7E8soJwP+7vHEjJjai6Y9FM9ol6E3dLbdn8VwKzfzdc+3tR\nhDQXQoo5icU5+wu8y7Kf3N9Wd0P1PZd+XkubeohooDMhrYaxnyFcUMdCyWrHXCTds+J6yTOH\nZLzc5pWLCWk5hCRGSAtMYk3yLqpogmfp+6t2hOTJOjrNu7EVfSsAAAxxSURBVHZrv9rmXzNC\nWgsjIj+lMiTnPCJHeQRq6SdyPUED+SulGEvLw9shGb8hQlqRLJ/i86yyf2VzlVmS/2HaTkYf\n1bFc4yG7ztpbkews3fN1VxE2lC+DbFQOa9cPIa0YIRGS1p2GlPOPgmQzV71mZ59iZ2fU3FDM\nVOfihusuOxVR8E/AtQ8kGVGQv4FESFt0Cck/W9UNqbmp+le+6bWuD/1yZacaFeMefJfwzrmh\nnzlzlcB+hqYj3iG7JUZI7nHcljMf8pDeVvbhea7rZv5ASLN2lLoVNSHx+UgbQ0gGQlrHJDYr\nvI3kbyZVe8PzHRiH3+W9f//+PeST9EatZoVJLtAyautu1nPRAxmd02bvCNtI22PttTMPNIXk\nG0iHQxbQz8t9f/7OdTx681gTh9l2Feaq9+rZdD3GuMuhfkRYHhlbomlxwl3xbULaHkJqQUhL\nTGLz8t9kvmZxSq13mqX5RVDr3UnHYxZS+rPilHQMvumgfKB8Hg7Pyb3zdx7SVxbTV/6fTJZP\ns1fxmme9FLcjb/SEdEusAzHFX/OsozIcI6RyO+lYeKvufSr7dN5p6zZjncPUsuQqPhAg/295\nFeZ87jv3Hy0atRSKpXqWrfkT0i075osjcwnUGZK9Tvd2bJQ/Xt6n6aj8K+1M9rNY/0zL/1Zr\nXD2ri+em3vnfmidHSLeMkKIhpFuXh/Oz+tclJON7TUh1Fmfz/HE7pENzkxGSf+ZSGtqGCDEn\nWuw32WhFKSHdiWzjqPwyEFJafsRwNUPn3zvknGXSJaYypLqBLBr76mDDdvkZORrHNYuzdoXP\nPBZCuguENDdCujc/my/NSIo5ulq/OmYRFTFd/sefsbMAmjPMP6pLwI7Zo11EWu0AHHfEaIUI\n6a7lFWWdGEuH/JSiw5ATidK8pYEbRK5mcVftY9/igqhGSMhVM3SsiR3Cq46bRUjIEdI0hIT4\njk1JhaUHNB0hIa5qz1yZUnoLFaWEhLGmrIXVe7lvY3XOREgYh5CCCAmD5DsGsn0RV8//b9aR\nV0IipPuTneZq7GG7KgHjNIZtH3oNIyT0I6RehIR+P+2Srthf3Zy/cIvrdSkhYZDiTUz5EdTy\nyE/nhR885plAN9kRIWGU+hBqzzVUHEfbfONbDCFhDEJqQUgYxQjpZ/9PV4yIbuOMIA8hYax6\nidS8W7DPjVeUEhLGI6QAQsJo1crd4I6qkm7kTO8QQsIVxu1uuOWdDBVCwswONx9RhpAwM0LS\n3WWFkwCUCAky4R0Jt74sKhASZEIh3fw6XYmQIENIc99lhZNADPeSESFhPvl7J5YeRCyEhLkQ\nkvwuK5wEZnc3q3UZQgIECAkQiBrSvz+PSeZx/2+uSQCLiBjS+0PS+DHLJHCHrE9nqj74Of36\n+rr8N/sizigihrRPdn9f8q9en3fJfo5J4A6VHxiYaT5CPeuo/J84o4gY0i55qb9+SXZzTAJ3\n6O5CSpK2f8gmgfuTfwj0JZn6A6GLlC63fH/lip7GGX8PlkjYtDyigpFR5rv6vyKnmcXdRnp+\nzb9iGwkqX19tIZW+v28tpPSHsdfu4X2WSeDe3GNI6b99fhxp9/iH40hQ8DLyOrqEVKQ080g4\nswEbZkcU6igPqdrtkLEOO+kQEm7cd5NSEdvHR1tM39dPhZBw4247JI4jIZJvo6SypXBI3xM6\nWlFIiUkxCSAeVu0AAUICBAgJEOCNfYAAb+wDBHhjHyDA2ygAAd7YBwiwRAIEeGMfIMAb+wAB\n3tgHCHBmAyBASIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBIgAAhAQKEBAgQEiBASIAAIQEC\nhAQIEBIgQEiAACEBAoQECBASIEBIgAAhAQKEBAgQEiBASIAAIQEChAQIEBIgQEiAACEBAoQE\nCBASIEBIgAAhAQKEBAgQEiBASIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBIgAAhAQKEBAgQ\nEiBASIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBIgAAhAQKEBAgQEiBASIAAIQEChAQIEBIg\nQEiAACEBAoQECBASIEBIgAAhAQKEBAgQEiBASIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBI\ngAAhAQKEBAgQEiBASIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBIgAAhAQKEBAgQEiBASIAA\nIQEC8UN6ekiSx+dZJwHEFjGkJL/jjyS3n2USwEJih7RP9u9p+rpPnuaYBLCQ2CHtkvfs6/fk\nYY5JAAuJHVKSGP+QTwJYSOyQflUh7eaYBLCQqCE9/nl6Tv5evnzfd+9tICRsTNSQCvmXu/c5\nJgEsJOZxpJeXp6fHx3yXw76zI0LC1nBmAyBASIAAIQECS4XEcSTcFEICBFi1AwQICRBYT0iJ\naZ5JAHOJGtK/P495Jo/7f3NNAlhExJDeH4xFzo9ZJgEsJGJI+2T39yX/6vV5x0mruCkRQ9ol\nL/XXL7yNAjcl+jUbQv+QTQJYCEskQCDuNtLza/4V20i4NTF3f/8w9to98MY+3JK4x5H2+XGk\n3eMfjiPhtqznzIbIkwCUCAkQWCKk/lPpCAkbQ0iAACEBAoQECBASIEBIgAC7vwEBQgIECAkQ\nICRAgJAAAUICBAgJECAkQICQAAFCAgQICRAgJECAkAABQgIECAkQICRAgJAAAUICBAgJECAk\nQICQAAFCAgQICRAgJECAkAABQgIECAkQICRAgJAAAUICBAgJECAkQICQAAFCAgQICRBYaUjA\nxlwxl+vDWWQawzASz2oGsvGRENIyVjOS1Qxk4yMhpGWsZiSrGcjGR0JIy1jNSFYzkI2PhJCW\nsZqRrGYgGx8JIS1jNSNZzUA2PhJCWsZqRrKagWx8JIS0jNWMZDUD2fhICGkZqxnJagay8ZEQ\n0jJWM5LVDGTjIyGkZaxmJKsZyMZHQkjLWM1IVjOQjY9kPaMHNoyQAAFCAgQICRAgJECAkAAB\nQgIECAkQICRAgJAAAUICBAgJECAkQICQAAFCAgQICRCYLaT9Ltnt37tuiMWb8NPDWkZy8W+J\nP2XeQF5+Jcmv1xWM5H2x2eQyX9i/inEjmevX+CO/qP9Dxw2xeBPe5zfs4v+yQi/B+26BkLyB\nPK/lJXndFSNZoukX+0MoRs6wM/0a/yW7l/Rll/xrvSEWb8Ivya/37K/Pr8VHknm85iNE5APZ\nXW54f0z2i4/kVz6GffxfTpoNw/xVjJ1hZ/o17pPny//+Tf603hCLN+HH4inHn4FDL8Hfqz6L\nRz2Qv/ns+57sFh9JstQv5/KX9Yc11bEz7EwjfkyyhfNL8th6QyxtE47/uwqM5NX57S00kF/J\nS/RBhEdSrujGTzq9/C2xfhVjZ9iZfo3eX5bF/tS0TPg9+bGCkfxIXhcIyRvIQ5L+2eVrvEuP\n5E+5ahd9xSV9ceaRsTPsvYb0lC+5Fx7Jn+TvEmsxgV/OY76Jv/xI0qdsb8PuKfpInGEQUv9I\ncq+76OuY/kjy9YZ1hJTtbPgVfzkQ+tuSib9AcoZBSP0jybzvoq/Yhdaosv3N6wgp20Z6jX9w\nwhvJU7Zqd0l6kUXSGkPaucPwboglOOEfCxzP8kbyK1+7XCAk7yVZ7K+cN5KHJNtQe1/ieKPz\n/MfOsLPutXt199q9LrXXzprw68OPJQ74uSOZ8mH00oEsd0TAG8lyu7/dqY6dYWca8Z/8r+1z\nc4jPuyEWf8LP8XfYBUeyWEgtv5zX+K+LN5JiObDAEa2M9YsYO8Pe35kNC8wvLSPJreHMhsvW\n0Xu2ZfJ38ZHsk+zstn38v7eZNZ7ZcFnZzeSzbDE+44a43JH8Wmg5EHhN7K8WHMiftfxyyjPc\nlvlTV/0qrpph5/o1FmfxGuMybojLHclSK1SB18T+asmBPP9Yxy+nPOd6gZGkbkgjZ9hFtuqA\nW0NIgAAhAQKEBAgQEiBASIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBIgAAhAQKEBAgQEiBA\nSIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBIgAAhAQKEBAgQEiBASIAAIQEChAQIEBIgQEiA\nACEBAoQECBASIEBIgAAhAQKEBAgQEiBASIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBIgAAh\nAQKEBAgQEiBASIAAIQEChAQIEBIgQEiAACEBAoQECBASIEBIgAAhAQKEBAgQEiBASIAAIQEC\nhAQIEBIgQEiAACEBAoQECBASIEBIgAAhAQKEBAgQEiBASIAAIQEC/wMDOhyptmRNVwAAAABJ\nRU5ErkJggg==", "text/plain": [ "plot without title" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# We'll reverse the matrix\n", "radiation_matrix <- radiation_slice[[1]]\n", "radiation_matrix <- radiation_matrix[,length(latitudes):1]\n", "\n", "# And also the latitudes values\n", "latitudes = (nc_file %>% activate(latitude) %>% hyper_tibble)$latitude\n", "latitudes <- latitudes[length(latitudes):1]\n", "\n", "# Calculate the lat/lon ratio\n", "lat_lon_ratio = length(latitudes)/length(longitudes)\n", "\n", "# And visualise them\n", "image(radiation_matrix, asp=lat_lon_ratio)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "Let's identify the coordinates of the maximum value" ] }, { "cell_type": "code", "execution_count": 208, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\t\n", "\n", "
rowcol
378113
\n" ], "text/latex": [ "\\begin{tabular}{ll}\n", " row & col\\\\\n", "\\hline\n", "\t 378 & 113\\\\\n", "\\end{tabular}\n" ], "text/markdown": [ "\n", "| row | col |\n", "|---|---|\n", "| 378 | 113 |\n", "\n" ], "text/plain": [ " row col\n", "[1,] 378 113" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "max_rad <- max(radiation_matrix, na.rm=TRUE)\n", "max_rad_idx <- which(radiation_matrix == max_rad, arr.ind=TRUE)\n", "\n", "max_rad_idx" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "We can also convert them into longitude and latitude, from this we can determine that at this specific time the maximum solar radiation falling on indonesia was on the island of Bacan." ] }, { "cell_type": "code", "execution_count": 209, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[1] 127.7 -0.8\n" ] } ], "source": [ "max_lon <- longitudes[max_rad_idx[1]]\n", "max_lat <- latitudes[max_rad_idx[2]]\n", "\n", "print(c(max_lon, max_lat))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "#### Questions\n", "\n", "Identify the location with the lowest solar radiation" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "### Solar Radiation to Power Conversion\n", "\n", "Let's convert from solar radiation to solar power, we'll use the standard conversion equation (further details in [Filik et al](https://www.e3s-conferences.org/articles/e3sconf/pdf/2017/10/e3sconf_asee2017_00046.pdf)).\n", "\n", "$$E(t) = A.r(t).g(t).p(t)$$\n", "\n", "* $A$: Total solar panel area in $m^{2}$\n", "* $r$: Panel efficiency (yield) parameter\n", "* $G$: Solar radiation on tilted panels\n", "* $p$: Performance ratio of the overall system\n", "\n", "We'll base our calculations on a 10 $m^{2}$ solar panel with 17.5% efficiency and a performace ratio of 60%. \n", "\n", "Inspecting the attributes of the solar radiation grid tells us the units are in joules per $m^{2}$, assuming that the solar radiation stays constant for the entire hour we can view the units in terms of watt-hours which will make calculations easier." ] }, { "cell_type": "code", "execution_count": 223, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\t\n", "\n", "
idnamevariablevalue
0 scale_factor ssrd 336.1979
1 add_offset ssrd 12895679
2 _FillValue ssrd -32767
3 missing_valuessrd -32767
4 units ssrd J m**-2
5 long_name ssrd Surface solar radiation downwards
6 standard_name ssrd surface_downwelling_shortwave_flux_in_air
\n" ], "text/latex": [ "\\begin{tabular}{r|llll}\n", " id & name & variable & value\\\\\n", "\\hline\n", "\t 0 & scale\\_factor & ssrd & 336.1979 \\\\\n", "\t 1 & add\\_offset & ssrd & 12895679 \\\\\n", "\t 2 & \\_FillValue & ssrd & -32767 \\\\\n", "\t 3 & missing\\_value & ssrd & -32767 \\\\\n", "\t 4 & units & ssrd & J m**-2 \\\\\n", "\t 5 & long\\_name & ssrd & Surface solar radiation downwards\\\\\n", "\t 6 & standard\\_name & ssrd & surface\\_downwelling\\_shortwave\\_flux\\_in\\_air\\\\\n", "\\end{tabular}\n" ], "text/markdown": [ "\n", "| id | name | variable | value |\n", "|---|---|---|---|\n", "| 0 | scale_factor | ssrd | 336.1979 |\n", "| 1 | add_offset | ssrd | 12895679 |\n", "| 2 | _FillValue | ssrd | -32767 |\n", "| 3 | missing_value | ssrd | -32767 |\n", "| 4 | units | ssrd | J m**-2 |\n", "| 5 | long_name | ssrd | Surface solar radiation downwards |\n", "| 6 | standard_name | ssrd | surface_downwelling_shortwave_flux_in_air |\n", "\n" ], "text/plain": [ " id name variable value \n", "1 0 scale_factor ssrd 336.1979 \n", "2 1 add_offset ssrd 12895679 \n", "3 2 _FillValue ssrd -32767 \n", "4 3 missing_value ssrd -32767 \n", "5 4 units ssrd J m**-2 \n", "6 5 long_name ssrd Surface solar radiation downwards \n", "7 6 standard_name ssrd surface_downwelling_shortwave_flux_in_air" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ncmeta::nc_atts(nc_filepath, 'ssrd')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "We'll create a conversion function that calculates the kWh output for a theoretical solar panel, in this case using data from between 8 and 9 AM at the location with highest solar radiation." ] }, { "cell_type": "code", "execution_count": 212, "metadata": {}, "outputs": [ { "data": { "text/html": [ "6.02215665332479" ], "text/latex": [ "6.02215665332479" ], "text/markdown": [ "6.02215665332479" ], "text/plain": [ "[1] 6.022157" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "radiation_to_power <- function(radiation, area=10, yield=0.175, pr=0.6, hours=1){\n", " kWh <- radiation * area * yield * pr * (hours/3600) / 1000\n", " return(kWh)\n", "}\n", "\n", "max_kWh <- radiation_to_power(max_rad)\n", "\n", "max_kWh" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "#### Questions\n", "\n", "Calculate the difference between this and the location with lowest solar radiation" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "### Wind Speed to Power Conversion\n", "\n", "Wind vectors are a useful way to communicate speed and direction, in particular they allow both direction to be expressed in terms of a cyclical variable (unlike degrees which has a discontinuity) - this is often a more useful feature for ML models. However in vector form we cant convert to power output, for that we'll have to calculate the wind speed as an intermediate.\n", "\n", "To calculate u and v from wind speed and direction we can use:\n", "* $u = speed \\times cos(direction)$\n", "* $v = speed \\times sin(direction)$\n", "\n", "To convert back to wind speed we can use:\n", "* $speed = \\sqrt{u^{2}+v^{2}}$\n", "\n", "
\n", "\n", "We'll calculate wind speed and then visualise it" ] }, { "cell_type": "code", "execution_count": 229, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAMAAADKOT/pAAAAVFBMVEUAAABNTU1oaGh8fHx9\nACWMjIyampqiBwanp6eysrK9vb3DIgDHx8fQ0NDZ2dnhPADh4eHp6entYgDw8PDyhAD1oQD3\nujz40HT75Jr/9Lf//8j///+x2akoAAAACXBIWXMAABJ0AAASdAHeZh94AAAgAElEQVR4nO2d\niXbaOABFTTtpadq0CQED4f//c/CuzfvzBveeM1NisCXbutYuRzcAGE20dAQAHgFEAhCASAAC\nEAlAACIBCEAkAAGIBCAAkQAEIBKAAEQCEIBIAAIQCUAAIgEIQCQAAYgEIACRAAQgEoAARAIQ\ngEgAAhAJQAAiAQhAJAABiAQgAJEABCASgABEAhCASAACEAlAACIBCEAkAAGIBCAAkQAEIBKA\nAEQCEIBIAAIQCUAAIgEIQCQAAYgEIACRAAQgEoAARAIQgEgAAhAJQAAiAQhAJAABiAQgAJEA\nBCASgABEAhCASAACEAlAACIBCEAkAAGIBCAAkQAEIBKAAEQCEIBIAAIQCUAAIgEIQCQAAYgE\nIACRAAQgEoAARAIQgEgAAhAJQAAiAQhAJAABiAQgAJEABCASgABEAhCASAACEAlAACIBCEAk\nAAGIBCAAkQAEIBKAAEQCEIBIAAIQCUAAIgEIQCQAAYgEIACRAAQgEoAARAIQgEgAAhAJQAAi\nAQhAJAABiAQgAJEABCASgIAZRIoANsaAVK4XZ4EgAJQgEoAARAIQgEgAAhAJQAAiAQhAJAAB\niAQgAJEABCASgABEAhCASAACEAlAACIBCEAkAAGIBCAAkQAEIBKAAEQCEIBIAAIQCUAAIgEI\nQCQAAYgEIACRAAQgEoAARAIQgEgAAhAJQAAiAQhAJAABiAQgAJEABCASgABEAhCASAACEAlA\nACIBCEAkAAHzi/T6LYpe3iYNAmBuZhQpSnf8HqXsJwkCYCHmFmkf7T9vt4999DpFEAALMbdI\nu+gz+fwZfZsiCICFmFukKDL+kAcBsBBzi/SjEGk3RRAACzGrSC+/Xt+iP/ePn/vm1gZEgo0x\nq0gZ6cfd5xRBACzEnP1I7++vry8vaZPDvtEjRIKtwcgGAAGIBCAAkQAELCUS/UjwUKxHpMhE\nEQTAfFC0AxCASAACEAlAACIBCEAkAAGIBCBggUGrHVq4EQk2xowivSISPCyzjv7efZ86CIBl\nmLWO9N6yeJAgCIBFmLex4TV6nzoIgCWg1Q5AACIBCEAkAAGIBCAAkQAEIBKAAEQCEIBIAAIQ\nCUAAIgEIQCQAAYgEIACRAAQgEoAARAIQgEgAAhAJQAAiAQhAJAABiAQgAJEABCASgABEAhCA\nSAACEAlAACIBCEAkAAGIBCAAkQAEIBKAAEQCEIBIAAIQCUAAIgEIQCQAAYgEIACRAAQgEoAA\nRAIQgEgAAhAJQAAiAQhAJAABiAQgAJEABCASgABEAhCASAACEAlAACIBCEAkAAGIBCAAkQAE\nIBKAAEQCEIBIAAIQCUAAIgEIQCQAAYgEIACRAAQgEoAARAIQgEgAAhAJQAAiAQhAJAABiAQg\nAJEABCASgABEAhCASAACEAlAACIBCEAkAAGIBCAAkQAEIBKAAEQCEIBIAAIQCUAAIgEIQCQA\nAYgEIACRAAQgEoAARAIQgEgAAhAJQAAiAQhAJAABiAQgAJEABCASgABEAhCASAACEAlAACIB\nCEAkAAGIBCAAkQAEIBKAAEQCEDCrSH9/vUQJL/u/UwUBsAgzivT5Lar4PkkQAAsxo0j7aPfn\nPf308baL9lMEAbAQM4q0i97Lz+/RboogABZiRpGiqO4PWRAAC0GOBCBg3jrS20f6iToSPBpz\nNn9/N1rtvn1OEgTAMszbj7RP+5F2L7/oR4LHgpENAAIQCUAAIgEIWEok+pHgoViPSJGJIggN\ncRyfz+dLxTX53/kcLx0xHdnppWd5P69RJ5Yepu7L6/Xrek0v3zX/N/l5/uXxOCbc5aFoV889\nUSWpyhUpSQRjE9yquJjcz+w05lANX95F+ioEyoK6/zc8rHWBSPUg0oBDNXyJSKN3WWEQPYg9\nke6pIbds81wdkjPruGuTNGG+7lxLjc7pE6r3QdYJIrVji5QX8M8PYtI1SduWSFMm7TwwKwd8\njEwJkToQKNulKSCOty7T15fhUf65p0h9xSsCvO+HSFODSDOBSCpmnY/UuYV79SJdr0Uxf9Mi\nuRpVIp06tjhcsqJgH5nKQLPW8v6xXiUzivS6gEgjGqAM8mpx3od0tU2ShLAQhkhfpUhZm2TH\nrHZYjaoIq0fLxtqZs2j3vmte8kQQhItapKvZcJdVluMxzcULUyNSVuCa/ilxpfl7GO/N0/kU\nQTggUiOIJGPexoZXY7b5REFMQTlKyGi4NbpDNmvSxRcpL3Fdim7aacPfdsHY4tFb7bJuv5EH\nqUS6+CQViq5V81VxPx2vrcFsfSgyp+lisM3rFgaROoBI07DN6xbmsUWKi2brsQcqKklBkVKV\nNjV6+Vw1Pl49k8y/+3crPSuPLZIxHnP4MY4ndyZFQKT7w3UzJiUnY+Y5IZGqz1mDysY7nqfn\n0UWq+nyGH+TUULQ752xLJLvh0VWp+FyV7ryO58dpJFCBSO0gEiK18uAiGQWYJLUYW25pua/r\ncc5F4qvSnS1SWk1SxXp6TIuCIlWzWYuztPanoOfy4CKlk8msdii7SaqbSadTXD7Frb2d9gZZ\nrCfHPgGnuaE6w/IxUe1JbSkMInUAkao9ESnMo4tkd5b4LnU8TLNIl1ykw0EW70koHxz2KRRn\nZFqUnVUcb3+myEw8h0gNJnXscMxFslqFHZHi08pFyjwKXw8/N6L/qBePLlLacFf12nsp6NLR\npGaR8iF3x+M/YcTVBIXxVSq/Wzq+2wKREAmRBDy+SEnprmYkWaZSV5HsZHjxTFq1SF8O6RnU\niFRcpqXjvCmeQKSyW/bqVw+yWZodVOog0jleq0ieJsUZ1OVJ5U/oeO3KU4h0Mxt4/SypS+ku\njltESkakrVIk+xFQnyOFMqWHWZpkehAJkRBJwHOIdPM6HU2fupgUx+eL2dYQMCkZb7c+kc49\nRMp+Vq1Ncl2iDXxD3domTyJSpZJMJL+WdBfp8Pu3PuqjyKceuSLVdK2l3xT9tUtE93g8bnOy\nHyIhEiIJeBqRbkUjlZ+SLp1KMPH5mg2JDmuUizQ6mv+NPoJFvmZLwKQ6jb5ui5ToSijaCZko\nCLui1CtHuqVp0s7GAo0No0X67+fYI9gUHvkq1Xh0zRaI1MaiD4gkZKog7CRlWtHh5RKxOUO7\nRqSR8fvvv/9+/v6tc8lSv6tIstCfCkRCJEQS8Iwi+Ukp9aCtQHPyplIUCuUuKlaX+vn73z9Z\ne8Wlj0jVZlXwT8VziVRN9HOSUm5S875x7C5jpRfpdhdp/EFS2kcAOReh+HPZPth1T0WpBZEQ\nCZEEPJtIKVe/mnTJVWou3uULTroqSUUaf4iMprHdtSKlZzKfSP9kT43FeUqR7OHQZt7SbtI5\nYFIl0noep3Z+5Hv0ZahTsEA013PBRvKcItkqGSLlS3nX75cV75w2h/WJ5EYxqFFgEO/QN4cN\nZy1XbCyIhEiIJOBZRfJLd4ZIDZWdOA6U7lYnkjmcIVyw863Kr8Lcq6iu5IqN5plFCudJ5xaV\nTl7jXfpHsmf3ZYunHdrqrYvfwaQvppeP4mlFshcDsbKkbD23OpVOeZ5UabQ+kXpZZPk0abwe\nGURCJEQS8MQi3cxpFaUTl2JN/No2h2QB40qkiyFSHK9g6PKlZtpEN5EwaSDPLFKmg5mMSo/y\nilJoFOohIFIxarX1tcw/xbMkTM5GK8gQjQqVWGJ1CIiESIgk4JlFumXloDT9lCJVJp3D77M8\nHI/xuRoIPpVI//WbKlvEoV9znffTK6W7gTyFSA1vNj+fr+Ub6r7ct4fV1ZNO+Rv8bJFams31\nnG7FW9urF9z2Einw23lN2tQbrJt5CpEaTLIm61kDHMoXA/k3+5i33LkipW/AnM8k4xXRRkw6\ni3QNSkeWNAxEQiREEvAcIt2ySkRgc7aCakCl5tLdqSzd2b+PZ31xX+5N1e5RrufdyaRNqbT2\nFV+fRqQk2QWypTaR6l9pWYy5szOxura+KTCqaUb0R4u0SpPqa7krAZEQCZEEPJFIt/R2OFvi\ndN3HplrSub5N22wuM0p3s1WTAiLVuNHPpPsXM51BR5Jxwog0gNn6kW5J4TtPjldzxqjpUt6A\ncDx6WU0+O8nJk7oPuRvD6RTKkVQirStPyh5YS8eikacVyXrAFc1dZgIz9MiWvLv75BsSm7lS\nOUpvnrKduxRLb4/M+X1rF2nZdZTbQaQEREKkkTyrSE6B2++CufrtDlU30cFa/eZs/3CmJnB/\ncTChSNNHvyvp9b6LtHQ82nhWkQI442vCJqV50/F4MGdIl7lSmSNpRQo3XlR1pLEiBbYroz+G\n1nWkVwMiGVjvbbhaA1INpe6312nhjvNcaRKRkkEUoe2ZSGGTeogU3q6L/ijOm/EIkUwQCZGG\ngkgmdlIMiJTrEugqMkt/SQO4TKVMUX97aGGw3iLVC7aKxobVd8NWIJKDKVKNSTV9g7Ghkkyk\nOM6C9L8pmwunMGl6kbpkNedLhxXZVwIiOSASIg0BkSyqdNkkUrFEuDvUIV09Uji6IR1jnrpi\nVNGydHU+W1MoBrpU+7P0CONPoJY09u2/uhpPLfqRBjDTmg2BTdX0pC+nXzZg0smrC8WZSyKR\nYnNOux2ZmxGZ4SY1izShStlFbP3VtcyNO/1+SZ5WpLpSg5EsG0QyxtXZu98zkTjvbBofR789\nwbKnRiRDhcEi5V9OkguUnW7tZ5/nSHnei0hrDAKREEnJ04qUFy/8rWaKdBJv4Y/5R6DSLJtG\nEWzhNrY0iVQ/srurSPlzpO0i9i4A5k+Bjmd/vuUzgRFpADOttNoi0vXLzQZcnLGvLmPftNA0\nnM6IWtCj+nF0vURqSfHl3JPuJ3VNWka7WJGdfRFK6F6tCURyNyISIg3g2UXy7k8g3daLlI8Z\nOqXcxSnNSRrGk6GthzEuHV2RqtQ9h0hVGA1xrMK9dWqh7nLMnOTsb2VJkBmyQ5ipH2moScbG\n26UcFZTsnHuTNorf/6cQyXXk6i7SYMWqoyiVlS0WpSHUj3gLZ5f1p1RFussFSNUpZgH3vXoz\n89QiZX3n7rZQYSogUmFS3gWbzZxJvcnsOaSfxoh0steECIt0tdocGvMfU5OWrxMXiqMHs4LL\nLTSxsNg9fzydq5fA38zMq0/zRP89lgGREAmRBDy5SIF1hUyRvuzFUByTrIl/WfEnNacSaRSH\nk71IUVVHchNwGamuGrWKdGtLutdscEXwOOalcUTr3jKRlpQvuUhjL+UMPLtI92qIu8keJhRW\nyTfpLtL93h9HZUE2pUhZnleK5CfQfhlS+w/bYuZmQYEAjO4uM8xrh6PnnMo2wbUtDhbk2UUK\nDEO2lgOvTLLzJOvNl1li7zSguQ9F0S7vAy5fPB4ue3YVqV6Arh7dWtv6qqsU+Lb7BbiUZUUa\nG4aBSAmIhEgjWXA+UsIlVE26uLgm3c7ipTqO9vqT16yH2PaojF7lfbtGzb9rj1i+/60q5CV1\nKiPR520wnkj9PLrle9+u9CMNZGGRnG7QGpHKH6Uf87XF21/a151jOR/JCWy4SF1+1xKrzOiu\n59DP0MDe5WdEGgAiJSASIo1kaZHs6XTXepNKsplnWpEOtkjFiAovYgNEavphS6yyePQ6kUES\npaVkK6BVryiESAbG/Ic0U3Ie90ZyDtt0rn/t7NAYxX6eJBGp9rdtMcoa4oWn2Mh23rGOSAYn\n88WVVTNZZ5HSF5vrRfJXZ6gVq94gX6Twr1ujNKtIl9W31hUgkgEiIdJQEKkBJ5XW2eMU7owD\njJ5unotkLbxlhtdVJHfebNNvW6J0GTApdgjJdey3FteiziFSQSAj8R74pjJ1IlnHHDvrvFgr\nL5Ad1Znka2JvaRGpLUaXWTOJLiZlK3vl12OGOIVBpAJEQqQRIFKJPyyhqQm8RiT7AIlIo4p3\n6cJ2DSJ5EawVydWqTqW2GCX7hb8Yc541xHUimc+nfK3J7Mos1tmESCWB8T0tU5MCejW9u3lA\njOKzUUXKwgxbVPzZYSxds0ltcaoVSUqRk3dbkbUY03VDpPmD6IiZ6jxtfLnsxobx63LFjka2\nxYUWvlgdRRqQJzWK1MXEbucdF6szdy2uhaKVn6QkRm0gUjOIFLggtV8i0sS7rDAIj5oBKNcA\ngfpJIZK173iR7JENvkhudJoECokU2qM5Rs0ijT1fg8ulWPV74AHsuyOMWRBEaqVOpNDCCY5I\niuAbRKpC7pAR1Yjk79V+NbpHvrBhAMWuA0cJ+XdsYDw6gkjt1FkTWEtYHvbxePMLd9diIp8l\nUqM5TSI5ezbG55pNRuoQ82rVsgGnnTT8XNM5TrdheVIxe8taLiA9VF2D+si3AyJSO4hUgUg1\nIFI7+Yq5fgL03ZK/Pfh4PLm1JE/ltupRYJpqvUiNsSl/3PIr2/vhZ1+E1JLE01kv2cdqpZjA\nYplZT1PwxQflS4CHRhWROmAk16rxLuTRZeRjzedwON5ib1ZSoKbUdUk7XzNz1/YLUfvDr+o3\nVQznGQhhdDaZ1yV0qmnEam/RiNgiUgcQqboQiBQGkUxqPLBLQW0ipSU8YaSyGefVu2RrRKpb\nY66DSkUia4uJGYgz3bxKqcXjJv9iruFvadzShSvN9fQCJjkrEtldFJlJQ56FiGTSSaSgSnat\nQDpX9nCvJ2WjhbJEUoagEqmqi7dQ3yphPGMC305PtS5ybXZk3Lh6kZKLWlOHagGRTGqz/eDa\nq1fj/rhFr3IR4/H8+/cv+edUvOTFzPosmeuSTheRuow5re9wcmKRp9U0gxh+3v2KWelCEi1P\nlfxquC45DCucI5IJIjWBSA0gksm5rrYZWnvVuj9OnSVP65JI/c7/PZ1KkfxJSvUitWuUv32i\nFW+3fHsVgUKk7CJmhs61jpaxBkCdSIbo8gF4iGRxrhu3f61NsoY67jBWbQrKRbrUijS05a5b\nuqrb2+hdcgxLvhy1RGrfRrRLOU2sNXtGpGlBpFoQqRFEcqgRyU6ztkOGO7ZI2pbf47FcCryc\nXHF1GSRSxwjU7e4Fm//+ehszb7YYn9C/FJbv6UfRvEI9jxm3reuOSBb5nNTQV0URvFYkaw3U\nkEjZi2WP3eae//7929nSSaT+OnVPqfVHCHpUijQwRyqX4B8go1VPqi5JdW16Hi/pyWtWCZFc\n4rjm6VOUq2yRjPTsz1Bypyed7KXzenEwREpUsls3BpvUI5UqjCxH4Haxq0vc6guOxXUwL0dL\nhlTO7Tf6ktLkkPUvnRvyJURyQaQ6EGlVIr1+i6KXt0mDGEddSj9XIlXrlprN3l6idlJLpmg8\n0KRUpGrQXXhF8GJTD5G610EUByorPqlLwwYRmEc71/UyVPWkwHOl6ZD5c6oYmGKPvM8cC0V6\nRpGidMfvUcp+kiCmJe/EMZ504cmyVbKunpenU1wxLHhbJIMsHpZI3fMkjUi38KIJFcWFyCTP\n3jKRXJ74dK8zJtXGjlVHg3wN2trpRcWNCuTQLed8MUYIBq53UP+5RdpH+8/b7WMfvU4RxLRY\nvnQSqXpcliLVTYfpQsNSKBdrMOsUIjUfpm1vQyT/y64NMKHD1hftzCzJKfL2aAkMXupAtjS3\nSLvoM/n8GX2bIohpQSREWo1IUWT8IQ9CiX93q/K96YqvkbVypFO4GzULs6wJN+hUxmc6kUIH\nHzXmZvxqSyGcDmvzsiRfd+3lq73I9iNxbpF+FCLtpghCSVAkc263dX9qRXIqw/HwCc3H4+nc\n0aQeIhUjTNvQGekzeon0Gu5X6WatCp4OpS1VatzzVgxzCV3b8ilZTZWZVaSXX69v0Z/7x899\nc2sDIgUjhEj9eFiRMtKPu88pglDij0EIvlnWdsq71G4BYoRIWf/6UmU72YFmx6uStovkeeOL\nlHUiVpWlOfuR3t9fX19e0iaHfaNHqxAphKVJnUcX9++i9yFt+B5TTeouUt2wzfEiNa+GPOi0\nbqKlNJspfMqj2vjb6ha6GtnPqOK7+/1gZEMvAiOu63IkO10b9pyHZ0pGz7tKpG4lOyPtWLs6\nOd86s6SU7JKnS+UVItU+QrxHY51IRsstIvUCkRAJkYbglDnOoVlArknhO1DaY7QYWCF1qHJX\nIpX3tE4jX6SaOTqdLoPxe+N4eVgbEClwVWpFumYVoIZHlDM7+rqkSOvvR8rxR/RcjDF3PUS6\nFGNLSg/sw2Y9tqfGrsnT6RR491gXkdy8o9rcegWcPTxr7Cd09ys7OadT8XA6B67I7XYLxvaS\nd7Y3XturVU+6rkqkyEQRhAY/n7DfMd6llFdsOZetqqkKdjhxB5GSlsRUpuYxDqHnpi9Rnqu0\nXoHwflUiXKtIxWASY/yqdxb+XtmFsz0yn5EFjkcU7dpAJERCJAn+QLDA7NSANs5vLplKxSFS\nqYxKUjEWr7Ge9C+LT2FS1nHoGlXV4+yEH/ZhWGeslQS9Y3Y4ZC/MC3krXrvcYYhPnA2Xv9ac\nfRpbdx/jdhUPq+Cj0ryEiNSJw8Hdcj57fbMha2yVylRuNt5V9a8iR2rvUbnHJ8mVsj/i/ECW\nR+X9NzQKpSVn9d5arH3S/5vfekfucsjOxPlzyz6H4oLe2rqgrOsQwt0h9CisLqFrUvXTWUX6\n++slrQG97P9OFcQEIBIirUqkz29Ga8L3SYKYCE+luGW0kFnScrZVjd5OT9KpamTqidmMdzZH\niVm3O5gIuojkyJf93/jeb2cfcA637B1mNUU2X4fyojbO/0737SySO9m45vHjX86ZRdpHuz/v\n6aePt936B60aeCKlywd30qg0ydhQH9CweTmGSEWGl6Sti+d6FaE8PXScSRAQydjTT2z9TXKy\n82Jz1thvrq/l+ZSdddOxw3v6IvnXq7aFJsCMIu2i9/Lz+/qnUTQSpwtxNwtkrJplI47LKXsk\nBzp58xWYDaPLT7kO3Yt2RbIq9iwD8heIGHISZkG0jFTWB3QuVqpL/jKS961M0y1n0eKFH4XG\nny8vktU7tJkO2TCIhEjkSCLc0UIhjWyTst/JI3LMq1ZeDauoPVRVpjSeVTJJpol2WIslj7Y1\nqq5MvH6SG3IOZk0mqMbFmqTurt/f+ES4uE0VtWKYG+o0Mn5v/jVvHentI/20tTpSGDt1WvKc\nyw9mG9o0Gt2rcO0/SVdm856i6crcHUQqstFrVdsoziScQAechL2GWNfrlAffdmGNh1jYkU7j\nqWyP7L3nFKlYiSvl2+on9rWCSDdEWkSk29992o+0e/m1pX6kBopWcEcjF7MIuFxk3Yb4+6a+\nwxAMCauinZfSBscvKFI2LKo9Wk0/KOqzt3xuZq0kHUSq+ZKRDeMwXqJXJ9LlsgaRjBypaiW/\n9kv2Zkt+1arhJ6qBUazmUJlb0wUkW7oFmhtNilcdFm325/rJWsaToqdhiDQO4wURYZGsdoZr\nty7QaTB6tKoWt35iVw0pdnnQT2BDIzncwvqvTqfYbi89h0p4XsHX6x1rMAyRxoJICYiESCM5\nxVYDd9gie1j4UlEtImVu6hmb6iydFoqrNbTva21zZfObZGw5X9zWjUISx5eOM/YRaSTH/GXj\nVRK7BClFWi6qgfX7+y+eXHjk1v6d58W6PPJEyq7Dxajf+ZlNcwaESBNgvCTiXI1mNZOW8Wm5\naEqWvDoHM6QU36T1GeXSWJjrpRIijQeRUhCpJ4jkcIqt4l1ZigstL7igSBLO4TpSglW4y3uo\n1uVRINJ2d7kjByLNzekU2w10ZhvQI3lUvcAu/I2R9FafG5XEpkyeSZ00QiQNiIRI/cNFpACn\nuGwHd8rdj1W0K8uvzvZQEWmRCA7BXr/WXiGoThz3C0TScKraHAx78o+mVktHdCRlrutsN0cU\nFtu2drrpiL50CU7jfiHS/NhtDqHM6X5flo7kKKrTcs+j9MjYtjWRcg6HkzHM2MydrAXOEGky\nEAmR+oBIdZycAUOXfMiJecGXjuMo6kUqu5Hyv4qpP1s0qXpZwaUaR9Rc0EMkMcfEJVskp2i9\ndAzHYFT+Ag+EzKTid9e0E2nDp5uvDm0sa4ZI84FIxe8QqRVEaqQoYBsiXY33Ny8dvVGUL0Wq\nq+tl5/dVirRp/v07HIpByeYSSuZzsSzxIZKcfBlWo0ei0OohROrQaJKe59bPtSB7J5Ux+Mn3\nCJGmIi7ml9vXe+siVdN16jKk6ncPI9Itff+HvYqhc0+TLYg0BYiESO0gUjtx1dH/SCJV6ab2\nB+k/+RtaN366JsablsIg0gj+/ftX91W5ck1i0oOIZFUUgr/IxhIay8dl57uKEYa1t6obZR+h\nI1DxNyKN4F+9Sfl1959cs0ZQizW8PahSJVKxmOTMUaxnpEe3/BUkgddiZdsQaQSI5IBIfUCk\nbiTXverEzIfXLB2pURjrXd5aHgmre2aERKovmdfgvYKk6ntHpMlI+h9uWeorXwqxbZFuyas5\nO6wVntB3ybwJufsSLDocD/5bTZvxihmINAOItBIQafMUqS5fu9pZnhFmIVx+Ow56Ze/9hpbV\nRESaj6N7q/qvyNiL//6b9vgWPZ/ma+NehW1+i3OIYz7lDJFmpvU9ChtmfDvYPNQ0rWZTjvoe\n7HQyZpyZy3Eg0qQg0vIg0oMSn9MpmMWfkgVQl+D30hEYRyrSkP3Cq7sj0swkj7ScTKa4f5V3\nHWxapEO6flpw6eUW3OUEEGkZEGkVINLWqRYjTF1KfUpc2l5d6ufSERjDKbvyQx5hp+C7exBp\nZuLYuPpV3pTKtK28adMiJc+zeNDj63CM44BKiDQzwXJBIVMmVPvdPRxW0GT2CCIN6Qw7HI/+\nCoaINDeItAoQaevcRboGVDKLeae2zqfj8bjxYQWLk78HdNBlPOSluzMizU+egxQzWi7FP/7w\nx+xJ2cDpNKxwDyXpykBJi93AjP2YvcfnjEgLcM9mjMvvz/gzXpPZ3FGYDhI7JbnSYRVFvA1y\nz1PGPYnyTgxEWgBEWg+ItF0O+RK4cdF46quUkq8mWX+gpD+w7H06pfWpvtNqnp6GNQI6UNaR\nEGlxnCnLzmsX7QH6iVZJrSgfAJt3rBtdUFnu9O/373SswdKAD8AAAA4SSURBVM//Kn7+/O2P\nP0i+mPl0H4tDPiwCkZYHkTYMIq2HpBJ0KetK/ouAS5Vy0bJeplNcvanWHGeUu5QUV34n/EzI\nRbp/WvpkH4+i1c5wCZEWw8xyLI3CtSab/CaaY4yS5ox/95wpI7Mpk2iqHCiQ3a2BGRpgjsfY\nMQmRlsIpHNg5UvjtIWVJMDHJVyuXKSl7JCW9pU9wMUa1JHTicDyeYrt4h0hLgUhTgUhzBrE4\n2TI2xkhiz6TCraKIZzrllfVylTY1fnwSZprdX3RkFLcPkZbjkLwQzmxAcMypNjpNEclP7PpV\n/vMnX+0ryYgmnyeZdObmHYJmnoRIy4FIYhBpgSDWRHlTnJKaOwQiUGuqGs7zJcafmKRWmCTw\nCUVK3ypbWGSV7RBpLRyzCqxNaFxeoIHc6nta+kSGI8hMTtPmSAfnJlWdeYi0IpLxeBnJX0mC\ncBfacDMnM58q8iirfHfZRmEviWXRITb8KP/SQtdkHuXvODdnYSLSKkEkRNLypCKl+K9IuN+m\nW9WNfnWHEnnVpktijzGFcOBqOfORzt3O56yOIknqkhgFMTqPfJUQaTNkt+3i15L8Xls7/0qa\nLZaOfCPJo6PKiTvx+3dofNKkXbHH08nMgpzaLCJtiFMcePuiK5KpVPnF0jFvZNgYjJ/+bkG5\nZFga2S7dy5OItCEQyQCR1hHEZnHuplGCy8tzRp3JqEvle3esjAsakSW1/n99hgxOOl/kZNSP\nCo1u1aKeR0TaHOkqQjmxe3ezSpSXXRUvR85/03D0pEY9tvUspe4Il2J8bUZjpSaZDtI9m5l0\n5lXVDVt4dMsqdtnXB0TaHIhUAyItEcTmSVvJ05JFckuNzvb7X9f0ZciGSbebXY+qDhMeC5Hr\nljabh1ejbH2TQ/rm6czerJ8o3Xa7FS30K29JdDAH2AVb6hHpkbDaj8/lwl9fmTjmoNfsJ0UF\nq9whMTL/WNan/VdlZE0exUJG4ajEqSiXWyll6no2LbjtJPq/bnxyilyxNmKI9MiUrXzlliKf\nqv6qtMqbK4qcIhfxci4zn7QYdjyawy2CRcU4W+I8seZW7Z2USZsim4h5zl6RvMFpVYj0yCDS\nbCDSo5PNzyj+sgtWRbdTVexLB0KY32ct6+lfaanmkIt0rSORp+Mg9GoXazjhJgewI9JTkGRM\n+UejoeFWddeaU52yX55MX8pM55CNOCv7r5yxSV/FkkftUbpm6yOZBpZ9YsIznwtEegoQaWoQ\n6dkwUmlhQbrZNOmSL6OS6HKzWstz4jhryTYbAcvWwW5Yg9cLk8ad2pIg0lOTWpQkX2s5MHdS\nUwNWa3rPoM1JIFvNiEoQCVKKxR/mDMxchGKWcCcEkSAFkcaBSDA/F3dl2aUjNB5Egnmx51R9\n3bZdNSpBJOjHmGHhl2IO1WZbuWtBJOgHIgVBJOjEtUj7g6c/+MuKSSO4MIgE7RjDFrLe2sHH\nKJHHcWEQCdpBpFYQCdqpun1SlQbMbrUtSg40RTwXBJGgA9cAQ/Z/pJ4jG0SCPph5Sg+VrrZJ\n08VvMRAJ+oBINSAS9GKYD9dquhEijdhlhUHAYMy6Tvddyi6kSeO2GIgEfUGkAIgEvammtXbe\npVxwaFvrQnYHkWAA/ZrAs4W7xq+CvGYQCSYmfuSMqASRYGIQSbfLCoMAUIJIICOc7TxyzagC\nkUBGSKSHbmEwQCSQgUhT77LCIGAOHr6NoQSRYCqeobGuBJFgKhBJvssKg4DJeRqJEhAJQAAi\nAQiYVaS/v16ihJf936mCAFiEGUX6/BZVfJ8kCHhCjkfjj/IlzofDMfvmME8sZhRpH+3+vKef\nPt520X6KIOAJMV6WXr1C/ZBIlE7cmOn96DOKtIvey8/v0W6KIOAJeTqRoqjuD1kQ8Hwkr4kp\n1ua/VMvAnk755vj+8dhyDI/fv/vuQY4EmybODco0Mt99m2dFcaLaafJ8ad460ttH+ok6EqjI\nXlp2KTQ6e2Nk8wxr6njM2fz93Wi1+/Y5SRDwbDyjSLe/+7Qfaffyi34kUJC/RrPBo9spL/tN\nHBNGNsCGqVoakkWKgusUxWXtKcuVJqouIRI8OPmrn41cq7btYUTnLSLBg/PYItGPBDNhi5Sa\nFBbp2Lu7yWA9IkUmiiAA5oOiHYAARAIQgEgAApjYByCAiX0AApjYByCAaRQAApjYByCAHAlA\nABP7AAQwsQ9AABP7AAQwsgFAACIBCEAkAAGIBCAAkQAEIBKAAEQCEIBIAAIQCUAAIgEIQCQA\nAYgEIACRAAQgEoAARAIQgEgAAhAJQAAiAQhAJAABiAQgAJEABCASgABEAhCASAACEAlAACIB\nCEAkAAGIBCAAkQAEIBKAAEQCEIBIAAIQCUAAIgEIQCQAAYgEIACRAAQgEoAARAIQgEgAAhAJ\nQAAiAQhAJAABiAQgAJEABCASgABEAhCASAACEAlAACIBCEAkAAGIBCAAkQAEIBKAAEQCEIBI\nAAIQCUAAIgEIQCQAAYgEIACRAAQgEoAARAIQgEgAAhAJQAAiAQhAJAABiAQgAJEABCASgABE\nAhCASAACEAlAACIBCEAkAAGIBCAAkQAEIBKAAEQCEIBIAAIQCUAAIgEIQCQAAYgEIACRAAQg\nEoAARAIQgEgAAhAJQAAiAQhAJAABiAQgAJEABCASgABEAhCASAACEAlAACIBCEAkAAHzi/T6\nLYpe3iYNAmBuZhQpSnf8HqXsJwkCYCHmFmkf7T9vt4999DpFEAALMbdIu+gz+fwZfZsiCICF\nmFukKDL+kAcBsBBzi/SjEGk3RRAACzGrSC+/Xt+iP/ePn/vm1gZEgo0xq0gZ6cfd5xRBACzE\nnP1I7++vry8vaZPDvtEjRIKtwcgGAAGIBCAAkQAELCUS/UjwUCASgACKdgACEAlAwHpEikym\nCQJgKmYV6e+vl1STl/3fqYIAWIQZRfr8ZmQ53ycJAmAhZhRpH+3+vKefPt52DFqFh2JGkXbR\ne/n5nWkU8FDMvmZD6A9ZEAALQY4EIGDeOtLbR/qJOhI8GnM2f383Wu2+MbEPHol5+5H2aT/S\n7uUX/UjwWKxnZMPMQQAoQSQAAUuI1D6UDpFgYyASgABEAhCASAACEAlAACIBCKD5G0AAIgEI\nQCQAAYgEIACRAAQgEoAARAIQgEgAAhAJQAAiAQhAJAABiAQgAJEABCASgABEAhCASAACEAlA\nACIBCEAkAAGIBCAAkQAEIBKAAEQCEIBIAAIQCUAAIgEIQCQAAYgEIACRAAQgEoAARAIQgEgA\nAhAJQAAiAQhAJAABKxUJYGMMSOV6cRYJoxvExGM1Edl4TBBpGVYTk9VEZOMxQaRlWE1MVhOR\njccEkZZhNTFZTUQ2HhNEWobVxGQ1Edl4TBBpGVYTk9VEZOMxQaRlWE1MVhORjccEkZZhNTFZ\nTUQ2HhNEWobVxGQ1Edl4TBBpGVYTk9VEZOMxQaRlWE1MVhORjccEkZZhNTFZTUQ2HpP1xB5g\nwyASgABEAhCASAACEAlAACIBCEAkAAGIBCAAkQAEIBKAAEQCEIBIAAIQCUAAIgEIQCQAAYgE\nIGAykfa7aLf/bNowF17Ar9/WEpM7f5d4lHkRef8RRT8+VhCTz8WSyT1d2LeiX0ymuo3f00X9\nvzVsmAsv4H26YTf/zQpdgs/dAiJ5EXlbyyX52GUxWcLpd/slFD0T7ES38W+0e7+976K/tRvm\nwgv4PfrxmTx9fiwek4SXIa8QkUdkd9/w+RLtF4/JjzQO+/lvzi2Jhnkr+ibYiW7jPnq7//9P\n9Kt2w1x4Ab9kpzx/Ag5dgj+D3sWjjsifNPl+RrvFYxItdXPuT9bvVqh9E+xEMX6Jksz5PXqp\n3TAXdQHPf68CMflw7t5CEfkRvc8eiXBM8oLu/Erf7s8S61b0TbAT3UbvybLYo6Ym4M/o+wpi\n8j36WEAkLyLfotuvXVriXTomv/Ki3ewFl9u7k0b6JthnFek1zbkXjsmv6M8SpZjAzXlJq/jL\nx+T2mrQ27F5nj4kTDURqj0nKx272MqYfk7TcsA6RksaGH/PnA6FnS8L8GZITDURqj0nC5272\ngl2oRJW0N69DpKSO9DF/54QXk9ekaHdXepEsaY0i7dxoeBvmIhjw9wX6s7yY/EhLlwuI5F2S\nxZ5yXky+RUlF7XOJ/kbn/Psm2Elb7T7cVruPpVrtrIA/vn1fosPPjcmYl9FLI7Jcj4AXk+Wa\nv91Q+ybYiWL8K33avlVdfN6GufADfpu/wS4Yk8VEqrk5H/NfFy8mWT6wQI9WgnUj+ibY5xvZ\nsEB6qYlJyhpGNtxrR59JzeTP4jHZR8notv38z9uENY5suBd2E9Ikm8XP2DAvbkx+LJQPBK6J\n/WnBiPxay83JR7gt86grbsWgBDvVbcxG8RrxMjbMixuTpQpUgWtif1oyIm/f13Fz8jHXC8Tk\n5orUM8EuUqsDeDQQCUAAIgEIQCQAAYgEIACRAAQgEoAARAIQgEgAAhAJQAAiAQhAJAABiAQg\nAJEABCASgABEAhCASAACEAlAACIBCEAkAAGIBCAAkQAEIBKAAEQCEIBIAAIQCUAAIgEIQCQA\nAYgEIACRAAQgEoAARAIQgEgAAhAJQAAiAQhAJAABiAQgAJEABCASgABEAhCASAACEAlAACIB\nCEAkAAGIBCAAkQAEIBKAAEQCEIBIAAIQCUAAIgEIQCQAAYgEIACRAAQgEoAARAIQgEgAAhAJ\nQAAiAQhAJAABiAQgAJEABCASgABEAhCASAACEAlAACIBCEAkAAGIBCAAkQAEIBKAAEQCEIBI\nAAIQCUAAIgEIQCQAAYgEIACRAAQgEoCA/wFyGrQjT00gVgAAAABJRU5ErkJggg==", "text/plain": [ "plot without title" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "u_slice <- nc_file %>% \n", " activate(u10) %>% \n", " hyper_filter(time=time<1043145) %>% \n", " hyper_array\n", "\n", "v_slice <- nc_file %>% \n", " activate(v10) %>% \n", " hyper_filter(time=time<1043145) %>% \n", " hyper_array\n", "\n", "u_matrix <- u_slice[[1]]\n", "v_matrix <- v_slice[[1]]\n", "\n", "speed_matrix <- sqrt(u_matrix**2 + v_matrix**2)\n", "speed_matrix <- speed_matrix[,length(latitudes):1]\n", "\n", "image(speed_matrix, asp=lat_lon_ratio)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "Again we'll use the maximum value of the resource as an example" ] }, { "cell_type": "code", "execution_count": 230, "metadata": {}, "outputs": [ { "data": { "text/html": [ "11.7539165645299" ], "text/latex": [ "11.7539165645299" ], "text/markdown": [ "11.7539165645299" ], "text/plain": [ "[1] 11.75392" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "max_speed <- max(speed_matrix, na.rm=TRUE)\n", "\n", "max_speed" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "We'll use the Helman equation to extrapolate from 10m up to 70m which is a more common height for wind turbines, to do this we'll use an estimated value for the hellman exponent - 1/7. Helman Equation:\n", "\n", "speed at turbine = data speed x (turbine height/data height)^hellman exponent\n", "\n", "If you are able to obtain surface roughness data for indonesia you can calculate a more accurate helman exponent using: \n", "\n", "hellman exponent = 1 / ln(turbine height/roughness)" ] }, { "cell_type": "code", "execution_count": 232, "metadata": {}, "outputs": [ { "data": { "text/html": [ "15.520685364153" ], "text/latex": [ "15.520685364153" ], "text/markdown": [ "15.520685364153" ], "text/plain": [ "[1] 15.52069" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "calc_turbine_speed <- function(speed, turbine_height=70, data_height=10, hellman_exponent=1/7){\n", " turbine_speed <- speed*(turbine_height/data_height)^hellman_exponent\n", " return(turbine_speed)\n", "}\n", "\n", "turbine_speed <- calc_turbine_speed(max_speed)\n", "\n", "turbine_speed" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "We now need to convert the wind speed into wind power, to do so we can use a power curve. In this example we'll use the power curve of the Gamsea G114, the turbine used in the Indonesian Sidrap wind farm. The data used for this can be found [here](https://www.thewindpower.net/turbine_en_599_gamesa_g114-2500.php)." ] }, { "cell_type": "code", "execution_count": 262, "metadata": {}, "outputs": [ { "data": { "text/html": [ "16: 2500" ], "text/latex": [ "\\textbf{16:} 2500" ], "text/markdown": [ "**16:** 2500" ], "text/plain": [ " 16 \n", "2500 " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "speed_to_power <- function(speed){\n", " pc <- c(0, 0, 0, 0, 0, 40, 102, 160, 257, 355, 471, 605, 750, 972, 1194, 1455, 1713, 1970, 2227, 2380, 2465, 2495, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2500, 2460, 2400, 2335, 2250)\n", " names(pc) <- seq(from=0, to=24, by=0.5)\n", " power <- pc[toString(round(speed, 0))]\n", " return(power)\n", "}\n", "\n", "power <- speed_to_power(turbine_speed)\n", "\n", "power" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "#### Questions\n", "\n", "Now try and calculate another layer that shows wind direction, how could this be useful in your analysis?\n", "For additional reading please look [here](http://colaweb.gmu.edu/dev/clim301/lectures/wind/wind-uv). You should also try to calculate power output for an offshore wind farm using a power curve of an offshore turbine." ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "### Financial Modelling\n", "\n", "Now that we can calcuate the power output for solar and wind farms we'll do some financial modelling with them.\n", "\n", "First we'll calculate the annual revenue of the wind farm, we'll use the number of turbines and capacity from the Sidrap wind farm. We'll assume a capacity factor of 30%, lifetime of 25 years and electricity price of $0.07/kWh." ] }, { "cell_type": "code", "execution_count": 298, "metadata": {}, "outputs": [ { "data": { "text/html": [ "13797000" ], "text/latex": [ "13797000" ], "text/markdown": [ "13797000" ], "text/plain": [ "[1] 13797000" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "turbine_capacity <- 2.5\n", "num_turbines <- 30\n", "capacity_factor <- 0.3\n", "dollars_per_kWh <- 0.07\n", "lifetime_yrs <- 25\n", "\n", "kWh <- turbine_capacity * num_turbines * capacity_factor * 1000 * 24 * 365\n", "annual_revenue <-dollars_per_kWh * kWh\n", "\n", "annual_revenue" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "We can then use this to calcute the Net Present Value (NPV).\n", "\n", "Here we'll assume the Sidrap capital expenditure (CAPEX) of 150,000,000. We'll assume an operating cost of 0 although in reality this would cost something. We'll assume a discount rate of 5%." ] }, { "cell_type": "code", "execution_count": 299, "metadata": {}, "outputs": [ { "data": { "text/html": [ "51597010" ], "text/latex": [ "51597010" ], "text/markdown": [ "51597010" ], "text/plain": [ "[1] 51597010" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "calc_NPV <- function(annual_revenue, i=0.05, lifetime_yrs=25, CAPEX=150000000, OPEX=0){\n", " costs <- rep(OPEX, lifetime_yrs)\n", " costs[1] <- costs[1] + CAPEX\n", " revenue <- rep(annual_revenue, lifetime_yrs)\n", "\n", " t <- seq(1, lifetime_yrs, 1)\n", " NPV <- sum((revenue - costs)/(1 + i)**t)\n", "\n", " return(round(NPV, 0))\n", "}\n", "\n", "calc_NPV(annual_revenue)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "#### Questions\n", "\n", "Using the assumptions provided in the coursework instructions you should try to estimate the net present value of an ammonia plant. Assume a lifetime of 25 years." ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# " ] } ], "metadata": { "kernelspec": { "display_name": "ESDA-Spatial", "language": "R", "name": "ir" }, "language_info": { "codemirror_mode": "r", "file_extension": ".r", "mimetype": "text/x-r-source", "name": "R", "pygments_lexer": "r", "version": "3.6.1" } }, "nbformat": 4, "nbformat_minor": 4 }