w3nmlshelmd.F90

Go to the documentation of this file.

#include "w3macros.h"
!/ ------------------------------------------------------------------- /
MODULE w3nmlshelmd
  !/
  !/                  +-----------------------------------+
  !/                  | 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_shel.nml for ww3_shel program
  !
  !/ ------------------------------------------------------------------- /
 
  ! module defaults
  IMPLICIT NONE
 
  PUBLIC
 
  ! domain structure
  TYPE nml_domain_t
    INTEGER                     :: iostyp
    CHARACTER(15)               :: start
    CHARACTER(15)               :: stop
    CHARACTER(1)                :: outffile
  END TYPE nml_domain_t
 
 
 
  ! input structure
  TYPE nml_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_forcing_t
  !
  TYPE nml_assim_t
    CHARACTER(13)               :: mean
    CHARACTER(13)               :: spec1d
    CHARACTER(13)               :: spec2d
  END TYPE nml_assim_t
  !
  TYPE nml_input_t
    TYPE(nml_forcing_t)         :: forcing
    TYPE(nml_assim_t)           :: assim
  END TYPE nml_input_t
 
  ! output type structure
  TYPE nml_field_t
    CHARACTER(1024)             :: list
  END TYPE nml_field_t
  !
  TYPE nml_point_t
    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(15)               :: 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
    TYPE(nml_output_time_t)         :: coupling
  END TYPE nml_output_date_t
 
 
  ! homogeneous input structure
  TYPE nml_homog_count_t
    INTEGER                     :: n_ic1
    INTEGER                     :: n_ic2
    INTEGER                     :: n_ic3
    INTEGER                     :: n_ic4
    INTEGER                     :: n_ic5
    INTEGER                     :: n_mdn
    INTEGER                     :: n_mth
    INTEGER                     :: n_mvs
    INTEGER                     :: n_lev
    INTEGER                     :: n_cur
    INTEGER                     :: n_wnd
    INTEGER                     :: n_ice
    INTEGER                     :: n_tau
    INTEGER                     :: n_rho
    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 w3nmlshel (MPI_COMM, NDSI, INFILE, NML_DOMAIN,            &
       NML_INPUT, NML_OUTPUT_TYPE, NML_OUTPUT_DATE,   &
       NML_HOMOG_COUNT, NML_HOMOG_INPUT, IERR)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         18-Apr-2018 |
    !/                  +-----------------------------------+
    !/
 
    !  1. Purpose :
    !
    !     Reads all the namelist to define the single 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
    !      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_DOMAIN_NML
    !      REPORT_DOMAIN_NML
    !      READ_INPUT_NML
    !      REPORT_INPUT_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
    !     ----------------------------------------------------------------
    !      WW3_SHEL  Prog.   N/A    Single grid main program.
    !     ----------------------------------------------------------------
    !
    !  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_t), INTENT(INOUT)          :: NML_INPUT
    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, 'W3NMLSHEL')
#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 namelist
    CALL read_domain_nml (ndsi, nml_domain)
    IF ( improc .EQ. nmplog ) CALL report_domain_nml (nml_domain)
 
    ! read input namelist
    CALL read_input_nml (ndsi, nml_input)
    IF ( improc .EQ. nmplog ) CALL report_input_nml (nml_input)
 
    ! read output type namelist
    CALL read_output_type_nml (ndsi, nml_output_type)
    IF ( improc .EQ. nmplog ) CALL report_output_type_nml (nml_output_type)
 
    ! read output date namelist
    CALL read_output_date_nml (ndsi, nml_output_date)
    IF ( improc .EQ. nmplog ) CALL report_output_date_nml (nml_output_date)
 
    ! read homogeneous namelist
    CALL read_homogeneous_nml (ndsi, nml_homog_count, nml_homog_input)
    IF ( improc .EQ. nmplog ) CALL report_homogeneous_nml (nml_homog_count, nml_homog_input)
 
    ! close namelist files
    CLOSE (ndsi)
    IF ( nmplog .EQ. improc ) CLOSE (ndsn)
 
  END SUBROUTINE w3nmlshel
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE read_domain_nml (NDSI, NML_DOMAIN)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         18-Apr-2018 |
    !/                  +-----------------------------------+
    !/
 
    !  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
    !     ----------------------------------------------------------------
    !      W3NMLSHEL 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 domain structure
    domain%IOSTYP = 1
    domain%START  = '19680606 000000'
    domain%STOP   = '19680607 000000'
 
    ! read domain 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%IOSTYP.LT.0.OR.domain%IOSTYP.GT.3) THEN
      WRITE (mdse,'(A,I3)') 'ERROR: BAD IOSTYP INPUT: ',domain%IOSTYP
      CALL extcde (2)
    END IF
 
    ! save namelist
    nml_domain = domain
 
  END SUBROUTINE read_domain_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
 
  SUBROUTINE read_input_nml (NDSI, NML_INPUT)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         22-Mar-2021 |
    !/                  +-----------------------------------+
    !/
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NDSI              Int.
    !      NML_INPUT         Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLSHEL 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_input_t), INTENT(INOUT)   :: NML_INPUT
 
    ! locals
    INTEGER                                :: IERR
    TYPE(nml_input_t) :: INPUT
    namelist /input_nml/ input
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
    ierr = 0
#ifdef W3_S
    CALL strace (ient, 'READ_INPUT_NML')
#endif
 
 
    ! set default values for input structure
    input%FORCING%WATER_LEVELS  = 'F'
    input%FORCING%CURRENTS      = 'F'
    input%FORCING%WINDS         = 'F'
    input%FORCING%ATM_MOMENTUM  = 'F'
    input%FORCING%AIR_DENSITY   = 'F'
    input%FORCING%ICE_CONC      = 'F'
    input%FORCING%ICE_PARAM1    = 'F'
    input%FORCING%ICE_PARAM2    = 'F'
    input%FORCING%ICE_PARAM3    = 'F'
    input%FORCING%ICE_PARAM4    = 'F'
    input%FORCING%ICE_PARAM5    = 'F'
    input%FORCING%MUD_DENSITY   = 'F'
    input%FORCING%MUD_THICKNESS = 'F'
    input%FORCING%MUD_VISCOSITY = 'F'
    input%ASSIM%MEAN            = 'F'
    input%ASSIM%SPEC1D          = 'F'
    input%ASSIM%SPEC2D          = 'F'
 
    ! read input namelist
    rewind(ndsi)
    READ (ndsi, nml=input_nml, iostat=ierr, iomsg=msg)
    IF (ierr.GT.0) THEN
      WRITE (mdse,'(A,/A)') &
           'ERROR: READ_INPUT_NML: namelist read error', &
           'ERROR: '//trim(msg)
      CALL extcde (3)
    END IF
 
    ! save namelist
    nml_input = input
 
  END SUBROUTINE read_input_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE read_output_type_nml (NDSI, NML_OUTPUT_TYPE)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         25-Sep-2020 |
    !/                  +-----------------------------------+
    !/
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NDSI              Int.
    !      NML_OUTPUT_TYPE   Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLSHEL 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_output_type_t), INTENT(INOUT) :: NML_OUTPUT_TYPE
 
    ! locals
    INTEGER                                :: IERR
    TYPE(nml_output_type_t) :: TYPE
    namelist /output_type_nml/ TYPE
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
    ierr = 0
#ifdef W3_S
    CALL strace (ient, 'READ_OUTPUT_TYPE_NML')
#endif
 
    ! set default values for output type structure
    type%FIELD%LIST   = 'unset'
    type%POINT%FILE   = 'points.list'
    type%TRACK%FORMAT = .true.
    type%PARTITION%X0 = 0
    type%PARTITION%XN = 0
    type%PARTITION%NX = 0
    type%PARTITION%Y0 = 0
    type%PARTITION%YN = 0
    type%PARTITION%NY = 0
    type%PARTITION%FORMAT = .true.
#ifdef W3_COU
    type%COUPLING%SENT      = 'unset'
    type%COUPLING%RECEIVED  = 'unset'
    type%COUPLING%COUPLET0  = .false.
#endif
    type%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 (4)
    END IF
 
    ! save namelist
    nml_output_type = TYPE
 
  END SUBROUTINE read_output_type_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE read_output_date_nml (NDSI, NML_OUTPUT_DATE)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         18-Apr-2018 |
    !/                  +-----------------------------------+
    !/
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NDSI              Int.
    !      NML_OUTPUT_DATE   Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLSHEL 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_output_date_t), INTENT(INOUT) :: NML_OUTPUT_DATE
 
    ! locals
    INTEGER                                :: IERR
    TYPE(nml_output_date_t) :: DATE
    namelist /output_date_nml/ date
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
    ierr = 0
#ifdef W3_S
    CALL strace (ient, 'READ_OUTPUT_DATE_NML')
#endif
 
    ! set default values for output_date structure
    date%FIELD%START = '19680606 000000'
    date%FIELD%STRIDE = '0'
    date%FIELD%STOP = '19680607 000000'
    date%FIELD%OUTFFILE = '0'
    date%POINT%OUTFFILE = '0'
    date%POINT%START = '19680606 000000'
    date%POINT%STRIDE = '0'
    date%POINT%STOP = '19680607 000000'
    date%TRACK%START = '19680606 000000'
    date%TRACK%STRIDE = '0'
    date%TRACK%STOP = '19680607 000000'
    date%RESTART%START = '19680606 000000'
    date%RESTART%STRIDE = '0'
    date%RESTART%STOP = '19680607 000000'
    date%RESTART2%START = '19680606 000000'
    date%RESTART2%STRIDE = '0'
    date%RESTART2%STOP = '19680607 000000'
    date%BOUNDARY%START = '19680606 000000'
    date%BOUNDARY%STRIDE = '0'
    date%BOUNDARY%STOP = '19680607 000000'
    date%PARTITION%START = '19680606 000000'
    date%PARTITION%STRIDE = '0'
    date%PARTITION%STOP = '19680607 000000'
    date%COUPLING%START = '19680606 000000'
    date%COUPLING%STRIDE = '0'
    date%COUPLING%STOP = '19680607 000000'
 
 
    ! 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 (5)
    END IF
 
    ! save namelist
    nml_output_date = date
 
  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 :         22-Mar-2021 |
    !/                  +-----------------------------------+
    !/
    !/
    !  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
    !     ----------------------------------------------------------------
    !      W3NMLSHEL 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_IC1   = 0
    homog_count%N_IC2   = 0
    homog_count%N_IC3   = 0
    homog_count%N_IC4   = 0
    homog_count%N_IC5   = 0
    homog_count%N_MDN   = 0
    homog_count%N_MTH   = 0
    homog_count%N_MVS   = 0
    homog_count%N_LEV   = 0
    homog_count%N_CUR   = 0
    homog_count%N_WND   = 0
    homog_count%N_ICE   = 0
    homog_count%N_TAU   = 0
    homog_count%N_RHO   = 0
    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 (6)
    END IF
 
    ! allocate the total count of homogeneous input
    homog_count%N_TOT = homog_count%N_IC1 + homog_count%N_IC2 + homog_count%N_IC3 + homog_count%N_IC4 + homog_count%N_IC5 + &
         homog_count%N_MDN + homog_count%N_MTH + homog_count%N_MVS + homog_count%N_LEV + homog_count%N_CUR + &
         homog_count%N_WND + homog_count%N_ICE + homog_count%N_TAU + homog_count%N_RHO + 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 (7)
    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 :         18-Apr-2018 |
    !/                  +-----------------------------------+
    !/
    !/
    !  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
    !     ----------------------------------------------------------------
    !      W3NMLSHEL 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),'IOSTYP = ', nml_domain%IOSTYP
    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_nml (NML_INPUT)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         22-Mar-2021 |
    !/                  +-----------------------------------+
    !/
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NML_INPUT         Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLSHEL 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_input_t), INTENT(IN) :: NML_INPUT
 
    ! locals
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
#ifdef W3_S
    CALL strace (ient, 'REPORT_INPUT_NML')
#endif
 
    WRITE (msg,'(A)') 'INPUT GRID % :'
    WRITE (ndsn,'(A)')
    WRITE (ndsn,10) trim(msg),'FORCING % WATER_LEVELS   = ', nml_input%FORCING%WATER_LEVELS
    WRITE (ndsn,10) trim(msg),'FORCING % CURRENTS       = ', nml_input%FORCING%CURRENTS
    WRITE (ndsn,10) trim(msg),'FORCING % WINDS          = ', nml_input%FORCING%WINDS
    WRITE (ndsn,10) trim(msg),'FORCING % ATM_MOMENTUM   = ', nml_input%FORCING%ATM_MOMENTUM
    WRITE (ndsn,10) trim(msg),'FORCING % AIR_DENSITY    = ', nml_input%FORCING%AIR_DENSITY
    WRITE (ndsn,10) trim(msg),'FORCING % ICE_CONC       = ', nml_input%FORCING%ICE_CONC
    WRITE (ndsn,10) trim(msg),'FORCING % ICE_PARAM1     = ', nml_input%FORCING%ICE_PARAM1
    WRITE (ndsn,10) trim(msg),'FORCING % ICE_PARAM2     = ', nml_input%FORCING%ICE_PARAM2
    WRITE (ndsn,10) trim(msg),'FORCING % ICE_PARAM3     = ', nml_input%FORCING%ICE_PARAM3
    WRITE (ndsn,10) trim(msg),'FORCING % ICE_PARAM4     = ', nml_input%FORCING%ICE_PARAM4
    WRITE (ndsn,10) trim(msg),'FORCING % ICE_PARAM5     = ', nml_input%FORCING%ICE_PARAM5
    WRITE (ndsn,10) trim(msg),'FORCING % MUD_DENSITY    = ', nml_input%FORCING%MUD_DENSITY
    WRITE (ndsn,10) trim(msg),'FORCING % MUD_THICKNESS  = ', nml_input%FORCING%MUD_THICKNESS
    WRITE (ndsn,10) trim(msg),'FORCING % MUD_VISCOSITY  = ', nml_input%FORCING%MUD_VISCOSITY
    WRITE (ndsn,10) trim(msg),'ASSIM % MEAN             = ', nml_input%ASSIM%MEAN
    WRITE (ndsn,10) trim(msg),'ASSIM % SPEC1D           = ', nml_input%ASSIM%SPEC1D
    WRITE (ndsn,10) trim(msg),'ASSIM % SPEC2D           = ', nml_input%ASSIM%SPEC2D
 
 
10  FORMAT (a,2x,a,a)
13  FORMAT (a,2x,a,l1)
 
  END SUBROUTINE report_input_nml
 
  !/ ------------------------------------------------------------------- /
 
 
 
 
 
 
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE report_output_type_nml (NML_OUTPUT_TYPE)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         25-Sep-2020 |
    !/                  +-----------------------------------+
    !/
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NML_OUTPUT_TYPE   Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLSHEL 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_output_type_t), INTENT(IN) :: NML_OUTPUT_TYPE
 
    ! locals
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
#ifdef W3_S
    CALL strace (ient, 'REPORT_OUTPUT_TYPE_NML')
#endif
 
    WRITE (msg,'(A)') 'OUTPUT TYPE % '
    WRITE (ndsn,'(A)')
    WRITE (ndsn,10) trim(msg),'FIELD % LIST         = ', trim(nml_output_type%FIELD%LIST)
    WRITE (ndsn,10) trim(msg),'POINT % FILE         = ', trim(nml_output_type%POINT%FILE)
    WRITE (ndsn,13) trim(msg),'TRACK % FORMAT       = ', nml_output_type%TRACK%FORMAT
    WRITE (ndsn,11) trim(msg),'PARTITION % X0       = ', nml_output_type%PARTITION%X0
    WRITE (ndsn,11) trim(msg),'PARTITION % XN       = ', nml_output_type%PARTITION%XN
    WRITE (ndsn,11) trim(msg),'PARTITION % NX       = ', nml_output_type%PARTITION%NX
    WRITE (ndsn,11) trim(msg),'PARTITION % Y0       = ', nml_output_type%PARTITION%Y0
    WRITE (ndsn,11) trim(msg),'PARTITION % YN       = ', nml_output_type%PARTITION%YN
    WRITE (ndsn,11) trim(msg),'PARTITION % NY       = ', nml_output_type%PARTITION%NY
    WRITE (ndsn,13) trim(msg),'PARTITION % FORMAT   = ', nml_output_type%PARTITION%FORMAT
#ifdef W3_COU
    WRITE (ndsn,10) trim(msg),'COUPLING % SENT         = ', trim(nml_output_type%COUPLING%SENT)
    WRITE (ndsn,10) trim(msg),'COUPLING % RECEIVED     = ', trim(nml_output_type%COUPLING%RECEIVED)
    WRITE (ndsn,13) trim(msg),'COUPLING % COUPLET0     = ', nml_output_type%COUPLING%COUPLET0
#endif
    WRITE (ndsn,10) trim(msg),'RESTART % EXTRA      = ', trim(nml_output_type%RESTART%EXTRA)
 
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 (NML_OUTPUT_DATE)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           M. Accensi              |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         18-Apr-2018 |
    !/                  +-----------------------------------+
    !/
    !/
    !  1. Purpose :
    !
    !
    !  2. Method :
    !
    !     See source term routines.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !      NML_OUTPUT_DATE   Type.
    !     ----------------------------------------------------------------
    !
    !  4. Subroutines used :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      STRACE    Subr. W3SERVMD SUBROUTINE tracing.
    !     ----------------------------------------------------------------
    !
    !  5. Called by :
    !
    !      Name      Type  Module   Description
    !     ----------------------------------------------------------------
    !      W3NMLSHEL 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_output_date_t), INTENT(IN) :: NML_OUTPUT_DATE
 
    ! locals
#ifdef W3_S
    INTEGER, SAVE                           :: IENT = 0
#endif
 
#ifdef W3_S
    CALL strace (ient, 'REPORT_OUTPUT_DATE_NML')
#endif
 
    WRITE (msg,'(A)') 'OUTPUT DATE MODEL GRID % '
    WRITE (ndsn,'(A)')
    WRITE (ndsn,10) trim(msg),'FIELD % START        = ', trim(nml_output_date%FIELD%START)
    WRITE (ndsn,10) trim(msg),'FIELD % STRIDE       = ', trim(nml_output_date%FIELD%STRIDE)
    WRITE (ndsn,10) trim(msg),'FIELD % STOP         = ', trim(nml_output_date%FIELD%STOP)
    WRITE (ndsn,10) trim(msg),'POINT % START        = ', trim(nml_output_date%POINT%START)
    WRITE (ndsn,10) trim(msg),'POINT % STRIDE       = ', trim(nml_output_date%POINT%STRIDE)
    WRITE (ndsn,10) trim(msg),'POINT % STOP         = ', trim(nml_output_date%POINT%STOP)
    WRITE (ndsn,10) trim(msg),'TRACK % START        = ', trim(nml_output_date%TRACK%START)
    WRITE (ndsn,10) trim(msg),'TRACK % STRIDE       = ', trim(nml_output_date%TRACK%STRIDE)
    WRITE (ndsn,10) trim(msg),'TRACK % STOP         = ', trim(nml_output_date%TRACK%STOP)
    WRITE (ndsn,10) trim(msg),'RESTART % START      = ', trim(nml_output_date%RESTART%START)
    WRITE (ndsn,10) trim(msg),'RESTART % STRIDE     = ', trim(nml_output_date%RESTART%STRIDE)
    WRITE (ndsn,10) trim(msg),'RESTART % STOP       = ', trim(nml_output_date%RESTART%STOP)
    WRITE (ndsn,10) trim(msg),'RESTART2 % START      = ', trim(nml_output_date%RESTART2%START)
    WRITE (ndsn,10) trim(msg),'RESTART2 % STRIDE     = ', trim(nml_output_date%RESTART2%STRIDE)
    WRITE (ndsn,10) trim(msg),'RESTART2 % STOP       = ', trim(nml_output_date%RESTART2%STOP)
    WRITE (ndsn,10) trim(msg),'BOUNDARY % START     = ', trim(nml_output_date%BOUNDARY%START)
    WRITE (ndsn,10) trim(msg),'BOUNDARY % STRIDE    = ', trim(nml_output_date%BOUNDARY%STRIDE)
    WRITE (ndsn,10) trim(msg),'BOUNDARY % STOP      = ', trim(nml_output_date%BOUNDARY%STOP)
    WRITE (ndsn,10) trim(msg),'PARTITION % START    = ', trim(nml_output_date%PARTITION%START)
    WRITE (ndsn,10) trim(msg),'PARTITION % STRIDE   = ', trim(nml_output_date%PARTITION%STRIDE)
    WRITE (ndsn,10) trim(msg),'PARTITION % STOP     = ', trim(nml_output_date%PARTITION%STOP)
#ifdef W3_COU
    WRITE (ndsn,10) trim(msg),'COUPLING % START    = ', trim(nml_output_date%COUPLING%START)
    WRITE (ndsn,10) trim(msg),'COUPLING % STRIDE   = ', trim(nml_output_date%COUPLING%STRIDE)
    WRITE (ndsn,10) trim(msg),'COUPLING % STOP     = ', trim(nml_output_date%COUPLING%STOP)
#endif
 
 
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 :         22-Mar-2021 |
    !/                  +-----------------------------------+
    !/
    !/
    !  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
    !     ----------------------------------------------------------------
    !      W3NMLSHEL 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_IC1       = ', nml_homog_count%N_IC1
    WRITE (ndsn,11) trim(msg),'N_IC2       = ', nml_homog_count%N_IC2
    WRITE (ndsn,11) trim(msg),'N_IC3       = ', nml_homog_count%N_IC3
    WRITE (ndsn,11) trim(msg),'N_IC4       = ', nml_homog_count%N_IC4
    WRITE (ndsn,11) trim(msg),'N_IC5       = ', nml_homog_count%N_IC5
    WRITE (ndsn,11) trim(msg),'N_MDN       = ', nml_homog_count%N_MDN
    WRITE (ndsn,11) trim(msg),'N_MTH       = ', nml_homog_count%N_MTH
    WRITE (ndsn,11) trim(msg),'N_MVS       = ', nml_homog_count%N_MVS
    WRITE (ndsn,11) trim(msg),'N_LEV       = ', nml_homog_count%N_LEV
    WRITE (ndsn,11) trim(msg),'N_CUR       = ', nml_homog_count%N_CUR
    WRITE (ndsn,11) trim(msg),'N_WND       = ', nml_homog_count%N_WND
    WRITE (ndsn,11) trim(msg),'N_ICE       = ', nml_homog_count%N_ICE
    WRITE (ndsn,11) trim(msg),'N_TAU       = ', nml_homog_count%N_TAU
    WRITE (ndsn,11) trim(msg),'N_RHO       = ', nml_homog_count%N_RHO
    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 w3nmlshelmd
 
!/ ------------------------------------------------------------------- /