diff --git a/model/src/w3iopomd.F90 b/model/src/w3iopomd.F90
index 0c15ea8c6..246ab4b71 100644
--- a/model/src/w3iopomd.F90
+++ b/model/src/w3iopomd.F90
@@ -118,6 +118,96 @@ MODULE W3IOPOMD
   CHARACTER(LEN=10), PARAMETER, PRIVATE :: VEROPT = '2021-04-06'
   CHARACTER(LEN=31), PARAMETER, PRIVATE ::                        &
        IDSTR = 'WAVEWATCH III POINT OUTPUT FILE'
+
+  !> Dimension name for the netCDF point output file, for NOPTS, the
+  !> Number of Output Points.
+  character(*), parameter, private :: DNAME_NOPTS = 'NOPTS'
+  
+  !> Dimension name for the netCDF point output file, for NSPEC.
+  character(*), parameter, private :: DNAME_NSPEC = 'NSPEC'
+  
+  !> Dimension name for the netCDF point output file, for VSIZE. This
+  !> is for the vector size for points, which is 2.
+  character(*), parameter, private :: DNAME_VSIZE = 'VSIZE'
+  
+  !> Dimension name for the netCDF point output file, for
+  !> NAMELEN. This is the length of the PTNME strings, which contains
+  !> the names of the points.
+  character(*), parameter, private :: DNAME_NAMELEN = 'NAMELEN'
+  
+  !> Dimension name for the netCDF point output file, for GRDIDLEN,
+  !> this is the length of the GRDID character array.
+  character(*), parameter, private :: DNAME_GRDIDLEN = 'GRDIDLEN'
+
+  !> Dimension name for the netCDF point output file, for TIME.
+  character(*), parameter, private :: DNAME_TIME = 'TIME'
+
+  !> Variable name for the netCDF point output file, for NK.
+  character(*), parameter, private :: VNAME_NK = 'NK'
+
+  !> Variable name for the netCDF point output file, for MTH.
+  character(*), parameter, private :: VNAME_MTH = 'MTH'
+
+  !> Variable name for the netCDF point output file, for PTLOC.
+  character(*), parameter, private :: VNAME_PTLOC = 'PTLOC'
+
+  !> Variable name for the netCDF point output file, for PTNME.
+  character(*), parameter, private :: VNAME_PTNME = 'PTNME'
+
+  !> Variable name for the netCDF point output file, for IW.
+  character(*), parameter, private :: VNAME_IW = 'IW'
+
+  !> Variable name for the netCDF point output file, for II.
+  character(*), parameter, private :: VNAME_II = 'II'
+
+  !> Variable name for the netCDF point output file, for IL.
+  character(*), parameter, private :: VNAME_IL = 'IL'
+
+  !> Variable name for the netCDF point output file, for DPO.
+  character(*), parameter, private :: VNAME_DPO = 'DPO'
+
+  !> Variable name for the netCDF point output file, for WAO.
+  character(*), parameter, private :: VNAME_WAO = 'WAO'
+
+  !> Variable name for the netCDF point output file, for WDO.
+  character(*), parameter, private :: VNAME_WDO = 'WDO'
+
+  !> Variable name for the netCDF point output file, for TAUAO.
+  character(*), parameter, private :: VNAME_TAUAO = 'TAUAO'
+
+  !> Variable name for the netCDF point output file, for TAIDO.
+  character(*), parameter, private :: VNAME_TAIDO = 'TAIDO'
+
+  !> Variable name for the netCDF point output file, for DAIRO.
+  character(*), parameter, private :: VNAME_DAIRO = 'DAIRO'
+
+  !> Variable name for the netCDF point output file, for ZET_SETO.
+  character(*), parameter, private :: VNAME_ZET_SETO = 'ZET_SETO'
+
+  !> Variable name for the netCDF point output file, for ASO.
+  character(*), parameter, private :: VNAME_ASO = 'ASO'
+
+  !> Variable name for the netCDF point output file, for CAO.
+  character(*), parameter, private :: VNAME_CAO = 'CAO'
+
+  !> Variable name for the netCDF point output file, for CDO.
+  character(*), parameter, private :: VNAME_CDO = 'CDO'
+
+  !> Variable name for the netCDF point output file, for ICEO.
+  character(*), parameter, private :: VNAME_ICEO = 'ICEO'
+
+  !> Variable name for the netCDF point output file, for ICEHO.
+  character(*), parameter, private :: VNAME_ICEHO = 'ICEHO'
+
+  !> Variable name for the netCDF point output file, for ICEFO.
+  character(*), parameter, private :: VNAME_ICEFO = 'ICEFO'
+
+  !> Variable name for the netCDF point output file, for GRDID.
+  character(*), parameter, private :: VNAME_GRDID = 'GRDID'
+
+  !> Variable name for the netCDF point output file, for SPCO.
+  character(*), parameter, private :: VNAME_SPCO = 'SPCO'
+  
   !/
 CONTAINS
   !/ ------------------------------------------------------------------- /
@@ -1024,6 +1114,457 @@ SUBROUTINE W3IOPE ( A )
     !/
   END SUBROUTINE W3IOPE
 
+  !> Read point output in netCDF format.
+  !>
+  !> @param[out] IOTST Test indictor for reading.
+  !> @param[in] IMOD Model number for W3GDAT etc.
+  !> @param[in] filename Name of file to read.
+  !> @param[inout] ncerr Error code, 0 for success, netCDF error code
+  !> otherwise.
+  !>
+  !> @author Edward Hartnett  @date 1-Nov-2023
+  !>
+  SUBROUTINE W3IOPON_READ(IOTST, IMOD, filename, ncerr)
+    use netcdf
+    USE W3GDATMD, ONLY: NTH, NK, NSPEC, FILEXT
+    USE W3ODATMD, ONLY: NDST, NDSE, IPASS => IPASS2, NOPTS, IPTINT, &
+         IL, IW, II, PTLOC, PTIFAC, DPO, WAO, WDO,   &
+         ASO, CAO, CDO, SPCO, PTNME, O2INIT, FNMPRE, &
+         GRDID, ICEO, ICEHO, ICEFO
+#ifdef W3_FLX5
+    USE W3ODATMD, ONLY: TAUAO, TAUDO, DAIRO
+#endif
+#ifdef W3_SETUP
+    USE W3ODATMD, ONLY: ZET_SETO
+#endif
+    IMPLICIT NONE
+
+    INTEGER, INTENT(OUT)          :: IOTST
+    INTEGER, INTENT(IN), OPTIONAL :: IMOD
+    character(*), intent(in) :: filename
+    integer, intent(inout) :: ncerr
+    integer :: fh
+    integer :: d_nopts, d_nspec, d_vsize, d_namelen, d_grdidlen    
+    integer :: d_nopts_len, d_nspec_len, d_vsize_len, d_namelen_len, d_grdidlen_len    
+    integer :: v_idtst, v_vertst, v_nk, v_mth, v_ptloc, v_ptnme
+    integer :: v_iw, v_ii, v_il, v_dpo, v_wao, v_wdo, v_tauao
+    integer :: v_taido, v_dairo, v_zet_seto, v_aso, v_cao, v_cdo, v_iceo
+    integer :: v_iceho, v_icefo, v_grdid, v_spco
+
+    print *,filename
+    IOTST = 0
+    
+    ! Open the netCDF file.
+    ncerr = nf90_open(filename, NF90_NOWRITE, fh)
+    if (ncerr .ne. 0) return
+
+    ! Read the dimension information for NOPTS.
+    ncerr = nf90_inq_dimid(fh, DNAME_NOPTS, d_nopts)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inquire_dimension(fh, d_nopts, len = d_nopts_len)
+    if (ncerr .ne. 0) return
+
+    ! Read the dimension information for NSPEC.
+    ncerr = nf90_inq_dimid(fh, DNAME_NSPEC, d_nspec)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inquire_dimension(fh, d_nspec, len = d_nspec_len)
+    if (ncerr .ne. 0) return
+
+    ! Read the dimension information for VSIZE.
+    ncerr = nf90_inq_dimid(fh, DNAME_VSIZE, d_vsize)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inquire_dimension(fh, d_vsize, len = d_vsize_len)
+    if (ncerr .ne. 0) return
+
+    ! Read the dimension information for NAMELEN.
+    ncerr = nf90_inq_dimid(fh, DNAME_NAMELEN, d_namelen)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inquire_dimension(fh, d_namelen, len = d_namelen_len)
+    if (ncerr .ne. 0) return
+
+    ! Read the dimension information for GRDIDLEN.
+    ncerr = nf90_inq_dimid(fh, DNAME_GRDIDLEN, d_grdidlen)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inquire_dimension(fh, d_grdidlen, len = d_grdidlen_len)
+    if (ncerr .ne. 0) return
+
+    ! Read scalar variables.
+    ncerr = nf90_inq_varid(fh, VNAME_NK, v_nk)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_nk, NK)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_MTH, v_mth)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_mth, NTH)
+    if (ncerr .ne. 0) return
+
+    ! Read vars with nopts as a dimension.
+    ncerr = nf90_inq_varid(fh, VNAME_PTLOC, v_ptloc)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_ptloc, PTLOC)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_PTNME, v_ptnme)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_ptnme, PTNME)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_IW, v_iw)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_iw, IW)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_II, v_ii)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_ii, II)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_IL, v_il)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_il, IL)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_DPO, v_dpo)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_dpo, DPO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_WAO, v_wao)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_wao, WAO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_WDO, v_wdo)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_wdo, WDO)
+    if (ncerr .ne. 0) return
+#ifdef W3_FLX5
+    ncerr = nf90_inq_varid(fh, VNAME_TAUAO, v_tauao)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_tauao, TAUAO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_TAIDO, v_taido)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_taido, TAIDO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_DAIRO, v_dairo)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_dairo, DAIRO)
+    if (ncerr .ne. 0) return
+#endif    
+#ifdef W3_SETUP
+    ncerr = nf90_inq_varid(fh, ZET_SETO, v_zet_seto)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_zet_seto, ZET_SETO)
+    if (ncerr .ne. 0) return
+#endif    
+    ncerr = nf90_inq_varid(fh, VNAME_ASO, v_aso)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_aso, ASO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_CAO, v_cao)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_cao, CAO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_CDO, v_cdo)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_cdo, CDO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_ICEO, v_iceo)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_iceo, ICEO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_ICEHO, v_iceho)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_iceho, ICEHO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_ICEFO, v_icefo)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_icefo, ICEFO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_GRDID, v_grdid)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_grdid, GRDID)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_inq_varid(fh, VNAME_SPCO, v_spco)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_get_var(fh, v_spco, SPCO)
+    if (ncerr .ne. 0) return
+    
+    ! Close the file.
+    ncerr = nf90_close(fh)
+    if (ncerr .ne. 0) return
+    
+  END SUBROUTINE W3IOPON_READ
+
+  !/ ------------------------------------------------------------------- /
+  !>
+  !> @brief Write point output in netCDF format.
+  !>
+  !> @param[in] IMOD Model number for W3GDAT etc.
+  !> @param[in] filename Name of file to write.
+  !> @param[in] timestep_only Will be 0 if whole model run should be
+  !> written, 1 if only one timestep should be written.
+  !> @param[inout] ncerr Error code, 0 for success, netCDF error code
+  !> otherwise.
+  !>
+  !> @author Edward Hartnett  @date 1-Nov-2023
+  !>
+  SUBROUTINE W3IOPON_WRITE(timestep_only, IMOD, filename, ncerr)
+    use netcdf
+    USE W3GDATMD, ONLY: NTH, NK, NSPEC, FILEXT
+    USE W3ODATMD, ONLY: NDST, NDSE, IPASS => IPASS2, NOPTS, IPTINT, &
+         IL, IW, II, PTLOC, PTIFAC, DPO, WAO, WDO,   &
+         ASO, CAO, CDO, SPCO, PTNME, O2INIT, FNMPRE, &
+         GRDID, ICEO, ICEHO, ICEFO
+#ifdef W3_FLX5
+    USE W3ODATMD, ONLY: TAUAO, TAUDO, DAIRO
+#endif
+#ifdef W3_SETUP
+    USE W3ODATMD, ONLY: ZET_SETO
+#endif
+    
+    IMPLICIT NONE
+    integer, intent(in) :: timestep_only ! 1 if only timestep should be written.
+    INTEGER, INTENT(IN) :: IMOD
+    character(*), intent(in) :: filename
+    integer, intent(inout) :: ncerr
+    integer :: fh, ndim, nvar, fmt
+    integer :: d_nopts, d_nspec, d_vsize, d_namelen, d_grdidlen, d_time
+    integer :: v_idtst, v_vertst, v_nk, v_mth, v_ptloc, v_ptnme
+    integer :: v_iw, v_ii, v_il, v_dpo, v_wao, v_wdo, v_tauao
+    integer :: v_taido, v_dairo, v_zet_seto, v_aso, v_cao, v_cdo, v_iceo
+    integer :: v_iceho, v_icefo, v_grdid, v_spco
+    character (len = *), parameter :: FILE_NAME = "f90tst_nc4.nc"
+    CHARACTER(LEN=31), PARAMETER :: IDSTR = 'WAVEWATCH III POINT OUTPUT FILE'
+    CHARACTER(LEN=10), PARAMETER :: VEROPT = '2021-04-06'
+
+    print *, 'WRITE ',filename, len(filename)
+
+    ! Create the netCDF file.
+    ncerr = nf90_create(filename, NF90_NETCDF4, fh)
+    if (ncerr .ne. 0) return
+
+    ! Define dimensions.
+    ncerr = nf90_def_dim(fh, DNAME_NOPTS, NOPTS, d_nopts)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_dim(fh, DNAME_NSPEC, NSPEC, d_nspec)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_dim(fh, DNAME_VSIZE, 2, d_vsize)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_dim(fh, DNAME_NAMELEN, 40, d_namelen)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_dim(fh, DNAME_GRDIDLEN, 13, d_grdidlen)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_dim(fh, DNAME_TIME, NF90_UNLIMITED, d_time)
+    if (ncerr .ne. 0) return
+
+    ! Define global attributes.
+    ncerr = nf90_put_att(fh, NF90_GLOBAL, 'title', IDSTR)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_att(fh, NF90_GLOBAL, 'version', VEROPT)
+    if (ncerr .ne. 0) return
+    
+    ! Define scalar variables.
+    ncerr = nf90_def_var(fh, VNAME_NK, NF90_INT, v_nk)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_MTH, NF90_INT, v_mth)
+    if (ncerr .ne. 0) return
+
+    ! Define vars with nopts as a dimension.
+    ncerr = nf90_def_var(fh, VNAME_PTLOC, NF90_INT, (/d_vsize, d_nopts/), v_ptloc)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_PTNME, NF90_CHAR, (/d_namelen, d_nopts/), v_ptnme)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_IW, NF90_INT, (/d_nopts, d_time/), v_iw)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_II, NF90_INT, (/d_nopts, d_time/), v_ii)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_IL, NF90_INT, (/d_nopts, d_time/), v_il)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_DPO, NF90_INT, (/d_nopts, d_time/), v_dpo)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_WAO, NF90_INT, (/d_nopts, d_time/), v_wao)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_WDO, NF90_INT, (/d_nopts, d_time/), v_wdo)
+    if (ncerr .ne. 0) return
+#ifdef W3_FLX5
+    ncerr = nf90_def_var(fh, VNAME_TAUAO, NF90_INT, (/d_nopts, d_time/), v_tauao)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_TAIDO, NF90_INT, (/d_nopts, d_time/), v_taido)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_DAIRO, NF90_INT, (/d_nopts, d_time/), v_dairo)
+    if (ncerr .ne. 0) return
+#endif    
+#ifdef W3_SETUP
+    ncerr = nf90_def_var(fh, VNAME_ZET_SETO, NF90_INT, (/d_nopts, d_time/), v_zet_seto)
+    if (ncerr .ne. 0) return
+#endif    
+    ncerr = nf90_def_var(fh, VNAME_ASO, NF90_INT, (/d_nopts, d_time/), v_aso)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_CAO, NF90_INT, (/d_nopts, d_time/), v_cao)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_CDO, NF90_INT, (/d_nopts, d_time/), v_cdo)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_ICEO, NF90_INT, (/d_nopts, d_time/), v_iceo)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_ICEHO, NF90_INT, (/d_nopts, d_time/), v_iceho)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_ICEFO, NF90_INT, (/d_nopts, d_time/), v_icefo)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_GRDID, NF90_CHAR, (/d_grdidlen, d_nopts, d_time/), v_grdid)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_def_var(fh, VNAME_SPCO, NF90_INT, (/d_nspec, d_nopts, d_time/), v_spco)
+    if (ncerr .ne. 0) return
+
+    ! Write the scalar data.
+    ncerr = nf90_put_var(fh, v_nk, NK)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_mth, NTH)
+    if (ncerr .ne. 0) return
+
+    ! Write the data with NOPTS as a dimension.
+    ncerr = nf90_put_var(fh, v_ptloc, PTLOC)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_ptnme, PTNME)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_iw, IW)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_ii, II)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_il, IL)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_dpo, DPO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_wao, WAO)
+    if (ncerr .ne. 0) return
+#ifdef W3_FLX5
+    ncerr = nf90_put_var(fh, v_tauao, TAUAO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_taido, TAIDO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_dairo, DAIRO)
+    if (ncerr .ne. 0) return
+#endif
+#ifdef W3_SETUP
+    ncerr = nf90_put_var(fh, v_zet_seto, ZET_SETO)
+    if (ncerr .ne. 0) return
+#endif
+    ncerr = nf90_put_var(fh, v_aso, ASO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_cao, CAO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_iceo, ICEO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_iceho, ICEHO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_icefo, ICEFO)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_grdid, GRDID)
+    if (ncerr .ne. 0) return
+    ncerr = nf90_put_var(fh, v_spco, SPCO)
+    if (ncerr .ne. 0) return
+
+    ! Close the file.
+    ncerr = nf90_close(fh)
+    if (ncerr .ne. 0) return
+
+  END SUBROUTINE W3IOPON_WRITE
+
+  !> Read/write point output in netCDF format.
+  !>
+  !> @param[in] INXOUT String indicating read/write. Must be 'READ' or
+  !> 'WRITE'.
+  !> @param[in] NDSOP File unit number.
+  !> @param[out] IOTST Error code:
+  !> - 0 No error.
+  !> - -1 Unexpected end of file when reading.
+  !> @param[in] IMOD Model number for W3GDAT etc.
+#ifdef W3_ASCII
+  !> @param[in] NDSOA File unit number for ASCII output.
+#endif
+  !>
+  !> @author Edward Hartnett  @date 1-Nov-2023
+  SUBROUTINE W3IOPON ( INXOUT, NDSOP, IOTST, IMOD &
+#ifdef W3_ASCII
+      ,NDSOA &
+#endif
+      )
+    USE W3GDATMD, ONLY: W3SETG
+    USE W3WDATMD, ONLY: W3SETW
+    USE W3ODATMD, ONLY: W3SETO, W3DMO2
+    !/
+    USE W3GDATMD, ONLY: NTH, NK, NSPEC, FILEXT
+    USE W3WDATMD, ONLY: TIME
+    USE W3ODATMD, ONLY: NDST, NDSE, IPASS => IPASS2, NOPTS, IPTINT, &
+         IW, PTLOC, PTIFAC, WAO, WDO,   &
+         SPCO, PTNME, O2INIT, FNMPRE
+#ifdef W3_FLX5
+    USE W3ODATMD, ONLY: TAUAO, TAUDO, DAIRO
+#endif
+    USE W3ODATMD, ONLY :  OFILES
+    !/
+#ifdef W3_SETUP
+    USE W3ODATMD, ONLY: ZET_SETO
+#endif
+    !/
+    USE W3SERVMD, ONLY: EXTCDE
+#ifdef W3_S
+    USE W3SERVMD, ONLY: STRACE
+#endif
+    use netcdf
+    IMPLICIT NONE
+
+    CHARACTER, INTENT(IN)         :: INXOUT*(*)
+    INTEGER, INTENT(IN)           :: NDSOP
+    INTEGER, INTENT(OUT)          :: IOTST
+    INTEGER, INTENT(IN), OPTIONAL :: IMOD
+#ifdef W3_ASCII
+    INTEGER, INTENT(IN), OPTIONAL :: NDSOA
+#endif
+
+    INTEGER :: IGRD, MK, MTH
+    character(len = 124) :: filename
+    integer :: ncerr
+
+    ! Optimistically assume success.
+    IOTST = 0
+
+    ! Has a model number been specified?
+    IF (PRESENT(IMOD)) THEN
+      IGRD = IMOD
+    ELSE
+      IGRD = 1
+    END IF
+
+    CALL W3SETO(IGRD, NDSE, NDST)
+    CALL W3SETG(IGRD, NDSE, NDST)
+    CALL W3SETW(IGRD, NDSE, NDST)
+
+    ! INXOUT must be 'READ' or 'WRITE'.
+    IF (INXOUT .NE. 'READ' .AND. INXOUT .NE. 'WRITE') THEN
+      WRITE (NDSE, 900) INXOUT
+      CALL EXTCDE(1)
+   END IF
+
+   ! Determine filename.
+   filename = ''
+   filename = transfer(FNMPRE(:LEN_TRIM(FNMPRE))//'out_pnt_nc.'//FILEXT(:LEN_TRIM(FILEXT)), filename)
+   !print *, filename
+   
+   ! Do a read or a write of the point file.
+   IF (INXOUT .EQ. 'READ') THEN
+      CALL W3IOPON_READ(IOTST, IMOD, filename, ncerr)
+   ELSE
+      CALL W3IOPON_WRITE(OFILES(2), IMOD, filename, ncerr)
+   ENDIF
+   if (ncerr .ne. 0) then
+      print *, nf90_strerror(ncerr)
+      CALL EXTCDE(21)      
+   endif
+   
+    !/
+    !/ End of W3IOPON ----------------------------------------------------- /
+    !/
+
+900 FORMAT (/' *** WAVEWATCH III ERROR IN W3IOPO :'/                &
+         '     ILEGAL INXOUT VALUE: ',A/)
+  END SUBROUTINE W3IOPON
+ 
+  !/ ------------------------------------------------------------------- /
   !> Read or write point output.
   !>
   !> This subroutine can either read or write the point output file,
@@ -1495,7 +2036,7 @@ SUBROUTINE W3IOPO ( INXOUT, NDSOP, IOTST, IMOD &
 #endif
              'ASO(I), CAO(I), CDO(I), ICEO(I), ICEHO(I):',        &
              ASO(I), CAO(I), CDO(I), ICEO(I), ICEHO(I),        &
-             'ICEFO(I), GRDID(I), (SPCO(J,I),J=1,NSPEC):',       & 
+             'ICEFO(I), GRDID(I), (SPCO(J,I),J=1,NSPEC):',       &
              ICEFO(I), GRDID(I), (SPCO(J,I),J=1,NSPEC)
 #endif
       ELSE
diff --git a/regtests/unittests/test_io.F90 b/regtests/unittests/test_io.F90
index 86e333679..ce89d5df6 100644
--- a/regtests/unittests/test_io.F90
+++ b/regtests/unittests/test_io.F90
@@ -97,9 +97,9 @@ program test_io
   print *, 'testing writing the WW3 binary point file in netCDF...'
 
   ! Write in netCDF.
-  ! call w3iopon('WRITE', ndsop, iotest)
-  ! if (iotest .ne. 0) stop 100
-  ! print *, 'OK!'
+  call w3iopon('WRITE', ndsop, iotest)
+  if (iotest .ne. 0) stop 100
+  print *, 'OK!'
   
   ! print *, 'testing reading the WW3 binary point file in netCDF...'
   ! call w3iopon('READ', ndsop, iotest)