## Analisis Hidrodinamika dan Transformasi
Gelombang bagi Re-design Breakwater di Pelabuhan Klungkung Bali

## ALIF FARHAN DARMAWAN, YESSI NIRWANA KURNIADI
Jurusan Teknik Sipil Institut Teknologi Nasional Email: aliffarhand27@gmail.com

## ABSTRAK
Pelabuhan  Gunaksa  Klungkung,  Bali  memerlukan  breakwateruntuk  melidungi pelabuhan  dari  gelombang  tinggi.  Analisis  saat  ini  menunjukan  bahwa breakwater  yang  ada  tidak  dapat  mereduksi  tinggi  gelombang  yang  terjadi. Tujuan  dari  penelitian  ini  adalah  untuk  membuktikan  bahwa  breakwater  yang ada saat ini tidak dapat mereduksi tinggi gelombang yang terjadi dengan analisis hidrodinamika  dan  transformasi  gelombang.  Pemodelan  hidrodinamika  di Pelabuhan  Gunaksa  dilakukan  dengan  cara  mensimulasi  gelombang  dan melakukan  kalibrasi  data  pasang  surut  pemodelan  dengan  pengamatan.  Hasil dari  simulasi  ini  menunjukkan  bahwa  saat  pasang  purnama  arus  bergerak daribarat daya ke timur laut dengan elevasi muka air tertinggi 0,9 m dan tinggi gelombang  sebesar  0,56–0,64  m.  Pada  kondisi  pasang  perbani  arus  bergerak dari  timur  laut  ke  barat  daya  dengan  elevasi  muka  air  terendah  -1,02  m  dan tinggi gelombang sebesar 0,48 – 0,56 m.Tinggi gelombang yang terjadi melebihi 0,5m, oleh sebab itu Pelabuhan Gunaksa membutuhkan re-design breakwater.
Kata kunci : pemodelan, hidrodinamika, transformasi gelombang

## ABSTRACT
Gunaksa Harbour  Klungkung, Bali needs breakwater to protect the harbour from high waves. Former analysis shows that the existing breakwater cannot reduce the height  waves.  The  purpose  of  this  research  is  to  proves  that  the  existing breakwater  cannot  reduce  the  height  waves  by  hydryodynamic  conditions  and wave  transformation..  Hydrodynamic  modeling  of  Gunaksa  Harbour,  simulating waves  and  calibrated  against  data  from  a  substantial  field.  Modeling  results indicated that atspring tide, the current waves from southwest to northeast with the highest water level 0,9 m with the wave height is 0,56- 0,64 m, while at neap tide, the current waves from northeast to southwest with the lowest water level- 1,02 m with the wave height is 0,48 – 0,56 m. The wave height requirements for harbour  is  0,5  m,  however  the  wave  height  is  higher  than  0,5  m,  therefore Gunanksa Harbour could need re-design breakwater.
Keywords : modeling, hydrodynamic, wave transformation

## 1 . PENDAHULUAN
Kabupaten Klungkung merupakan salah satu dari sembilan kabupaten/kota di Propinsi Bali. Kabupaten  Klungkung  dibagi  dalam  tiga wilayah  kecamatan  di  daratan  Pulau  Bali  dan  satu wilayah  kecamatan  di  Kepulauan  Nusa  Penida.  Pertumbuhan  dan  perkembangannya, Kecamatan  Nusa  Penida  cukup  tertinggal  jauh  dibandingkan  dengan  tiga  kecamatan  yang berada  di  Klungkung  daratan.  Kendala  utamanya  adalah  masalah  aksesibilitas.  Pemerintah Kabupaten Klungkung membangun Pelabuhan Angkutan Sungai dan Penyeberangan (ASDP) yang memiliki rute dari Padang Bai sampai ke Lembar yang dimulai pada Oktober 2003 dan selesai pada bulan Juli 2006.
Pelabuhan ASDP Klungkung ini memiliki tinggi gelombang yang dominan terutama dari arah timur, tenggara, timur, yaitu 1,1 − 1,4 m. Kriteria tinggi gelombang pada pelabuhan memiliki tinggi gelombang maksimum 0,5 m. Tinggi gelombang dominan yang terjadi melebihi 0,5 m yang  mengakibatkan  kapal  tidak  dapat  berlabuh  dengan  baik.  Pelabuhan  ASDP  ini memerlukan breakwater untuk menjaga area sekitar pelabuhan dari gelombang.
Analisis  diperlukan  untuk  mengetahui  bahwa breakwateryang  ada  mampu  melindungi pelabuhan  dari  gelombang.  Penelitian  dilakukan  dengan  analisis  hidrodinamika  dan transformasi gelombang.

## 2. TINJAUAN PUSTAKA

## 2.1 Pemodelan Hidrodinamika 2DH

## 1. Flow Model FM
MIKE 21 memiliki kemampuan untuk memodelkan fenomena-fenomena yang terkait dengan sungai, danau, teluk, pantai, dan laut. MIKE 21 Flow Model FM adalah satu sistem modeling berbasis  pada  satu  pendekatan flexible  mesh  yang  dikembangkan  untuk  aplikasi  didalam oceanografi,  rekayasa  pantai,  dan  muara  sungai.  Flow  Model  FMitu  sendiri  terdiri  atas beberapa  modul,  diantaranya Hydrodinamic  Module.  Hydrodinamic  Module  adalah  model matematika  untuk  menghitung  perilaku hidrodinamika air  terhadap  berbagai  macam  fungsi gaya. Persamaan 1 , Persamaan 2 , dan Persamaan 3 berikut menggambarkan aliran dan perbedaan muka air.
𝜕𝜁 𝜕𝑡 + 𝜕𝑝 𝜕𝑥 + 𝜕𝑞 𝜕𝑦 = 𝜕𝑑 𝜕𝑡 …(1) 𝜕𝑝 𝜕𝑡 + 𝜕 𝜕𝑥 ( 𝑝 2 ℎ ) + 𝜕 𝜕𝑦 ( 𝑝𝑞 ℎ ) + 𝑔ℎ 𝜕𝜁 𝜕𝑥 + 𝑔𝑝√𝑝 2 +𝑞 2 𝐶 2 .ℎ 2 − 1 𝜌 𝑤 [ 𝜕 𝜕𝑥 (ℎ𝜏 𝑥𝑥 ) + 𝜕 𝜕𝑦 (ℎ𝜏 𝑥𝑦 )] − Ω 𝑞 − 𝑓𝑉𝑉 𝑥 + ℎ 𝜌 𝑤 𝜕 𝜕𝑥 (𝑝 𝑎 ) = 0 …(2)
𝜕𝑞 𝜕𝑡 + 𝜕 𝜕𝑦 ( 𝑞 2 ℎ ) + 𝜕 𝜕𝑥 ( 𝑝𝑞 ℎ ) + 𝑔ℎ 𝜕𝜁 𝜕𝑦 + 𝑔𝑝√𝑝 2 +𝑞 2 𝐶 2 .ℎ 2 − 1 𝜌 𝑤 [ 𝜕 𝜕𝑦 (ℎ𝜏 𝑥𝑥 ) + 𝜕 𝜕𝑥 (ℎ𝜏 𝑥𝑦 )] − Ω 𝑝 − 𝑓𝑉𝑉 𝑦 + ℎ 𝜌 𝑤 𝜕 𝜕𝑥𝑦 (𝑝 𝑎 ) = 0 …(3)
dengan: ℎ(𝑥, 𝑦, 𝑡) = kedalaman air [m], 𝑑(𝑥, 𝑦, 𝑡) = kedalaman air dalam berbagai waktu [m], 𝜁(𝑥, 𝑦, 𝑡) = elevasi permukaan [m],
𝑝, 𝑞(𝑥, 𝑦, 𝑡) = flux density dalam arah 𝑥 dan 𝑦 [m 3 /s/m] = (uh,vh); (u,v) = depth average velocity dalam arah 𝑥 dan 𝑦 ,
𝐶(𝑥, 𝑦) = tahanan Chezy [m 1/2 /s],

## Analisis Hidrodinamika dan Transformasi Gelombang bagi Re-design Breakwater di Pelabuhan Klungkung Bali
𝑔 = kecepatan gravitasi [m/s 2 ], 𝑓(𝑉) = faktor gesekan angin,
𝑉, 𝑉 𝑥 , 𝑉 𝑦 (𝑥, 𝑦, 𝑡) = kecepatan angin dalam arah 𝑥 dan 𝑦 [m/s],
Ω(𝑥, 𝑦) = parameter Coriolis [s -1 ], 𝑝 𝑎 (𝑥, 𝑦, 𝑡) = tekanan atmosfer [kg/m/s 2 ], 𝜌 𝑤 = berat jenis air [kg/m 3 ], (𝑥, 𝑦) = kordinat ruang [m], 𝑡 = waktu [s], 𝜏 𝑥𝑥 , 𝜏 𝑥𝑦 , 𝜏 𝑦𝑦 = komponen effectiveshear stress.

## 2. Spectral Waves FM
Prediksi  gelombang  dan  analisis  gelombang  dalam  skala  regional  dan  skala  lokal  serta perhitungan  transportasi  sedimen  yang  sebagian  besar  ditentukan  oleh  kondisi  gelombang dan wave-induced  currents.  Wave-incuded  current    disebabkan  oleh    gradient  radiation stressed    yang  terjadi  pada  surf  zone.  Koordinat  kartesius  keseimbangan  gaya  gelombang dapat dilihat pada Persamaan 4 dan koordinat spherical gaya gelombang dapat dilihat pada Persamaan 5 .
𝜕𝑁 𝜕𝑡 + ∇ ∙ (𝑣⃗𝑁) = 𝑆 𝜎 …(4)
𝑁̂ = 𝑁𝑅 2 cos ∅ = 𝐸𝑅 2 cos ∅ 𝜎 …(5)
dengan:
𝑁 (𝑥⃗, 𝜎, 𝜃, 𝑡) =  rapat gaya, 𝑡 =  waktu [detik], 𝑥⃗(𝑥, 𝑦) =  kordinat kartesian, 𝑣⃗(𝑐 𝑥 , 𝑐 𝑦 , 𝑐 𝜎 , 𝑐 𝜃 ) =  kecepatan propagansi grup gelombang empat dimensi [m/s],
𝑆 = source,
𝑥⃗(𝜙, 𝜆) =  kordinat spherical, dimana 𝜙 = 𝑙𝑎𝑡𝑖𝑡𝑢𝑑𝑒 dan 𝜆 = 𝑙𝑜𝑛𝑔𝑖𝑡𝑢𝑑𝑒 , 𝐸 =  kecepatan propagansi grup gelombang empat dimensi [m/s],
𝑅 = source.

## 2.2 Analisis Hidrodinamika dan Transformasi Gelombang
Berdasarkan  penelitian  Fadli  (2014)  kondisi  hidrodinamika  dilihat  pada  dua  kondisi,  yaitu pada  kondisi  pasang  purnama  dan  kondisi  pasang  perbani.  Pada  kondisi  pasang  purnama akan menghasilkan nilai pasang surut tertinggi, hal ini disebabkan karena posisi bulan, bumi, dan  matahari  yang  segaris  sehingga  gaya  tarik  gravitasi  dalam  keadaan  maksimum, sedangkan pada kondisi pasang perbani akan menghasilkan nilai pasang surut terendah, hal ini disebabkan karena posisi bumi, bulan, dan matahari yang membentuk sudut 90 o sehingga gaya  tarik  gravitasinya  minimum.  Pada  pemodelan  transformasi  gelombang  dapat  dilihat elevasi muka air di mulut kolam pelabuhan.

## 3. METODOLOGI PENELITIAN

## 3.1   Lokasi Penelitian
Penelitian  dilakukan  pada  kawasan  pelabuhan  Gunaksa  Kabupaten  Klungkungsecara geografis terletak di lintang -8,573311 o dan bujur 115,433595 o ( Gambar 1 ).
Gambar 1. Lokasi penelitian (Sumber: Google Earth, 2018)

## 3.2   Alur Penelitian
Penelitian ini dilakukan untuk menentukan rumusan masalah yang terjadi pada breakwater di Pelabuhan Gunaksa. Data-data yang diperlukan dalam penelitian adalah data batimetri, data pasang surut, data angin, dan data gelombang. Data-data tersebut diperlukan dalam analisis hidrodinamika  dan  transformasi  gelombang.  Validasi  data  pasang  surut  dan  pemodelannya mengikuti  fase  dan  amplitudo  pasang  surut  pengamatan.  Pola  hidrodinamika  disajikan berupa  elevasi  muka  air  dan  arah  arus.  Transformasi  gelombang  dimodelkan  dengan bantuan software  MIKE  21  modul  Spectral  Waves  FM.  Hasil  dari  pemodelan  transformasi gelombang berupa arah arus dan tinggi gelombang.

## 4. PENGOLAHAN DATA DAN PEMBAHASAN

## 4.1   Pengolahan Data Validasi
Berikut ini merupakan hasil pengolahan data pasang surut ( Gambar 2) .
Gambar 2. Kurva pasang surut hasil pengamatan
Gambar 2 menunjukkan bahwa pada rentang 8 Juli hingga 22 Juli tinggi muka air terbesar 1,1  m  terjadi  pada  tanggal  12  Juli  dan 13  Juli,  sedangkan  tinggi  muka  air  terendah-1,1  m terjadi pada tanggal 12 Juli dan 13 Juli.Jenis pasang surutnya termasuk jenis pasang surut campuran semi diurnal. Hal ini ditandai dengan terjadi satu kali pasang dan satu kali surut pada  hari  yang  sama  namun  beberapa  kalijuga  terjadi  dua  kali  pasang  dan  dua  kali  surut pada hari yang sama.

## 4.2   Validasi Data
Berikut ini  merupakan parameter yang digunakan Tabel 1 dan validasi data pasang surut pengamatan dan pemodelan ( Gambar 3) .

## Analisis Hidrodinamika dan Transformasi Gelombang bagi Re-design Breakwater di Pelabuhan Klungkung Bali
Tabel 1. Parameter yang Digunakan dalam Setup PemodelanHidrodinamika Parameter Diterapkan dalam Pemodelan Mesh and Bathymetry  Hasil interpolasi mesh Simulation Period No. of time step =360 Time step interval =3600 s Simulation start date =7/8/2006 Simulation end date =7/22/2006 Eddy Viscosity Constant Value = 0,28 Bed Resistance Manning = 37 m 1/3 /s Boundary Condition
Hasil prediksi tiap boundary

## 4.3   Simulasi Pemodelan
Dalam  pemodelandilakukan  kalibrasi  pemodelan  agar  mendapat  hasil  yang  valid  antara pengamatan dengan pemodelan. Percobaan pemodelan dengan berbagai nilai manning atau Bed Resistence  menghasilkan rekapitulasi yang dapat dilihat pada Tabel 2 .
Tabel 2. Rekapitulasi Hasil Kalibrasi Bed Resistence Error Method Elevasi Muka Air 28 RMSE 0,2957 Korelasi(%) 81,4605 30 RMSE 0,2982 Korelasi(%) 81,2799 32 RMSE 0,2960 Korelasi(%) 81,4301 35 RMSE 0,2967 Korelasi(%) 81,4632 37 RMSE 0,2954 Korelasi(%) 81,4853
Berdasarkan  hasil  tersebut  dipilih  nilai manning  =  37  m 1/3 karena  memiliki  nilai  RMSE  dan nilai korelasi yang cukup tinggi dibandingkan yang lainnya. Hasil kalibrasi dapat dilihat pada Gambar 3 .
Gambar 3. Grafik perbandingan pasang surut pengamatan dan pemodelan
Gambar  3 menunjukan  bahwa  pasang  surut  pengamatan  dan  pemodelan  hampir  sama, namun pada beberapa bagian time step terdapat elevasi muka air yang berbeda.

## 4.4  Hasil Analisis Hidrodinamika
Berikut  ini  adalah  hasil  dari  simulasi  hidrodinamika  pada  saat  angin  barat  (Desember  s.d Februari)  dan  angin  timur  (Juni  s.d  Agustus)  dengan  kondisi  pasang  purnama  dan  kondisi pasang perbani.

## 1. Kondisi pasang purnama pada saat angin barat (Desember s.d Februari)
Gambar 4. Kondisi pasang purnama saat angin barat pada tanggal 12 Juli 2006
Gambar  4 menunjukan  kondisi  pasang  purnama  saat  angin  barat  yang  menunjukan beberapa kondisi yaitu padasaat surut menuju pasang yang terjadi pada tanggal 11 Juli 2006 pukul 23.00, pasang tertinggi yang terjadi pada tanggal 12 Juli 2006 pukul 02.00  , pasang menuju surut yang terjadi pada tanggal 12 Juli 2006 pukul 05.00, dan surut terendah yang terjadi  pada  tanggal  12  Juli  2006  pukul  08.00,  arus  di  sekitar  pelabuhan  ASDP  Klungkung Bali  bergerak  dari  arah  barat  daya  ke  arah  timur  laut.  Elevasi  muka  air  pada  saat  surut menuju  pasang  sebesar  0,54  m,  pada  saat  pasang  tertinggi  sebesar  0,99  m,  pada  saat pasang menuju surut sebesar 0,55 m, dan pada saat surut terendah sebesar-0,46 m.
U

## Analisis Hidrodinamika dan Transformasi Gelombang bagi Re-design Breakwater di Pelabuhan Klungkung Bali

## 2. Kondisi pasang perbani pada saat angin barat (Desember s.d Februari)

## Gambar 5. Kondisi pasang perbani saat angin baratpada tanggal 12 Juli 2006
Gambar 5 menunjukan kondisi pasang perbanisaat angin barat yang masing-masing terjadi pada tanggal 12 Juli 2006 dengan beberapa kondisi, yaitu surut menuju pasang yang terjadi pada  pukul  08.00,  pasang  tertinggi  yang  terjadi  pada  pukul  11.00,  pasang  menuju  surut yang  terjadi  pada  pukul  14.00,  dan  surut  terendah  yang  terjadi  pada  pukul  17.00,  arus  di sekitar  pelabuhan  ASDP  Klungkung  Bali  bergerak  dari  arah  timur  laut  ke  arah  barat  daya. Elevasi  muka  air  pada  saat  surut  menuju  pasang  sebesar  -0,46  m,  pada  saat  pasang tertinggi sebesar -1,02 m, pada saat pasang menuju surut sebesar -0,71 m, dan pada saat surut terendah sebesar -0,2 m.

## 3. Kondisi pasang purnama pada saat angin timur (Desember s.d Februari)
Gambar 6. Kondisi pasang purnama saat angin timurpada tanggal 12 Juli 2006
Gambar 6 menunjukan kondisi pasang purnama pada saat angin timur yang masing-masing terjadi pada tanggal 12 Juli 2006 dengan beberapa kondisi yaitu surut menuju pasang yang
U
U
terjadi pada  pukul 00.00 , pasang tertinggi yang  terjadi pada pukul 03.00, pasang menuju surut yang terjadi pada pukul 06.00, dan surut terendah yang terjadi pada pukul 09.00, arus di sekitar pelabuhan ASDP Klungkung Bali bergerak dari arah barat daya ke arah timur laut . Elevasi muka air pada saat surut menuju pasang sebesar 0,76 m, pada saat pasang tertinggi sebesar  0,94  m,  pada  saat  pasang  menuju  surut  sebesar  0,23  m,  dan  pada  saat  surut terendah sebesar -0,75 m.

## 4. Kondisi pasang perbani pada saat angin timur (Juni s.d Agustus)
Gambar 7. Kondisi pasang perbani saat angin timurpada tanggal 12 Juli 2006
Gambar 7 menunjukan kondisi pasang perbani pada saat angin timur yang masing-masing terjadi pada tanggal 12 Juli 2006 dengan beberapa kondisi yaitu surut menuju pasang yang terjadi  pada  pukul  08.00,  pasang  tertinggi  yang  terjadi  pada  pukul  11.00,  pasang  menuju surut yang terjadi pada pukul 14.00, dan surut terendah yang terjadi pada pukul 17.00, arus di sekitar pelabuhan ASDP Klungkung Bali bergerak dari arah timur laut ke arah barat daya. Elevasi  muka  air  pada  saat  surut  menuju  pasang  sebesar  -0,46  m,  pada  saat  pasang tertinggi sebesar -1,02 m, pada saat pasang menuju surut sebesar -0,71 m, dan pada saat surut terendah sebesar -0,20 m.

## 4.5  Hasil Analisis Transformasi Gelombang
Berikut  ini  adalah  hasil  dari  simulasi  transformasi  gelombang  pada  saat  angin  barat (Desember s.d Februari) dan angin timur (Juni s.d Agustus) dengan kondisi pasang purnama dan kondisi pasang perbani.

## U

## Analisis Hidrodinamika dan Transformasi Gelombang bagi Re-design Breakwater di Pelabuhan Klungkung Bali
1. Kondisi pasang purnama pada saat angin barat (Desember s.d Februari)

## Gambar 8. Transformasi Gelombang pada pasang purnama
saat angin barat pada tanggal 12 Juli 2006

## 2. Kondisi pasang perbani pada saat angin barat (Desember s.d Februari)
Gambar 9. Transformasi gelombang dengan kondisi pasang perbani saat angin barat pada tanggal 12 Juli 2006
U
U
3. Kondisi pasang purnama pada saat angin timur (Juni s.d Agustus)

## Gambar 10. Transformasi gelombang dengan kondisi pasang purnama saat angin timur pada tanggal 12 Juli 2006

## 4. Kondisi pasang perbani pada saat angin timur (Juni s.d Agustus)
Gambar 11. Transformasi gelombang dengan kondisi pasang perbani saat angin timur pada tanggal 12 Juli 2006
Gambar  8  s.d  Gambar  11 menunjukan  transformasi  gelombang  yang  terjadi  di  sekitar mulut  kolam  pelabuhan.  Pada  saat  pasang  purnama  baik  pada  saat  angin  barat  maupun angin  timur  arus  bergerak  dari  arah  barat  daya  ke  arah  timur  laut  yang  memiliki  tinggi gelombang  sebesar  0,56  m − 0,64  m,  sedangkan  pada  kondisi  pasang  perbani  baik  saat angin  barat  maupun  angin  timur  arus  bergerak  dari  arah  timur  laut  ke  arah  barat  daya memiliki tinggi gelombang sebesar 0,48 m − 0,56 m.
U
U

## Analisis Hidrodinamika dan Transformasi Gelombang bagi Re-design Breakwater di Pelabuhan Klungkung Bali

## 5. KESIMPULAN

## 5.1 Kesimpulan
1. Hasil dari validasi pasang surut pemodelan memiliki tingkat kesesuaian terhadap kondisi real atau kondisi di lapangan yang dilihat dari nilai RMSE ( Root Mean Square Error) dan nilai korelasinnya. Pada penelitian ini digunakan Bed Resistencedengan manningsebesar 37  m 1/3 dengan  RMSE  sebesar  0,2954  m  atau  sekitar  29  cm  dan  korelasi  sebesar 81,4853%.
2. Kondisi  arus  di  perairan  Pelabuhan  Klungkung  Bali  pada  saat  pasang  purnama,  arus bergerak meninggalkan pelabuhan dari arah barat daya ke arah timur laut yang memiliki elevasi muka air tertinggi 0,9 m, sedangkan pada kondisi pasang perbani, arus bergerak meninggalkan  pelabuhan  dari  arah  timur  laut  ke arah  barat  daya  yang  memiliki elevasi muka air terendah -1,02 m.
3. Tinggi gelombang terbesar yang terjadi di sekitar mulut pelabuhan pada kondisi pasang purnama baik saat angin barat maupun angin timur memiliki tinggi gelombang sebesar 0,56 m – 0,64 m, sedangkan pada kondisi pasang perbani baik saat angin barat maupun angin timur memiliki tinggi gelombang sebesar 0,48 m − 0,56 m.
4. Dengan  tinggi  gelombang  yang  terjadi  di  sekitar  mulut  kolam  pelabuhan  saat  pasang purnama sebesar 0,56 m − 0,64 m dan saat pasang perbani sebesar 0,48 m − 0,56 m; tinggi gelombang tersebut melebihi dari syarat tinggi gelombang untuk pelabuhan yaitu 0,5 m; oleh karena itu Pelabuhan ASDP Klungkung membutuhkan re-design breakwater.

## UCAPAN TERIMA KASIH
Ucapan terima kasih kepada PT LAPI ITB yang telah memberikan data-data yang dibutuhkan pada penelitian ini.

## DAFTAR RUJUKAN
Fadli,  M.,  Radjawane, I.  dan Susanna.  (2014).  Pemodelan  Hidrodinamika  di  Perairan  Teluk Ambon. Oseanografi , ___ (___), 13-17. PT LAPI ITB. (2006). Laporan Akhir Pembangunan Pelabuhan ASDP Gunaksa. Bali Triatmodjo, B. (2016). Teknik Pelabuhan. Yogyakarta: Beta Offset. Triatmodjo, B. (2016). Teknik Pantai. Yogyakarta: Beta Offset.