{ "cells": [ { "cell_type": "markdown", "id": "035efd04-ead5-45e0-8b2c-1f0d58290ffd", "metadata": {}, "source": [ "# GOES Total Precipitable Water\n", "`ABI-L2-TPWF` is not an RGB prodcut, it is a Level 2 *derived* product." ] }, { "cell_type": "code", "execution_count": 1, "id": "26100f36-2d6c-40fd-ae6f-5a8174415d93", "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "from toolbox.cartopy_tools import (\n", " common_features,\n", ") # <-- This is from https://github.com/blaylockbk/Carpenter_Workshop\n", "import matplotlib.pyplot as plt\n", "from goes2go.data import goes_nearesttime, goes_latest\n", "from goes2go.tools import abi_crs" ] }, { "cell_type": "markdown", "id": "a341b10b-b735-4c3f-97d8-2c193478d986", "metadata": {}, "source": [ "# TPW CONUS" ] }, { "cell_type": "code", "execution_count": 2, "id": "bc701b1c-c9b5-4e2f-9a04-ae255e16a69e", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " _______________________________\n", " | Satellite: noaa-goes16 |\n", " | Product: ABI-L2-TPWC |\n", " | Domain: C |\n", "📦 Finished downloading [1] files to [/p/cwfs/blaylock/data/noaa-goes16/ABI-L2-TPWC].\n", "📚 Finished reading [1] files into xarray.Dataset.\n" ] }, { "data": { "text/html": [ "
<xarray.Dataset>\n",
"Dimensions: (y: 300, x: 500, number_of_time_bounds: 2, number_of_image_bounds: 2, number_of_LZA_bounds: 2, number_of_SZA_bounds: 2, number_of_lat_bounds: 2, sounding_emissive_bands: 7)\n",
"Coordinates: (12/15)\n",
" t datetime64[ns] 2020-12-...\n",
" * y (y) float32 0.1281 ... ...\n",
" * x (x) float32 -0.1012 ......\n",
" y_image float32 0.08624\n",
" x_image float32 -0.03136\n",
" retrieval_local_zenith_angle float32 80.0\n",
" ... ...\n",
" sounding_emissive_wavelengths (sounding_emissive_bands) float32 ...\n",
" sounding_emissive_band_ids (sounding_emissive_bands) int8 ...\n",
" dataset_name <U72 'OR_ABI-L2-TPWC-M6...\n",
" date_created <U22 '2020-12-14T00:05:...\n",
" time_coverage_start <U22 '2020-12-14T00:01:...\n",
" time_coverage_end <U22 '2020-12-14T00:03:...\n",
"Dimensions without coordinates: number_of_time_bounds, number_of_image_bounds, number_of_LZA_bounds, number_of_SZA_bounds, number_of_lat_bounds, sounding_emissive_bands\n",
"Data variables: (12/30)\n",
" TPW (y, x) float32 ...\n",
" DQF_Overall (y, x) float32 ...\n",
" DQF_Retrieval (y, x) float32 ...\n",
" DQF_SkinTemp (y, x) float32 ...\n",
" time_bounds (number_of_time_bounds) datetime64[ns] ...\n",
" goes_imager_projection int32 -2147483647\n",
" ... ...\n",
" percent_uncorrectable_L0_errors float64 0.0\n",
" percent_uncorrectable_GRB_errors float64 0.0\n",
" total_attempted_retrievals float64 8.576e+04\n",
" mean_obs_modeled_diff_sounding_emissive_bands (sounding_emissive_bands) float32 ...\n",
" std_dev_obs_modeled_diff_sounding_emissive_bands (sounding_emissive_bands) float32 ...\n",
" filename <U108 'noaa-goes16/ABI-...\n",
"Attributes: (12/26)\n",
" naming_authority: gov.nesdis.noaa\n",
" Conventions: CF-1.7\n",
" Metadata_Conventions: Unidata Dataset Discovery v1.0\n",
" standard_name_vocabulary: CF Standard Name Table (v35, 20 July 2016)\n",
" institution: DOC/NOAA/NESDIS > U.S. Department of Commerce,...\n",
" project: GOES\n",
" ... ...\n",
" processing_level: National Aeronautics and Space Administration ...\n",
" cdm_data_type: Image\n",
" timeline_id: ABI Mode 6\n",
" production_data_source: Realtime\n",
" id: 07a697ed-a150-4c8f-b3a9-a00af27a7b20\n",
" path: ['noaa-goes16/ABI-L2-TPWC/2020/349/00/OR_ABI-L...array('2020-12-14T00:02:26.722499968', dtype='datetime64[ns]')array([0.1281 , 0.12782, 0.12754, ..., 0.04494, 0.04466, 0.04438],\n",
" dtype=float32)array([-0.10122, -0.10094, -0.10066, ..., 0.03794, 0.03822, 0.0385 ],\n",
" dtype=float32)array(0.08624, dtype=float32)
array(-0.03136, dtype=float32)
array(80., dtype=float32)
array(70., dtype=float32)
array(180., dtype=float32)
array(70., dtype=float32)
array([ 6.17, 6.93, 7.34, 10.33, 11.19, 12.27, 13.27], dtype=float32)
array([ 8, 9, 10, 13, 14, 15, 16], dtype=int8)
array('OR_ABI-L2-TPWC-M6_G16_s20203490001080_e20203490003453_c20203490005368.nc',\n",
" dtype='<U72')array('2020-12-14T00:05:36.8Z', dtype='<U22')array('2020-12-14T00:01:08.0Z', dtype='<U22')array('2020-12-14T00:03:45.3Z', dtype='<U22')[150000 values with dtype=float32]
[150000 values with dtype=float32]
[150000 values with dtype=float32]
[150000 values with dtype=float32]
array(['2020-12-14T00:01:08.080396032', '2020-12-14T00:03:45.364604032'],\n",
" dtype='datetime64[ns]')array(-2147483647, dtype=int32)
array([0.12824, 0.04424], dtype=float32)
array([-0.10136, 0.03864], dtype=float32)
array(0.)
array(-75.199997)
array(35786.023438)
array(9.96921e+36, dtype=float32)
array(0.)
array(0.722243)
array(48.375961)
array(23.279434)
array(12.999699)
array(-2147483647, dtype=int32)
array(-2147483647, dtype=int32)
array(-2147483647, dtype=int32)
array([ 0., 80.], dtype=float32)
array([ 0., 70.], dtype=float32)
array([ 0., 180.], dtype=float32)
array([-70., 70.], dtype=float32)
array(0.)
array(0.)
array(85759.)
array([-0.337454, 0.485801, -0.59802 , -0.02091 , -0.048114, -0.43329 ,\n",
" -1.010593], dtype=float32)array([1.025441, 0.959408, 0.794338, 0.719077, 0.788357, 0.851414, 0.647595],\n",
" dtype=float32)array('noaa-goes16/ABI-L2-TPWC/2020/349/00/OR_ABI-L2-TPWC-M6_G16_s20203490001080_e20203490003453_c20203490005368.nc',\n",
" dtype='<U108')<xarray.Dataset>\n",
"Dimensions: (y: 1086, x: 1086, number_of_time_bounds: 2, number_of_image_bounds: 2, number_of_LZA_bounds: 2, number_of_SZA_bounds: 2, number_of_lat_bounds: 2, sounding_emissive_bands: 7)\n",
"Coordinates: (12/15)\n",
" t datetime64[ns] 2020-12-...\n",
" * y (y) float32 0.1519 ... ...\n",
" * x (x) float32 -0.1519 ......\n",
" y_image float32 0.0\n",
" x_image float32 0.0\n",
" retrieval_local_zenith_angle float32 80.0\n",
" ... ...\n",
" sounding_emissive_wavelengths (sounding_emissive_bands) float32 ...\n",
" sounding_emissive_band_ids (sounding_emissive_bands) int8 ...\n",
" dataset_name <U72 'OR_ABI-L2-TPWF-M6...\n",
" date_created <U22 '2020-12-14T00:10:...\n",
" time_coverage_start <U22 '2020-12-14T00:00:...\n",
" time_coverage_end <U22 '2020-12-14T00:09:...\n",
"Dimensions without coordinates: number_of_time_bounds, number_of_image_bounds, number_of_LZA_bounds, number_of_SZA_bounds, number_of_lat_bounds, sounding_emissive_bands\n",
"Data variables: (12/30)\n",
" TPW (y, x) float32 ...\n",
" DQF_Overall (y, x) float32 ...\n",
" DQF_Retrieval (y, x) float32 ...\n",
" DQF_SkinTemp (y, x) float32 ...\n",
" time_bounds (number_of_time_bounds) datetime64[ns] ...\n",
" goes_imager_projection int32 -2147483647\n",
" ... ...\n",
" percent_uncorrectable_L0_errors float64 0.0\n",
" percent_uncorrectable_GRB_errors float64 0.0\n",
" total_attempted_retrievals float64 4.077e+05\n",
" mean_obs_modeled_diff_sounding_emissive_bands (sounding_emissive_bands) float32 ...\n",
" std_dev_obs_modeled_diff_sounding_emissive_bands (sounding_emissive_bands) float32 ...\n",
" filename <U108 'noaa-goes16/ABI-...\n",
"Attributes: (12/26)\n",
" naming_authority: gov.nesdis.noaa\n",
" Conventions: CF-1.7\n",
" Metadata_Conventions: Unidata Dataset Discovery v1.0\n",
" standard_name_vocabulary: CF Standard Name Table (v35, 20 July 2016)\n",
" institution: DOC/NOAA/NESDIS > U.S. Department of Commerce,...\n",
" project: GOES\n",
" ... ...\n",
" processing_level: National Aeronautics and Space Administration ...\n",
" cdm_data_type: Image\n",
" timeline_id: ABI Mode 6\n",
" production_data_source: Realtime\n",
" id: e4eab11a-4da8-454d-903c-07472bed0229\n",
" path: ['noaa-goes16/ABI-L2-TPWF/2020/349/00/OR_ABI-L...array('2020-12-14T00:04:56.709048064', dtype='datetime64[ns]')array([ 0.1519 , 0.15162, 0.15134, ..., -0.15134, -0.15162, -0.1519 ],\n",
" dtype=float32)array([-0.1519 , -0.15162, -0.15134, ..., 0.15134, 0.15162, 0.1519 ],\n",
" dtype=float32)array(0., dtype=float32)
array(0., dtype=float32)
array(80., dtype=float32)
array(70., dtype=float32)
array(180., dtype=float32)
array(70., dtype=float32)
array([ 6.17, 6.93, 7.34, 10.33, 11.19, 12.27, 13.27], dtype=float32)
array([ 8, 9, 10, 13, 14, 15, 16], dtype=int8)
array('OR_ABI-L2-TPWF-M6_G16_s20203490000113_e20203490009421_c20203490010545.nc',\n",
" dtype='<U72')array('2020-12-14T00:10:54.5Z', dtype='<U22')array('2020-12-14T00:00:11.3Z', dtype='<U22')array('2020-12-14T00:09:42.1Z', dtype='<U22')[1179396 values with dtype=float32]
[1179396 values with dtype=float32]
[1179396 values with dtype=float32]
[1179396 values with dtype=float32]
array(['2020-12-14T00:00:11.309766016', '2020-12-14T00:09:42.108329088'],\n",
" dtype='datetime64[ns]')array(-2147483647, dtype=int32)
array([ 0.151872, -0.151872], dtype=float32)
array([-0.151872, 0.151872], dtype=float32)
array(0.)
array(-75.199997)
array(35786.023438)
array(9.96921e+36, dtype=float32)
array(0.)
array(0.669217)
array(65.253807)
array(27.692272)
array(13.492225)
array(-2147483647, dtype=int32)
array(-2147483647, dtype=int32)
array(-2147483647, dtype=int32)
array([ 0., 80.], dtype=float32)
array([ 0., 70.], dtype=float32)
array([ 0., 180.], dtype=float32)
array([-70., 70.], dtype=float32)
array(0.)
array(0.)
array(407682.)
array([ 0.052707, 0.743839, -0.450023, -0.318673, -0.310954, -0.676081,\n",
" -1.089596], dtype=float32)array([1.171488, 1.118285, 0.913385, 0.767696, 0.853697, 0.914093, 0.71145 ],\n",
" dtype=float32)array('noaa-goes16/ABI-L2-TPWF/2020/349/00/OR_ABI-L2-TPWF-M6_G16_s20203490000113_e20203490009421_c20203490010545.nc',\n",
" dtype='<U108')