w3nmlmultimd.F90

Go to the documentation of this file.

#include "w3macros.h"
!/ ------------------------------------------------------------------- /
MODULE w3nmlmultimd
  !/
  !/                  +-----------------------------------+
  !/                  | WAVEWATCH III           NOAA/NCEP |
  !/                  |           M. Accensi              |
  !/                  |                                   |
  !/                  |                        FORTRAN 90 |
  !/                  | Last update :         22-Mar-2021 |
  !/                  +-----------------------------------+
  !/
  !/    For updates see subroutines.
  !/
  !  1. Purpose :
  !
  !     Manages namelists from configuration file ww3_multi.nml for ww3_multi program
  !
  !/ ------------------------------------------------------------------- /
 
  ! module defaults
  IMPLICIT NONE
 
  PUBLIC
 
  ! domain definition structure
  TYPE nml_domain_t
    INTEGER                     :: nrinp
    INTEGER                     :: nrgrd
    LOGICAL                     :: unipts
    INTEGER                     :: iostyp
    LOGICAL                     :: upproc
    LOGICAL                     :: pshare
    LOGICAL                     :: flghg1
    LOGICAL                     :: flghg2
    CHARACTER(15)               :: start
    CHARACTER(15)               :: stop
  END TYPE nml_domain_t
 
 
 
  ! model grid data structure
  TYPE nml_model_forcing_t
    CHARACTER(13)               :: water_levels
    CHARACTER(13)               :: currents
    CHARACTER(13)               :: winds
    CHARACTER(13)               :: atm_momentum
    CHARACTER(13)               :: air_density
    CHARACTER(13)               :: ice_conc
    CHARACTER(13)               :: ice_param1
    CHARACTER(13)               :: ice_param2
    CHARACTER(13)               :: ice_param3
    CHARACTER(13)               :: ice_param4
    CHARACTER(13)               :: ice_param5
    CHARACTER(13)               :: mud_density
    CHARACTER(13)               :: mud_thickness
    CHARACTER(13)               :: mud_viscosity
  END TYPE nml_model_forcing_t
  !
  TYPE nml_model_assim_t
    CHARACTER(13)               :: mean
    CHARACTER(13)               :: spec1d
    CHARACTER(13)               :: spec2d
  END TYPE nml_model_assim_t
  !
  TYPE nml_model_resource_t
    INTEGER                     :: rank_id
    INTEGER                     :: group_id
    REAL(4)                     :: comm_frac(2)
    LOGICAL                     :: bound_flag
  END TYPE nml_model_resource_t
  !
  TYPE nml_model_grid_t
    CHARACTER(13)               :: name
    TYPE(nml_model_forcing_t)   :: forcing
    TYPE(nml_model_assim_t)     :: assim
    TYPE(nml_model_resource_t)  :: resource
  END TYPE nml_model_grid_t
 
 
 
  ! input grid data structure
  TYPE nml_input_forcing_t
    LOGICAL                     :: water_levels
    LOGICAL                     :: currents
    LOGICAL                     :: winds
    LOGICAL                     :: atm_momentum
    LOGICAL                     :: air_density
    LOGICAL                     :: ice_conc
    LOGICAL                     :: ice_param1
    LOGICAL                     :: ice_param2
    LOGICAL                     :: ice_param3
    LOGICAL                     :: ice_param4
    LOGICAL                     :: ice_param5
    LOGICAL                     :: mud_density
    LOGICAL                     :: mud_thickness
    LOGICAL                     :: mud_viscosity
  END TYPE nml_input_forcing_t
  !
  TYPE nml_input_assim_t
    LOGICAL                     :: mean
    LOGICAL                     :: spec1d
    LOGICAL                     :: spec2d
  END TYPE nml_input_assim_t
  !
  TYPE nml_input_grid_t
    CHARACTER(13)               :: name
    TYPE(nml_input_forcing_t)   :: forcing
    TYPE(nml_input_assim_t)     :: assim
  END TYPE nml_input_grid_t
 
 
 
  ! output type structure
  TYPE nml_field_t
    CHARACTER(1024)             :: list
  END TYPE nml_field_t
  !
  TYPE nml_point_t
    CHARACTER(13)               :: name
    CHARACTER(64)               :: file
  END TYPE nml_point_t
  !
  TYPE nml_track_t
    LOGICAL                     :: format
  END TYPE nml_track_t
  !
  TYPE nml_partition_t
    INTEGER                     :: x0
    INTEGER                     :: xn
    INTEGER                     :: nx
    INTEGER                     :: y0
    INTEGER                     :: yn
    INTEGER                     :: ny
    LOGICAL                     :: format
  END TYPE nml_partition_t
  !
#ifdef W3_COU
  TYPE nml_coupling_t
    CHARACTER(1024)             :: sent
    CHARACTER(1024)             :: received
    LOGICAL                     :: couplet0
  END TYPE nml_coupling_t
#endif
  !
  TYPE nml_restart_t
    CHARACTER(1024)             :: extra
  END TYPE nml_restart_t
  !
  TYPE nml_output_type_t
    TYPE(nml_point_t)           :: point
    TYPE(nml_field_t)           :: field
    TYPE(nml_track_t)           :: track
    TYPE(nml_partition_t)       :: partition
#ifdef W3_COU
    TYPE(nml_coupling_t)       :: coupling
#endif
    TYPE(nml_restart_t)         :: restart
  END TYPE nml_output_type_t
 
 
 
  ! output date structure
  TYPE nml_output_time_t
    CHARACTER(15)               :: start
    CHARACTER(15)               :: stride
    CHARACTER(15)               :: stop
    CHARACTER(1)                :: outffile
    !
  END TYPE nml_output_time_t
  !
  TYPE nml_output_date_t
    TYPE(nml_output_time_t)     :: field
    TYPE(nml_output_time_t)     :: point
    TYPE(nml_output_time_t)     :: track
    TYPE(nml_output_time_t)     :: restart
    TYPE(nml_output_time_t)     :: restart2
    TYPE(nml_output_time_t)     :: boundary
    TYPE(nml_output_time_t)     :: partition
#ifdef W3_COU
    TYPE(nml_output_time_t)     :: coupling
#endif
  END TYPE nml_output_date_t
 
 
 
  ! homogeneous input structure
  TYPE nml_homog_count_t
    INTEGER                     :: n_mov
    INTEGER                     :: n_tot
  END TYPE nml_homog_count_t
  !
  TYPE nml_homog_input_t
    CHARACTER(15)               :: name
    CHARACTER(15)               :: date
    REAL                        :: value1
    REAL                        :: value2
    REAL                        :: value3
  END TYPE nml_homog_input_t
 
 
  ! miscellaneous
  CHARACTER(256)                :: msg
  INTEGER                       :: ndsn
 
 
 
CONTAINS
  !/ ------------------------------------------------------------------- /
  SUBROUTINE w3nmlmultidef (MPI_COMM, NDSI, INFILE, NML_DOMAIN, IERR)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |           F. Ardhuin              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         15-May-2018 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
    !
    !  1. Purpose :
    !
    !     Reads the domain definition namelist to define the number of
    !     model and forcing grids
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      MPI_COMM   Int.  Public   Communicator used in the wave MODEL.
    !      NDSI       Int.
    !      INFILE     Char.
    !      NML_DOMAIN type.
    !      IERR       Int.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !      READ_DOMAIN_NML
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      WMINITNML Subr.   N/A    Model configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
    USE wmmdatmd, ONLY: mdse, improc, nmplog
#ifdef W3_MPI
    USE wmmdatmd, ONLY: mpi_comm_mwave
#endif
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    INTEGER, INTENT(IN)                       :: MPI_COMM, NDSI
    CHARACTER*(*), INTENT(IN)                 :: INFILE
    TYPE(nml_domain_t), INTENT(OUT)           :: NML_DOMAIN
    INTEGER, INTENT(OUT)                      :: IERR
 
    ! locals
#ifdef W3_MPI
    INTEGER                                 :: IERR_MPI
#endif
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
    ierr = 0
#ifdef W3_S
    CALL strace (ient, 'W3NMLMULTIDEF')
#endif
 
#ifdef W3_MPI
    mpi_comm_mwave = mpi_comm
    CALL mpi_comm_rank ( mpi_comm_mwave, improc, ierr_mpi )
    improc = improc + 1
#endif
 
    ! open namelist log file
    IF ( nmplog .EQ. improc ) THEN
      ndsn = 3
      OPEN (ndsn, file=trim(infile)//'.log', form='formatted', iostat=ierr)
      IF (ierr.NE.0) THEN
        WRITE (mdse,'(A)') 'ERROR: open full nml file '//trim(infile)//'.log failed'
        RETURN
      END IF
    END IF
 
    ! open input file
    OPEN (ndsi, file=trim(infile), form='formatted', status='old', iostat=ierr)
    IF (ierr.NE.0) THEN
      WRITE (mdse,'(A)') 'ERROR: open input file '//trim(infile)//' failed'
      RETURN
    END IF
 
    ! read domain def namelist
    CALL read_domain_nml (ndsi, nml_domain)
    IF ( nmplog .EQ. improc ) CALL report_domain_nml (nml_domain)
 
    ! close namelist files
    CLOSE (ndsi)
    IF ( nmplog .EQ. improc ) CLOSE (ndsn)
 
  END SUBROUTINE w3nmlmultidef
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE w3nmlmulticonf (MPI_COMM, NDSI, INFILE, NML_DOMAIN,       &
       NML_INPUT_GRID, NML_MODEL_GRID,           &
       NML_OUTPUT_TYPE, NML_OUTPUT_DATE,         &
       NML_HOMOG_COUNT, NML_HOMOG_INPUT, IERR)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |           F. Ardhuin              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         15-May-2018 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
 
    !  1. Purpose :
    !
    !     Reads all the namelist to define the multi grid
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      MPI_COMM  Int.  Public   Communicator used in the wave MODEL.
    !      NDSI
    !      INFILE
    !      NML_DOMAIN
    !      NML_INPUT_GRID
    !      NML_MODEL_GRID
    !      NML_OUTPUT_TYPE
    !      NML_OUTPUT_DATE
    !      NML_HOMOG_COUNT
    !      NML_HOMOG_INPUT
    !      IERR
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !      READ_INPUT_GRID_NML
    !      REPORT_INPUT_GRID_NML
    !      READ_MODEL_GRID_NML
    !      REPORT_MODEL_GRID_NML
    !      READ_OUTPUT_TYPE_NML
    !      REPORT_OUTPUT_TYPE_NML
    !      READ_OUTPUT_DATE_NML
    !      REPORT_OUTPUT_DATE_NML
    !      READ_HOMOGENEOUS_NML
    !      REPORT_HOMOGENEOUS_NML
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      WMINITNML Subr.   N/A    Model configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
 
    USE wmmdatmd, ONLY: mdse, improc, nmplog
#ifdef W3_MPI
    USE wmmdatmd, ONLY: mpi_comm_mwave
#endif
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    INTEGER, INTENT(IN)                       :: MPI_COMM, NDSI
    CHARACTER*(*), INTENT(IN)                 :: INFILE
    TYPE(nml_domain_t), INTENT(INOUT)         :: NML_DOMAIN
    TYPE(nml_input_grid_t), INTENT(INOUT)     :: NML_INPUT_GRID(:)
    TYPE(nml_model_grid_t), INTENT(INOUT)     :: NML_MODEL_GRID(:)
    TYPE(nml_output_type_t), INTENT(INOUT)    :: NML_OUTPUT_TYPE(:)
    TYPE(nml_output_date_t), INTENT(INOUT)    :: NML_OUTPUT_DATE(:)
    TYPE(nml_homog_count_t), INTENT(INOUT)    :: NML_HOMOG_COUNT
    TYPE(nml_homog_input_t), ALLOCATABLE, INTENT(INOUT)    :: NML_HOMOG_INPUT(:)
    INTEGER, INTENT(OUT)                      :: IERR
 
    ! locals
#ifdef W3_MPI
    INTEGER                                 :: IERR_MPI
#endif
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
    ierr = 0
#ifdef W3_S
    CALL strace (ient, 'W3NMLMULTICONF')
#endif
 
#ifdef W3_MPI
    mpi_comm_mwave = mpi_comm
    CALL mpi_comm_rank ( mpi_comm_mwave, improc, ierr_mpi )
    improc = improc + 1
#endif
 
    ! open namelist log file
    IF ( nmplog .EQ. improc ) THEN
      ndsn = 3
      OPEN (ndsn, file=trim(infile)//'.log', form='formatted', status='old', iostat=ierr)
      IF (ierr.NE.0) THEN
        WRITE (mdse,'(A)') 'ERROR: open full nml file '//trim(infile)//'.log failed'
        RETURN
      END IF
    END IF
 
    ! open input file
    OPEN (ndsi, file=trim(infile), form='formatted', status='old', iostat=ierr)
    IF (ierr.NE.0) THEN
      WRITE (mdse,'(A)') 'ERROR: open INPUT file '//trim(infile)//' failed'
      RETURN
    END IF
 
    ! read input grid namelist
    CALL read_input_grid_nml (ndsi, nml_domain%NRINP, nml_input_grid)
    IF ( nmplog .EQ. improc ) CALL report_input_grid_nml (nml_domain%NRINP, nml_input_grid)
 
    ! read model grid namelist
    CALL read_model_grid_nml (ndsi, nml_domain%NRGRD, nml_model_grid)
    IF ( nmplog .EQ. improc ) CALL report_model_grid_nml (nml_domain%NRGRD, nml_model_grid)
 
    ! read output type namelist
    CALL read_output_type_nml (ndsi, nml_domain%NRGRD, nml_output_type)
    IF ( nmplog .EQ. improc ) CALL report_output_type_nml (nml_domain%NRGRD, nml_output_type)
 
    ! read output date namelist
    CALL read_output_date_nml (ndsi, nml_domain%NRGRD, nml_output_date)
    IF ( nmplog .EQ. improc ) CALL report_output_date_nml (nml_domain%NRGRD,nml_output_date)
 
    ! read homogeneous namelist
    CALL read_homogeneous_nml (ndsi, nml_homog_count, nml_homog_input)
    IF ( nmplog .EQ. improc ) CALL report_homogeneous_nml (nml_homog_count, nml_homog_input)
 
 
    ! close namelist files
    CLOSE (ndsi)
    IF ( nmplog .EQ. improc ) CLOSE (ndsn)
 
  END SUBROUTINE w3nmlmulticonf
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE read_domain_nml (NDSI, NML_DOMAIN)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |           F. Ardhuin              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         15-May-2018 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
 
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NDSI            Int.
    !      NML_DOMAIN      Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLMULTIDEF Subr.   N/A    Namelist configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
 
    USE wmmdatmd, ONLY: mdse
    USE w3servmd, ONLY: extcde
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    INTEGER, INTENT(IN)               :: NDSI
    TYPE(nml_domain_t), INTENT(OUT)   :: NML_DOMAIN
 
    ! locals
    INTEGER                   :: IERR
    TYPE(nml_domain_t) :: DOMAIN
    namelist /domain_nml/ domain
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
    ierr = 0
#ifdef W3_S
    CALL strace (ient, 'READ_DOMAIN_NML')
#endif
 
    ! set default values for model definition data
    domain%NRINP  = 0
    domain%NRGRD  = 1
    domain%UNIPTS = .false.
    domain%IOSTYP = 1
    domain%UPPROC = .false.
    domain%PSHARE = .false.
    domain%FLGHG1 = .false.
    domain%FLGHG2 = .false.
    domain%START  = '19680606 000000'
    domain%STOP   = '19680607 000000'
 
    ! read model def namelist
    rewind(ndsi)
    READ (ndsi, nml=domain_nml, iostat=ierr, iomsg=msg)
    IF (ierr.GT.0) THEN
      WRITE (mdse,'(A,/A)') &
           'ERROR: READ_DOMAIN_NML: namelist read error', &
           'ERROR: '//trim(msg)
      CALL extcde (1)
    END IF
 
    ! set/check return values
    IF (domain%NRINP.LT.0) THEN
      WRITE (mdse,'(A,I3)') 'ERROR: BAD NRINP INPUT: ',domain%NRINP
      CALL extcde (2)
    END IF
    IF (domain%NRGRD.LE.0) THEN
      WRITE (mdse,'(A,I3)') 'ERROR: BAD NRGRD INPUT: ',domain%NRGRD
      CALL extcde (3)
    END IF
    IF (domain%IOSTYP.LT.0.or.domain%IOSTYP.GT.3) THEN
      WRITE (mdse,'(A,I3)') 'ERROR: BAD IOSTYP INPUT: ',domain%IOSTYP
      CALL extcde (4)
    END IF
 
    ! save namelist
    nml_domain = domain
 
  END SUBROUTINE read_domain_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE read_input_grid_nml (NDSI, NRINP, NML_INPUT_GRID)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         22-Mar-2021 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
    !/    22-Mar-2021 : Update namelist, new input fields   ( version 7.13 )
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NDSI              Int.
    !      NRINP             Int.
    !      NML_INPUT_GRID    Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLMULTICONF Subr.   N/A    Namelist configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
 
    USE wmmdatmd, ONLY: mdse
    USE w3servmd, ONLY: extcde
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    INTEGER, INTENT(IN)               :: NDSI, NRINP
    TYPE(nml_input_grid_t), INTENT(INOUT) :: NML_INPUT_GRID(NRINP)
 
    ! locals
    INTEGER                   :: IERR, I
    INTEGER, PARAMETER :: MAX_NRINP = 99
    TYPE(nml_input_grid_t) :: INPUT(MAX_NRINP)
    namelist /input_grid_nml/ input
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
    ierr = 0
#ifdef W3_S
    CALL strace (ient, 'READ_INPUT_GRID_NML')
#endif
 
    ! test NRINP
    IF (nrinp.GT.max_nrinp) THEN
      WRITE (mdse,'(A,/A,I6,/A,I6)') &
           'ERROR: READ_INPUT_GRID_NML: NRINP > MAX_NRINP',      &
           'ERROR: READ_INPUT_GRID_NML:     NRINP = ',    nrinp, &
           'ERROR: READ_INPUT_GRID_NML: MAX_NRINP = ',max_nrinp
      CALL extcde (6)
    END IF
 
    ! if no input grids, then exit
    IF (nrinp.EQ.0) RETURN
 
    ! set default values for grid INPUT data
    DO i = 1,nrinp
      input(i)%NAME                  = 'unset'
      input(i)%FORCING%WATER_LEVELS  = .false.
      input(i)%FORCING%CURRENTS      = .false.
      input(i)%FORCING%WINDS         = .false.
      input(i)%FORCING%ATM_MOMENTUM  = .false.
      input(i)%FORCING%AIR_DENSITY   = .false.
      input(i)%FORCING%ICE_CONC      = .false.
      input(i)%FORCING%ICE_PARAM1    = .false.
      input(i)%FORCING%ICE_PARAM2    = .false.
      input(i)%FORCING%ICE_PARAM3    = .false.
      input(i)%FORCING%ICE_PARAM4    = .false.
      input(i)%FORCING%ICE_PARAM5    = .false.
      input(i)%FORCING%MUD_DENSITY   = .false.
      input(i)%FORCING%MUD_THICKNESS = .false.
      input(i)%FORCING%MUD_VISCOSITY = .false.
      input(i)%ASSIM%MEAN            = .false.
      input(i)%ASSIM%SPEC1D          = .false.
      input(i)%ASSIM%SPEC2D          = .false.
    END DO
 
    ! read input grid namelist
    rewind(ndsi)
    READ (ndsi, nml=input_grid_nml, iostat=ierr, iomsg=msg)
    IF (ierr.GT.0) THEN
      WRITE (mdse,'(A,/A)') &
           'ERROR: READ_INPUT_GRID_NML: namelist read error', &
           'ERROR: '//trim(msg)
      CALL extcde (7)
    END IF
 
    ! set/check RETURN values
    DO i = 1,nrinp
      IF (trim(input(i)%NAME).EQ.'unset') THEN
        WRITE (mdse,10) 'ERROR: READ_INPUT_GRID_NML: REQUIRED INPUT(',i,')%NAME NOT SET'
        CALL extcde (8)
      END IF
    END DO
 
    ! save namelist
    nml_input_grid = input(1:nrinp)
 
10  FORMAT (a,i0,a)
 
  END SUBROUTINE read_input_grid_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE read_model_grid_nml (NDSI, NRGRD, NML_MODEL_GRID)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         22-Mar-2021 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
    !/    22-Mar-2021 : Update namelist, new input fields   ( version 7.13 )
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NDSI              Int.
    !      NRGRD             Int.
    !      NML_MODEL_GRID    Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLMULTICONF Subr.   N/A    Namelist configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
 
    USE wmmdatmd, ONLY: mdse
    USE w3servmd, ONLY: extcde
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    INTEGER, INTENT(IN)               :: NDSI, NRGRD
    TYPE(nml_model_grid_t), INTENT(INOUT) :: NML_MODEL_GRID(NRGRD)
 
    ! locals
    INTEGER                   :: IERR, I
    INTEGER, PARAMETER :: MAX_NRGRD = 99
    TYPE(nml_model_grid_t) :: MODEL(MAX_NRGRD)
    namelist /model_grid_nml/ model
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
    ierr = 0
#ifdef W3_S
    CALL strace (ient, 'READ_MODEL_GRID_NML')
#endif
 
    ! test NRGRD
    IF (nrgrd.GT.max_nrgrd) THEN
      WRITE (mdse,'(A,/A,I6,/A,I6)') &
           'ERROR: READ_MODEL_GRID_NML: NRGRD > MAX_NRGRD',      &
           'ERROR: READ_MODEL_GRID_NML:     NRGRD = ',    nrgrd, &
           'ERROR: READ_MODEL_GRID_NML: MAX_NRGRD = ',max_nrgrd
      CALL extcde (9)
    END IF
 
    ! if no model grids, then exit
    IF (nrgrd.EQ.0) RETURN
 
    ! set default values for model input data
    DO i = 1,nrgrd
      model(i)%NAME                  = 'unset'
      model(i)%FORCING%WATER_LEVELS  = 'no'
      model(i)%FORCING%CURRENTS      = 'no'
      model(i)%FORCING%WINDS         = 'no'
      model(i)%FORCING%ATM_MOMENTUM  = 'no'
      model(i)%FORCING%AIR_DENSITY   = 'no'
      model(i)%FORCING%ICE_CONC      = 'no'
      model(i)%FORCING%ICE_PARAM1    = 'no'
      model(i)%FORCING%ICE_PARAM2    = 'no'
      model(i)%FORCING%ICE_PARAM3    = 'no'
      model(i)%FORCING%ICE_PARAM4    = 'no'
      model(i)%FORCING%ICE_PARAM5    = 'no'
      model(i)%FORCING%MUD_DENSITY   = 'no'
      model(i)%FORCING%MUD_THICKNESS = 'no'
      model(i)%FORCING%MUD_VISCOSITY = 'no'
      model(i)%ASSIM%MEAN            = 'no'
      model(i)%ASSIM%SPEC1D          = 'no'
      model(i)%ASSIM%SPEC2D          = 'no'
      model(i)%RESOURCE%RANK_ID      = i
      model(i)%RESOURCE%GROUP_ID     = 1
      model(i)%RESOURCE%COMM_FRAC    = (/0.00,1.00/)
      model(i)%RESOURCE%BOUND_FLAG   = .false.
    END DO
 
    ! read model namelist
    rewind(ndsi)
    READ (ndsi, nml=model_grid_nml, iostat=ierr, iomsg=msg)
    IF (ierr.NE.0) THEN
      WRITE (mdse,'(A,/A)') &
           'ERROR: READ_MODEL_GRID_NML: namelist read error', &
           'ERROR: '//trim(msg)
      CALL extcde (10)
    END IF
 
    ! set/check return values
    DO i = 1,nrgrd
      IF (trim(model(i)%NAME).EQ.'unset') THEN
        WRITE (mdse,10) 'ERROR: READ_MODEL_GRID_NML: required MODEL(',i,')%NAME not set'
        CALL extcde (11)
      END IF
    END DO
 
    ! save namelist
    nml_model_grid = model(1:nrgrd)
 
10  FORMAT (a,i0,a)
 
  END SUBROUTINE read_model_grid_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE read_output_type_nml (NDSI, NRGRD, NML_OUTPUT_TYPE)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         25-Sep-2020 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
    !/    25-Sep-2020 : Update namelist                     ( version 7.10 )
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NDSI              Int.
    !      NRGRD             Int.
    !      NML_OUTPUT_TYPE   Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLMULTICONF Subr.   N/A    Namelist configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
 
    USE wmmdatmd, ONLY: mdse
    USE w3servmd, ONLY: extcde
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    INTEGER, INTENT(IN)                :: NDSI, NRGRD
    TYPE(nml_output_type_t), INTENT(INOUT) :: NML_OUTPUT_TYPE(NRGRD)
 
    ! locals
    INTEGER                   :: IERR, I
    INTEGER, PARAMETER        :: MAX_NRGRD = 99
    TYPE(nml_output_type_t)   :: ALLTYPE
    TYPE(nml_output_type_t)   :: ITYPE(MAX_NRGRD)
    namelist /output_type_nml/ alltype, itype
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
    ierr = 0
#ifdef W3_S
    CALL strace (ient, 'READ_OUTPUT_TYPE_NML')
#endif
 
    ! if no model grids, then exit
    IF (nrgrd.EQ.0) RETURN
 
    ! set default values for output type data
    DO i=1,nrgrd
      itype(i)%FIELD%LIST   = 'unset'
      itype(i)%POINT%NAME   = 'unset'
      itype(i)%POINT%FILE   = 'points.list'
      itype(i)%TRACK%FORMAT = .true.
      itype(i)%PARTITION%X0 = 0
      itype(i)%PARTITION%XN = 0
      itype(i)%PARTITION%NX = 0
      itype(i)%PARTITION%Y0 = 0
      itype(i)%PARTITION%YN = 0
      itype(i)%PARTITION%NY = 0
      itype(i)%PARTITION%FORMAT = .true.
#ifdef W3_COU
      itype(i)%COUPLING%SENT      = 'unset'
      itype(i)%COUPLING%RECEIVED  = 'unset'
      itype(i)%COUPLING%COUPLET0  = .false.
#endif
      itype(i)%RESTART%EXTRA = 'unset'
    END DO
    alltype%FIELD%LIST   = 'unset'
    alltype%POINT%NAME   = 'unset'
    alltype%POINT%FILE   = 'points.list'
    alltype%TRACK%FORMAT = .true.
    alltype%PARTITION%X0 = 0
    alltype%PARTITION%XN = 0
    alltype%PARTITION%NX = 0
    alltype%PARTITION%Y0 = 0
    alltype%PARTITION%YN = 0
    alltype%PARTITION%NY = 0
    alltype%PARTITION%FORMAT = .true.
#ifdef W3_COU
    alltype%COUPLING%SENT      = 'unset'
    alltype%COUPLING%RECEIVED  = 'unset'
    alltype%COUPLING%COUPLET0  = .false.
#endif
    alltype%RESTART%EXTRA = 'unset'
 
 
    ! read OUTPUT_TYPE namelist
    rewind(ndsi)
    READ (ndsi, nml=output_type_nml, iostat=ierr, iomsg=msg)
    IF (ierr.GT.0) THEN
      WRITE (mdse,'(A,/A)') &
           'ERROR: READ_OUTPUT_TYPE_NML: namelist read error', &
           'ERROR: '//trim(msg)
      CALL extcde (12)
    END IF
 
    ! propagate general type on all grids
    DO i=1,nrgrd
      itype(i) = alltype
    END DO
 
    ! read OUTPUT_TYPE namelist
    rewind(ndsi)
    READ (ndsi, nml=output_type_nml, iostat=ierr, iomsg=msg)
    IF (ierr.GT.0) THEN
      WRITE (mdse,'(A,/A)') &
           'ERROR: READ_OUTPUT_TYPE_NML: namelist read error', &
           'ERROR: '//trim(msg)
      CALL extcde (13)
    END IF
 
    ! save namelist
    nml_output_type = itype(1:nrgrd)
 
 
  END SUBROUTINE read_output_type_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE read_output_date_nml (NDSI, NRGRD, NML_OUTPUT_DATE)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         15-May-2018 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NDSI              Int.
    !      NRGRD             Int.
    !      NML_OUTPUT_DATE   Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLMULTICONF Subr.   N/A    Namelist configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
 
    USE wmmdatmd, ONLY: mdse
    USE w3servmd, ONLY: extcde
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    INTEGER, INTENT(IN)                :: NDSI, NRGRD
    TYPE(nml_output_date_t), INTENT(INOUT) :: NML_OUTPUT_DATE(NRGRD)
 
    ! locals
    INTEGER                   :: IERR, I
    INTEGER, PARAMETER        :: MAX_NRGRD = 99
    TYPE(nml_output_date_t)   :: ALLDATE
    TYPE(nml_output_date_t)   :: IDATE(MAX_NRGRD)
    namelist /output_date_nml/ alldate, idate
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
    ierr = 0
#ifdef W3_S
    CALL strace (ient, 'READ_OUTPUT_DATE_NML')
#endif
 
    ! if no model grids, then exit
    IF (nrgrd.EQ.0) RETURN
 
    ! set default values for OUTPUT_DATE input data
    DO i=1,nrgrd
      idate(i)%FIELD%START = '19680606 000000'
      idate(i)%FIELD%STRIDE = '0'
      idate(i)%FIELD%STOP = '19680607 000000'
      idate(i)%FIELD%OUTFFILE = '0'
      idate(i)%POINT%OUTFFILE = '0'
      idate(i)%POINT%START = '19680606 000000'
      idate(i)%POINT%STRIDE = '0'
      idate(i)%POINT%STOP = '19680607 000000'
      idate(i)%TRACK%START = '19680606 000000'
      idate(i)%TRACK%STRIDE = '0'
      idate(i)%TRACK%STOP = '19680607 000000'
      idate(i)%RESTART%START = '19680606 000000'
      idate(i)%RESTART%STRIDE = '0'
      idate(i)%RESTART%STOP = '19680607 000000'
      idate(i)%RESTART2%START = '19680606 000000'
      idate(i)%RESTART2%STRIDE = '0'
      idate(i)%RESTART2%STOP = '19680607 000000'
      idate(i)%BOUNDARY%START = '19680606 000000'
      idate(i)%BOUNDARY%STRIDE = '0'
      idate(i)%BOUNDARY%STOP = '19680607 000000'
      idate(i)%PARTITION%START = '19680606 000000'
      idate(i)%PARTITION%STRIDE = '0'
      idate(i)%PARTITION%STOP = '19680607 000000'
#ifdef W3_COU
      idate(i)%COUPLING%START = '19680606 000000'
      idate(i)%COUPLING%STRIDE = '0'
      idate(i)%COUPLING%STOP = '19680607 000000'
#endif
    END DO
    alldate%FIELD%START = '19680606 000000'
    alldate%FIELD%STRIDE = '0'
    alldate%FIELD%STOP = '19680607 000000'
    alldate%FIELD%OUTFFILE = '0'
    alldate%POINT%START = '19680606 000000'
    alldate%POINT%STRIDE = '0'
    alldate%POINT%STOP = '19680607 000000'
    alldate%POINT%OUTFFILE = '0'
    alldate%TRACK%START = '19680606 000000'
    alldate%TRACK%STRIDE = '0'
    alldate%TRACK%STOP = '19680607 000000'
    alldate%RESTART%START = '19680606 000000'
    alldate%RESTART%STRIDE = '0'
    alldate%RESTART%STOP   = '19680607 000000'
    alldate%RESTART2%START  = '19680606 000000'
    alldate%RESTART2%STRIDE = '0'
    alldate%RESTART2%STOP   = '19680607 000000'
    alldate%BOUNDARY%START = '19680606 000000'
    alldate%BOUNDARY%STRIDE = '0'
    alldate%BOUNDARY%STOP = '19680607 000000'
    alldate%PARTITION%START = '19680606 000000'
    alldate%PARTITION%STRIDE = '0'
    alldate%PARTITION%STOP = '19680607 000000'
#ifdef W3_COU
    alldate%COUPLING%START = '19680606 000000'
    alldate%COUPLING%STRIDE = '0'
    alldate%COUPLING%STOP = '19680607 000000'
#endif
 
 
    ! read OUTPUT_DATE namelist
    rewind(ndsi)
    READ (ndsi, nml=output_date_nml, iostat=ierr, iomsg=msg)
    IF (ierr.GT.0) THEN
      WRITE (mdse,'(A,/A)') &
           'ERROR: READ_OUTPUT_DATE_NML: namelist read error', &
           'ERROR: '//trim(msg)
      CALL extcde (14)
    END IF
 
    ! propagate general date to all grids
    DO i=1,nrgrd
      idate(i) = alldate
    END DO
 
    ! read OUTPUT_DATE namelist
    rewind(ndsi)
    READ (ndsi, nml=output_date_nml, iostat=ierr, iomsg=msg)
    IF (ierr.GT.0) THEN
      WRITE (mdse,'(A,/A)') &
           'ERROR: READ_OUTPUT_DATE_NML: namelist read error', &
           'ERROR: '//trim(msg)
      CALL extcde (15)
    END IF
 
    ! save namelist
    nml_output_date = idate(1:nrgrd)
 
 
  END SUBROUTINE read_output_date_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE read_homogeneous_nml (NDSI, NML_HOMOG_COUNT, NML_HOMOG_INPUT)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         15-May-2018 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NDSI               Int.
    !      NML_HOMOG_COUNT    Type.
    !      NML_HOMOG_INPUT    Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLMULTICONF Subr.   N/A    Namelist configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
 
    USE wmmdatmd, ONLY: mdse
    USE w3servmd, ONLY: extcde
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    INTEGER, INTENT(IN)                      :: NDSI
    TYPE(nml_homog_count_t), INTENT(OUT)  :: NML_HOMOG_COUNT
    TYPE(nml_homog_input_t), ALLOCATABLE, INTENT(OUT)   :: NML_HOMOG_INPUT(:)
 
    ! locals
    INTEGER                   :: IERR, I
    TYPE(nml_homog_count_t)  :: HOMOG_COUNT
    namelist /homog_count_nml/  homog_count
    TYPE(nml_homog_input_t), ALLOCATABLE   :: HOMOG_INPUT(:)
    namelist /homog_input_nml/   homog_input
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
    ierr = 0
#ifdef W3_S
    CALL strace (ient, 'READ_HOMOGENEOUS_NML')
#endif
 
 
    ! set default values for homogeneous number structure
    homog_count%N_MOV   = 0
    homog_count%N_TOT   = 0
 
 
    ! read homogeneous count namelist
    rewind(ndsi)
    READ (ndsi, nml=homog_count_nml, iostat=ierr, iomsg=msg)
    IF (ierr.GT.0) THEN
      WRITE (mdse,'(A,/A)') &
           'ERROR: READ_HOMOGENEOUS_NML: namelist HOMOG_COUNT read error', &
           'ERROR: '//trim(msg)
      CALL extcde (16)
    END IF
 
    ! allocate the total count of homogeneous input
    homog_count%N_TOT = homog_count%N_MOV
    ALLOCATE(homog_input(homog_count%N_TOT))
    ALLOCATE(nml_homog_input(homog_count%N_TOT))
 
    ! set default values for homogeneous structure
    IF (homog_count%N_TOT .NE. 0 ) THEN
      DO i=1,homog_count%N_TOT
        homog_input(i)%NAME      = 'unset'
        homog_input(i)%DATE      = '19680606 000000'
        homog_input(i)%VALUE1    = 0.
        homog_input(i)%VALUE2    = 0.
        homog_input(i)%VALUE3    = 0.
      END DO
    END IF
 
    ! read homogeneous input namelist
    rewind(ndsi)
    READ (ndsi, nml=homog_input_nml, iostat=ierr, iomsg=msg)
    IF (ierr.GT.0) THEN
      WRITE (mdse,'(A,/A)') &
           'ERROR: READ_HOMOGENEOUS_NML: namelist HOMOG_INPUT_NML read error', &
           'ERROR: '//trim(msg)
      CALL extcde (17)
    END IF
 
    ! save namelist
    nml_homog_count = homog_count
    nml_homog_input = homog_input
 
 
  END SUBROUTINE read_homogeneous_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE report_domain_nml (NML_DOMAIN)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         15-May-2018 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NML_DOMAIN       Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLMULTICONF Subr.   N/A    Namelist configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
 
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    TYPE(nml_domain_t), INTENT(IN) :: NML_DOMAIN
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
#ifdef W3_S
    CALL strace (ient, 'REPORT_DOMAIN_NML')
#endif
 
    WRITE (msg,'(A)') 'DOMAIN % '
    WRITE (ndsn,'(A)')
    WRITE (ndsn,11) trim(msg),'NRINP  = ', nml_domain%NRINP
    WRITE (ndsn,11) trim(msg),'NRGRD  = ', nml_domain%NRGRD
    WRITE (ndsn,13) trim(msg),'UNIPTS = ', nml_domain%UNIPTS
    WRITE (ndsn,11) trim(msg),'IOSTYP = ', nml_domain%IOSTYP
    WRITE (ndsn,13) trim(msg),'UPPROC = ', nml_domain%UPPROC
    WRITE (ndsn,13) trim(msg),'PSHARE = ', nml_domain%PSHARE
    WRITE (ndsn,13) trim(msg),'FLGHG1 = ', nml_domain%FLGHG1
    WRITE (ndsn,13) trim(msg),'FLGHG2 = ', nml_domain%FLGHG2
    WRITE (ndsn,10) trim(msg),'START  = ', trim(nml_domain%START)
    WRITE (ndsn,10) trim(msg),'STOP   = ', trim(nml_domain%STOP)
 
10  FORMAT (a,2x,a,a)
11  FORMAT (a,2x,a,i8)
13  FORMAT (a,2x,a,l1)
 
  END SUBROUTINE report_domain_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE report_input_grid_nml (NRINP, NML_INPUT_GRID)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         22-Mar-2021 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
    !/    22-Mar-2021 : Update namelist, new input fields   ( version 7.13 )
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NRINP             Int.
    !      NML_INPUT_GRID    Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLMULTICONF Subr.   N/A    Namelist configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
 
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    INTEGER, INTENT(IN)            :: NRINP
    TYPE(nml_input_grid_t), INTENT(IN) :: NML_INPUT_GRID(NRINP)
 
    ! locals
    INTEGER              :: I
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
#ifdef W3_S
    CALL strace (ient, 'REPORT_INPUT_GRID_NML')
#endif
 
    DO i = 1,nrinp
      WRITE (msg,'(A,I0.2,A)') 'INPUT GRID ',i,' % '
      WRITE (ndsn,'(A)')
      WRITE (ndsn,10) trim(msg),'NAME                     = ', trim(nml_input_grid(i)%NAME)
      WRITE (ndsn,13) trim(msg),'FORCING % WATER_LEVELS   = ', nml_input_grid(i)%FORCING%WATER_LEVELS
      WRITE (ndsn,13) trim(msg),'FORCING % CURRENTS       = ', nml_input_grid(i)%FORCING%CURRENTS
      WRITE (ndsn,13) trim(msg),'FORCING % WINDS          = ', nml_input_grid(i)%FORCING%WINDS
      WRITE (ndsn,13) trim(msg),'FORCING % ATM_MOMENTUM   = ', nml_input_grid(i)%FORCING%ATM_MOMENTUM
      WRITE (ndsn,13) trim(msg),'FORCING % AIR_DENSITY    = ', nml_input_grid(i)%FORCING%AIR_DENSITY
      WRITE (ndsn,13) trim(msg),'FORCING % ICE_CONC       = ', nml_input_grid(i)%FORCING%ICE_CONC
      WRITE (ndsn,13) trim(msg),'FORCING % ICE_PARAM1     = ', nml_input_grid(i)%FORCING%ICE_PARAM1
      WRITE (ndsn,13) trim(msg),'FORCING % ICE_PARAM2     = ', nml_input_grid(i)%FORCING%ICE_PARAM2
      WRITE (ndsn,13) trim(msg),'FORCING % ICE_PARAM3     = ', nml_input_grid(i)%FORCING%ICE_PARAM3
      WRITE (ndsn,13) trim(msg),'FORCING % ICE_PARAM4     = ', nml_input_grid(i)%FORCING%ICE_PARAM4
      WRITE (ndsn,13) trim(msg),'FORCING % ICE_PARAM5     = ', nml_input_grid(i)%FORCING%ICE_PARAM5
      WRITE (ndsn,13) trim(msg),'FORCING % MUD_DENSITY    = ', nml_input_grid(i)%FORCING%MUD_DENSITY
      WRITE (ndsn,13) trim(msg),'FORCING % MUD_THICKNESS  = ', nml_input_grid(i)%FORCING%MUD_THICKNESS
      WRITE (ndsn,13) trim(msg),'FORCING % MUD_VISCOSITY  = ', nml_input_grid(i)%FORCING%MUD_VISCOSITY
      WRITE (ndsn,13) trim(msg),'ASSIM % MEAN             = ', nml_input_grid(i)%ASSIM%MEAN
      WRITE (ndsn,13) trim(msg),'ASSIM % SPEC1D           = ', nml_input_grid(i)%ASSIM%SPEC1D
      WRITE (ndsn,13) trim(msg),'ASSIM % SPEC2D           = ', nml_input_grid(i)%ASSIM%SPEC2D
    END DO
    WRITE (*,*)
 
10  FORMAT (a,2x,a,a)
13  FORMAT (a,2x,a,l1)
 
  END SUBROUTINE report_input_grid_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE report_model_grid_nml (NRGRD, NML_MODEL_GRID)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         22-Mar-2021 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
    !/    22-Mar-2021 : Update namelist, new input fields   ( version 7.13 )
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NRGRD             Int.
    !      NML_MODEL_GRID    Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLMULTICONF Subr.   N/A    Namelist configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
 
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    INTEGER, INTENT(IN)            :: NRGRD
    TYPE(nml_model_grid_t), INTENT(IN) :: NML_MODEL_GRID(NRGRD)
 
    ! locals
    INTEGER              :: I
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
#ifdef W3_S
    CALL strace (ient, 'REPORT_MODEL_GRID_NML')
#endif
 
    DO i = 1,nrgrd
      WRITE (msg,'(A,I0.4,A)') 'MODEL GRID ',i,' % '
      WRITE (ndsn,'(A)')
      WRITE (ndsn,10) trim(msg),'NAME                     = ', trim(nml_model_grid(i)%NAME)
      WRITE (ndsn,10) trim(msg),'FORCING % WATER_LEVELS   = ', trim(nml_model_grid(i)%FORCING%WATER_LEVELS)
      WRITE (ndsn,10) trim(msg),'FORCING % CURRENTS       = ', trim(nml_model_grid(i)%FORCING%CURRENTS)
      WRITE (ndsn,10) trim(msg),'FORCING % WINDS          = ', trim(nml_model_grid(i)%FORCING%WINDS)
      WRITE (ndsn,10) trim(msg),'FORCING % ATM_MOMENTUM   = ', trim(nml_model_grid(i)%FORCING%ATM_MOMENTUM)
      WRITE (ndsn,10) trim(msg),'FORCING % AIR_DENSITY    = ', trim(nml_model_grid(i)%FORCING%AIR_DENSITY)
      WRITE (ndsn,10) trim(msg),'FORCING % ICE_CONC       = ', trim(nml_model_grid(i)%FORCING%ICE_CONC)
      WRITE (ndsn,10) trim(msg),'FORCING % ICE_PARAM1     = ', trim(nml_model_grid(i)%FORCING%ICE_PARAM1)
      WRITE (ndsn,10) trim(msg),'FORCING % ICE_PARAM2     = ', trim(nml_model_grid(i)%FORCING%ICE_PARAM2)
      WRITE (ndsn,10) trim(msg),'FORCING % ICE_PARAM3     = ', trim(nml_model_grid(i)%FORCING%ICE_PARAM3)
      WRITE (ndsn,10) trim(msg),'FORCING % ICE_PARAM4     = ', trim(nml_model_grid(i)%FORCING%ICE_PARAM4)
      WRITE (ndsn,10) trim(msg),'FORCING % ICE_PARAM5     = ', trim(nml_model_grid(i)%FORCING%ICE_PARAM5)
      WRITE (ndsn,10) trim(msg),'FORCING % MUD_DENSITY    = ', trim(nml_model_grid(i)%FORCING%MUD_DENSITY)
      WRITE (ndsn,10) trim(msg),'FORCING % MUD_THICKNESS  = ', trim(nml_model_grid(i)%FORCING%MUD_THICKNESS)
      WRITE (ndsn,10) trim(msg),'FORCING % MUD_VISCOSITY  = ', trim(nml_model_grid(i)%FORCING%MUD_VISCOSITY)
      WRITE (ndsn,10) trim(msg),'ASSIM % MEAN             = ', trim(nml_model_grid(i)%ASSIM%MEAN)
      WRITE (ndsn,10) trim(msg),'ASSIM % SPEC1D           = ', trim(nml_model_grid(i)%ASSIM%SPEC1D)
      WRITE (ndsn,10) trim(msg),'ASSIM % SPEC2D           = ', trim(nml_model_grid(i)%ASSIM%SPEC2D)
      WRITE (ndsn,11) trim(msg),'RESOURCE % RANK_ID       = ', nml_model_grid(i)%RESOURCE%RANK_ID
      WRITE (ndsn,11) trim(msg),'RESOURCE % GROUP_ID      = ', nml_model_grid(i)%RESOURCE%GROUP_ID
      WRITE (ndsn,12) trim(msg),'RESOURCE % COMM_FRAC     = ', nml_model_grid(i)%RESOURCE%COMM_FRAC(1), &
           nml_model_grid(i)%RESOURCE%COMM_FRAC(2)
      WRITE (ndsn,13) trim(msg),'RESOURCE % BOUND_FLAG    = ', nml_model_grid(i)%RESOURCE%BOUND_FLAG
    END DO
    WRITE (ndsn,'(A)')
 
10  FORMAT (a,2x,a,a)
11  FORMAT (a,2x,a,i8)
12  FORMAT (a,2x,a,'(',f5.2,',',f5.2,' )')
13  FORMAT (a,2x,a,l1)
 
  END SUBROUTINE report_model_grid_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE report_output_type_nml (NRGRD, NML_OUTPUT_TYPE)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         25-Sep-2020 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
    !/    25-Sep-2020 : Update namelist                     ( version 7.10 )
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NRGRD             Int.
    !      NML_OUTPUT_TYPE   Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLMULTICONF Subr.   N/A    Namelist configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
 
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    INTEGER, INTENT(IN)             :: NRGRD
    TYPE(nml_output_type_t), INTENT(IN) :: NML_OUTPUT_TYPE(NRGRD)
 
    ! locals
    INTEGER              :: I
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
#ifdef W3_S
    CALL strace (ient, 'REPORT_OUTPUT_TYPE_NML')
#endif
 
    DO i=1,nrgrd
      WRITE (msg,'(A,I1,A)') 'OUTPUT TYPE MODEL GRID',i, ' % '
      WRITE (ndsn,'(A)')
      WRITE (ndsn,10) trim(msg),'FIELD % LIST         = ', trim(nml_output_type(i)%FIELD%LIST)
      WRITE (ndsn,10) trim(msg),'POINT % NAME         = ', trim(nml_output_type(i)%POINT%NAME)
      WRITE (ndsn,10) trim(msg),'POINT % FILE         = ', trim(nml_output_type(i)%POINT%FILE)
      WRITE (ndsn,13) trim(msg),'TRACK % FORMAT       = ', nml_output_type(i)%TRACK%FORMAT
      WRITE (ndsn,11) trim(msg),'PARTITION % X0       = ', nml_output_type(i)%PARTITION%X0
      WRITE (ndsn,11) trim(msg),'PARTITION % XN       = ', nml_output_type(i)%PARTITION%XN
      WRITE (ndsn,11) trim(msg),'PARTITION % NX       = ', nml_output_type(i)%PARTITION%NX
      WRITE (ndsn,11) trim(msg),'PARTITION % Y0       = ', nml_output_type(i)%PARTITION%Y0
      WRITE (ndsn,11) trim(msg),'PARTITION % YN       = ', nml_output_type(i)%PARTITION%YN
      WRITE (ndsn,11) trim(msg),'PARTITION % NY       = ', nml_output_type(i)%PARTITION%NY
      WRITE (ndsn,13) trim(msg),'PARTITION % FORMAT   = ', nml_output_type(i)%PARTITION%FORMAT
#ifdef W3_COU
      WRITE (ndsn,10) trim(msg),'COUPLING % SENT         = ', trim(nml_output_type(i)%COUPLING%SENT)
      WRITE (ndsn,10) trim(msg),'COUPLING % RECEIVED     = ', trim(nml_output_type(i)%COUPLING%RECEIVED)
      WRITE (ndsn,13) trim(msg),'COUPLING % COUPLET0     = ', nml_output_type(i)%COUPLING%COUPLET0
#endif
      WRITE (ndsn,10) trim(msg),'RESTART % EXTRA      = ', trim(nml_output_type(i)%RESTART%EXTRA)
    END DO
    WRITE (ndsn,'(A)')
 
10  FORMAT (a,2x,a,a)
11  FORMAT (a,2x,a,i8)
13  FORMAT (a,2x,a,l1)
 
  END SUBROUTINE report_output_type_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE report_output_date_nml (NRGRD, NML_OUTPUT_DATE)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         15-May-2018 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NRGRD             Int.
    !      NML_OUTPUT_DATE   Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLMULTICONF Subr.   N/A    Namelist configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
 
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    INTEGER, INTENT(IN)             :: NRGRD
    TYPE(nml_output_date_t), INTENT(IN) :: NML_OUTPUT_DATE(NRGRD)
 
    ! locals
    INTEGER              :: I
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
#ifdef W3_S
    CALL strace (ient, 'REPORT_OUTPUT_DATE_NML')
#endif
 
    DO i=1,nrgrd
      WRITE (msg,'(A,I1,A)') 'OUTPUT DATE MODEL GRID',i, ' % '
      WRITE (ndsn,'(A)')
      WRITE (ndsn,10) trim(msg),'FIELD % START        = ', trim(nml_output_date(i)%FIELD%START)
      WRITE (ndsn,10) trim(msg),'FIELD % STRIDE       = ', trim(nml_output_date(i)%FIELD%STRIDE)
      WRITE (ndsn,10) trim(msg),'FIELD % STOP         = ', trim(nml_output_date(i)%FIELD%STOP)
      WRITE (ndsn,10) trim(msg),'POINT % START        = ', trim(nml_output_date(i)%POINT%START)
      WRITE (ndsn,10) trim(msg),'POINT % STRIDE       = ', trim(nml_output_date(i)%POINT%STRIDE)
      WRITE (ndsn,10) trim(msg),'POINT % STOP         = ', trim(nml_output_date(i)%POINT%STOP)
      WRITE (ndsn,10) trim(msg),'TRACK % START        = ', trim(nml_output_date(i)%TRACK%START)
      WRITE (ndsn,10) trim(msg),'TRACK % STRIDE       = ', trim(nml_output_date(i)%TRACK%STRIDE)
      WRITE (ndsn,10) trim(msg),'TRACK % STOP         = ', trim(nml_output_date(i)%TRACK%STOP)
      WRITE (ndsn,10) trim(msg),'RESTART % START      = ', trim(nml_output_date(i)%RESTART%START)
      WRITE (ndsn,10) trim(msg),'RESTART % STRIDE     = ', trim(nml_output_date(i)%RESTART%STRIDE)
      WRITE (ndsn,10) trim(msg),'RESTART % STOP       = ', trim(nml_output_date(i)%RESTART%STOP)
      WRITE (ndsn,10) trim(msg),'RESTART2 % START     = ', trim(nml_output_date(i)%RESTART2%START)
      WRITE (ndsn,10) trim(msg),'RESTART2 % STRIDE    = ', trim(nml_output_date(i)%RESTART2%STRIDE)
      WRITE (ndsn,10) trim(msg),'RESTART2 % STOP      = ', trim(nml_output_date(i)%RESTART2%STOP)
      WRITE (ndsn,10) trim(msg),'BOUNDARY % START     = ', trim(nml_output_date(i)%BOUNDARY%START)
      WRITE (ndsn,10) trim(msg),'BOUNDARY % STRIDE    = ', trim(nml_output_date(i)%BOUNDARY%STRIDE)
      WRITE (ndsn,10) trim(msg),'BOUNDARY % STOP      = ', trim(nml_output_date(i)%BOUNDARY%STOP)
      WRITE (ndsn,10) trim(msg),'PARTITION % START    = ', trim(nml_output_date(i)%PARTITION%START)
      WRITE (ndsn,10) trim(msg),'PARTITION % STRIDE   = ', trim(nml_output_date(i)%PARTITION%STRIDE)
      WRITE (ndsn,10) trim(msg),'PARTITION % STOP     = ', trim(nml_output_date(i)%PARTITION%STOP)
#ifdef W3_COU
      WRITE (ndsn,10) trim(msg),'COUPLING % START    = ', trim(nml_output_date(i)%COUPLING%START)
      WRITE (ndsn,10) trim(msg),'COUPLING % STRIDE   = ', trim(nml_output_date(i)%COUPLING%STRIDE)
      WRITE (ndsn,10) trim(msg),'COUPLING % STOP     = ', trim(nml_output_date(i)%COUPLING%STOP)
#endif
    END DO
    WRITE (ndsn,'(A)')
 
10  FORMAT (a,2x,a,a)
 
  END SUBROUTINE report_output_date_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE report_homogeneous_nml (NML_HOMOG_COUNT, NML_HOMOG_INPUT)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         15-May-2018 |
    !/                  +-----------------------------------+
    !/
    !/    09-Aug-2016 : Adding comments                     ( version 5.12 )
    !/    15-May-2018 : Update namelist                     ( version 6.05 )
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NML_HOMOG_COUNT    Type.
    !      NML_HOMOG_INPUT    Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLMULTICONF Subr.   N/A    Namelist configuration routine.
    !     ----------------------------------------------------------------
    !
    !  6. Error messages :
    !
    !     None.
    !
    !  7. Remarks :
    !
    !  8. Structure :
    !
    !     See source code.
    !
    !  9. Switches :
    !
    !     !/MPI  Uses MPI communications
    !
    ! 10. Source code :
    !
    !/ ------------------------------------------------------------------- /
 
#ifdef W3_S
    USE w3servmd, ONLY: strace
#endif
 
    IMPLICIT NONE
 
    TYPE(nml_homog_count_t), INTENT(IN)  :: NML_HOMOG_COUNT
    TYPE(nml_homog_input_t), INTENT(IN) :: NML_HOMOG_INPUT(NML_HOMOG_COUNT%N_TOT)
 
    ! locals
    INTEGER              :: I
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
#ifdef W3_S
    CALL strace (ient, 'REPORT_HOMOGENEOUS_NML')
#endif
 
    WRITE (msg,'(A)') 'HOMOG_COUNT % '
    WRITE (ndsn,'(A)')
    WRITE (ndsn,11) trim(msg),'N_MOV       = ', nml_homog_count%N_MOV
 
    IF (nml_homog_count%N_TOT .NE. 0) THEN
      DO i=1,nml_homog_count%N_TOT
        WRITE (msg,'(A,I5,A)') 'HOMOG_INPUT(',i,') % '
        WRITE (ndsn,'(A)')
        WRITE (ndsn,10) trim(msg),'NAME      = ', trim(nml_homog_input(i)%NAME)
        WRITE (ndsn,10) trim(msg),'DATE      = ', trim(nml_homog_input(i)%DATE)
        WRITE (ndsn,14) trim(msg),'VALUE1    = ', nml_homog_input(i)%VALUE1
        WRITE (ndsn,14) trim(msg),'VALUE2    = ', nml_homog_input(i)%VALUE2
        WRITE (ndsn,14) trim(msg),'VALUE3    = ', nml_homog_input(i)%VALUE3
        WRITE (ndsn,'(A)')
      END DO
    END IF
 
 
10  FORMAT (a,2x,a,a)
11  FORMAT (a,2x,a,i8)
14  FORMAT (a,2x,a,f8.2)
 
 
  END SUBROUTINE report_homogeneous_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
END MODULE w3nmlmultimd
 
!/ ------------------------------------------------------------------- /