subroutine plotcoast(region,ifnum)

c - See DRU-13, p. 4-6, 5/20,2019

c - Cannibalized from NADCON 5.01 to VERTCON 3.0

c - 2018 03 26 : To put in new region names for VERTCON3, specifically
c - Remove: prvi , guamcnmi 
c - Add:  pr, vi, guam, cnmi

c - 2016 01 07
c     Forced the Alaska region to plot the
c     islands of St. George, St. Matthew and St. Paul
c     (St. Lawrence is already plotted), as well
c     as 35 missing Aleutian Islands

c - 2015 09 23:  Added four new regions:
c      St. Lawrence Island, Alaska
c      St. Matthew Island, Alaska
c      St. George Island, Alaska
c      St. Paul Island, Alaska

      character*10 region

c - Write GMT-based commands to create a shoreline
c - based on region.
c - Use GMT-default coastline for:
c      conus
c      alaska
c - Use Dru's custome coastline for:
c      hawaii
c      prvi
c      as
c      guamcnmi
c      stlawrence
c      stmatthew 
c      stgeorge  
c      stpaul    


      if    (trim(region).eq.'conus')then
        write(ifnum,100) 
      elseif(trim(region).eq.'alaska')then
        write(ifnum,109) 
      elseif(trim(region).eq.'hawaii')then
        write(ifnum,101) 
c     elseif(trim(region).eq.'prvi')then
      elseif(trim(region).eq.'pr' .or. trim(region).eq.'vi')then
        write(ifnum,102) 
c     elseif(trim(region).eq.'guamcnmi')then
      elseif(trim(region).eq.'guam' .or. trim(region).eq.'cnmi')then
        write(ifnum,103) 
      elseif(trim(region).eq.'as')then
        write(ifnum,104) 
      elseif(trim(region).eq.'stlawrence')then
        write(ifnum,105) 
      elseif(trim(region).eq.'stmatthew')then
        write(ifnum,106) 
      elseif(trim(region).eq.'stgeorge')then
        write(ifnum,107) 
      elseif(trim(region).eq.'stpaul')then
        write(ifnum,108) 
     
      else
        write(6,200)region
        stop 10001
      endif
      return

  200 format('FATAL in plotcoast.f: Unknown region: ',a10)

c - CONUS 
  100 format(
     *'pscoast -Df -R -JM -W0.25p -N1 -N2 -A1200 -O >>  plot.ps')

c - HAWAII
  101 format(
     *'psxy /home/dru/Boundaries/HI_Hawaii.gmt -R -JM -O -K',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/HI_Kahoolawe.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/HI_Kauai.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/HI_Lanai.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/HI_Maui.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/HI_Molokai.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/HI_Nihau.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/HI_Oahu.gmt -R -JM -O ',
     *' >> plot.ps')

c - PRVI
  102 format(
     *'psxy /home/dru/Boundaries/VQ_StCroix.gmt -R -JM -O -K',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/VQ_StJohn.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/VQ_StThomas.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/PR_Culebra.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/PR_Desecheo.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/PR_Mona.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/PR_Viequez.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/PR_PuertoRico.gmt -R -JM -O ',
     *' >> plot.ps')
 
c - GUAM/CNMI
  103 format(
     *'psxy /home/dru/Boundaries/GQ_Guam.gmt -R -JM -O -K',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_Rota.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_Saipan.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_Tinian.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_Anatahan.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_Sarigan.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_Guguan.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_Alamagan.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_Pagan.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_Agrihan.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_Asuncion.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_MaugWest.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_MaugEast.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_MaugNorth.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/CQ_Pajaros.gmt -R -JM -O',
     *' >> plot.ps')

c - AS  
  104 format(
     *'psxy /home/dru/Boundaries/AS_Aunuu.gmt -R -JM -O -K',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AS_Nuutele.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AS_Ofu.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AS_Olosega.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AS_Rose.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AS_Swains.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AS_Tau.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AS_Tutila.gmt -R -JM -O ',
     *' >> plot.ps')

c - St. Lawrence Island, Alaska
  105 format(
     *'psxy /home/dru/Boundaries/AK_StLawrence.gmt -R -JM -O -K',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_NorthPunuk.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_MiddlePunuk.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_SouthPunuk.gmt -R -JM -O ',
     *' >> plot.ps')

c - St. Matthew Island, Alaska
  106 format(
     *'psxy /home/dru/Boundaries/AK_StMatthew.gmt -R -JM -O -K',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Hall.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Pinnacle.gmt -R -JM -O ',
     *' >> plot.ps')

c - St. George Island, Alaska
  107 format(
     *'psxy /home/dru/Boundaries/AK_StGeorge.gmt -R -JM -O',
     *' >> plot.ps')

c - St. Paul Island, Alaska
  108 format(
     *'psxy /home/dru/Boundaries/AK_StPaul.gmt -R -JM -O -K',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Otter.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Walrus.gmt -R -JM -O ',
     *' >> plot.ps')

c - ALASKA
  109 format(
     *'pscoast -Df -R -JM -W0.25p -N1 -N2 -A1200 -O -K ',
     *' >>  plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_NorthPunuk.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_MiddlePunuk.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_SouthPunuk.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_StMatthew.gmt -R -JM -O -K',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Hall.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Pinnacle.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_StGeorge.gmt -R -JM -O -K',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_StPaul.gmt -R -JM -O -K',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Otter.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Walrus.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Adak.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Agattu.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Amchitka.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Amlia.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Amukta.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Andreanof.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Atka.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Attu.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Buldir.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Carlisle.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Chagulak.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Chugul.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Fenimore.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_FourMountains.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_GreatSitkin.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Herbert.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Igitkin.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Ikiginak.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Kagalaska.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Kagami.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Kanaga.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Kasatochi.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Kiska.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Koniuji.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Oglodak.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Seguam.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Semichi.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Semisopochnoi.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Shemya.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Tagalak.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Tanaga.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Ulak.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Uliaga.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Umak.gmt -R -JM -O -K ',
     *' >> plot.ps',/,
     *'psxy /home/dru/Boundaries/AK_Yunaska.gmt -R -JM -O ',
     *' >> plot.ps')

      end