adding test_io3

regtests/unittests/CMakeLists.txt

@@ -42,6 +42,7 @@ copy_test_data_2(ww3_grid.inp ww3_grid.inp)
 unit_test(test_io_points_bin)
 unit_test(test_io)
 unit_test(test_io2)
+unit_test(test_io3)
 unit_test(test_io_restart_bin)
 
 

regtests/unittests/test_io2.F90

@@ -1,7 +1,7 @@
 ! This is a test for model IO for WW3.
 !
 ! Ed Hartnett 10/14/23
-program test_io
+program test_io2
   use w3iopomd
   use w3gdatmd
   use w3wdatmd
@@ -19,7 +19,7 @@ program test_io
   real :: expected_loc_1
   integer :: write_test_file
 
-  print *, 'Testing WW3 netCDF point file code.'
+  print *, 'Testing WW3 netCDF point file code some more.'
 
   ! These are mysterious but have to be called or else the IPASS
   ! variable does not exist and w3iopo() crashes.
@@ -119,5 +119,5 @@ program test_io
   print *, 'OK!'
 
   print *, 'SUCCESS!'
-end program test_io
+end program test_io2
 

regtests/unittests/test_io3.F90

@@ -0,0 +1,109 @@
+! This is a test for model IO for WW3.
+!
+! Ed Hartnett 10/14/23
+program test_io3
+  use w3iopomd
+  use w3gdatmd
+  use w3wdatmd
+  use w3odatmd
+  use w3iogrmd
+  use w3adatmd  
+  implicit none
+
+  integer, target :: i, j, k, l
+  integer :: ndsop, iotest, imod, ndstst, ierr, ndsbul, ndsm
+  integer :: ndstrc, ntrace
+  real :: m2km
+  character*7 expected_ptnme
+  character*6 my_fmt
+  real :: expected_loc_1
+  integer :: write_test_file
+
+  print *, 'Testing WW3 netCDF point file code even more.'
+
+  ! These are mysterious but have to be called or else the IPASS
+  ! variable does not exist and w3iopo() crashes.
+  call w3nmod(1, 6, 6)
+  call w3setg(1, 6, 6)
+  call w3ndat(6, 6)
+  call w3setw(1, 6, 6)
+  call w3nout(6, 6)
+  call w3seto(1, 6, 6)
+
+  ndsm   = 20
+  ndsop  = 20
+  ndsbul = 0
+  ndstrc =  6
+  ntrace = 10
+  imod = 1
+
+  ! Create a point output file needed for this test.
+  print *, 'Creating point output test file for this test...'
+  if (write_test_file() .ne. 0) stop 1
+
+  ! 2.  Read model definition file.
+  CALL W3IOGR('READ', NDSM)
+  WRITE (NDSO,920) GNAME
+920 FORMAT ('  Grid name : ',A/)  
+
+  ! Read the file out_pnt.ww3 from the model/tests/data directory.
+  print *, 'OK!'
+  print *, 'Reading the point output test file for this test...'
+  call w3iopo('READ', ndsop, iotest)
+  if (iotest .ne. 0) stop 10
+  close(ndsop)
+
+  ! ! Make sure we got the values we expected.
+  ! if (nopts .ne. 11) stop 11
+  ! expected_loc_1 = 0.0
+  ! do i = 1, nopts
+  !    ! Check ptnme and ptloc arrays.
+  !    print *, ptnme(i), ptloc(1, i), ptloc(2, i)
+  !    if (i .lt. 10) then
+  !       my_fmt = '(a,i1)'
+  !    else
+  !       my_fmt = '(a,i2)'
+  !    endif
+  !    write(fmt = my_fmt, unit=expected_ptnme) 'Point', i
+  !    if (ptnme(i) .ne. expected_ptnme) stop 20
+  !    print *, expected_loc_1
+  !    if (ptloc(1, i) .ne. expected_loc_1) stop 21
+  !    expected_loc_1 = expected_loc_1 + 5000.0
+  !    if (ptloc(2, i) .ne. 0) stop 22
+  ! end do
+
+  print *, 'OK!'
+  print *, 'initializing some data...'
+  ipass2 = 0
+  do i = 1, nopts
+     do j = 1, nspec
+        spco(j, i) = 0.0
+     end do
+  end do
+
+  print *, 'OK!'
+  print *, 'testing writing the WW3 binary point file in netCDF...'
+
+  ! Write in netCDF.
+  ofiles(2) = 1
+  call w3iopon('WRITE', ndsop, iotest, imod)
+  if (iotest .ne. 0) stop 100
+  print *, 'OK!'
+
+  ! Another timestep in netCDF.
+  call w3iopon('WRITE', ndsop, iotest, imod)
+  if (iotest .ne. 0) stop 100
+  print *, 'OK!'
+
+  print *, 'testing reading the WW3 binary point file in netCDF...'
+  ipass2 = 0
+  call w3iopon('READ', ndsop, iotest)
+  print *, iotest
+  if (iotest .ne. 0) stop 100
+  print *, 'OK!'
+
+  print *, 'OK!'
+
+  print *, 'SUCCESS!'
+end program test_io3
+