10/11/24, 13:10 FOR_COR_PC1_JATENG_Lismasari.

ipynb - Colab
PM

!pip install xarray

!pip install matplotlib
!pip install eofs
!pip install cartopy
!pip install Basemap
!pip install netcdf4

https://colab.research.google.com/drive/1jxGDhOcV030lRCBtUqWsOryXlxLG9scW#printMode=true 1/5
10/11/24, 13:10 FOR_COR_PC1_JATENG_Lismasari.ipynb - Colab
PM
Requirement already satisfied: xarray in /usr/local/lib/python3.10/dist-packages (2024.10.0)
Requirement already satisfied: numpy>=1.24 in /usr/local/lib/python3.10/dist-packages (from xarray) (1.26.4)
Requirement already satisfied: packaging>=23.1 in /usr/local/lib/python3.10/dist-packages (from xarray) (24.1)
Requirement already satisfied: pandas>=2.1 in /usr/local/lib/python3.10/dist-packages (from xarray) (2.2.2)
Requirement already satisfied: python-dateutil>=2.8.2 in /usr/local/lib/python3.10/dist-packages (from pandas>=2.1->xarray) (2.8.2)
Requirement already satisfied: pytz>=2020.1 in /usr/local/lib/python3.10/dist-packages (from pandas>=2.1->xarray) (2024.2)
Requirement already satisfied: tzdata>=2022.7 in /usr/local/lib/python3.10/dist-packages (from pandas>=2.1->xarray) (2024.2)
Requirement already satisfied: six>=1.5 in /usr/local/lib/python3.10/dist-packages (from python-dateutil>=2.8.2->pandas>=2.1->xarray)
Requirement already satisfied: matplotlib in /usr/local/lib/python3.10/dist-packages (3.8.0)
Requirement already satisfied: contourpy>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (1.3.0)
Requirement already satisfied: cycler>=0.10 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (0.12.1)
Requirement already satisfied: fonttools>=4.22.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (4.54.1)
Requirement already satisfied: kiwisolver>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (1.4.7)
Requirement already satisfied: numpy<2,>=1.21 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (1.26.4)
Requirement already satisfied: packaging>=20.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (24.1)
Requirement already satisfied: pillow>=6.2.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (10.4.0)
Requirement already satisfied: pyparsing>=2.3.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (3.2.0)
Requirement already satisfied: python-dateutil>=2.7 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (2.8.2)
Requirement already satisfied: six>=1.5 in /usr/local/lib/python3.10/dist-packages (from python-dateutil>=2.7->matplotlib) (1.16.0)
Collecting eofs
Downloading eofs-1.4.1-py3-none-any.whl.metadata (3.6 kB)
Requirement already satisfied: numpy in /usr/local/lib/python3.10/dist-packages (from eofs) (1.26.4)
Downloading eofs-1.4.1-py3-none-any.whl (1.1 MB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 1.1/1.1 MB 12.0 MB/s eta 0:00:00
Installing collected packages: eofs
Successfully installed eofs-1.4.1
Collecting cartopy
Downloading Cartopy-0.24.1-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl.metadata (7.9 kB)
Requirement already satisfied: numpy>=1.23 in /usr/local/lib/python3.10/dist-packages (from cartopy) (1.26.4)
Requirement already satisfied: matplotlib>=3.6 in /usr/local/lib/python3.10/dist-packages (from cartopy) (3.8.0)
Requirement already satisfied: shapely>=1.8 in /usr/local/lib/python3.10/dist-packages (from cartopy) (2.0.6)
Requirement already satisfied: packaging>=21 in /usr/local/lib/python3.10/dist-packages (from cartopy) (24.1)
Requirement already satisfied: pyshp>=2.3 in /usr/local/lib/python3.10/dist-packages (from cartopy) (2.3.1)
Requirement already satisfied: pyproj>=3.3.1 in /usr/local/lib/python3.10/dist-packages (from cartopy) (3.7.0)
Requirement already satisfied: contourpy>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (1.3.0)
Requirement already satisfied: cycler>=0.10 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (0.12.1)
Requirement already satisfied: fonttools>=4.22.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (4.54.1)
Requirement already satisfied: kiwisolver>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (1.4.7)
Requirement already satisfied: pillow>=6.2.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (10.4.0)
Requirement already satisfied: pyparsing>=2.3.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (3.2.0)
Requirement already satisfied: python-dateutil>=2.7 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (2.8.2
Requirement already satisfied: certifi in /usr/local/lib/python3.10/dist-packages (from pyproj>=3.3.1->cartopy) (2024.8.30)
Requirement already satisfied: six>=1.5 in /usr/local/lib/python3.10/dist-packages (from python-dateutil>=2.7->matplotlib>=3.6->carto
Downloading Cartopy-0.24.1-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (11.7 MB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 11.7/11.7 MB 77.3 MB/s eta 0:00:00
Installing collected packages: cartopy
Successfully installed cartopy-0.24.1
Collecting Basemap
Downloading basemap-1.4.1-cp310-cp310-manylinux1_x86_64.whl.metadata (9.1 kB)
Collecting basemap-data<1.4,>=1.3.2 (from Basemap)
Downloading basemap_data-1.3.2-py2.py3-none-any.whl.metadata (2.7 kB)
Requirement already satisfied: pyshp<2.4,>=1.2 in /usr/local/lib/python3.10/dist-packages (from Basemap) (2.3.1)
Requirement already satisfied: matplotlib<3.9,>=1.5 in /usr/local/lib/python3.10/dist-packages (from Basemap) (3.8.0)
Collecting pyproj<3.7.0,>=1.9.3 (from Basemap)
Downloading pyproj-3.6.1-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl.metadata (31 kB)
Collecting packaging<24.0,>=16.0 (from Basemap)
Downloading packaging-23.2-py3-none-any.whl.metadata (3.2 kB)
Requirement already satisfied: numpy<1.27,>=1.21 in /usr/local/lib/python3.10/dist-packages (from Basemap) (1.26.4)
Requirement already satisfied: contourpy>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (1.3.
Requirement already satisfied: cycler>=0.10 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (0.12.1)
Requirement already satisfied: fonttools>=4.22.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (4.5
Requirement already satisfied: kiwisolver>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (1.4
Requirement already satisfied: pillow>=6.2.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (10.4.0)
Requirement already satisfied: pyparsing>=2.3.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (3.2.
Requirement already satisfied: python-dateutil>=2.7 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (
Requirement already satisfied: certifi in /usr/local/lib/python3.10/dist-packages (from pyproj<3.7.0,>=1.9.3->Basemap) (2024.8.30)
Requirement already satisfied: six>=1.5 in /usr/local/lib/python3.10/dist-packages (from python-dateutil>=2.7->matplotlib<3.9,>=1.5->
Downloading basemap-1.4.1-cp310-cp310-manylinux1_x86_64.whl (935 kB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 936.0/936.0 kB 12.4 MB/s eta 0:00:00
Downloading basemap_data-1.3.2-py2.py3-none-any.whl (30.5 MB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 30.5/30.5 MB 54.8 MB/s eta 0:00:00
Downloading packaging-23.2-py3-none-any.whl (53 kB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 53.0/53.0 kB 5.7 MB/s eta 0:00:00
Downloading pyproj-3.6.1-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (8.3 MB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 8.3/8.3 MB 101.9 MB/s eta 0:00:00
Installing collected packages: pyproj, packaging, basemap-data, Basemap
Attempting uninstall: pyproj
Found existing installation: pyproj 3.7.0
Uninstalling pyproj-3.7.0:
Successfully uninstalled pyproj-3.7.0
Attempting uninstall: packaging
Found existing installation: packaging 24.1
Uninstalling packaging-24.1:
Successfully uninstalled packaging-24.1

https://colab.research.google.com/drive/1jxGDhOcV030lRCBtUqWsOryXlxLG9scW#printMode=true 2/5
10/11/24, 13:10 FOR_COR_PC1_JATENG_Lismasari.ipynb - Colab
PM Successfully installed Basemap-1.4.1 basemap-data-1.3.2 packaging-23.2 pyproj-3.6.1
Collecting netcdf4
Downloading netCDF4-1.7.2-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl.metadata (1.8 kB)
Collecting cftime (from netcdf4)
Downloading cftime-1.6.4.post1-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl.metadata (8.7 kB)
Requirement already satisfied: certifi in /usr/local/lib/python3.10/dist-packages (from netcdf4) (2024.8.30)
Requirement already satisfied: numpy in /usr/local/lib/python3.10/dist-packages (from netcdf4) (1.26.4)
Downloading netCDF4-1.7.2-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (9.1 MB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 9.1/9.1 MB 42.0 MB/s eta 0:00:00
Downloading cftime-1.6.4.post1-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (1.3 MB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 1.3/1.3 MB 54.3 MB/s eta 0:00:00

import numpy as np
import numpy.ma as ma
import xarray as xr
import matplotlib.pyplot as plt
import matplotlib as mpl
import matplotlib.colors as colors
import matplotlib.patches as mpatches
from mpl_toolkits.basemap import Basemap
from numpy.fft import fft, fftfreq, ifft
import pandas as pd

from scipy.fftpack import *

#import cartopy.crs as ccrs
#import cartopy.feature as cfeature

#from scipy import signal

#from eofs.standard import Eof

plt.rcParams['font.family'] = 'deJavu Sans'

plt.rcParams['font.size'] = 35 #22
plt.rcParams['font.weight'] = 'bold' #22
#plt.rcParams['figure.figsize'] = [35,10]
import matplotlib.ticker as ticker

from google.colab import drive

drive.mount('/content/drive')

Drive already mounted at /content/drive; to attempt to forcibly remount, call drive.mount("/content/drive", force_remount=True).

pc_1=pd.read_excel('/content/drive/MyDrive/Fisika Komputasi/data_curah_hujan_jawa_tengah.xlsx')
pc_1_nc=pc_1[0:480] #48-16
#northern=northern_bmkg[120:480] #90-19
pc_1

Tahun pc1 pc2 pc3

0 01-01-1993 8888.0 0.1 7.7

1 01-02-1993 13.0 3.2 126.6

2 01-03-1993 0.0 0.0 22.0

3 01-04-1993 3.4 6.3 1.9

4 01-05-1993 0.8 0.5 0.0

... ... ... ... ...

366 01-08-2023 0.0 0.0 0.0

367 01-09-2023 NaN 0.0 0.0

368 01-10-2023 0.0 0.0 0.0

369 01-11-2023 0.0 0.0 7.0

370 01-12-2023 2.0 25.4 0.5

371 rows × 4 columns

ds_sst_slp_wind_era5=xr.open_dataset("/content/drive/MyDrive/Fisika Komputasi/FULL_LONG-1940-2022-ver1-5-XX.nc")

t_start='1993-01-01'
t_stop='2023-12-01'

min_lon = 40

https://colab.research.google.com/drive/1jxGDhOcV030lRCBtUqWsOryXlxLG9scW#printMode=true 3/5
10/11/24, 13:10 FOR_COR_PC1_JATENG_Lismasari.ipynb - Colab
PM max_lon = 250
min_lat = -35
max_lat = 35

#ds_uwind_sel = ds_uwind.sel(time=slice(t_start, t_stop), lat=slice(max_lat,min_lat), lon=slice(min_lon,max_lon))

#ds_vwind_sel = ds_vwind.sel(time=slice(t_start, t_stop), lat=slice(max_lat,min_lat), lon=slice(min_lon,max_lon))
#ds_slp_sel = ds_slp.sel(time=slice(t_start, t_stop), lat=slice(max_lat,min_lat), lon=slice(min_lon,max_lon))
ds_era5_sel = ds_sst_slp_wind_era5.sel(time=slice(t_start, t_stop), latitude=slice(max_lat,min_lat), longitude=slice(min_lon,max_lon))
#ds_era5_sel

ds_era5_sel_climo = ds_era5_sel.groupby('time.month').mean('time')
ds_era5_sel_anoms = (ds_era5_sel.groupby('time.month') - ds_era5_sel_climo)

ds_era5_sel

xarray.Dataset

► Dimensions: (time: 360, latitude: 281, longitude: 841)

▼ Coordinates:
longitude (longitude) float32 40.0 40.25 40.5 ... 249.8 250.0
latitude (latitude) float32 35.0 34.75 34.5 ... -34.75 -35.0
time (time) datetime64[ns] 1993-01-01 ... 2022-12-01

▼ Data variables:
u10 (time, latitude, longitude) float32 ...
v10 (time, latitude, longitude) float32 ...
msl (time, latitude, longitude) float32 ...
sst (time, latitude, longitude) float32 ...

► Indexes: (3)

▼ Attributes:
Conventions : CF-1.6
history : 2023-03-10 06:55:34 GMT by grib_to_netcdf-2.25.1: /opt/ecmwf/mars-client/bin/grib_to
_netcdf.bin -S param -o /cache/data6/adaptor.mars.internal-1678431328.9136436-787
5-2-7e6d2ff1-c136-4379-91b1-0200a1d8d906.nc /cache/tmp/7e6d2ff1-c136-4379-91b1
-0200a1d8d906-adaptor.mars.internal-1678431297.9416697-7875-3-tmp.grib

print(pc_1.columns)

Index(['Tahun', 'pc1', 'pc2', 'pc3'], dtype='object')

rr_fcor=[]
list_tp=pc_1['pc1'].to_list()
for k in range(len(list_tp)):
rr_fcor.append([[list_tp[k] for i in range(len(ds_era5_sel.longitude))] for j in range(len(ds_era5_sel.latitude))])

print(len(rr_fcor), len(rr_fcor[0]), len(rr_fcor[0][0]))

371 281 841

#rr_fcor

rr_fcorx=[[list_tp for i in range(len(ds_era5_sel.latitude))] for j in range(len(ds_era5_sel.longitude))]

print(len(rr_fcorx), len(rr_fcorx[0]), len(rr_fcorx[0][0]))

841 281 371

for ii in range(2):
for jj in range(3):
# for kk in range(len(northern.tp)):
rr_fcorx[jj][ii]=(pc_1.pc1.to_list())

print(len(rr_fcorx), len(rr_fcorx[0]), len(rr_fcorx[0][0]))

841 281 371

https://colab.research.google.com/drive/1jxGDhOcV030lRCBtUqWsOryXlxLG9scW#printMode=true 4/5
10/11/24, 13:10 FOR_COR_PC1_JATENG_Lismasari.ipynb - Colab
PM rr_fcorxx=[[np.nan for i in range(2)] for j in range(3)]
for ii in range(2):
for jj in range(3):
# for kk in range(len(northern.tp)):
rr_fcorxx[jj][ii]=(pc_1.pc1.to_list())

print(len(rr_fcorxx), len(rr_fcorxx[0]), len(rr_fcorxx[0][0]))

3 2 371

#print(len(northern.years), len(ds_era5_sel.longitude), len(ds_era5_sel.latitude))

#len(ds_era5_sel.longitude)

for i in range(len(ds_era5_sel.latitude)):
for j in range(len(ds_era5_sel.longitude)):
#for k in range(len(northern.tp)):
rr_fcorxx[jj][ii]=(pc_1.pc1.to_list())

print(len(rr_fcor), len(rr_fcor[0]), len(rr_fcor[0][0]))

371 281 841

sst_rr_xar=xr.DataArray(
data = rr_fcor,
dims = ['time','latitude','longitude'],
coords = {'time':pd.to_datetime(pc_1.Tahun),'latitude': ds_era5_sel.latitude,'longitude': ds_era5_sel.longitude},
attrs = {'rr': 'rr_stasiun data', 'units': 'mm'}
)

sst_rr_xar

https://colab.research.google.com/drive/1jxGDhOcV030lRCBtUqWsOryXlxLG9scW#printMode=true 5/5