Peramalan Gelombang di Perairan Kabupaten Indramayu dengan Pemodelan Numerik SWAN 41.01A

## FAJAR WIDHIARNO, YATI MULIATI
Jurusan Teknik Sipil, Institut Teknologi Nasional-Bandung e-mail: fajar_widhiarno@yahoo.co.id

## ABSTRAK
Karakteristik gelombang merupakan hal yang sangat penting dalam perencanaan bangunan pantai. Melakukan pengukuran gelombang dalam jangka waktu yang panjang membutuhkan biaya yang sangat mahal. Diperlukan alternatif lain untuk meramalkan  gelombang  yaitu  dengan  menggunakan  data  angin.  Belum  ada lembaga  yang  melakukan  peramalan  gelombang  secara  berkelanjutan  di Indonesia. Tujuan dari penelitian ini adalah untuk mengetahui tinggi gelombang signifikan (𝐻 𝑠 ) , periode gelombang signifikan (𝑇 𝑠 ) dan arah gelombang dengan menggunakan  metode  pemodelan  numerik  SWAN  41.01A  yang  selanjutnya divalidasi  dengan  menggunakan  metode  SPM.  Hasil  dari  peramalan  gelombang didapat nilai (𝐻 𝑠 ) setinggi 0,4-0,5 m dengan (𝑇 𝑠 ) selama 2-3 detik dengan musim timuran  yang  mendominasi  arah  angin,  tinggi  gelombang  maksimum  setinggi 0,8-0,9 m dengan periode 2-2,5 detik, dan adanya pembelokan arah gelombang diakibatkan oleh pendangkalan kedalaman laut.
Kata Kunci: peramalan gelombang, SWAN, tinggi gelombang signifikan, periode signifikan, arah gelombang.

## ABSTRACT
Characteristic  of  wave  is  the  important  thing  in  seashore  building  construction planning. The cost for the wave forecast in long duration is very expensive. Wind data is an alternative method for the wave forecast. There is no institution yet to forecast  waves  sustainably  in  indonesia.  The  purpose  of  this  research  is  to determine  the  significant  wave  height (𝐻 𝑠 ) ,  significant  wave  period (𝑇 𝑠 ) and direction of the waves using SWAN 41.01A numerical modelling method which is then  validated  using  SPM  method.  The  results  of  forecast  obtains  values  wave height (𝐻 𝑠 ) 0.4-0.5  m  with (𝑇 𝑠 ) during  2-3  second  with  easterlies  season dominate the wind direction, the maximum wave height 0.8-0.9 m with a period of 2-2.5 second, and bending the direction of the waves caused by silting ocean depths.
Keywords: wave forecast, SWAN, significant wave height, significant period, wave direction.

## 1. PENDAHULUAN
Gelombang laut dapat terjadi diakibatkan oleh berbagai macam faktor salah satunya adalah gelombang  yang  disebabkan  oleh  angin  (gelombang  angin).  Gelombang  merupakan  hal yang  penting  dalam  merencanakan  suatu  bangunan  di  pesisir  pantai.  Di  Indonesia  belum ada  instansi  atau  lembaga  yang  melakukan  peramalan  gelombang  secara  berkelanjutan. Diperlukan  biaya  yang  sangat  mahal  untuk  melakukan  pengukuran  gelombang  dalam jangka waktu yang panjang maka dari itu perlu alternatif lain untuk meramalkan gelombang yaitu menggunakan data angin. Salah satu metode yang dapat meramalkan gelombang dari data  angin  adalah Simulating  Wave  Nearshore  (SWAN)  yang  dikembangkan  oleh  DELFT University dari Belanda.

## 2. TINJAUAN PUSTAKA

## 2.1 Gelombang
Menurut  Triatmodjo  (1999)  gelombang  adalah  pergerakan  naik  dan  turunnya  air  dengan arah  tegak  lurus  permukaan  air  laut  yang  membentuk  kurva/grafik  sinusoidal.  Gelombang yang  terjadi  di  lautan  dapat  diklasifikasikan  menjadi  beberapa  macam  berdasarkan pembangkitnya, salah satunya pembangkitan gelombang laut yang terjadi disebabkan oleh angin. Mekanisme terjadinya gelombang oleh angin terjadi ketika angin yang berhembus di atas  permukaan  air  akan  memindahkan  energinya  ke  air.  Kecepatan  angin  akan menimbulkan  tegangan  pada  permukaan  laut,  sehingga  permukaan  air  yang  semulanya tenang  akan  terganggu  dan  timbul  riak  gelombang  kecil  di  atas  permukaan  air.  Apabila kecepatan  angin  bertambah,  riak  tersebut  menjadi  semakin  besar,  dan  apabila  angin berhembus  terus  akhirnya  akan  terbentuk  gelombang.  Semakin  lama  dan  semakin  kuat angin berhembus maka semakin besar gelombang yang terbentuk.

## 2.2 Angin
Menurut  Triatmodjo  (1999)  angin  merupakan  sirkulasi  yang  kurang  lebih  sejajar  dengan permukaan  bumi.  Angin  terjadi  akibat  adanya  perubahan  ataupun  perbedaan  suhu  antara suatu  tempat  dengan  tempat  yang  lain.  Salah  satu  contoh  yang  dapat  diambil  adalah perubahan  suhu  yang  terjadi  antara  daratan  dan  lautan.  Daratan  cenderung  lebih  cepat menerima dan melepaskan panas. Oleh sebab itu, maka siang hari terjadi angin laut yang diakibatkan oleh naiknya udara daratan yang digantikan oleh udara dari laut. Pada malam hari terjadi sebaliknya, yaitu terjadi angin darat yang diakibatkan oleh naiknya udara di laut dan digantikan oleh udara dari darat.

## 2.3 Peramalan Gelombang
Peramalan gelombang dapat dilakukan dengan berbagai macam metode untuk mendapatkan tinggi  dan  periode  gelombang  signifikan  diantaranya,  Metode Groen  &  Dorrestein,  Metode Shore  Protection  Manual  (SPM)  dan  Permalan  Gelombang  dengan  Pemodelan  Numerik SWAN.  Dalam  penelitain  ini  untuk  meramalkan  gelombang  menggunakan  metode  SWAN ( Simulating  Waves  Nearshore)  adalah  model  numerik  gelombang  generasi  ketiga  yang digunakan untuk memperkirakan parameter gelombang di daerah perairan skala kecil seperti pantai, danau, dan estuari dengan menggunakan data angin, kondisi dasar (batimetri), dan arus  serta  menggunakan  persamaan  dasar  kesetimbangan  aksi  ( action  balance)  dengan source dan sink (Booij, dkk., 1999 dan The SWAN Team, 2006). Model  SWAN  dibangun  berdasarkan  persamaan  kesetimbangan  gerak  gelombang  yang mendeskripsikan perubahan spektrum gelombang menggunakan persamaan kesetimbangan energi gelombang yang menjelaskan evolusi spektrum gelombang laut. Persamaan 1 dan
Persamaan  2 menjelaskan  tentang  keseimbangan  gerak  gelombang  dalam  koordinat Kartesian (Holthuijsen, 2007), adalah:
     total y x S N C N C N C y N C x N t                ... (1) br ds b ds w ds nl nl in total S S S S S S S , , , 4 3      
… (2)
dengan :
N t   = perubahan densitas aksi terhadap waktu,
N C x x 
 = penjalaran energi gelombang dalam ruang waktu dalam arah 𝑋 , N C y y   = penjalaran energi gelombang dalam ruang waktu dalam arah 𝑌 ,
N C     = pergeseran frekuensi karena variasi kedalaman,
N C   
 = refraksi akibat variasi kedalaman. in S = merupakan nilai input angin,
3 nl S dan 4 nl S
= merupakan suku nonlinear,
w ds S , , b ds S , dan br ds S , = merupakan suku disipasi gelombang.
Dalam  Adiprabowo  (2015)  dijelaskan  bahwa  evolusi  spektrum  gelombang  ini  tedapat  tiga fungsi  sumber  yang  mempresentasikan  input  angin,  interaksi  non-linear  dan  disipasi. Spektrum  gelombang  merupakan  distribusi  energi  terhadap  periode  atau  frekuensi amplitudo.  Bentuk  spektrum  dipengaruhi  oleh  durasi  angin  bertiup,  kecepatan  angin,  arah angin,  dan fetch  sesuai  dengan  energi  dan  periode  dari  masing-masing  gelombang  pada suatu  kondisi  perairan.  Pada  penelitian  ini  spektrum  yang  digunakan  yaitu  spektrum JONSWAP seperti tertulis dalam Persamaan 3 yaitu:
      f e f g f E fp f    4 25 . 1 5 4 2
2         … (3)
dengan:
  f E = densitas variansi (m 2 /s),
 = konstanta tak berdimensi,  = 0,0081,
g = percepatan gravitasi (m 2 /s),
f = frekuensi (Hz),
fp = frekuensi puncak dari Pierson-Moskowitsz,
 = faktor penguat spektrum.
Input  angin  yang  digunakan  untuk  mejalankan  model  SWAN  adalah  kecepatan  angin  pada ketinggian  10  m  dari  permukaan  laut (𝑈 10 ) sedangkan  untuk  keperluan  perhitungan  input
tersebut  diubah  menjadi  kecepatan  gesekan  (fiction  velocity)  dengan  perumusan (Holthuijsen, 2007) seperti pada Persamaan 4 berikut:
𝑢 ∗2 =   𝐶 𝐷 𝑈 10 2
… (4)
dengan 𝐶 𝐷 adalah drag coeficient (Wu, 1982 di dalam Holthuijsen, 2007) ditulis dalam Persamaan 5 berikut:
𝐶 𝐷 (𝑈 10 ) = 1.2875 ∗ 10
untuk 𝑈 10 < 7.5 m/s 𝐶 𝐷 (𝑈 10 ) = (0.8 + 0.065 𝑠 𝑚 ∗ 𝑈 10 ) untuk 𝑈 10 ≥ 7.5 m/s … (5)
Dalam  Aisyiah  (2011)  dijelaskan  bahwa  disipasi  energi  gelombang  di  dalam  model  SWAN dipresentasikan  oleh  suku sink (𝑆 𝑑𝑠 ) .  Suku  tersebut  merupakan  penjumlahan  dari  ketiga kontributor  lainnya,  yaitu  disipasi whitecapping (𝑆 𝑑𝑠,𝑤 ) ,  disipasi  gesekan  dasar (𝑆 𝑑𝑠,𝑏 ) ,  dan disipasi gelombang pecah yang disebabkan kedalaman (𝑆 𝑑𝑠,𝑏𝑟 ) . Disipasi whitecapping (𝑆 𝑑𝑠,𝑤 ) dikontrol oleh kecuraman gelombang berdasarkan model pulse- based (Hasselmann, 1974 dalam SWAN UM, 2008), yang diadaptasi oleh kelompok WAMDI (1988) pada Persamaan 6 berikut:
         , ~ ~ , , E k k
S w ds   
… (6)
dengan:
𝑆 𝑑𝑠,𝑤 (𝜎, 𝜃) = densitas disipasi whitecapping,
Γ
= koefisien kemiringan,
𝜎̃
= frekuensi rata-rata, 𝑘 = bilangan gelombang,
𝑘̃ = bilangan gelombang rata-rata,
    , E = energi densitas.
dimana 𝛤 merupakan koefisien kemiringan yang diestimasi dengan menutup keseimbangan energi pada keadaan fully developed sea yang bergantung pada input angin. Disipasi gesekan dasar (𝑆 𝑑𝑠,𝑏 ) adalah mekanisme dominan dari disipasi di continental shelf seas  dengan  dasar  berpasir  (Bertotti  dan  Cavaleri,  1994  dalam  SWAN  UM,  2008)  yang dinyatakan dalam Persamaan 7 berikut:
     , ) ( sinh 2 2 2 , E kd g C S b b ds          
… (7)
dengan:
𝑆 𝑑𝑠,𝑏 = disipasi gesekan dasar, 𝐶 𝑏 = koefisien gesekan dasar,
𝜎
= frekuensi rata-rata, 𝑘 = bilangan gelombang, 𝑑
= bilangan gelombang rata-rata.
Dimana C b adalah  koefisien  gesekan  dasar  yang  bergantung  pada  gerakan  orbital  yang ditunjukan oleh U rms .
Disipasi  gelombang  pecah  karena  kedalaman  menurut  Elderberky  dan  Battjes  (SWAN  UM, 2008) untuk memformulasikan versi spektrum dari model bore-based (Battjes dan Janssen, 1978  dalam  SWAN  UM,  2008)  tetapi  tidak  merubah  bentuk  spektrumnya  dengan Persamaan 8 yang digunakan dalam SWAN adalah:
        , , , E E D S tot tot br ds 
… (8)
dengan:
𝑆 𝑑𝑠,𝑏𝑟 (𝜎, 𝜃) = densitas disipasi gelombang pecah karena kedalaman,
𝐷 𝑡𝑜𝑡 = rata-rata disipasi energi per unit area horizontal, 𝐸 𝑡𝑜𝑡
= energi gelombang total.

## 3. METODOLOGI PENELITIAN

## 3.1 Deskripsi Penelitian
Penelitian ini menggunakan sebuah software yang bernama SWAN dengan memasukan data angin dan data batimetri sebagai input datanya. SWAN ini memerlukan software tambahan untuk membantu proses running yaitu menggunakan Notepad++ untuk membuat script dan MATLAB untuk memplot hasil dari running SWAN. Pertama kali yang harus dilakukan yaitu menyiapkan  data  angin  dan  batimetri,  selanjutnya  membuat  sebuah script  yang  berisi sebuah  parameter  dan  kondisi  yang  akan  dimodelkan,  setelah  itu running  menggunakan SWAN  dan  hasil  dari  SWAN  diplot  menggunakan  MATLAB  agar  dapat  dilihat  hasil pemodelannya.

## 4. ANALISIS DATA DAN PEMBAHASAN

## 4.1 Pengolahan Data Angin
Dalam  penelitian ini  data  angin  yang  digunakan  yaitu  data  angin  yang  berbentuk  stationer yang  didapat  dari  stasiun  klimatologi  Penggung  di  bandara  Penggung  Cirebon.  Data  ini digunakan  karena  di  daerah  Kabupaten  Indramayu  tidak  tersedianya  stasiun  yang menyediakan data angin. Data angin yang digunakan yaitu data angin jam-jaman sepanjang 10  tahun  dari  tahun  1998-2007.  Berdasarkan  penelitian  yang  dilakukan  oleh  PPN  (2009) Kejawanan Cirebon didapat data sebagai berikut.
Tabel 1. Tabel Windrose Arah (Knot) <5 5-10 10-15 15-20 >20 Total <5 5-10 10-15 15-20 >20 Total Utara 2.918 1.433 96 13 0 4.460 3,33 1,63 0,11 0,01 0,00 5,09 Timur Laut 3.247 2.487 283 28 2 6.047 3,70 2,84 0,32 0,03 0,00 6,90 Timur 5.565 4.905 886 106 16 11.478 6,35 5,60 1,01 0,12 0,02 13,10 Tenggara 1.041 893 340 63 7 2.344 1,19 1,02 0,39 0,07 0,01 2,67 Selatan 1.390 1.671 1.058 739 111 4.969 1,59 1,91 1,21 0,84 0,13 5,67 Barat Daya 308 111 103 87 15 624 0,35 0,13 0,12 0,10 0,02 0,71 Barat 3.078 1.451 129 26 2 4.686 3,51 1,66 0,15 0,03 0,00 5,35 Barat Laut 1.209 294 22 9 1 1.535 1,38 0,34 0,03 0,01 0,00 1,75 Sumber : Laporan Rencana Pengembangan Pelabuhan Perikanan Nusantara  Kejawanan Cirebon (PPN Kejawanan Cirebon),  2009.
Jumlah Jam Persentase

## Gambar 1. Windrose
Dari  hasil  peramalan  yang  dilakukan  oleh  PPN  Kejawanan  arah  angin  dominan  terjadi  dari arah timur dengan kecepatan 10-15 knot atau 5-7,5 m/s.

## 4.2 Pengolahan Data Batimetri
Dalam penelitian ini data batimetri yang digunakan yaitu dari General Bathymetric Chart of the  Oceans  (GEBCO).  Berdasarkan  GEBCO  (2014)  untuk  mendapatkan  data  batimetri pertama masuk ke browser dengan link https://www.bodc.ac.uk/data/online_delivery/gebco/ lalu  login.  Kedua  masukan  koordinat  pembatas  area  ke  dalam  kolom  yang  sudah  tersedia atau  dengan  cara  zoom  ke  area  yang  ditinjau  kemudian  tekan  kombinasi  tombol  shift+klik kiri  pada  mouse  lalu  drag  area  yang  akan  ditinjau  seperti  pada Gambar  2 .  Ketiga  pilih format data yang akan diunduh seperti pada Gambar 3 . Keempat klik add data to basket dan  klik view  basket  maka  tampilan  untuk  mengunduh  akan  seperti Gambar  4 lalu  klik download.  Setelah  diunduh  hasil  dari  GEBCO  berformat  ASC  bila  dibuka  dengan menggunakan Notepad++ akan seperti Gambar 5 .
Data  GEBCO  ini  dapat  langsung  digunakan  kedalam software  SWAN  tanpa  harus  ada pengolahan  terlebih  dahulu.  Untuk  jenis  format  data  tidak  hanya  ASC  saja  yang  dapat digunakan namun dapat juga menggunakan format .txt, .bot dan .pts.
Gambar 3. Format data output Gambar 4. Unduh data GEBCO
Gambar 5. Output data batimetri

## 4.3 Membuat Script
Dalam penelitian ini untuk membuat sebuah script menggunakan Notepad++ karena mudah pemakaiannya  dan  bentuknya  pun  rapi  dan  dilengkapi  banyak toolbar untuk  memudahkan pengoprasian, tidak seperti Notepad biasa tampilannya yang terlalu sederhana. Dalam membuat sebuah script harus mengikuti ketentuan-ketentuan yang telah ada didalam software  SWAN  atau  dapat  menggunakan  SWAN  User  Manual  dalam  SWAN  (1993-2015) sebagai panduan. Pada Gambar 6 adalah script yang digunakan pada penelitian ini.

## Gambar 6. Script penelitian
Pada  pembuatan script  ini  yang  perlu  diperhatikan  saat  memasukan  data  batimetri  pada bagian  ke-4  dan  ke-3  dimana  pada  bagian  ke-3  atau domain  model  harus  lebih  kecil  dari pada  bagian  ke-4  atau bathymetric  input.  Untuk  non-stationer  tidak  dapat  menggunakan script di atas karena untuk non-stationer memerlukan data hari, bulan dan tahun pada data angin  yang  akan  digunakan  demikian  juga  dengan  semua  fungsi  yang  ada  di  atas  akan berubah harus disesuaikan dengan kondisi non-stationer.
Setelah pembuatan script selesai save as dimana folder data batimetri dan angin disimpan usahakan  nama file  batimetri  sama  dengan  script  ini  lalu  ditambahakan  dengan  .swn  lalu save  misalnya  Indramayu.swn,  dengan  menambahkan  .swn  pada  nama  file  otomatis  file tersebut berubah menjadi swn yang selanjutnya digunakan untuk running di software SWAN.

## 4.4 Mengoperasikan SWAN
Untuk  melakukan  proses  running  SWAN  Pertama  buka  aplikasi  SWAN  lalu  ketik  “cd  lokasi penyimpanan folder swn” enter maka akan seperti Gambar 7 . Selanjutnya ketik “SWANRUN nama  file  dengan  format  swn” enter  maka  SWAN  akan  running  seperti  pada Gambar  8 . Setelah running  selesai  hasil  filenya  akan  secara  otomatis  tersimpan  di  folder  yang  sama dengan format .mat yang selanjutnya diplot dengan menggunakan MATLAB.

## Peramalan Gelombang di Perairan Kabupaten Indramayu dengan Pemodelan Numerik SWAN 41.01A
Hasil  dari running SWAN ini bila dibuka  menggunakan MATLAB akan  didapatkan data nilai DEPTH  (kedalaman  laut),  DIR  (arah  gelombang), 𝐻 𝑠𝑖𝑔 (tinggi  gelombang  signifikan), 𝐻 𝑢 (angin arah 𝑈 ), 𝐻 𝑣 (angin arah 𝑉 ), 𝑇𝑀 01 (periode gelombang signifikan), 𝑇𝑀𝑀 10 (periode energi gelombang), 𝑋 𝑝 (koordinat 𝑋 ) dan 𝑌 𝑝 (koordinat 𝑌 ).
Gambar 7. Lokasi penyimpanan folder

## Gambar 8. Running SWAN
4.5 Hasil Analisis Hasil  dari running  SWAN  selanjutnya  diplot  dengan  menggunakan  MATLAB  yang  akan didapatkan  sebuah  pemodelan  gelombang  dengan  gambar  dua  dimensi  yang menggambarkan  terjadinya  gelombang  di  perairan  Kabupaten  Indramayu  seperti  pada Gambar 9, Gambar 10 dan Gambar 11 .
Gambar 9. Batimetri
Gambar 10. H signifikan (𝑯 𝒔 )

## Peramalan Gelombang di Perairan Kabupaten Indramayu dengan Pemodelan Numerik SWAN 41.01A
Dari  hasil  pemodelan  yang  dilakukan  dengan  menggunakan software  SWAN  41.01A  dapat ditarik beberapa hasil analisis diantaranya:
1. Gelombang  maksimum  yang  setinggi  0,8-0,9  m  dengan  periode  selama  2-2,5  detik terjadi di letak geografis 107,7º BT dan -6,2º LS.
2. Pada  jarak  sejauh  44,528  km  dari  titik  awal  lokasi  penelitian  (-6,5º  LS  s/d  -6,1º  LS), tinggi  gelombang (𝐻 𝑠 ) yang  paling  dominan  setinggi  0,4-0,5  m  dengan  periode (𝑇 𝑠 ) selama 2-3 detik.
3. Pada  letak  geografis  108,2º-108,4º  BT  dan  6,2º-6,3º  LS  mengalami  pembelokan  arah yang disebabkan oleh pendangkalan kedalaman laut.
4. Dilihat  dari  arah  angin  yang  didominasi  dari  arah  timur  hal  ini  disebabkan  angin  yang bertiup  paling  besar  ketika  musim  timuran,  dimana  tekanan  udara  di  Asia  rendah  dan tekanan  udara  di  Australia  tinggi  yang  mengakibatkan  angin  bertiup  dari  Australia  ke Asia.
5. Validasi model SWAN ini dilakukan dengan cara membandingkan dengan metode SPM, dimana  gelombang  maksimum  hasil  SWAN  sebesar  0,8-0,9  m  dan  hasil  SPM  sebesar 0,875 m, hal ini menunjukan bahwa model SWAN dapat digunakan.
6. Metode  SPM  juga  menunjukan  gelombang  di  perairan  dalam,  untuk  itu  coba  diperiksa lokasi  gelombang  maksimum  pada  hasil  SWAN  yaitu  pada  kedalaman (𝑑)   =  30 m dengan menggunakan perhitungan berikut:
𝐿 0 =  1,56 𝑇 2 =  1,56 (2,5) 2 =  9,75 m 𝑑/𝐿 0 =  30/9,75  =  3,076 Plot  hasil 𝑑/𝐿 0 pada  tabel  Wiegel  akan  didapat  nilai 𝑑/𝐿 ,  karena  nilai 𝑑/𝐿 0 ≥  1 maka nilai 𝑑/𝐿 = 1 , mengingat ketentuan perairan dalam 𝑑/𝐿 > 0,5 , maka pada lokasi 𝑑 = 30 m tersebut benar sudah termasuk perairan dalam.

## 5. KESIMPULAN

## 5.1 Kesimpulan
Berdasarkan  hasil  penelitian  yang  dilakukan  untuk  memodelkan  gelombang  yang  terjadi  di perairan  Kabupaten  Indramayu  dengan  menggunakan  metode Simulating  Wave  Nearshore (SWAN)  dan  menggunakan  data  angin  tahun  1998-2007  dapat  diambil  beberapa  simpulan sebagai berikut:
1. Gelombang  maksimum  yang  terjadi  setinggi  0,8-0,9  m  dengan  periode  selama  2-2,5 detik.
2. Tinggi  gelombang (𝐻 𝑠 ) yang  paling  dominan  setinggi  0,4-0,5  m  dengan  periode (𝑇 𝑠 ) selama 2-3 detik.
3. Terjadinya pembelokan arah gelombang yang disebabkan oleh pendangkalan kedalaman laut.
4. Arah  angin  yang  paling  mendominasi  besarnya  angin  di  perairan  Kabupaten  Indramayu terjadi pada saat musim timuran.

## DAFTAR RUJUKAN
Asiyah, Z. (2011). Kajian Gelombang Laut di Perairan Selatan Provinsi Banten Menggunakan Model SWAN. Tugas Akhir, Program Studi Oseanografi, ITB.
Adiprabowo,  S.  R.  (2015).  Studi  Karakteristik  Gelombang  di  Perairan  Selat  Sunda.  Tugas Akhir, Program Studi Oseanografi, ITB.
GEBCO. (2014). Data  Batimetri  Online.  Dipetik Maret 18, 2016, dari https://www.bodc.ac.uk/data/online_delivery/gebco/.
PPN. (2009). Laporan Rencana Pengembangan Pelabuhan Perikanan Nusantara Kejawanan. Cirebon: PPN Kejawanan.
SWAN.  (1993-2015). SWAN  Scientific  and  Technical  Documentation.  Netherlands:  Delft University of Technology.
SWAN. (1993-2015). SWAN User Manual. Netherlands: Delft University of Technology. Triatmodjo, B. (1999). Teknik Pantai. Yogyakarta: Penerbit Beta Offset.