{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "c6654f9c", "metadata": {}, "outputs": [], "source": [ "#conda create -n cptec python=3.10\n", "\n", "#conda activate cptec\n", "#conda install -c conda-forge xarray matplotlib pycurl cfgrib netcdf4 pynio xarray matplotlib pycurl cfgrib netcdf4 pynio dask ipywidgets scipy mpi4py esmpy\n", "\n", "#conda install -c anaconda ipykernel\n", "\n", "# Instalar o pacote wesmf\n", "#conda install -c conda-forge xesmf\n", "#pip install jupyter" ] }, { "cell_type": "code", "execution_count": 1, "id": "c1b51b57", "metadata": {}, "outputs": [], "source": [ "import cptecmodel.CPTEC_ETA as ETA\n", "import cptecmodel.CPTEC_WRF as WRF\n", "import cptecmodel.CPTEC_BAM as BAM\n", "# Adicionar o Pacote UTIL\n", "import cptecmodel.CPTEC_UTIL as cpt\n" ] }, { "cell_type": "code", "execution_count": 2, "id": "b3d63701", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "#### Regional (ams_08km) #####\n", "\n", "Forecast data available for reading between 20230110 and 20230120.\n", "\n", "Surface variables: t2m, u10m, v10m, slp, psfc, precip\n", " terrain, sbcape, sbcin, pw.\n", "Level variables: t, u, v, rh, g, omega.\n", "\n", "levels (hPa): 1020 1000 950 925 900 850 800 750 700 650 600 \n", " 550 500 450 400 350 300 250 200 150 100 50.\n", "\n", "Frequency: hourly frequency [0,1,2,...,22,23].\n", "\n", "\n", "#### US National Weather Service - NCEP (WMC) (ams_07km) #####\n", "\n", "Forecast data available for reading between 20230110 and 20230120.\n", "\n", "Surface variables: t2m, u10m, v10m, slp, psfc, precip\n", " terrain, sbcape, sbcin, pw.\n", "Level variables: t, u, v, rh, g, omega.\n", "\n", "levels (hPa): 1000 975 950 925 900 875 850 825 800 775 750 700 650\n", " 600 550 500 450 400 350 300 250 200 150 100 50.\n", "\n", "Frequency: hourly frequency [0,1,2,...,22,23].\n", "\n" ] } ], "source": [ "# Inicializar os Modelos\n", "eta = ETA.model()\n", "wrf = WRF.model()" ] }, { "cell_type": "code", "execution_count": 4, "id": "45b0119d", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Help on function regrid in module CPTEC_UTIL:\n", "\n", "regrid(source, target)\n", " A função regrid muda a grade do arquivo fonte para a grade do arquivo destino\n", " \n", " Parametros\n", " ------------------------------------------------------------------------------------------------------------ \n", " source : Xarray da grade fonte\n", " target : Xarray da grade destino\n", " ------------------------------------------------------------------------------------------------------------ \n", " \n", " regrid(source,target)\n", " \n", " ------------------------------------------------------------------------------------------------------------ \n", " \n", " Retorna um Xarray contendo com os dados do Fonte na grade de Destino\n", " \n", " ------------------------------------------------------------------------------------------------------------\n", "\n" ] } ], "source": [ "help(cpt.regrid)" ] }, { "cell_type": "code", "execution_count": 5, "id": "e2174c04", "metadata": {}, "outputs": [], "source": [ "# Ler Dados do Eta\n", "data_eta = eta.load(date='2022111500', var='t2m', steps=3)" ] }, { "cell_type": "code", "execution_count": 6, "id": "61f385d5", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n",
       "Dimensions:    (time: 4, latitude: 931, longitude: 875)\n",
       "Coordinates:\n",
       "  * time       (time) datetime64[ns] 2022-11-15 ... 2022-11-15T03:00:00\n",
       "  * latitude   (latitude) float64 -55.0 -54.92 -54.84 ... 19.24 19.32 19.4\n",
       "  * longitude  (longitude) float64 270.0 270.1 270.2 270.2 ... 339.8 339.8 339.9\n",
       "Data variables:\n",
       "    t2m        (time, latitude, longitude) float32 5.288 5.288 ... 24.99 25.01\n",
       "Attributes:\n",
       "    center:   National Institute for Space Research - INPE\n",
       "    model:    Eta Regional (ams_08km)
" ], "text/plain": [ "\n", "Dimensions: (time: 4, latitude: 931, longitude: 875)\n", "Coordinates:\n", " * time (time) datetime64[ns] 2022-11-15 ... 2022-11-15T03:00:00\n", " * latitude (latitude) float64 -55.0 -54.92 -54.84 ... 19.24 19.32 19.4\n", " * longitude (longitude) float64 270.0 270.1 270.2 270.2 ... 339.8 339.8 339.9\n", "Data variables:\n", " t2m (time, latitude, longitude) float32 5.288 5.288 ... 24.99 25.01\n", "Attributes:\n", " center: National Institute for Space Research - INPE\n", " model: Eta Regional (ams_08km)" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "data_eta" ] }, { "cell_type": "code", "execution_count": 7, "id": "c80d5128", "metadata": {}, "outputs": [], "source": [ "# Ler dados do WRF\n", "data_wrf = wrf.load(date='2022111500', var='t2m', steps=3)" ] }, { "cell_type": "code", "execution_count": 8, "id": "6671029e", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n",
       "Dimensions:    (time: 4, latitude: 1081, longitude: 1019)\n",
       "Coordinates:\n",
       "  * time       (time) datetime64[ns] 2022-11-15 ... 2022-11-15T03:00:00\n",
       "  * latitude   (latitude) float64 -57.9 -57.83 -57.76 ... 17.56 17.63 17.7\n",
       "  * longitude  (longitude) float64 269.3 269.4 269.5 269.5 ... 340.4 340.5 340.6\n",
       "Data variables:\n",
       "    t2m        (time, latitude, longitude) float32 4.065 4.085 ... 25.82 25.86\n",
       "Attributes:\n",
       "    center:   National Institute for Space Research - INPE\n",
       "    model:    The Weather Research and Forecasting (ams_07km)
" ], "text/plain": [ "\n", "Dimensions: (time: 4, latitude: 1081, longitude: 1019)\n", "Coordinates:\n", " * time (time) datetime64[ns] 2022-11-15 ... 2022-11-15T03:00:00\n", " * latitude (latitude) float64 -57.9 -57.83 -57.76 ... 17.56 17.63 17.7\n", " * longitude (longitude) float64 269.3 269.4 269.5 269.5 ... 340.4 340.5 340.6\n", "Data variables:\n", " t2m (time, latitude, longitude) float32 4.065 4.085 ... 25.82 25.86\n", "Attributes:\n", " center: National Institute for Space Research - INPE\n", " model: The Weather Research and Forecasting (ams_07km)" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "data_wrf" ] }, { "cell_type": "code", "execution_count": 9, "id": "15fc36ad", "metadata": {}, "outputs": [], "source": [ "# Mudar a grade do Eta no padrao WRF \n", "out = cpt.regrid(data_eta, data_wrf)" ] }, { "cell_type": "code", "execution_count": 10, "id": "34e19064", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n",
       "Dimensions:    (time: 4, latitude: 1081, longitude: 1019)\n",
       "Coordinates:\n",
       "  * time       (time) datetime64[ns] 2022-11-15 ... 2022-11-15T03:00:00\n",
       "  * latitude   (latitude) float64 -57.9 -57.83 -57.76 ... 17.56 17.63 17.7\n",
       "  * longitude  (longitude) float64 269.3 269.4 269.5 269.5 ... 340.4 340.5 340.6\n",
       "Data variables:\n",
       "    t2m        (time, latitude, longitude) float32 0.0 0.0 0.0 ... 0.0 0.0 0.0\n",
       "Attributes:\n",
       "    regrid_method:  bilinear
" ], "text/plain": [ "\n", "Dimensions: (time: 4, latitude: 1081, longitude: 1019)\n", "Coordinates:\n", " * time (time) datetime64[ns] 2022-11-15 ... 2022-11-15T03:00:00\n", " * latitude (latitude) float64 -57.9 -57.83 -57.76 ... 17.56 17.63 17.7\n", " * longitude (longitude) float64 269.3 269.4 269.5 269.5 ... 340.4 340.5 340.6\n", "Data variables:\n", " t2m (time, latitude, longitude) float32 0.0 0.0 0.0 ... 0.0 0.0 0.0\n", "Attributes:\n", " regrid_method: bilinear" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "out" ] }, { "cell_type": "code", "execution_count": 11, "id": "f53db79a", "metadata": {}, "outputs": [], "source": [ "# Mudar a grade do WRF no padrao do Eta\n", "out2 = cpt.regrid(data_wrf, data_eta)" ] }, { "cell_type": "code", "execution_count": 12, "id": "d2c12606", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n",
       "Dimensions:    (time: 4, latitude: 931, longitude: 875)\n",
       "Coordinates:\n",
       "  * time       (time) datetime64[ns] 2022-11-15 ... 2022-11-15T03:00:00\n",
       "  * latitude   (latitude) float64 -55.0 -54.92 -54.84 ... 19.24 19.32 19.4\n",
       "  * longitude  (longitude) float64 270.0 270.1 270.2 270.2 ... 339.8 339.8 339.9\n",
       "Data variables:\n",
       "    t2m        (time, latitude, longitude) float32 5.273 5.25 5.216 ... 0.0 0.0\n",
       "Attributes:\n",
       "    regrid_method:  bilinear
" ], "text/plain": [ "\n", "Dimensions: (time: 4, latitude: 931, longitude: 875)\n", "Coordinates:\n", " * time (time) datetime64[ns] 2022-11-15 ... 2022-11-15T03:00:00\n", " * latitude (latitude) float64 -55.0 -54.92 -54.84 ... 19.24 19.32 19.4\n", " * longitude (longitude) float64 270.0 270.1 270.2 270.2 ... 339.8 339.8 339.9\n", "Data variables:\n", " t2m (time, latitude, longitude) float32 5.273 5.25 5.216 ... 0.0 0.0\n", "Attributes:\n", " regrid_method: bilinear" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "out2" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.8" } }, "nbformat": 4, "nbformat_minor": 5 }