{"metadata":{"kernelspec":{"language":"python","display_name":"Python 3","name":"python3"},"language_info":{"name":"python","version":"3.10.13","mimetype":"text/x-python","codemirror_mode":{"name":"ipython","version":3},"pygments_lexer":"ipython3","nbconvert_exporter":"python","file_extension":".py"},"kaggle":{"accelerator":"none","dataSources":[{"sourceId":56537,"databundleVersionId":8015876,"sourceType":"competition"},{"sourceId":8681881,"sourceType":"datasetVersion","datasetId":5204766},{"sourceId":8682000,"sourceType":"datasetVersion","datasetId":5204846}],"dockerImageVersionId":30732,"isInternetEnabled":true,"language":"python","sourceType":"notebook","isGpuEnabled":false}},"nbformat_minor":4,"nbformat":4,"cells":[{"cell_type":"code","source":"!pip install netCDF4\nimport netCDF4 as nc\n\nds = nc.Dataset(\"/kaggle/input/climsim-info/ClimSim_low-res_grid-info.nc\")\n\nds","metadata":{"execution":{"iopub.status.busy":"2024-06-13T11:21:21.711518Z","iopub.execute_input":"2024-06-13T11:21:21.712763Z","iopub.status.idle":"2024-06-13T11:21:36.701072Z","shell.execute_reply.started":"2024-06-13T11:21:21.712717Z","shell.execute_reply":"2024-06-13T11:21:36.699836Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"import polars as pl\n\ndf = pl.read_csv(\"/kaggle/input/leap-atmospheric-physics-ai-climsim/train.csv\", \n                    columns=[\"state_t_59\"], n_rows=384)\ndf = df.to_pandas()\ndf[\"lon\"] = ds[\"lon\"][:]\ndf[\"lat\"] = ds[\"lat\"][:]\ndf","metadata":{"execution":{"iopub.status.busy":"2024-06-13T11:21:41.533224Z","iopub.execute_input":"2024-06-13T11:21:41.533633Z","iopub.status.idle":"2024-06-13T11:21:42.753618Z","shell.execute_reply.started":"2024-06-13T11:21:41.533599Z","shell.execute_reply":"2024-06-13T11:21:42.75233Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"import matplotlib.pyplot as plt\nimport cartopy.crs as ccrs\nimport cartopy.feature as cfeature\n\nfig, ax = plt.subplots(nrows=1,ncols=1,\n                        subplot_kw={'projection': ccrs.Robinson(central_longitude=179.5)},\n                        figsize=(8.,3.))\nax.set_global()\nx = df['lon']\ny = df['lat']\nz = df[\"state_t_59\"]\n# cmap = \"turbo\"\nh = ax.tricontourf(x, y, z, transform=ccrs.PlateCarree(),\n#                    levels = clevels, \n#                    cmap=cmap, \n                   extend='both')\nh.cmap.set_under('silver')\nax.coastlines()\nfig.subplots_adjust(right=0.9)\ncbar_ax = fig.add_axes([0.925, 0.145, 0.020, 0.71])\nfig.colorbar(h, cax=cbar_ax)\n# fig.set_facecolor('w')","metadata":{"execution":{"iopub.status.busy":"2024-06-13T11:21:47.628987Z","iopub.execute_input":"2024-06-13T11:21:47.629505Z","iopub.status.idle":"2024-06-13T11:21:49.818424Z","shell.execute_reply.started":"2024-06-13T11:21:47.629474Z","shell.execute_reply":"2024-06-13T11:21:49.816848Z"},"trusted":true},"execution_count":null,"outputs":[]}]}