{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "fb5ca344",
"metadata": {
"ExecuteTime": {
"end_time": "2022-11-28T02:35:25.994025Z",
"start_time": "2022-11-28T02:35:25.953Z"
}
},
"outputs": [],
"source": [
"# This is not working\n",
"# The working solution is through xarray open_mfdataset, extend time as dimension, and assign time value"
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "f9c66c19",
"metadata": {},
"outputs": [],
"source": [
"library('ncdf4')\n",
"library('stringi')\n",
"library('stringr')"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "251affe0",
"metadata": {},
"outputs": [],
"source": [
"setwd('/bog/amuttaqin/Datasets/TRMM/')"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "eb986b8c",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"'/bog/amuttaqin/Datasets/TRMM'"
],
"text/latex": [
"'/bog/amuttaqin/Datasets/TRMM'"
],
"text/markdown": [
"'/bog/amuttaqin/Datasets/TRMM'"
],
"text/plain": [
"[1] \"/bog/amuttaqin/Datasets/TRMM\""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"getwd()"
]
},
{
"cell_type": "code",
"execution_count": 33,
"id": "4d89d444",
"metadata": {},
"outputs": [],
"source": [
"nc.files <- list.files(pattern= \"3B42.20121001.+nc4\")"
]
},
{
"cell_type": "code",
"execution_count": 34,
"id": "3278cb85",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"- '3B42.20121001.00.7.HDF.nc4'
- '3B42.20121001.03.7.HDF.nc4'
- '3B42.20121001.06.7.HDF.nc4'
- '3B42.20121001.09.7.HDF.nc4'
- '3B42.20121001.12.7.HDF.nc4'
- '3B42.20121001.15.7.HDF.nc4'
- '3B42.20121001.18.7.HDF.nc4'
- '3B42.20121001.21.7.HDF.nc4'
\n"
],
"text/latex": [
"\\begin{enumerate*}\n",
"\\item '3B42.20121001.00.7.HDF.nc4'\n",
"\\item '3B42.20121001.03.7.HDF.nc4'\n",
"\\item '3B42.20121001.06.7.HDF.nc4'\n",
"\\item '3B42.20121001.09.7.HDF.nc4'\n",
"\\item '3B42.20121001.12.7.HDF.nc4'\n",
"\\item '3B42.20121001.15.7.HDF.nc4'\n",
"\\item '3B42.20121001.18.7.HDF.nc4'\n",
"\\item '3B42.20121001.21.7.HDF.nc4'\n",
"\\end{enumerate*}\n"
],
"text/markdown": [
"1. '3B42.20121001.00.7.HDF.nc4'\n",
"2. '3B42.20121001.03.7.HDF.nc4'\n",
"3. '3B42.20121001.06.7.HDF.nc4'\n",
"4. '3B42.20121001.09.7.HDF.nc4'\n",
"5. '3B42.20121001.12.7.HDF.nc4'\n",
"6. '3B42.20121001.15.7.HDF.nc4'\n",
"7. '3B42.20121001.18.7.HDF.nc4'\n",
"8. '3B42.20121001.21.7.HDF.nc4'\n",
"\n",
"\n"
],
"text/plain": [
"[1] \"3B42.20121001.00.7.HDF.nc4\" \"3B42.20121001.03.7.HDF.nc4\"\n",
"[3] \"3B42.20121001.06.7.HDF.nc4\" \"3B42.20121001.09.7.HDF.nc4\"\n",
"[5] \"3B42.20121001.12.7.HDF.nc4\" \"3B42.20121001.15.7.HDF.nc4\"\n",
"[7] \"3B42.20121001.18.7.HDF.nc4\" \"3B42.20121001.21.7.HDF.nc4\""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"nc.files"
]
},
{
"cell_type": "code",
"execution_count": 35,
"id": "5c9556db",
"metadata": {},
"outputs": [],
"source": [
"ncin <- nc_open(nc.files[1])"
]
},
{
"cell_type": "code",
"execution_count": 36,
"id": "65321538",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"File 3B42.20121001.00.7.HDF.nc4 (NC_FORMAT_NETCDF4_CLASSIC):\n",
"\n",
" 1 variables (excluding dimension variables):\n",
" float precipitation[nlat,nlon] (Chunking: [201,560]) (Compression: level 4)\n",
" units: mm/hr\n",
" coordinates: nlon nlat\n",
" _FillValue: -9999.900390625\n",
"\n",
" 2 dimensions:\n",
" nlon Size:560 \n",
" long_name: longitude\n",
" standard_name: longitude\n",
" units: degrees_east\n",
" nlat Size:201 \n",
" long_name: latitude\n",
" standard_name: latitude\n",
" units: degrees_north\n",
"\n",
" 5 global attributes:\n",
" Grid.GridHeader: BinMethod=ARITHMETIC_MEAN;\n",
"Registration=CENTER;\n",
"LatitudeResolution=0.25;\n",
"LongitudeResolution=0.25;\n",
"NorthBoundingCoordinate=50;\n",
"SouthBoundingCoordinate=-50;\n",
"EastBoundingCoordinate=180;\n",
"WestBoundingCoordinate=-180;\n",
"Origin=SOUTHWEST;\n",
"\n",
" FileHeader: AlgorithmID=3B42;\n",
"AlgorithmVersion=3B42_7.0;\n",
"FileName=3B42.20121001.00.7.HDF;\n",
"GenerationDateTime=2013-01-17T13:42:19.000Z;\n",
"StartGranuleDateTime=2012-09-30T22:30:00.000Z;\n",
"StopGranuleDateTime=2012-10-01T01:29:59.999Z;\n",
"GranuleNumber=;\n",
"NumberOfSwaths=0;\n",
"NumberOfGrids=1;\n",
"GranuleStart=;\n",
"TimeInterval=3_HOUR;\n",
"ProcessingSystem=PPS;\n",
"ProductVersion=7;\n",
"MissingData=;\n",
"\n",
" FileInfo: DataFormatVersion=m;\n",
"TKCodeBuildVersion=1;\n",
"MetadataVersion=m;\n",
"FormatPackage=HDF Version 4.2 Release 4, January 25, 2009;\n",
"BlueprintFilename=TRMM.V7.3B42.blueprint.xml;\n",
"BlueprintVersion=BV_13;\n",
"TKIOVersion=1.6;\n",
"MetadataStyle=PVL;\n",
"EndianType=LITTLE_ENDIAN;\n",
"\n",
" GridHeader: BinMethod=ARITHMETIC_MEAN;\n",
"Registration=CENTER;\n",
"LatitudeResolution=0.25;\n",
"LongitudeResolution=0.25;\n",
"NorthBoundingCoordinate=50;\n",
"SouthBoundingCoordinate=-50;\n",
"EastBoundingCoordinate=180;\n",
"WestBoundingCoordinate=-180;\n",
"Origin=SOUTHWEST;\n",
"\n",
" history: 2022-11-25 05:33:50 GMT Hyrax-1.16.3 https://disc2.gesdisc.eosdis.nasa.gov/opendap/TRMM_L3/TRMM_3B42.7/2012/274/3B42.20121001.00.7.HDF.nc4?precipitation[880:1439][99:299],nlon[880:1439],nlat[99:299]"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"ncin"
]
},
{
"cell_type": "code",
"execution_count": 37,
"id": "73c0ef10",
"metadata": {},
"outputs": [],
"source": [
"lon <- ncvar_get(nc=ncin, \"nlon\")"
]
},
{
"cell_type": "code",
"execution_count": 38,
"id": "a668f86e",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"- 40.125
- 40.375
- 40.625
- 40.875
- 41.125
- 41.375
- 41.625
- 41.875
- 42.125
- 42.375
- 42.625
- 42.875
- 43.125
- 43.375
- 43.625
- 43.875
- 44.125
- 44.375
- 44.625
- 44.875
- 45.125
- 45.375
- 45.625
- 45.875
- 46.125
- 46.375
- 46.625
- 46.875
- 47.125
- 47.375
- 47.625
- 47.875
- 48.125
- 48.375
- 48.625
- 48.875
- 49.125
- 49.375
- 49.625
- 49.875
- 50.125
- 50.375
- 50.625
- 50.875
- 51.125
- 51.375
- 51.625
- 51.875
- 52.125
- 52.375
- 52.625
- 52.875
- 53.125
- 53.375
- 53.625
- 53.875
- 54.125
- 54.375
- 54.625
- 54.875
- 55.125
- 55.375
- 55.625
- 55.875
- 56.125
- 56.375
- 56.625
- 56.875
- 57.125
- 57.375
- 57.625
- 57.875
- 58.125
- 58.375
- 58.625
- 58.875
- 59.125
- 59.375
- 59.625
- 59.875
- 60.125
- 60.375
- 60.625
- 60.875
- 61.125
- 61.375
- 61.625
- 61.875
- 62.125
- 62.375
- 62.625
- 62.875
- 63.125
- 63.375
- 63.625
- 63.875
- 64.125
- 64.375
- 64.625
- 64.875
- 65.125
- 65.375
- 65.625
- 65.875
- 66.125
- 66.375
- 66.625
- 66.875
- 67.125
- 67.375
- 67.625
- 67.875
- 68.125
- 68.375
- 68.625
- 68.875
- 69.125
- 69.375
- 69.625
- 69.875
- 70.125
- 70.375
- 70.625
- 70.875
- 71.125
- 71.375
- 71.625
- 71.875
- 72.125
- 72.375
- 72.625
- 72.875
- 73.125
- 73.375
- 73.625
- 73.875
- 74.125
- 74.375
- 74.625
- 74.875
- 75.125
- 75.375
- 75.625
- 75.875
- 76.125
- 76.375
- 76.625
- 76.875
- 77.125
- 77.375
- 77.625
- 77.875
- 78.125
- 78.375
- 78.625
- 78.875
- 79.125
- 79.375
- 79.625
- 79.875
- 80.125
- 80.375
- 80.625
- 80.875
- 81.125
- 81.375
- 81.625
- 81.875
- 82.125
- 82.375
- 82.625
- 82.875
- 83.125
- 83.375
- 83.625
- 83.875
- 84.125
- 84.375
- 84.625
- 84.875
- 85.125
- 85.375
- 85.625
- 85.875
- 86.125
- 86.375
- 86.625
- 86.875
- 87.125
- 87.375
- 87.625
- 87.875
- 88.125
- 88.375
- 88.625
- 88.875
- 89.125
- 89.375
- 89.625
- 89.875
- ⋯
- 130.125
- 130.375
- 130.625
- 130.875
- 131.125
- 131.375
- 131.625
- 131.875
- 132.125
- 132.375
- 132.625
- 132.875
- 133.125
- 133.375
- 133.625
- 133.875
- 134.125
- 134.375
- 134.625
- 134.875
- 135.125
- 135.375
- 135.625
- 135.875
- 136.125
- 136.375
- 136.625
- 136.875
- 137.125
- 137.375
- 137.625
- 137.875
- 138.125
- 138.375
- 138.625
- 138.875
- 139.125
- 139.375
- 139.625
- 139.875
- 140.125
- 140.375
- 140.625
- 140.875
- 141.125
- 141.375
- 141.625
- 141.875
- 142.125
- 142.375
- 142.625
- 142.875
- 143.125
- 143.375
- 143.625
- 143.875
- 144.125
- 144.375
- 144.625
- 144.875
- 145.125
- 145.375
- 145.625
- 145.875
- 146.125
- 146.375
- 146.625
- 146.875
- 147.125
- 147.375
- 147.625
- 147.875
- 148.125
- 148.375
- 148.625
- 148.875
- 149.125
- 149.375
- 149.625
- 149.875
- 150.125
- 150.375
- 150.625
- 150.875
- 151.125
- 151.375
- 151.625
- 151.875
- 152.125
- 152.375
- 152.625
- 152.875
- 153.125
- 153.375
- 153.625
- 153.875
- 154.125
- 154.375
- 154.625
- 154.875
- 155.125
- 155.375
- 155.625
- 155.875
- 156.125
- 156.375
- 156.625
- 156.875
- 157.125
- 157.375
- 157.625
- 157.875
- 158.125
- 158.375
- 158.625
- 158.875
- 159.125
- 159.375
- 159.625
- 159.875
- 160.125
- 160.375
- 160.625
- 160.875
- 161.125
- 161.375
- 161.625
- 161.875
- 162.125
- 162.375
- 162.625
- 162.875
- 163.125
- 163.375
- 163.625
- 163.875
- 164.125
- 164.375
- 164.625
- 164.875
- 165.125
- 165.375
- 165.625
- 165.875
- 166.125
- 166.375
- 166.625
- 166.875
- 167.125
- 167.375
- 167.625
- 167.875
- 168.125
- 168.375
- 168.625
- 168.875
- 169.125
- 169.375
- 169.625
- 169.875
- 170.125
- 170.375
- 170.625
- 170.875
- 171.125
- 171.375
- 171.625
- 171.875
- 172.125
- 172.375
- 172.625
- 172.875
- 173.125
- 173.375
- 173.625
- 173.875
- 174.125
- 174.375
- 174.625
- 174.875
- 175.125
- 175.375
- 175.625
- 175.875
- 176.125
- 176.375
- 176.625
- 176.875
- 177.125
- 177.375
- 177.625
- 177.875
- 178.125
- 178.375
- 178.625
- 178.875
- 179.125
- 179.375
- 179.625
- 179.875
\n"
],
"text/latex": [
"\\begin{enumerate*}\n",
"\\item 40.125\n",
"\\item 40.375\n",
"\\item 40.625\n",
"\\item 40.875\n",
"\\item 41.125\n",
"\\item 41.375\n",
"\\item 41.625\n",
"\\item 41.875\n",
"\\item 42.125\n",
"\\item 42.375\n",
"\\item 42.625\n",
"\\item 42.875\n",
"\\item 43.125\n",
"\\item 43.375\n",
"\\item 43.625\n",
"\\item 43.875\n",
"\\item 44.125\n",
"\\item 44.375\n",
"\\item 44.625\n",
"\\item 44.875\n",
"\\item 45.125\n",
"\\item 45.375\n",
"\\item 45.625\n",
"\\item 45.875\n",
"\\item 46.125\n",
"\\item 46.375\n",
"\\item 46.625\n",
"\\item 46.875\n",
"\\item 47.125\n",
"\\item 47.375\n",
"\\item 47.625\n",
"\\item 47.875\n",
"\\item 48.125\n",
"\\item 48.375\n",
"\\item 48.625\n",
"\\item 48.875\n",
"\\item 49.125\n",
"\\item 49.375\n",
"\\item 49.625\n",
"\\item 49.875\n",
"\\item 50.125\n",
"\\item 50.375\n",
"\\item 50.625\n",
"\\item 50.875\n",
"\\item 51.125\n",
"\\item 51.375\n",
"\\item 51.625\n",
"\\item 51.875\n",
"\\item 52.125\n",
"\\item 52.375\n",
"\\item 52.625\n",
"\\item 52.875\n",
"\\item 53.125\n",
"\\item 53.375\n",
"\\item 53.625\n",
"\\item 53.875\n",
"\\item 54.125\n",
"\\item 54.375\n",
"\\item 54.625\n",
"\\item 54.875\n",
"\\item 55.125\n",
"\\item 55.375\n",
"\\item 55.625\n",
"\\item 55.875\n",
"\\item 56.125\n",
"\\item 56.375\n",
"\\item 56.625\n",
"\\item 56.875\n",
"\\item 57.125\n",
"\\item 57.375\n",
"\\item 57.625\n",
"\\item 57.875\n",
"\\item 58.125\n",
"\\item 58.375\n",
"\\item 58.625\n",
"\\item 58.875\n",
"\\item 59.125\n",
"\\item 59.375\n",
"\\item 59.625\n",
"\\item 59.875\n",
"\\item 60.125\n",
"\\item 60.375\n",
"\\item 60.625\n",
"\\item 60.875\n",
"\\item 61.125\n",
"\\item 61.375\n",
"\\item 61.625\n",
"\\item 61.875\n",
"\\item 62.125\n",
"\\item 62.375\n",
"\\item 62.625\n",
"\\item 62.875\n",
"\\item 63.125\n",
"\\item 63.375\n",
"\\item 63.625\n",
"\\item 63.875\n",
"\\item 64.125\n",
"\\item 64.375\n",
"\\item 64.625\n",
"\\item 64.875\n",
"\\item 65.125\n",
"\\item 65.375\n",
"\\item 65.625\n",
"\\item 65.875\n",
"\\item 66.125\n",
"\\item 66.375\n",
"\\item 66.625\n",
"\\item 66.875\n",
"\\item 67.125\n",
"\\item 67.375\n",
"\\item 67.625\n",
"\\item 67.875\n",
"\\item 68.125\n",
"\\item 68.375\n",
"\\item 68.625\n",
"\\item 68.875\n",
"\\item 69.125\n",
"\\item 69.375\n",
"\\item 69.625\n",
"\\item 69.875\n",
"\\item 70.125\n",
"\\item 70.375\n",
"\\item 70.625\n",
"\\item 70.875\n",
"\\item 71.125\n",
"\\item 71.375\n",
"\\item 71.625\n",
"\\item 71.875\n",
"\\item 72.125\n",
"\\item 72.375\n",
"\\item 72.625\n",
"\\item 72.875\n",
"\\item 73.125\n",
"\\item 73.375\n",
"\\item 73.625\n",
"\\item 73.875\n",
"\\item 74.125\n",
"\\item 74.375\n",
"\\item 74.625\n",
"\\item 74.875\n",
"\\item 75.125\n",
"\\item 75.375\n",
"\\item 75.625\n",
"\\item 75.875\n",
"\\item 76.125\n",
"\\item 76.375\n",
"\\item 76.625\n",
"\\item 76.875\n",
"\\item 77.125\n",
"\\item 77.375\n",
"\\item 77.625\n",
"\\item 77.875\n",
"\\item 78.125\n",
"\\item 78.375\n",
"\\item 78.625\n",
"\\item 78.875\n",
"\\item 79.125\n",
"\\item 79.375\n",
"\\item 79.625\n",
"\\item 79.875\n",
"\\item 80.125\n",
"\\item 80.375\n",
"\\item 80.625\n",
"\\item 80.875\n",
"\\item 81.125\n",
"\\item 81.375\n",
"\\item 81.625\n",
"\\item 81.875\n",
"\\item 82.125\n",
"\\item 82.375\n",
"\\item 82.625\n",
"\\item 82.875\n",
"\\item 83.125\n",
"\\item 83.375\n",
"\\item 83.625\n",
"\\item 83.875\n",
"\\item 84.125\n",
"\\item 84.375\n",
"\\item 84.625\n",
"\\item 84.875\n",
"\\item 85.125\n",
"\\item 85.375\n",
"\\item 85.625\n",
"\\item 85.875\n",
"\\item 86.125\n",
"\\item 86.375\n",
"\\item 86.625\n",
"\\item 86.875\n",
"\\item 87.125\n",
"\\item 87.375\n",
"\\item 87.625\n",
"\\item 87.875\n",
"\\item 88.125\n",
"\\item 88.375\n",
"\\item 88.625\n",
"\\item 88.875\n",
"\\item 89.125\n",
"\\item 89.375\n",
"\\item 89.625\n",
"\\item 89.875\n",
"\\item ⋯\n",
"\\item 130.125\n",
"\\item 130.375\n",
"\\item 130.625\n",
"\\item 130.875\n",
"\\item 131.125\n",
"\\item 131.375\n",
"\\item 131.625\n",
"\\item 131.875\n",
"\\item 132.125\n",
"\\item 132.375\n",
"\\item 132.625\n",
"\\item 132.875\n",
"\\item 133.125\n",
"\\item 133.375\n",
"\\item 133.625\n",
"\\item 133.875\n",
"\\item 134.125\n",
"\\item 134.375\n",
"\\item 134.625\n",
"\\item 134.875\n",
"\\item 135.125\n",
"\\item 135.375\n",
"\\item 135.625\n",
"\\item 135.875\n",
"\\item 136.125\n",
"\\item 136.375\n",
"\\item 136.625\n",
"\\item 136.875\n",
"\\item 137.125\n",
"\\item 137.375\n",
"\\item 137.625\n",
"\\item 137.875\n",
"\\item 138.125\n",
"\\item 138.375\n",
"\\item 138.625\n",
"\\item 138.875\n",
"\\item 139.125\n",
"\\item 139.375\n",
"\\item 139.625\n",
"\\item 139.875\n",
"\\item 140.125\n",
"\\item 140.375\n",
"\\item 140.625\n",
"\\item 140.875\n",
"\\item 141.125\n",
"\\item 141.375\n",
"\\item 141.625\n",
"\\item 141.875\n",
"\\item 142.125\n",
"\\item 142.375\n",
"\\item 142.625\n",
"\\item 142.875\n",
"\\item 143.125\n",
"\\item 143.375\n",
"\\item 143.625\n",
"\\item 143.875\n",
"\\item 144.125\n",
"\\item 144.375\n",
"\\item 144.625\n",
"\\item 144.875\n",
"\\item 145.125\n",
"\\item 145.375\n",
"\\item 145.625\n",
"\\item 145.875\n",
"\\item 146.125\n",
"\\item 146.375\n",
"\\item 146.625\n",
"\\item 146.875\n",
"\\item 147.125\n",
"\\item 147.375\n",
"\\item 147.625\n",
"\\item 147.875\n",
"\\item 148.125\n",
"\\item 148.375\n",
"\\item 148.625\n",
"\\item 148.875\n",
"\\item 149.125\n",
"\\item 149.375\n",
"\\item 149.625\n",
"\\item 149.875\n",
"\\item 150.125\n",
"\\item 150.375\n",
"\\item 150.625\n",
"\\item 150.875\n",
"\\item 151.125\n",
"\\item 151.375\n",
"\\item 151.625\n",
"\\item 151.875\n",
"\\item 152.125\n",
"\\item 152.375\n",
"\\item 152.625\n",
"\\item 152.875\n",
"\\item 153.125\n",
"\\item 153.375\n",
"\\item 153.625\n",
"\\item 153.875\n",
"\\item 154.125\n",
"\\item 154.375\n",
"\\item 154.625\n",
"\\item 154.875\n",
"\\item 155.125\n",
"\\item 155.375\n",
"\\item 155.625\n",
"\\item 155.875\n",
"\\item 156.125\n",
"\\item 156.375\n",
"\\item 156.625\n",
"\\item 156.875\n",
"\\item 157.125\n",
"\\item 157.375\n",
"\\item 157.625\n",
"\\item 157.875\n",
"\\item 158.125\n",
"\\item 158.375\n",
"\\item 158.625\n",
"\\item 158.875\n",
"\\item 159.125\n",
"\\item 159.375\n",
"\\item 159.625\n",
"\\item 159.875\n",
"\\item 160.125\n",
"\\item 160.375\n",
"\\item 160.625\n",
"\\item 160.875\n",
"\\item 161.125\n",
"\\item 161.375\n",
"\\item 161.625\n",
"\\item 161.875\n",
"\\item 162.125\n",
"\\item 162.375\n",
"\\item 162.625\n",
"\\item 162.875\n",
"\\item 163.125\n",
"\\item 163.375\n",
"\\item 163.625\n",
"\\item 163.875\n",
"\\item 164.125\n",
"\\item 164.375\n",
"\\item 164.625\n",
"\\item 164.875\n",
"\\item 165.125\n",
"\\item 165.375\n",
"\\item 165.625\n",
"\\item 165.875\n",
"\\item 166.125\n",
"\\item 166.375\n",
"\\item 166.625\n",
"\\item 166.875\n",
"\\item 167.125\n",
"\\item 167.375\n",
"\\item 167.625\n",
"\\item 167.875\n",
"\\item 168.125\n",
"\\item 168.375\n",
"\\item 168.625\n",
"\\item 168.875\n",
"\\item 169.125\n",
"\\item 169.375\n",
"\\item 169.625\n",
"\\item 169.875\n",
"\\item 170.125\n",
"\\item 170.375\n",
"\\item 170.625\n",
"\\item 170.875\n",
"\\item 171.125\n",
"\\item 171.375\n",
"\\item 171.625\n",
"\\item 171.875\n",
"\\item 172.125\n",
"\\item 172.375\n",
"\\item 172.625\n",
"\\item 172.875\n",
"\\item 173.125\n",
"\\item 173.375\n",
"\\item 173.625\n",
"\\item 173.875\n",
"\\item 174.125\n",
"\\item 174.375\n",
"\\item 174.625\n",
"\\item 174.875\n",
"\\item 175.125\n",
"\\item 175.375\n",
"\\item 175.625\n",
"\\item 175.875\n",
"\\item 176.125\n",
"\\item 176.375\n",
"\\item 176.625\n",
"\\item 176.875\n",
"\\item 177.125\n",
"\\item 177.375\n",
"\\item 177.625\n",
"\\item 177.875\n",
"\\item 178.125\n",
"\\item 178.375\n",
"\\item 178.625\n",
"\\item 178.875\n",
"\\item 179.125\n",
"\\item 179.375\n",
"\\item 179.625\n",
"\\item 179.875\n",
"\\end{enumerate*}\n"
],
"text/markdown": [
"1. 40.125\n",
"2. 40.375\n",
"3. 40.625\n",
"4. 40.875\n",
"5. 41.125\n",
"6. 41.375\n",
"7. 41.625\n",
"8. 41.875\n",
"9. 42.125\n",
"10. 42.375\n",
"11. 42.625\n",
"12. 42.875\n",
"13. 43.125\n",
"14. 43.375\n",
"15. 43.625\n",
"16. 43.875\n",
"17. 44.125\n",
"18. 44.375\n",
"19. 44.625\n",
"20. 44.875\n",
"21. 45.125\n",
"22. 45.375\n",
"23. 45.625\n",
"24. 45.875\n",
"25. 46.125\n",
"26. 46.375\n",
"27. 46.625\n",
"28. 46.875\n",
"29. 47.125\n",
"30. 47.375\n",
"31. 47.625\n",
"32. 47.875\n",
"33. 48.125\n",
"34. 48.375\n",
"35. 48.625\n",
"36. 48.875\n",
"37. 49.125\n",
"38. 49.375\n",
"39. 49.625\n",
"40. 49.875\n",
"41. 50.125\n",
"42. 50.375\n",
"43. 50.625\n",
"44. 50.875\n",
"45. 51.125\n",
"46. 51.375\n",
"47. 51.625\n",
"48. 51.875\n",
"49. 52.125\n",
"50. 52.375\n",
"51. 52.625\n",
"52. 52.875\n",
"53. 53.125\n",
"54. 53.375\n",
"55. 53.625\n",
"56. 53.875\n",
"57. 54.125\n",
"58. 54.375\n",
"59. 54.625\n",
"60. 54.875\n",
"61. 55.125\n",
"62. 55.375\n",
"63. 55.625\n",
"64. 55.875\n",
"65. 56.125\n",
"66. 56.375\n",
"67. 56.625\n",
"68. 56.875\n",
"69. 57.125\n",
"70. 57.375\n",
"71. 57.625\n",
"72. 57.875\n",
"73. 58.125\n",
"74. 58.375\n",
"75. 58.625\n",
"76. 58.875\n",
"77. 59.125\n",
"78. 59.375\n",
"79. 59.625\n",
"80. 59.875\n",
"81. 60.125\n",
"82. 60.375\n",
"83. 60.625\n",
"84. 60.875\n",
"85. 61.125\n",
"86. 61.375\n",
"87. 61.625\n",
"88. 61.875\n",
"89. 62.125\n",
"90. 62.375\n",
"91. 62.625\n",
"92. 62.875\n",
"93. 63.125\n",
"94. 63.375\n",
"95. 63.625\n",
"96. 63.875\n",
"97. 64.125\n",
"98. 64.375\n",
"99. 64.625\n",
"100. 64.875\n",
"101. 65.125\n",
"102. 65.375\n",
"103. 65.625\n",
"104. 65.875\n",
"105. 66.125\n",
"106. 66.375\n",
"107. 66.625\n",
"108. 66.875\n",
"109. 67.125\n",
"110. 67.375\n",
"111. 67.625\n",
"112. 67.875\n",
"113. 68.125\n",
"114. 68.375\n",
"115. 68.625\n",
"116. 68.875\n",
"117. 69.125\n",
"118. 69.375\n",
"119. 69.625\n",
"120. 69.875\n",
"121. 70.125\n",
"122. 70.375\n",
"123. 70.625\n",
"124. 70.875\n",
"125. 71.125\n",
"126. 71.375\n",
"127. 71.625\n",
"128. 71.875\n",
"129. 72.125\n",
"130. 72.375\n",
"131. 72.625\n",
"132. 72.875\n",
"133. 73.125\n",
"134. 73.375\n",
"135. 73.625\n",
"136. 73.875\n",
"137. 74.125\n",
"138. 74.375\n",
"139. 74.625\n",
"140. 74.875\n",
"141. 75.125\n",
"142. 75.375\n",
"143. 75.625\n",
"144. 75.875\n",
"145. 76.125\n",
"146. 76.375\n",
"147. 76.625\n",
"148. 76.875\n",
"149. 77.125\n",
"150. 77.375\n",
"151. 77.625\n",
"152. 77.875\n",
"153. 78.125\n",
"154. 78.375\n",
"155. 78.625\n",
"156. 78.875\n",
"157. 79.125\n",
"158. 79.375\n",
"159. 79.625\n",
"160. 79.875\n",
"161. 80.125\n",
"162. 80.375\n",
"163. 80.625\n",
"164. 80.875\n",
"165. 81.125\n",
"166. 81.375\n",
"167. 81.625\n",
"168. 81.875\n",
"169. 82.125\n",
"170. 82.375\n",
"171. 82.625\n",
"172. 82.875\n",
"173. 83.125\n",
"174. 83.375\n",
"175. 83.625\n",
"176. 83.875\n",
"177. 84.125\n",
"178. 84.375\n",
"179. 84.625\n",
"180. 84.875\n",
"181. 85.125\n",
"182. 85.375\n",
"183. 85.625\n",
"184. 85.875\n",
"185. 86.125\n",
"186. 86.375\n",
"187. 86.625\n",
"188. 86.875\n",
"189. 87.125\n",
"190. 87.375\n",
"191. 87.625\n",
"192. 87.875\n",
"193. 88.125\n",
"194. 88.375\n",
"195. 88.625\n",
"196. 88.875\n",
"197. 89.125\n",
"198. 89.375\n",
"199. 89.625\n",
"200. 89.875\n",
"201. ⋯\n",
"202. 130.125\n",
"203. 130.375\n",
"204. 130.625\n",
"205. 130.875\n",
"206. 131.125\n",
"207. 131.375\n",
"208. 131.625\n",
"209. 131.875\n",
"210. 132.125\n",
"211. 132.375\n",
"212. 132.625\n",
"213. 132.875\n",
"214. 133.125\n",
"215. 133.375\n",
"216. 133.625\n",
"217. 133.875\n",
"218. 134.125\n",
"219. 134.375\n",
"220. 134.625\n",
"221. 134.875\n",
"222. 135.125\n",
"223. 135.375\n",
"224. 135.625\n",
"225. 135.875\n",
"226. 136.125\n",
"227. 136.375\n",
"228. 136.625\n",
"229. 136.875\n",
"230. 137.125\n",
"231. 137.375\n",
"232. 137.625\n",
"233. 137.875\n",
"234. 138.125\n",
"235. 138.375\n",
"236. 138.625\n",
"237. 138.875\n",
"238. 139.125\n",
"239. 139.375\n",
"240. 139.625\n",
"241. 139.875\n",
"242. 140.125\n",
"243. 140.375\n",
"244. 140.625\n",
"245. 140.875\n",
"246. 141.125\n",
"247. 141.375\n",
"248. 141.625\n",
"249. 141.875\n",
"250. 142.125\n",
"251. 142.375\n",
"252. 142.625\n",
"253. 142.875\n",
"254. 143.125\n",
"255. 143.375\n",
"256. 143.625\n",
"257. 143.875\n",
"258. 144.125\n",
"259. 144.375\n",
"260. 144.625\n",
"261. 144.875\n",
"262. 145.125\n",
"263. 145.375\n",
"264. 145.625\n",
"265. 145.875\n",
"266. 146.125\n",
"267. 146.375\n",
"268. 146.625\n",
"269. 146.875\n",
"270. 147.125\n",
"271. 147.375\n",
"272. 147.625\n",
"273. 147.875\n",
"274. 148.125\n",
"275. 148.375\n",
"276. 148.625\n",
"277. 148.875\n",
"278. 149.125\n",
"279. 149.375\n",
"280. 149.625\n",
"281. 149.875\n",
"282. 150.125\n",
"283. 150.375\n",
"284. 150.625\n",
"285. 150.875\n",
"286. 151.125\n",
"287. 151.375\n",
"288. 151.625\n",
"289. 151.875\n",
"290. 152.125\n",
"291. 152.375\n",
"292. 152.625\n",
"293. 152.875\n",
"294. 153.125\n",
"295. 153.375\n",
"296. 153.625\n",
"297. 153.875\n",
"298. 154.125\n",
"299. 154.375\n",
"300. 154.625\n",
"301. 154.875\n",
"302. 155.125\n",
"303. 155.375\n",
"304. 155.625\n",
"305. 155.875\n",
"306. 156.125\n",
"307. 156.375\n",
"308. 156.625\n",
"309. 156.875\n",
"310. 157.125\n",
"311. 157.375\n",
"312. 157.625\n",
"313. 157.875\n",
"314. 158.125\n",
"315. 158.375\n",
"316. 158.625\n",
"317. 158.875\n",
"318. 159.125\n",
"319. 159.375\n",
"320. 159.625\n",
"321. 159.875\n",
"322. 160.125\n",
"323. 160.375\n",
"324. 160.625\n",
"325. 160.875\n",
"326. 161.125\n",
"327. 161.375\n",
"328. 161.625\n",
"329. 161.875\n",
"330. 162.125\n",
"331. 162.375\n",
"332. 162.625\n",
"333. 162.875\n",
"334. 163.125\n",
"335. 163.375\n",
"336. 163.625\n",
"337. 163.875\n",
"338. 164.125\n",
"339. 164.375\n",
"340. 164.625\n",
"341. 164.875\n",
"342. 165.125\n",
"343. 165.375\n",
"344. 165.625\n",
"345. 165.875\n",
"346. 166.125\n",
"347. 166.375\n",
"348. 166.625\n",
"349. 166.875\n",
"350. 167.125\n",
"351. 167.375\n",
"352. 167.625\n",
"353. 167.875\n",
"354. 168.125\n",
"355. 168.375\n",
"356. 168.625\n",
"357. 168.875\n",
"358. 169.125\n",
"359. 169.375\n",
"360. 169.625\n",
"361. 169.875\n",
"362. 170.125\n",
"363. 170.375\n",
"364. 170.625\n",
"365. 170.875\n",
"366. 171.125\n",
"367. 171.375\n",
"368. 171.625\n",
"369. 171.875\n",
"370. 172.125\n",
"371. 172.375\n",
"372. 172.625\n",
"373. 172.875\n",
"374. 173.125\n",
"375. 173.375\n",
"376. 173.625\n",
"377. 173.875\n",
"378. 174.125\n",
"379. 174.375\n",
"380. 174.625\n",
"381. 174.875\n",
"382. 175.125\n",
"383. 175.375\n",
"384. 175.625\n",
"385. 175.875\n",
"386. 176.125\n",
"387. 176.375\n",
"388. 176.625\n",
"389. 176.875\n",
"390. 177.125\n",
"391. 177.375\n",
"392. 177.625\n",
"393. 177.875\n",
"394. 178.125\n",
"395. 178.375\n",
"396. 178.625\n",
"397. 178.875\n",
"398. 179.125\n",
"399. 179.375\n",
"400. 179.625\n",
"401. 179.875\n",
"\n",
"\n"
],
"text/plain": [
" [1] 40.125 40.375 40.625 40.875 41.125 41.375 41.625 41.875 42.125\n",
" [10] 42.375 42.625 42.875 43.125 43.375 43.625 43.875 44.125 44.375\n",
" [19] 44.625 44.875 45.125 45.375 45.625 45.875 46.125 46.375 46.625\n",
" [28] 46.875 47.125 47.375 47.625 47.875 48.125 48.375 48.625 48.875\n",
" [37] 49.125 49.375 49.625 49.875 50.125 50.375 50.625 50.875 51.125\n",
" [46] 51.375 51.625 51.875 52.125 52.375 52.625 52.875 53.125 53.375\n",
" [55] 53.625 53.875 54.125 54.375 54.625 54.875 55.125 55.375 55.625\n",
" [64] 55.875 56.125 56.375 56.625 56.875 57.125 57.375 57.625 57.875\n",
" [73] 58.125 58.375 58.625 58.875 59.125 59.375 59.625 59.875 60.125\n",
" [82] 60.375 60.625 60.875 61.125 61.375 61.625 61.875 62.125 62.375\n",
" [91] 62.625 62.875 63.125 63.375 63.625 63.875 64.125 64.375 64.625\n",
"[100] 64.875 65.125 65.375 65.625 65.875 66.125 66.375 66.625 66.875\n",
"[109] 67.125 67.375 67.625 67.875 68.125 68.375 68.625 68.875 69.125\n",
"[118] 69.375 69.625 69.875 70.125 70.375 70.625 70.875 71.125 71.375\n",
"[127] 71.625 71.875 72.125 72.375 72.625 72.875 73.125 73.375 73.625\n",
"[136] 73.875 74.125 74.375 74.625 74.875 75.125 75.375 75.625 75.875\n",
"[145] 76.125 76.375 76.625 76.875 77.125 77.375 77.625 77.875 78.125\n",
"[154] 78.375 78.625 78.875 79.125 79.375 79.625 79.875 80.125 80.375\n",
"[163] 80.625 80.875 81.125 81.375 81.625 81.875 82.125 82.375 82.625\n",
"[172] 82.875 83.125 83.375 83.625 83.875 84.125 84.375 84.625 84.875\n",
"[181] 85.125 85.375 85.625 85.875 86.125 86.375 86.625 86.875 87.125\n",
"[190] 87.375 87.625 87.875 88.125 88.375 88.625 88.875 89.125 89.375\n",
"[199] 89.625 89.875 90.125 90.375 90.625 90.875 91.125 91.375 91.625\n",
"[208] 91.875 92.125 92.375 92.625 92.875 93.125 93.375 93.625 93.875\n",
"[217] 94.125 94.375 94.625 94.875 95.125 95.375 95.625 95.875 96.125\n",
"[226] 96.375 96.625 96.875 97.125 97.375 97.625 97.875 98.125 98.375\n",
"[235] 98.625 98.875 99.125 99.375 99.625 99.875 100.125 100.375 100.625\n",
"[244] 100.875 101.125 101.375 101.625 101.875 102.125 102.375 102.625 102.875\n",
"[253] 103.125 103.375 103.625 103.875 104.125 104.375 104.625 104.875 105.125\n",
"[262] 105.375 105.625 105.875 106.125 106.375 106.625 106.875 107.125 107.375\n",
"[271] 107.625 107.875 108.125 108.375 108.625 108.875 109.125 109.375 109.625\n",
"[280] 109.875 110.125 110.375 110.625 110.875 111.125 111.375 111.625 111.875\n",
"[289] 112.125 112.375 112.625 112.875 113.125 113.375 113.625 113.875 114.125\n",
"[298] 114.375 114.625 114.875 115.125 115.375 115.625 115.875 116.125 116.375\n",
"[307] 116.625 116.875 117.125 117.375 117.625 117.875 118.125 118.375 118.625\n",
"[316] 118.875 119.125 119.375 119.625 119.875 120.125 120.375 120.625 120.875\n",
"[325] 121.125 121.375 121.625 121.875 122.125 122.375 122.625 122.875 123.125\n",
"[334] 123.375 123.625 123.875 124.125 124.375 124.625 124.875 125.125 125.375\n",
"[343] 125.625 125.875 126.125 126.375 126.625 126.875 127.125 127.375 127.625\n",
"[352] 127.875 128.125 128.375 128.625 128.875 129.125 129.375 129.625 129.875\n",
"[361] 130.125 130.375 130.625 130.875 131.125 131.375 131.625 131.875 132.125\n",
"[370] 132.375 132.625 132.875 133.125 133.375 133.625 133.875 134.125 134.375\n",
"[379] 134.625 134.875 135.125 135.375 135.625 135.875 136.125 136.375 136.625\n",
"[388] 136.875 137.125 137.375 137.625 137.875 138.125 138.375 138.625 138.875\n",
"[397] 139.125 139.375 139.625 139.875 140.125 140.375 140.625 140.875 141.125\n",
"[406] 141.375 141.625 141.875 142.125 142.375 142.625 142.875 143.125 143.375\n",
"[415] 143.625 143.875 144.125 144.375 144.625 144.875 145.125 145.375 145.625\n",
"[424] 145.875 146.125 146.375 146.625 146.875 147.125 147.375 147.625 147.875\n",
"[433] 148.125 148.375 148.625 148.875 149.125 149.375 149.625 149.875 150.125\n",
"[442] 150.375 150.625 150.875 151.125 151.375 151.625 151.875 152.125 152.375\n",
"[451] 152.625 152.875 153.125 153.375 153.625 153.875 154.125 154.375 154.625\n",
"[460] 154.875 155.125 155.375 155.625 155.875 156.125 156.375 156.625 156.875\n",
"[469] 157.125 157.375 157.625 157.875 158.125 158.375 158.625 158.875 159.125\n",
"[478] 159.375 159.625 159.875 160.125 160.375 160.625 160.875 161.125 161.375\n",
"[487] 161.625 161.875 162.125 162.375 162.625 162.875 163.125 163.375 163.625\n",
"[496] 163.875 164.125 164.375 164.625 164.875 165.125 165.375 165.625 165.875\n",
"[505] 166.125 166.375 166.625 166.875 167.125 167.375 167.625 167.875 168.125\n",
"[514] 168.375 168.625 168.875 169.125 169.375 169.625 169.875 170.125 170.375\n",
"[523] 170.625 170.875 171.125 171.375 171.625 171.875 172.125 172.375 172.625\n",
"[532] 172.875 173.125 173.375 173.625 173.875 174.125 174.375 174.625 174.875\n",
"[541] 175.125 175.375 175.625 175.875 176.125 176.375 176.625 176.875 177.125\n",
"[550] 177.375 177.625 177.875 178.125 178.375 178.625 178.875 179.125 179.375\n",
"[559] 179.625 179.875"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"lon"
]
},
{
"cell_type": "code",
"execution_count": 39,
"id": "2e650502",
"metadata": {},
"outputs": [],
"source": [
"lat <- ncvar_get(nc=ncin, \"nlat\")"
]
},
{
"cell_type": "code",
"execution_count": 40,
"id": "121fa873",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"- -25.125
- -24.875
- -24.625
- -24.375
- -24.125
- -23.875
- -23.625
- -23.375
- -23.125
- -22.875
- -22.625
- -22.375
- -22.125
- -21.875
- -21.625
- -21.375
- -21.125
- -20.875
- -20.625
- -20.375
- -20.125
- -19.875
- -19.625
- -19.375
- -19.125
- -18.875
- -18.625
- -18.375
- -18.125
- -17.875
- -17.625
- -17.375
- -17.125
- -16.875
- -16.625
- -16.375
- -16.125
- -15.875
- -15.625
- -15.375
- -15.125
- -14.875
- -14.625
- -14.375
- -14.125
- -13.875
- -13.625
- -13.375
- -13.125
- -12.875
- -12.625
- -12.375
- -12.125
- -11.875
- -11.625
- -11.375
- -11.125
- -10.875
- -10.625
- -10.375
- -10.125
- -9.875
- -9.625
- -9.375
- -9.125
- -8.875
- -8.625
- -8.375
- -8.125
- -7.875
- -7.625
- -7.375
- -7.125
- -6.875
- -6.625
- -6.375
- -6.125
- -5.875
- -5.625
- -5.375
- -5.125
- -4.875
- -4.625
- -4.375
- -4.125
- -3.875
- -3.625
- -3.375
- -3.125
- -2.875
- -2.625
- -2.375
- -2.125
- -1.875
- -1.625
- -1.375
- -1.125
- -0.875
- -0.625
- -0.375
- -0.125
- 0.125
- 0.375
- 0.625
- 0.875
- 1.125
- 1.375
- 1.625
- 1.875
- 2.125
- 2.375
- 2.625
- 2.875
- 3.125
- 3.375
- 3.625
- 3.875
- 4.125
- 4.375
- 4.625
- 4.875
- 5.125
- 5.375
- 5.625
- 5.875
- 6.125
- 6.375
- 6.625
- 6.875
- 7.125
- 7.375
- 7.625
- 7.875
- 8.125
- 8.375
- 8.625
- 8.875
- 9.125
- 9.375
- 9.625
- 9.875
- 10.125
- 10.375
- 10.625
- 10.875
- 11.125
- 11.375
- 11.625
- 11.875
- 12.125
- 12.375
- 12.625
- 12.875
- 13.125
- 13.375
- 13.625
- 13.875
- 14.125
- 14.375
- 14.625
- 14.875
- 15.125
- 15.375
- 15.625
- 15.875
- 16.125
- 16.375
- 16.625
- 16.875
- 17.125
- 17.375
- 17.625
- 17.875
- 18.125
- 18.375
- 18.625
- 18.875
- 19.125
- 19.375
- 19.625
- 19.875
- 20.125
- 20.375
- 20.625
- 20.875
- 21.125
- 21.375
- 21.625
- 21.875
- 22.125
- 22.375
- 22.625
- 22.875
- 23.125
- 23.375
- 23.625
- 23.875
- 24.125
- 24.375
- 24.625
- 24.875
\n"
],
"text/latex": [
"\\begin{enumerate*}\n",
"\\item -25.125\n",
"\\item -24.875\n",
"\\item -24.625\n",
"\\item -24.375\n",
"\\item -24.125\n",
"\\item -23.875\n",
"\\item -23.625\n",
"\\item -23.375\n",
"\\item -23.125\n",
"\\item -22.875\n",
"\\item -22.625\n",
"\\item -22.375\n",
"\\item -22.125\n",
"\\item -21.875\n",
"\\item -21.625\n",
"\\item -21.375\n",
"\\item -21.125\n",
"\\item -20.875\n",
"\\item -20.625\n",
"\\item -20.375\n",
"\\item -20.125\n",
"\\item -19.875\n",
"\\item -19.625\n",
"\\item -19.375\n",
"\\item -19.125\n",
"\\item -18.875\n",
"\\item -18.625\n",
"\\item -18.375\n",
"\\item -18.125\n",
"\\item -17.875\n",
"\\item -17.625\n",
"\\item -17.375\n",
"\\item -17.125\n",
"\\item -16.875\n",
"\\item -16.625\n",
"\\item -16.375\n",
"\\item -16.125\n",
"\\item -15.875\n",
"\\item -15.625\n",
"\\item -15.375\n",
"\\item -15.125\n",
"\\item -14.875\n",
"\\item -14.625\n",
"\\item -14.375\n",
"\\item -14.125\n",
"\\item -13.875\n",
"\\item -13.625\n",
"\\item -13.375\n",
"\\item -13.125\n",
"\\item -12.875\n",
"\\item -12.625\n",
"\\item -12.375\n",
"\\item -12.125\n",
"\\item -11.875\n",
"\\item -11.625\n",
"\\item -11.375\n",
"\\item -11.125\n",
"\\item -10.875\n",
"\\item -10.625\n",
"\\item -10.375\n",
"\\item -10.125\n",
"\\item -9.875\n",
"\\item -9.625\n",
"\\item -9.375\n",
"\\item -9.125\n",
"\\item -8.875\n",
"\\item -8.625\n",
"\\item -8.375\n",
"\\item -8.125\n",
"\\item -7.875\n",
"\\item -7.625\n",
"\\item -7.375\n",
"\\item -7.125\n",
"\\item -6.875\n",
"\\item -6.625\n",
"\\item -6.375\n",
"\\item -6.125\n",
"\\item -5.875\n",
"\\item -5.625\n",
"\\item -5.375\n",
"\\item -5.125\n",
"\\item -4.875\n",
"\\item -4.625\n",
"\\item -4.375\n",
"\\item -4.125\n",
"\\item -3.875\n",
"\\item -3.625\n",
"\\item -3.375\n",
"\\item -3.125\n",
"\\item -2.875\n",
"\\item -2.625\n",
"\\item -2.375\n",
"\\item -2.125\n",
"\\item -1.875\n",
"\\item -1.625\n",
"\\item -1.375\n",
"\\item -1.125\n",
"\\item -0.875\n",
"\\item -0.625\n",
"\\item -0.375\n",
"\\item -0.125\n",
"\\item 0.125\n",
"\\item 0.375\n",
"\\item 0.625\n",
"\\item 0.875\n",
"\\item 1.125\n",
"\\item 1.375\n",
"\\item 1.625\n",
"\\item 1.875\n",
"\\item 2.125\n",
"\\item 2.375\n",
"\\item 2.625\n",
"\\item 2.875\n",
"\\item 3.125\n",
"\\item 3.375\n",
"\\item 3.625\n",
"\\item 3.875\n",
"\\item 4.125\n",
"\\item 4.375\n",
"\\item 4.625\n",
"\\item 4.875\n",
"\\item 5.125\n",
"\\item 5.375\n",
"\\item 5.625\n",
"\\item 5.875\n",
"\\item 6.125\n",
"\\item 6.375\n",
"\\item 6.625\n",
"\\item 6.875\n",
"\\item 7.125\n",
"\\item 7.375\n",
"\\item 7.625\n",
"\\item 7.875\n",
"\\item 8.125\n",
"\\item 8.375\n",
"\\item 8.625\n",
"\\item 8.875\n",
"\\item 9.125\n",
"\\item 9.375\n",
"\\item 9.625\n",
"\\item 9.875\n",
"\\item 10.125\n",
"\\item 10.375\n",
"\\item 10.625\n",
"\\item 10.875\n",
"\\item 11.125\n",
"\\item 11.375\n",
"\\item 11.625\n",
"\\item 11.875\n",
"\\item 12.125\n",
"\\item 12.375\n",
"\\item 12.625\n",
"\\item 12.875\n",
"\\item 13.125\n",
"\\item 13.375\n",
"\\item 13.625\n",
"\\item 13.875\n",
"\\item 14.125\n",
"\\item 14.375\n",
"\\item 14.625\n",
"\\item 14.875\n",
"\\item 15.125\n",
"\\item 15.375\n",
"\\item 15.625\n",
"\\item 15.875\n",
"\\item 16.125\n",
"\\item 16.375\n",
"\\item 16.625\n",
"\\item 16.875\n",
"\\item 17.125\n",
"\\item 17.375\n",
"\\item 17.625\n",
"\\item 17.875\n",
"\\item 18.125\n",
"\\item 18.375\n",
"\\item 18.625\n",
"\\item 18.875\n",
"\\item 19.125\n",
"\\item 19.375\n",
"\\item 19.625\n",
"\\item 19.875\n",
"\\item 20.125\n",
"\\item 20.375\n",
"\\item 20.625\n",
"\\item 20.875\n",
"\\item 21.125\n",
"\\item 21.375\n",
"\\item 21.625\n",
"\\item 21.875\n",
"\\item 22.125\n",
"\\item 22.375\n",
"\\item 22.625\n",
"\\item 22.875\n",
"\\item 23.125\n",
"\\item 23.375\n",
"\\item 23.625\n",
"\\item 23.875\n",
"\\item 24.125\n",
"\\item 24.375\n",
"\\item 24.625\n",
"\\item 24.875\n",
"\\end{enumerate*}\n"
],
"text/markdown": [
"1. -25.125\n",
"2. -24.875\n",
"3. -24.625\n",
"4. -24.375\n",
"5. -24.125\n",
"6. -23.875\n",
"7. -23.625\n",
"8. -23.375\n",
"9. -23.125\n",
"10. -22.875\n",
"11. -22.625\n",
"12. -22.375\n",
"13. -22.125\n",
"14. -21.875\n",
"15. -21.625\n",
"16. -21.375\n",
"17. -21.125\n",
"18. -20.875\n",
"19. -20.625\n",
"20. -20.375\n",
"21. -20.125\n",
"22. -19.875\n",
"23. -19.625\n",
"24. -19.375\n",
"25. -19.125\n",
"26. -18.875\n",
"27. -18.625\n",
"28. -18.375\n",
"29. -18.125\n",
"30. -17.875\n",
"31. -17.625\n",
"32. -17.375\n",
"33. -17.125\n",
"34. -16.875\n",
"35. -16.625\n",
"36. -16.375\n",
"37. -16.125\n",
"38. -15.875\n",
"39. -15.625\n",
"40. -15.375\n",
"41. -15.125\n",
"42. -14.875\n",
"43. -14.625\n",
"44. -14.375\n",
"45. -14.125\n",
"46. -13.875\n",
"47. -13.625\n",
"48. -13.375\n",
"49. -13.125\n",
"50. -12.875\n",
"51. -12.625\n",
"52. -12.375\n",
"53. -12.125\n",
"54. -11.875\n",
"55. -11.625\n",
"56. -11.375\n",
"57. -11.125\n",
"58. -10.875\n",
"59. -10.625\n",
"60. -10.375\n",
"61. -10.125\n",
"62. -9.875\n",
"63. -9.625\n",
"64. -9.375\n",
"65. -9.125\n",
"66. -8.875\n",
"67. -8.625\n",
"68. -8.375\n",
"69. -8.125\n",
"70. -7.875\n",
"71. -7.625\n",
"72. -7.375\n",
"73. -7.125\n",
"74. -6.875\n",
"75. -6.625\n",
"76. -6.375\n",
"77. -6.125\n",
"78. -5.875\n",
"79. -5.625\n",
"80. -5.375\n",
"81. -5.125\n",
"82. -4.875\n",
"83. -4.625\n",
"84. -4.375\n",
"85. -4.125\n",
"86. -3.875\n",
"87. -3.625\n",
"88. -3.375\n",
"89. -3.125\n",
"90. -2.875\n",
"91. -2.625\n",
"92. -2.375\n",
"93. -2.125\n",
"94. -1.875\n",
"95. -1.625\n",
"96. -1.375\n",
"97. -1.125\n",
"98. -0.875\n",
"99. -0.625\n",
"100. -0.375\n",
"101. -0.125\n",
"102. 0.125\n",
"103. 0.375\n",
"104. 0.625\n",
"105. 0.875\n",
"106. 1.125\n",
"107. 1.375\n",
"108. 1.625\n",
"109. 1.875\n",
"110. 2.125\n",
"111. 2.375\n",
"112. 2.625\n",
"113. 2.875\n",
"114. 3.125\n",
"115. 3.375\n",
"116. 3.625\n",
"117. 3.875\n",
"118. 4.125\n",
"119. 4.375\n",
"120. 4.625\n",
"121. 4.875\n",
"122. 5.125\n",
"123. 5.375\n",
"124. 5.625\n",
"125. 5.875\n",
"126. 6.125\n",
"127. 6.375\n",
"128. 6.625\n",
"129. 6.875\n",
"130. 7.125\n",
"131. 7.375\n",
"132. 7.625\n",
"133. 7.875\n",
"134. 8.125\n",
"135. 8.375\n",
"136. 8.625\n",
"137. 8.875\n",
"138. 9.125\n",
"139. 9.375\n",
"140. 9.625\n",
"141. 9.875\n",
"142. 10.125\n",
"143. 10.375\n",
"144. 10.625\n",
"145. 10.875\n",
"146. 11.125\n",
"147. 11.375\n",
"148. 11.625\n",
"149. 11.875\n",
"150. 12.125\n",
"151. 12.375\n",
"152. 12.625\n",
"153. 12.875\n",
"154. 13.125\n",
"155. 13.375\n",
"156. 13.625\n",
"157. 13.875\n",
"158. 14.125\n",
"159. 14.375\n",
"160. 14.625\n",
"161. 14.875\n",
"162. 15.125\n",
"163. 15.375\n",
"164. 15.625\n",
"165. 15.875\n",
"166. 16.125\n",
"167. 16.375\n",
"168. 16.625\n",
"169. 16.875\n",
"170. 17.125\n",
"171. 17.375\n",
"172. 17.625\n",
"173. 17.875\n",
"174. 18.125\n",
"175. 18.375\n",
"176. 18.625\n",
"177. 18.875\n",
"178. 19.125\n",
"179. 19.375\n",
"180. 19.625\n",
"181. 19.875\n",
"182. 20.125\n",
"183. 20.375\n",
"184. 20.625\n",
"185. 20.875\n",
"186. 21.125\n",
"187. 21.375\n",
"188. 21.625\n",
"189. 21.875\n",
"190. 22.125\n",
"191. 22.375\n",
"192. 22.625\n",
"193. 22.875\n",
"194. 23.125\n",
"195. 23.375\n",
"196. 23.625\n",
"197. 23.875\n",
"198. 24.125\n",
"199. 24.375\n",
"200. 24.625\n",
"201. 24.875\n",
"\n",
"\n"
],
"text/plain": [
" [1] -25.125 -24.875 -24.625 -24.375 -24.125 -23.875 -23.625 -23.375 -23.125\n",
" [10] -22.875 -22.625 -22.375 -22.125 -21.875 -21.625 -21.375 -21.125 -20.875\n",
" [19] -20.625 -20.375 -20.125 -19.875 -19.625 -19.375 -19.125 -18.875 -18.625\n",
" [28] -18.375 -18.125 -17.875 -17.625 -17.375 -17.125 -16.875 -16.625 -16.375\n",
" [37] -16.125 -15.875 -15.625 -15.375 -15.125 -14.875 -14.625 -14.375 -14.125\n",
" [46] -13.875 -13.625 -13.375 -13.125 -12.875 -12.625 -12.375 -12.125 -11.875\n",
" [55] -11.625 -11.375 -11.125 -10.875 -10.625 -10.375 -10.125 -9.875 -9.625\n",
" [64] -9.375 -9.125 -8.875 -8.625 -8.375 -8.125 -7.875 -7.625 -7.375\n",
" [73] -7.125 -6.875 -6.625 -6.375 -6.125 -5.875 -5.625 -5.375 -5.125\n",
" [82] -4.875 -4.625 -4.375 -4.125 -3.875 -3.625 -3.375 -3.125 -2.875\n",
" [91] -2.625 -2.375 -2.125 -1.875 -1.625 -1.375 -1.125 -0.875 -0.625\n",
"[100] -0.375 -0.125 0.125 0.375 0.625 0.875 1.125 1.375 1.625\n",
"[109] 1.875 2.125 2.375 2.625 2.875 3.125 3.375 3.625 3.875\n",
"[118] 4.125 4.375 4.625 4.875 5.125 5.375 5.625 5.875 6.125\n",
"[127] 6.375 6.625 6.875 7.125 7.375 7.625 7.875 8.125 8.375\n",
"[136] 8.625 8.875 9.125 9.375 9.625 9.875 10.125 10.375 10.625\n",
"[145] 10.875 11.125 11.375 11.625 11.875 12.125 12.375 12.625 12.875\n",
"[154] 13.125 13.375 13.625 13.875 14.125 14.375 14.625 14.875 15.125\n",
"[163] 15.375 15.625 15.875 16.125 16.375 16.625 16.875 17.125 17.375\n",
"[172] 17.625 17.875 18.125 18.375 18.625 18.875 19.125 19.375 19.625\n",
"[181] 19.875 20.125 20.375 20.625 20.875 21.125 21.375 21.625 21.875\n",
"[190] 22.125 22.375 22.625 22.875 23.125 23.375 23.625 23.875 24.125\n",
"[199] 24.375 24.625 24.875"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"lat"
]
},
{
"cell_type": "code",
"execution_count": 41,
"id": "cad73d79",
"metadata": {},
"outputs": [],
"source": [
"Time <- vector(length=length(nc.files))"
]
},
{
"cell_type": "code",
"execution_count": 42,
"id": "fcc6c09d",
"metadata": {},
"outputs": [],
"source": [
"nc_close(ncin)"
]
},
{
"cell_type": "code",
"execution_count": 95,
"id": "cb777d16",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[1] \"1\" \"20121001\"\n",
"[1] \"2\" \"20121001\"\n"
]
}
],
"source": [
"#for (i in 1:length(nc.files)){\n",
"for (i in 1:2){\n",
"#i=1\n",
" ncin <- nc_open(nc.files[i])\n",
" med <- ncvar_get(nc=ncin, \"precipitation\")\n",
" raw_date <- unlist(str_split(nc.files[i], \"[.]\"))[2]\n",
" raw_hour <- unlist(str_split(nc.files[i], \"[.]\"))[3]\n",
"\n",
" strtime <- paste(raw_date, raw_hour, sep=\" \")\n",
" Time[i] <- as.POSIXct(strtime, format=\"%Y%m%d %H\", tz=\"UTC\")\n",
" #Time[i] <- as.Date(strtime, format=\"%Y%m%d %H\", tz=\"UTC\")\n",
" # Time[i] <- as.Date(strptime(raw_time, \"%Y%m%d\"))\n",
" nc_close(ncin)\n",
"\n",
" dimLON <- ncdim_def('lon', \"degree\", longname='longitude', vals=lon)\n",
" dimLAT <- ncdim_def('lat', \"degree\", longname='latitude', vals=lat)\n",
" dimTime <- ncdim_def('time', \"hours since 2012-01-01 00:00:0.0\", longname='time', \n",
" calendar=\"standard\", vals=Time[i], unlim=TRUE)\n",
" fillvalue <- 1e32\n",
" dlname <- \"Precipitation\"\n",
" var_def <- ncvar_def(name=\"precip\", units=\"mm hr-1\", list(dimLAT,dimLON,dimTime),\n",
" fillvalue,dlname,prec=\"double\")\n",
" ncname <- paste0(\"precip_TRMM_\", raw_time,\"_time_.nc\")\n",
" ncpath <- \"\"\n",
" ncfname <- paste0(ncpath, ncname)\n",
" ncoutput <- nc_create(ncfname,list(var_def),force_v4=TRUE, verbose=F)\n",
" ncvar_put(ncoutput, var_def, med)\n",
"\n",
" ncatt_put(ncoutput,\"lon\",\"axis\",\"X\")\n",
" ncatt_put(ncoutput,\"lat\",\"axis\",\"Y\")\n",
" ncatt_put(ncoutput,\"time\",\"axis\",\"T\")\n",
"\n",
" nc_close(ncoutput)\n",
"\n",
" print(c(i,raw_time))\n",
"\n",
" rm(var_def, ncoutput, med)\n",
"}"
]
},
{
"cell_type": "code",
"execution_count": 96,
"id": "17f53057",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"1349049600"
],
"text/latex": [
"1349049600"
],
"text/markdown": [
"1349049600"
],
"text/plain": [
"[1] 1349049600"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Time[1]"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "eff36acc",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "24415ba6",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "R",
"language": "R",
"name": "ir"
},
"language_info": {
"codemirror_mode": "r",
"file_extension": ".r",
"mimetype": "text/x-r-source",
"name": "R",
"pygments_lexer": "r",
"version": "4.2.2"
},
"toc": {
"base_numbering": 1,
"nav_menu": {},
"number_sections": true,
"sideBar": true,
"skip_h1_title": false,
"title_cell": "Table of Contents",
"title_sidebar": "Contents",
"toc_cell": false,
"toc_position": {},
"toc_section_display": true,
"toc_window_display": false
},
"varInspector": {
"cols": {
"lenName": 16,
"lenType": 16,
"lenVar": 40
},
"kernels_config": {
"python": {
"delete_cmd_postfix": "",
"delete_cmd_prefix": "del ",
"library": "var_list.py",
"varRefreshCmd": "print(var_dic_list())"
},
"r": {
"delete_cmd_postfix": ") ",
"delete_cmd_prefix": "rm(",
"library": "var_list.r",
"varRefreshCmd": "cat(var_dic_list()) "
}
},
"types_to_exclude": [
"module",
"function",
"builtin_function_or_method",
"instance",
"_Feature"
],
"window_display": false
}
},
"nbformat": 4,
"nbformat_minor": 5
}