w3ounfmetamd.F90

Go to the documentation of this file.

 
!/ ------------------------------------------------------------------- /
MODULE w3ounfmetamd
  !/
  !/    02-Nov-2020 : Creation                            ( version 7.12 )
  !/    26-Jan-2021 : Added TP and alternative dir/mag    ( version 7.12 )
  !/                  metadata for directional fields.
  !/    16-Dec-2020 : Added user partition templates      ( version 7.12 )
  !/                  and coordinate reference system.
  !/                  Freeform meta data uses linked list.
  !/    02-Feb-2021 : Improved partitioned parameter      ( version 7.12 )
  !/                  template string implementation.
  !/    22-Mar-2021 : Adds extra coupling fields          ( version 7.13 )
  !/    02-Sep-2021 : Add coordinates attribute           ( version 7.12 )
  !/
  !/ ------------------------------------------------------------------- /
  !/
  USE netcdf
  USE constants, ONLY: tpiinv
  USE w3gdatmd, ONLY: sig, nk, gtype, ungtype
#ifdef W3_RTD
  USE w3gdatmd, ONLY : flagunr, polat, polon
#endif
#ifdef W3_SMC
  USE w3smcomd, ONLY : smcotype
#endif
  USE w3odatmd, ONLY: ptmeth, ptfcut, nogrp, noge, ngrpp,           &
       ndse, fnmpre, noswll
  USE w3servmd, ONLY: extcde, str_to_upper
 
  USE w3metamd
 
  IMPLICIT NONE
 
  PUBLIC
 
  LOGICAL, PRIVATE :: DEBUG = .false. 
 
  CHARACTER(LEN=*), PARAMETER :: fn_meta = "ounfmeta.inp"
 
  CHARACTER(LEN=*), PARAMETER :: ipart_token  = "<IPART>"
 
  CHARACTER(LEN=*), PARAMETER :: spart_token  = "<SPART>"
 
  CHARACTER(LEN=*), PARAMETER :: spart_token_ = "<SPART_>"
 
  TYPE meta_t
    REAL :: fsc
    REAL :: vmin
    REAL :: vmax = unsetr                   
    CHARACTER(LEN=24)  :: units = unsetc    
    CHARACTER(LEN=50)  :: ename = unsetc    
    CHARACTER(LEN=80)  :: varnm = unsetc, & 
         varnl = unsetc    
    CHARACTER(LEN=120) :: varns = unsetc, & 
         varng = unsetc, & 
         varnd = unsetc    
    CHARACTER(LEN=512) :: varnc = unsetc    
    TYPE(meta_list_t) :: extra
 
    ! For updating meta only:
    INTEGER :: ifi = 0, &           !< Group index to update
         ifj = 0, &           
         ifc = 1              
    CHARACTER(LEN=6) :: fldid = ''
  ENDTYPE meta_t
 
  TYPE field_t
    TYPE(meta_t), POINTER :: meta(:)
  END TYPE field_t
 
  TYPE group_t
    TYPE(field_t), ALLOCATABLE :: field(:)
  END TYPE group_t
 
  TYPE(group_t), ALLOCATABLE :: group(:)
 
  TYPE(meta_list_t) :: global_meta
 
  LOGICAL           :: fl_default_gbl_meta = .true.
 
  ! Storage for coordinate reference system (CRS) metadata:
  CHARACTER(LEN=128) :: crs_name = ''
  TYPE(meta_list_t)  :: crs_meta
  LOGICAL            :: crs_is_default = .false. 
 
  CHARACTER(LEN=256) :: coords_attr = ''
 
  TYPE part_tmpl_t
    CHARACTER(LEN=128)              :: tmpl
    CHARACTER(LEN=128), ALLOCATABLE :: part_text(:)
    INTEGER(KIND=2)                 :: np
    TYPE(part_tmpl_t), POINTER      :: next
  END TYPE part_tmpl_t
 
  TYPE(part_tmpl_t), POINTER :: part_tmpl
 
  INTEGER            :: ncvartype
  CHARACTER(LEN=30)  :: dircom
  CHARACTER(LEN=128) :: partcom
  CHARACTER(LEN=15)  :: snamep(5)
 
  LOGICAL, PRIVATE   :: vector
  LOGICAL            :: flrtd = .false. 
 
CONTAINS
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE init_meta(VEC)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         22-Mar-2021 |
    !/                  +-----------------------------------+
    !/
    !/    09-Nov-2020 : Creation                            ( version 7.12 )
    !/    22-Mar-2021 : Added vector flag                   ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Allocates space for the META_T arrays and sets some constants.
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
 
    LOGICAL, INTENT(IN), OPTIONAL  :: VEC
    !/
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    LOGICAL :: FLGNML
    INTEGER :: I, J
 
    vector = .true.
    IF(PRESENT(vec)) vector = vec
#ifdef W3_RTD
    ! Is the grid really rotated?
    IF ( polat < 90. ) flrtd = .true.
#endif
#if defined W3_SMC && defined W3_RTD
    ! SMC type 3/4 outputs are currently on standard pole grid only
    IF(smcotype .EQ. 3 .OR. smcotype .EQ. 4) flrtd = .false.
#endif
 
    ! 1. Allocate nested GROUP, FIELD structure:
    ALLOCATE(group(nogrp))
    DO i = 1,nogrp
      ALLOCATE(group(i)%FIELD(noge(i)))
      DO j = 1,noge(i)
        ALLOCATE(group(i)%FIELD(j)%META(3)) ! Hardcode to 3 components for the moment
      ENDDO
    ENDDO
 
    ! 1.1 Make sure partitioned template pointer is null (i.e. empty list)
    NULLIFY(part_tmpl)
 
    ! 2. Set direction convention:
    dircom = ""
#ifdef W3_RTD
    IF( flrtd ) THEN
      IF ( flagunr ) THEN
        dircom = 'True North'
      ELSE IF ( .NOT. flagunr ) THEN
        dircom = 'Rotated Pole Grid North'
      ENDIF
    ENDIF
#endif
 
    !  Set partitioning method comment and standard name templates:
    IF( ptmeth .LE. 3 ) THEN
      snamep(1) = 'wind'
      snamep(2) = 'primary swell'
      snamep(3) = 'secondary swell'
      snamep(4) = 'tertiary swell'
      snamep(5) = 'swell'
    ELSE
      snamep(1) = 'wind'
      snamep(2) = 'swell'
      snamep(3:5) = ''
    ENDIF
 
    IF ( ptmeth .EQ. 1 ) THEN
      partcom = "Wind sea and swells defined using topographic " //   &
           "partitions and partition wave-age cut-off "     //   &
           "(WWIII default scheme)"
    ELSE IF ( ptmeth .EQ. 2 ) THEN
      partcom = "Wind sea and swells defined using topographic " //   &
           "partitions and spectral wave-age cut-off"
    ELSE IF ( ptmeth .EQ. 3 ) THEN
      partcom = "Wave components defined using topographic "     //   &
           "partitions only"
    ELSE IF ( ptmeth .EQ. 4 ) THEN
      partcom = "Wind sea and swell defined using spectral "     //   &
           "wave-age cut-off"
    ELSE IF ( ptmeth .EQ. 5 ) THEN
      WRITE(partcom, '("Wave components defined using ",F5.3,'   //   &
           '"Hz spectral frequency cutoff")') ptfcut
    ELSE
      WRITE(partcom, '("PTM_",I1,"_Unknown")') ptmeth
    ENDIF
 
    ! 3. Set the default values for the OUNF netCDF meta data.
    CALL default_meta()
 
    ! Set the default coordiante reference system (if applicable)
    CALL default_crs_meta()
 
    ! If the ounfmeta.inp exists, read this in to override defaults:
    INQUIRE(file=trim(fnmpre)//"ounfmeta.inp", exist=flgnml)
    IF(flgnml) THEN
      CALL read_meta()
    ENDIF
 
  END SUBROUTINE init_meta
  !
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE teardown_meta()
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         09-Nov-2020 |
    !/                  +-----------------------------------+
    !/
    !/    09-Nov-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     De-allocates memory used for the META_T arrays
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
    !/
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    INTEGER :: I, J
 
    DO i = 1,nogrp
      DO j = 1,noge(i)
        DEALLOCATE(group(i)%FIELD(j)%META)
      ENDDO
      DEALLOCATE(group(i)%FIELD)
    ENDDO
    DEALLOCATE(group)
 
    CALL del_meta_list(global_meta)
    CALL del_meta_list(crs_meta)
 
  END SUBROUTINE teardown_meta
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
 
  SUBROUTINE next_line(NDMI, BUF, ILINE, EOF, NEW_SECTION)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         09-Nov-2020 |
    !/                  +-----------------------------------+
    !/
    !/    09-Nov-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Reads the next valid line from the user meta input file.
    !
    !  2. Method :
    !     Lines are repeatedly read in from the input file until
    !     a valid input line is reached. Blank lines and comment lines
    !     (starting with $) are skipped.
    !
    !     If the end of file is reached before any valid line is read
    !     then EOF is set to true.
    !
    !     If the next valid line is a new section marker (META or TEMPLATE)
    !     then the NEW_SECTION flag is set to true.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       NDMI    Int.  I  Unit number of input file
    !       BUF    Char.  O  Next input line read from file
    !       ILINE   Int. I/O Line number of file
    !       EOF    Bool.  O  True if end-of-file is reached.
    !       NEW_SECTION
    !              Bool.  O  True if new section marker found
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
 
    INTEGER, INTENT(IN)       :: NDMI
    CHARACTER(*), INTENT(OUT) :: BUF
    INTEGER, INTENT(INOUT)    :: ILINE
    LOGICAL, INTENT(OUT)      :: EOF
    LOGICAL, INTENT(OUT), OPTIONAL :: NEW_SECTION
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    INTEGER :: IERR
    CHARACTER(LEN=10) :: TEST
 
    eof = .false.
 
    ! Keep reading from file until we read a line that is not:
    !   - a blank line
    !   - a comment line (starting with $)
    !   - the end of the file
    DO
      READ(ndmi, '(A)', iostat=ierr, err=101, end=100) buf
 
      iline = iline + 1
 
      ! Remove any tab characters from buffer (replace with a space)
      CALL notabs(buf)
 
      ! Empty line?
      IF(trim(buf) == '') THEN
        IF(debug) WRITE(*,'(I5,1X,A20)') iline, '[blank line]'
        cycle
      ENDIF
 
      IF(trim(buf) == "$ DEBUG ON") THEN
        WRITE(*,'(I5,1X,A20)') iline, '[DEBUG ON]'
        debug = .true.
        cycle
      ENDIF
 
      IF(trim(buf) == "$ DEBUG OFF") THEN
        WRITE(*,'(I5,1X,A20)') iline, '[DEBUG OFF]'
        debug = .false.
        cycle
      ENDIF
 
      ! Read first token on line:
      READ(buf, *) test
 
      ! Check for comment:
      IF(test(1:1) == "$" .OR. trim(buf) == '') THEN
        IF(debug) WRITE(*,'(I5,1X,A20)') iline, '[comment line]'
        cycle
      ENDIF
 
      ! Check if is section header
      IF(PRESENT(new_section)) THEN
        CALL str_to_upper(test)
        SELECT CASE(test)
        CASE ("META", "TEMPLATE", "CRS")
          new_section = .true.
        CASE DEFAULT
          new_section = .false.
        END SELECT
      ENDIF
 
      ! Anything else can be considered the "next line"
      RETURN
    ENDDO
 
    !/    Escape locations
    !
    !     End of file
100 CONTINUE
    buf = ''
    eof = .true.
    RETURN
    !
    !     I/O Error
101 CONTINUE
    WRITE(ndse, 1000) fn_meta, iline, ierr
    CALL extcde(10)
    !
1000 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '/           &
         '     ERROR READING METADATA FILE'/                    &
         '     FILENAME = ', a /                                &
         '     LINE NO = ', i5 /                                &
         '     IOSTAT =',i5 /)
    !
  END SUBROUTINE next_line
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE notabs(STR)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         02-Nov-2020 |
    !/                  +-----------------------------------+
    !/
    !/    02-Nov-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Replaces tab characters in a string with a space.
    !
    !  2. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       STR    Char.  I/O  Character string to process
    !     ----------------------------------------------------------------
    !
    !  3. Remarks :
    !
    !     Assumes ASCII encoding! Tab character is ASCII value 9.
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
    CHARACTER(*), INTENT(INOUT) :: STR
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    INTEGER, PARAMETER :: ASCII_TAB = 9
    INTEGER :: SLEN
    INTEGER :: I
    !
    slen = len_trim(str)
 
    DO i=1,slen
      IF(ichar(str(i:i)) == ascii_tab) THEN
        str(i:i) = ' '
      ENDIF
    ENDDO
 
  END SUBROUTINE notabs
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  FUNCTION replace_char(STR, C, REP) RESULT(OSTR)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         02-Fev-2021 |
    !/                  +-----------------------------------+
    !/
    !/    02-Feb-2021 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Replaces single characters in a string. Returns a new string,
    !
    !  2. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       STR    CharArr  I  Character string to process
    !       C      Char     I  Character to search for
    !       REP    Char     I  Character to substitute
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
    CHARACTER(*)        :: str
    CHARACTER           :: c, rep
    CHARACTER(LEN(STR)) :: ostr
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    INTEGER :: i
 
    ostr = trim(str)
    DO
      i = index(trim(ostr), c)
      IF(i .LE. 0) EXIT
      ostr(i:i) = rep
    ENDDO
 
  END FUNCTION replace_char
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE read_meta()
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         26-Jan-2021 |
    !/                  +-----------------------------------+
    !/
    !/    09-Nov-2020 : Creation                            ( version 7.12 )
    !/    26-Jan-2021 : Added TP and alternative dir/mag    ( version 7.12 )
    !/                  metadata for directional fields.
    !/
    !
    !  1. Purpose :
    !
    !     Reads meta data entries from the ountmeta.inp file and update
    !     default values set via the DEFAULT_META subroutine.
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    INTEGER ::  IERR, I, NDMI
    CHARACTER(LEN=256) :: BUF
    TYPE(meta_t), POINTER :: PMETA
    CHARACTER(LEN=32) :: TEST, TESTU
 
    INTEGER :: IFI, IFJ, IFC
    INTEGER :: ILINE, MCNT
    LOGICAL :: EOF
 
    ndmi = 60
    iline = 0
    mcnt = 0
 
    OPEN(unit=ndmi, file=trim(fnmpre)//trim(fn_meta),                &
         status="OLD", iostat=ierr)
 
    IF(ierr .NE. 0) THEN
      WRITE(ndse, 5010) trim(fnmpre)//trim(fn_meta), ierr
      CALL extcde(10)
    ENDIF
 
    ! Loop over file, skipping comments or blank lines, until we find
    ! a META line.
    DO
      CALL next_line(ndmi, buf, iline, eof)
      IF(eof) EXIT
 
      ! Read first token on line:
      READ(buf, *) test
 
      ! A new meta-data section will start with the keyword "META"
      testu = test
      CALL str_to_upper(testu)
      IF(testu == "META") THEN
        mcnt = mcnt + 1
 
        IF(debug) WRITE(*,'(I5,1X,A20,1X,A)') iline, '[META header]', trim(buf)
 
        ! Get the IFI, IFJ, IFC values from the header:
        i = index(buf, trim(test)) + 4 ! Handles lower/mixed-case META keyword
        CALL decode_header(buf(i:), iline, ifi, ifj, ifc)
        IF(ifi .EQ. -1) THEN
          WRITE(ndse, 5011) trim(buf(i:)), trim(fn_meta), iline
          CALL extcde(10)
        ENDIF
 
        ! IFI = 999 is a section for the "global" meta data
        IF(ifi .EQ. 999) THEN
          CALL read_freeform_meta_list(ndmi, iline, global_meta)
          cycle
        ENDIF
 
        ! Error checking on size of IFJ, ICOMP, IPART.
        IF(ifi .LT. 1 .OR. ifi .GT. nogrp) THEN
          WRITE(ndse,5013) nogrp, trim(fn_meta), iline
          CALL extcde(1)
        ENDIF
        IF(ifj .LT. 1 .OR. ifj .GT. noge(ifi)) THEN
          WRITE(ndse,5014) noge(ifi), trim(fn_meta), iline
          CALL extcde(1)
        ENDIF
        IF(ifc .LT. 1 .OR. ifc .GT. 3) THEN
          WRITE(ndse,5015) trim(fn_meta), iline
          CALL extcde(1)
        ENDIF
 
        ! Select correct variable metadata entry:
        pmeta => group(ifi)%FIELD(ifj)%META(ifc)
 
        ! Update the metadata with values from file:
        CALL read_meta_pairs(ndmi, pmeta, iline)
 
      ELSE IF(testu == "TEMPLATE") THEN
        backspace(ndmi) ! We will reprocess this line
        CALL read_part_tmpl(ndmi, iline)
        cycle
 
      ELSE IF(testu == "CRS") THEN
        backspace(ndmi) ! We will reprocess this line
        CALL read_crs_meta(ndmi, iline)
        cycle
 
      ELSE
        ! Anything else is a syntax error
        WRITE(ndse, 5012) trim(fn_meta), iline, trim(buf)
        CALL extcde(10)
      ENDIF
    ENDDO
 
    CLOSE(ndmi)
    !WRITE(*, 5000) MCNT, N_GBLMETA, N_CRSMETA
    RETURN
    !
5000 FORMAT(/' Read in: ',i3,' variable metadata entries' /          &
         '     and: ',i3,' global meta data entries' /           &
         '     and: ',i3,' CRS meta data entries' /)
    !
5010 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '/           &
         '     ERROR OPENING METADATA FILE'/                    &
         '     FILENAME = ', a /                                &
         '     IOSTAT =', i5 /)
    !
5011 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '/           &
         '     UNKNOWN FIELD ID: ',a /                          &
         '     FILENAME = ', a /                                &
         '     LINE =', i5 /)
    !
5012 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '/           &
         '     SYNTAX ERROR ' /                                 &
         '     FILENAME = ', a /                                &
         '     LINE =', i5 /                                    &
         '  => ', a /)
    !
5013 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : ' /          &
         '     IFI value should be in range 1,',i2 /            &
         '     FILENAME = ', a /                                &
         '     LINE =', i5 /                                    &
         '  => ', a /)
    !
5014 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : ' /          &
         '     IFJ value should be in range 1,',i2 /            &
         '     FILENAME = ', a /                                &
         '     LINE =', i5 /                                    &
         '  => ', a /)
    !
5015 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : ' /          &
         '     IFC value should be in range 1,3' /              &
         '     FILENAME = ', a /                                &
         '     LINE =', i5 /                                    &
         '  => ', a /)
    !
  END SUBROUTINE read_meta
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE decode_header(BUF, ILINE, IFI, IFJ, IFC)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         02-Feb-2021 |
    !/                  +-----------------------------------+
    !/
    !/    09-Nov-2020 : Creation                            ( version 7.12 )
    !/    02-Feb-2021 : NODEFAULT option for Global meta    ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Decode the META header line.
    !
    !  2. Method:
    !
    !     The internal WW3 field can be specified either as an [IFI, IFJ]
    !     integer combination, or a field ID tag (such as "HS").
    !
    !     Both forms can also specify an optional component (IFC) integer
    !     value for multi-component fields (defaults to 1).
    !
    !     Field name ID tags are case-insensitive, HS == hs == Hs.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       BUF    Char.  I  Input header string (without leading META tag)
    !       ILINE   Int.  I  Line number (for error reporting)
    !       IFI     Int.  O  Output group number
    !       IFJ     Int.  O  Output field number
    !       IFC     Int.  O  Component number (defaults to 1)
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    USE w3iogomd, ONLY: w3fldtoij
 
    IMPLICIT NONE
 
    CHARACTER(*), INTENT(IN) :: BUF
    INTEGER, INTENT(IN) :: ILINE
    INTEGER, INTENT(OUT) :: IFI, IFJ, IFC
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    INTEGER :: IERR, I
    CHARACTER(LEN=10) :: FLD, OFLD, OPT
 
    ifi = 0
    ifj = 1
    ifc = 1
    fld = ''
 
    ! Is first value an int?
    READ(buf, *, iostat=ierr) ifi
    IF(ierr .EQ. 0) THEN
      ! Try reading 3 values:
      READ(buf, *, iostat=ierr) ifi, ifj, ifc
      IF(ierr .NE. 0) THEN
        ! Try just two values:
        READ(buf, *, iostat=ierr) ifi, ifj
      ENDIF
    ELSE
      ! Try reading field ID plus component
      READ(buf, *, iostat=ierr) fld, ifc
 
      IF(ierr .NE. 0) THEN
        ! Try just fldid
        READ(buf, *, iostat=ierr) fld
      ENDIF
    ENDIF
 
    IF(ierr .NE. 0) THEN
      WRITE(ndse, 6000) trim(fn_meta), iline, trim(buf)
      CALL extcde(10)
    ENDIf
 
    ofld = fld
    CALL str_to_upper(fld) ! field names are case insensitive
 
    ! If string value (FLDID), then decode into IFI,IFJ value:
    IF(fld /= '') THEN
      ! Special case for "global" attributes
      IF(trim(fld) == "GLOBAL") THEN
        IF(debug) WRITE(*,'(6X,A20,1X,A)') '[GLOBAL meta sec.]', trim(buf)
        ifi = 999  ! Marker for global section
 
        ! check for any options:
        i = index(buf, trim(ofld)) + len_trim(ofld)
        opt = adjustl(buf(i:))
        CALL str_to_upper(opt)
        SELECT CASE(trim(opt))
        CASE("")
          CONTINUE ! no option
        CASE("NODEFAULT")
          fl_default_gbl_meta = .false.
          IF(debug) WRITE(*,'(6X,A20,1X,A)') '[GLOBAL meta]', 'Defaults disabled'
        CASE DEFAULT
          WRITE(ndse, *) "Unknown GLOBAL extra option: [", trim(opt), "]"
        END SELECT
      ELSE
        IF(debug) WRITE(*,'(6X,A20,1X,A)') '[Decoding field ID]', trim(buf)
        CALL w3fldtoij(fld, ifi, ifj, 1, 1, 1)
      ENDIF
    ENDIF
 
    IF(debug) WRITE(*,'(6X,A20,1X,I4,2I2)') '[IFI, IFJ, IFC]', ifi,ifj,ifc
    !
6000 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '/           &
         '     SYNTAX ERROR IN SECTION HEADER. ' /              &
         '     FILENAME = ', a /                                &
         '     LINE NO =', i5 /                                 &
         '  => ', a /)
    !
  END SUBROUTINE decode_header
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE read_meta_pairs(NDMI, META, ILINE)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         09-Nov-2020 |
    !/                  +-----------------------------------+
    !/
    !/    09-Nov-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Reads in attribute name/value pairs and updates the relevant
    !     values in the META type.
    !
    !  2. Method:
    !
    !     Keeps looping over input lines in file until next META section
    !     or EOF is found. Splits meta pairs on the = character.
    !
    !     Note - the "extra" metadata pair can also provide a variable
    !     type ("c", "i", or "r"; for character, int or real respectively)
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       NDMI      Int.   I  Unit number of metadata input file
    !       META  Int.Ptr.   O  Pointer to META type
    !       ILINE     Int. I/O  Current line number in file
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
    INTEGER, INTENT(IN)                  :: NDMI
    TYPE(meta_t), INTENT(INOUT), POINTER :: META
    INTEGER, INTENT(INOUT)               :: ILINE
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !
    CHARACTER(LEN=256) :: BUF
    CHARACTER(LEN=128) :: ATTN, ATTV, TMP
    CHARACTER(LEN=16)  :: ATT_TYPE!, TEST
    INTEGER            :: I, IERR
    REAL               :: R
    LOGICAL            :: EOF, NEW
    TYPE(meta_pair_t)  :: EXTRA
 
    ! Keep reading lines until we hit EOF or anoter META keyword
    DO
      CALL next_line(ndmi, buf, iline, eof, new_section=new)
      IF(eof) THEN
        backspace(ndmi)
        RETURN
      ENDIF
 
      IF(new) THEN
        IF(debug) WRITE(*,'(I5,1X,A20)') iline, '[--end of section--]'
        iline = iline - 1
        backspace(ndmi)
        EXIT
      ENDIF
 
      IF(debug) WRITE(*,'(I5,1X,A20,1X,A)') iline, '[META pair]', trim(buf)
 
      ! Meta data should be formatted as "attr_name = attr_value"
      i = index(buf, "=")
      IF( i .LT. 1 ) THEN
        WRITE(ndse, 7000) fn_meta, iline, trim(buf)
        CALL extcde(10)
      ENDIF
 
      attn = adjustl(buf(1:i-1))
      attv = adjustl(buf(i+1:))
 
      ! Some compilers won't read an "empty" string unless quoted:
      IF(trim(attv) == '') THEN
        attv='""'
      ENDIF
 
      ierr = 0
      SELECT CASE(trim(attn))
        ! Character variables
        ! Note: Using internal reads will allow the use of quote marks in strings
      CASE("varnm")
        READ(attv, *, iostat=ierr) meta%VARNM
 
      CASE("ename")
        READ(attv, *, iostat=ierr) meta%ENAME
 
      CASE("standard_name", "varns")
        READ(attv, *, iostat=ierr) meta%VARNS
 
      CASE("long_name", "varnl")
        READ(attv, *, iostat=ierr) meta%VARNL
 
      CASE("globwave_name", "varng")
        READ(attv, *, iostat=ierr) meta%VARNG
 
      CASE("direction_reference", "dir_ref", "varnd")
        READ(attv, *, iostat=ierr) meta%VARND
 
      CASE("comment", "varnc")
        READ(attv, *, iostat=ierr) meta%VARNC
 
      CASE("units")
        READ(attv, *, iostat=ierr) meta%UNITS
 
        ! Real variables
      CASE("valid_min", "vmin")
        READ(attv, *, iostat=ierr) meta%VMIN
 
      CASE("valid_max", "vmax")
        READ(attv, *, iostat=ierr) meta%VMAX
 
      CASE("scale_factor", "fsc")
        READ(attv, *, iostat=ierr) meta%FSC
 
        ! Default case will be the "extra" meta data variable
      CASE DEFAULT
        tmp = attv
        CALL get_attval_type(tmp, iline, attv, att_type)
 
        IF(debug) THEN
          WRITE(*,'(I5,1X,A20,1X,A)') iline, '[META extra]', &
               trim(attn)//' = '//trim(attv)//' (type: '//trim(att_type)//")"
        ENDIF
 
        extra%ATTNAME = trim(attn)
        extra%ATTVAL = trim(attv)
        extra%TYPE = trim(att_type)
        CALL meta_list_append(meta%EXTRA, extra)
 
      END SELECT
 
      IF(ierr /= 0) THEN
        WRITE(ndse, 7002) fn_meta, iline, trim(buf)
        CALL extcde(10)
      ENDIF
 
    ENDDO
 
    RETURN
    !
7000 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '/           &
         '     SYNTAX ERROR IN METADATA FILE ' /                &
         '     SHOULD BE "attr_name = attr_value" ' /           &
         '     FILENAME = ', a /                                &
         '     LINE NO =', i5 /                                 &
         '  => ', a /)
    !
7002 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '/           &
         '     IO ERROR READING ATTRIBUTE' /                    &
         '     FILENAME = ', a /                                &
         '     LINE NO =', i5 /                                 &
         '  => ', a /)
    !
  END SUBROUTINE read_meta_pairs
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE get_attval_type(BUF, ILINE, ATTV, ATT_TYPE)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         09-Nov-2020 |
    !/                  +-----------------------------------+
    !/
    !/    09-Nov-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Gets the attribute value and optional variable type from
    !     the passed in string.
    !
    !  2. Method:
    !
    !     If two freeform values can be read from the input string, it is
    !     assumed to be a value and type, otherwise if only one value can
    !     be read the type is assumed to be "character".
    !
    !     It is important to quote strings if they contain spaces.
    !
    !     Valid types are "c" "r/f", and "i" for character, real/float and
    !     integer values.
 
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       BUF      Char.   I  Input string to process
    !       ILINE    Int.    I  Line number (for error reporting)
    !       ATTV     Char.   O  Attribute value
    !       ATT_TYPE Char.   O  Attribute type
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
 
    CHARACTER(*), INTENT(IN)  :: BUF
    INTEGER, INTENT(IN)       :: ILINE
    CHARACTER(*), INTENT(OUT) :: ATTV, ATT_TYPE
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !
    REAL :: R
    INTEGER :: I, IERR
 
    ! Get attribute and type (default to "c" if no type set)
    att_type = 'c'
    attv = ''
    READ(buf, *, iostat=ierr) attv, att_type
    IF(ierr /= 0) READ(buf, *, iostat=ierr) attv
 
    ! Check attr values are valid w.r.t. attr type
    SELECT CASE(trim(att_type))
 
    CASE("i")
      READ(attv, *, iostat=ierr) i
      IF(ierr .ne. 0) then
        WRITE(ndse, 8001) "INTEGER", trim(fn_meta), iline, trim(attv)
        CALL extcde(10)
      ENDIF
 
    CASE("r", "f")
      READ(attv, *, iostat=ierr) r
      IF(ierr .ne. 0) THEN
        WRITE(ndse, 8001) "REAL/FLOAT", trim(fn_meta), iline, trim(attv)
        CALL extcde(10)
      ENDIF
 
    CASE("c")
      ! Always ok.
 
    CASE DEFAULT
      WRITE(ndse, 8002) trim(fn_meta), iline, trim(buf)
      CALL extcde(10)
 
    END SELECT
    !
8001 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '/           &
         '     VALUE IS NOT A VALID ', a /                      &
         '     FILENAME = ', a /                                &
         '     LINE NO =', i5 /                                 &
         '  => ', a /)
    !
8002 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '/           &
         '     ATTRIBUTE TYPE SHOULD BE ONE OF [c,i,r] '/       &
         '     FILENAME = ', a /                                &
         '     LINE NO =', i5 /                                 &
         '  => ', a /)
    !
  END SUBROUTINE get_attval_type
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE read_freeform_meta_list(NDMI, ILINE, METALIST)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         16-Dec-2020 |
    !/                  +-----------------------------------+
    !/
    !/    16-Dec-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Reads in freeform attribute name/value pairs.
    !
    !  2. Method:
    !
    !     Keeps looping over input lines in file until next section
    !     or EOF is found. Splits meta pairs on the = character.
    !
    !     Freeform metadata pairs can also provide a variable type
    !     ("c", "i", or "r"; for character, int or real respectively).
    !     String values with spaces should be quoted.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       NDMI     Char.   I  Unit number of metadata input file
    !       ILINE     Int. I/O  Current line number in file
    !       METALIST Type. I/O  A META_LIST_T object to append to
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
    INTEGER, INTENT(IN)               :: NDMI
    INTEGER, INTENT(INOUT)            :: ILINE
    TYPE(meta_list_t), INTENT(INOUT)  :: METALIST
 
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    CHARACTER(LEN=256) :: BUF
    CHARACTER(LEN=128) :: ATTN, ATTV, TMP
    CHARACTER(LEN=16)  :: ATT_TYPE
    INTEGER            :: I, IERR
    REAL               :: R
    LOGICAL            :: EOF, NEW
    TYPE(meta_pair_t)  :: META
    !
    ! Keep reading lines until we hit EOF or anoter META keyword
    DO
      CALL next_line(ndmi, buf, iline, eof, new_section=new)
      IF(eof) THEN
        backspace(ndmi)
        RETURN
      ENDIF
 
      IF(new) THEN
        IF(debug) WRITE(*,'(I5,1X,A20)') iline, '[--end of section--]'
        iline = iline - 1
        backspace(ndmi)
        EXIT
      ENDIF
 
      ! Split attr name/value pair
      i = index(buf, "=")
      IF( i .LT. 1 ) THEN
        WRITE(ndse, 9000) trim(fn_meta), iline, trim(buf)
        CALL extcde(10)
      ENDIF
 
      attn = adjustl(buf(1:i-1))
      tmp = adjustl(buf(i+1:))
 
      ! Get type, if set:
      CALL get_attval_type(tmp, iline, attv, att_type)
 
      IF(debug) THEN
        WRITE(*,'(I5,1X,A20,1X,A)') iline, '[FREEFORM meta]', &
             trim(attn)//' = '//trim(attv)//' (type: '//trim(att_type)//")"
      ENDIF
 
 
      meta%ATTNAME = trim(attn)
      meta%ATTVAL = trim(attv)
      meta%TYPE = trim(att_type)
 
      CALL meta_list_append(metalist, meta)
 
    ENDDO
    !
9000 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '/           &
         '     SYNTAX ERROR IN METADATA FILE ' /                &
         '     SHOULD BE "attr_name = attr_value" ' /           &
         '     FILENAME = ', a /                                &
         '     LINE NO =', i5 /                                 &
         '  => ', a /)
    !
  END SUBROUTINE read_freeform_meta_list
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE read_crs_meta(NDMI, ILINE)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         07-Dec-2020 |
    !/                  +-----------------------------------+
    !/
    !/    07-Dec-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Reads in metadata for the coordinate reference system (CRS)
    !     scalar variable. The "grid_mapping_name" must be supplied as
    !     an attribute.
    !
    !  2. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       NDMI     Char.   I  Unit number of metadata input file
    !       ILINE     Int. I/O  Current line number in file
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
    INTEGER, INTENT(IN)     :: NDMI
    INTEGER, INTENT(INOUT)  :: ILINE
    !/
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    CHARACTER(LEN=128) :: BUF, PREV_NAME
    INTEGER            :: I, IERR
 
    prev_name = crs_name
 
    ! Re-read header line (we only want the second field)
    READ(ndmi, '(A)') buf
    READ(buf, *, iostat=ierr) crs_name, crs_name
    IF(ierr /= 0 ) THEN
      WRITE(ndse,1000)
      WRITE(ndse,2000) trim(fn_meta), iline, trim(buf)
      CALL extcde(10)
    ENDIF
    IF(debug) WRITE(*,'(I5,1X,A20,1X,A)') iline, '[CRS id]', trim(crs_name)
 
    IF(crs_meta%N .NE. 0) THEN
      IF(crs_is_default) THEN
        WRITE(ndse,1001) trim(prev_name)
        crs_is_default = .false.
      ELSE
        WRITE(ndse,1002) trim(prev_name)
      ENDIF
      WRITE(ndse,2000) trim(fn_meta), iline, trim(buf)
      CALL del_meta_list(crs_meta)
    ENDIF
 
    CALL read_freeform_meta_list(ndmi, iline, crs_meta)
 
    ! Check that "grid_mapping_name" is defined
    IF(.NOT.  meta_list_has_attr(crs_meta, "grid_mapping_name")) THEN
      WRITE(ndse, 1003)
      WRITE(ndse, 2000) trim(fn_meta), iline, ""
      CALL extcde(10)
    ENDIF
 
    RETURN
 
1000 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '/           &
         '     ERROR READING CRS HEADER - MISSING CRS NAME?' )
    !
1001 FORMAT (/' *** WARNING : USER DEFINED CRS SECTION WILL ' /      &
         '     OVERIDE DEFAULT CRS DEFINITION FOR GRID' /       &
         '     PREV CRS = ', a )
    !
1002 FORMAT (/' *** WARNING : DUPLICATE CRS SECTION WILL ' /         &
         '     OVERRIDE PREVIOUS CRS DEFINITION' /              &
         '     PREV CRS = ', a )
    !
1003 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '/           &
         '     CRS SECTION DOES NOT CONTAIN MANDATORY '/        &
         '     ATTRIBUTE "grid_mapping_name"' )
 
2000 FORMAT ( '     FILENAME = ', a /                                &
         '     LINE NO  = ', i5 /                               &
         '  => ', a /)
    !
 
  END SUBROUTINE read_crs_meta
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE default_crs_meta()
    IMPLICIT NONE
 
    TYPE(meta_pair_t) :: META
 
    IF(flrtd) THEN
#ifdef W3_RTD
      ! Rotated pole location
      crs_name = 'rotated_pole'
      CALL meta_list_append(crs_meta,                                 &
           'grid_mapping_name', 'rotated_latitude_longitude')
      CALL meta_list_append(crs_meta,                                 &
           'grid_north_pole_latitude', polat)
      CALL meta_list_append(crs_meta,                                 &
           'grid_north_pole_longitude', polon)
      crs_is_default = .true.
#endif
    ELSE IF(gtype .EQ. ungtype) THEN
      !       ! What do we want for unstructure grids?
    ELSE
      ! Lat/lon grid
      crs_name = 'crs'
      CALL meta_list_append(crs_meta,                                &
           'grid_mapping_name', 'latitude_longitude')
      ! TODO: Default to a spherical Earth?
      CALL meta_list_append(crs_meta,                                &
           'semi_major_axis', 6371000.0)
      CALL meta_list_append(crs_meta,                                &
           'inverse_flattening', 0.0)
    ENDIF
 
  END SUBROUTINE default_crs_meta
 
  !/ ------------------------------------------------------------------- /
 
  !/ ------------------------------------------------------------------- /
  FUNCTION getmeta(IFI, IFJ, ICOMP, IPART) RESULT(META)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         02-Nov-2020 |
    !/                  +-----------------------------------+
    !/
    !/    09-Nov-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Returns a META_T type containig the netCDF matadata for the
    !     requested field
    !
    !  2. Method :
    !
    !     A copy of the meta-data is returned, rather than a pointer. This
    !     is because in the case of paritioned parameters, the metadata
    !     will be updated with the partition number.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       IFI     Int.  I  Output group number
    !       IFJ     Int.  I  Output field number
    !       ICOMP   Int.  I  Component number (defaults to 1)
    !       IPART   Int.  I  Partition number (defaults to 1)
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
    INTEGER, INTENT(IN) :: ifi, ifj
    INTEGER, INTENT(IN), OPTIONAL :: icomp, ipart
    !/
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    INTEGER :: ifp, ifc
    TYPE(meta_t) :: meta ! Not pointer as we might need to modify it
 
    ifc = 1
    ifp = 1
    IF(PRESENT(icomp)) ifc = icomp
    IF(PRESENT(ipart)) ifp = ipart
 
    ! Error checking on size of IFJ, ICOMP, IPART.
    IF(ifi .LT. 1 .OR. ifi .GT. nogrp) THEN
      WRITE(ndse,1000) nogrp
      CALL extcde(1)
    ENDIF
    IF(ifj .LT. 1 .OR. ifj .GT. noge(ifi)) THEN
      WRITE(ndse,1001) noge(ifi)
      CALL extcde(1)
    ENDIF
    IF(ifc .LT. 1 .OR. ifc .GT. 3) THEN
      WRITE(ndse,1002)
      CALL extcde(1)
    ENDIF
 
    meta = meta_deep_copy(group(ifi)%FIELD(ifj)%META(ifc))
 
    ! For partitioned data, expand in the partition number:
    IF(ifi .EQ. 4) THEN
      CALL add_partno(meta, ifp)
    ENDIF
 
    RETURN
 
1000 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '            / &
         '     GETMETA: IFI value should be in range 1,',i2     / )
    !
1001 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '           / &
         '     GETMETA: IFJ value should be in range 1,',i2    / )
    !
1002 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '           / &
         '     GETMETA: IFC value should be in range 1,3'      / )
    !
  END FUNCTION getmeta
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE read_part_tmpl(NDMI, ILINE)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         04-Dec-2020 |
    !/                  +-----------------------------------+
    !/
    !/    04-Dec-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Reads in a TEMPLATE section from file.
    !     This section defines a list of text strings that will be used
    !     to replace a "placeholder string" when generating metadata for
    !     partitioned parameters.
    !
    !     Format of section is:
    !
    !       TEMPLATE <placeholder_string>
    !         Value for partition IPART=0
    !         Value for partition IPART=1
    !         Value for partition IPART=2
    !         ...
    !         Value for partition IPART=N
    !
    !  2. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       NDMI    Int.  I/O  Unit number
    !       ILINE   Int.  I/O  Line number
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
    INTEGER, INTENT(IN)               :: NDMI
    INTEGER, INTENT(INOUT)            :: ILINE
    !/
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    CHARACTER(LEN=256) :: BUF, ID
    INTEGER            :: IERR
    LOGICAL            :: EOF, NEW
    TYPE(part_tmpl_t), POINTER  :: P
 
    ! Re-read META line to get template string ID (the 2nd field)
    READ(ndmi, '(A)') buf
    READ(buf, *, iostat=ierr) id, id
    IF(ierr /= 0) THEN
      WRITE(ndse, 1000) fn_meta, iline, buf
      CALL extcde(10)
    ENDIF
    id = "<" // trim(id) // ">"
 
    IF(debug) WRITE(*,'(I5,1X,A20,1X,A)') iline, '[template id]', trim(id)
 
    ! Extend list of partition template types:
    IF(ASSOCIATED(part_tmpl)) THEN
      ! Got to end of list
      p => part_tmpl
      DO WHILE(ASSOCIATED(p%NEXT))
        p => p%NEXT
      ENDDO
      ALLOCATE(p%NEXT)
      p => p%NEXT
    ELSE
      ALLOCATE(part_tmpl)
      p => part_tmpl
    ENDIF
 
    ! Set template id and read template strings from file:
    p%TMPL = trim(id)
    ALLOCATE(p%PART_TEXT(0:noswll))
    NULLIFY(p%NEXT)
    p%NP = 0
 
    DO
      CALL next_line(ndmi, buf, iline, eof, new_section=new)
      IF(eof) THEN
        backspace(ndmi)
        RETURN
      ENDIF
 
      IF(new) THEN
        ! Start of new meta data entry
        IF(debug) WRITE(*,'(I5,1X,A20)') iline, '[--end of section--]'
        iline = iline - 1
        backspace(ndmi)
        EXIT
      ENDIF
 
      ! Check we have not exceeded NOSWLL
      IF(p%NP .GT. noswll) THEN
        WRITE(*,*) "Too many partition entries (NOSWLL=",noswll,"). Ignoring"
        cycle
      ENDIF
 
      ! Add string to array of partition text
      IF(debug) THEN
        WRITE(*,'(I5,1X,A20,1X,I1,1X,A)') iline, '[part template]', &
             p%NP, trim(buf)
      ENDIF
 
      p%PART_TEXT(p%NP) = trim(adjustl(buf)) ! Zero indexed
      p%NP = p%NP + 1
    ENDDO
 
    RETURN
    !
1000 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '/           &
         '     ERROR READING PART HEADER - MISSING TEMPLATE ID?'/ &
         '     FILENAME = ', a /                                &
         '     LINE NO =', i5 /                                 &
         '  => ', a /)
    !
  END SUBROUTINE read_part_tmpl
 
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE print_part_tmpl()
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         04-Dec-2020 |
    !/                  +-----------------------------------+
    !/
    !/    04-Dec-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Prints the patition templates to screen (for debug use).
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    TYPE(part_tmpl_t), POINTER :: P
    INTEGER :: I
 
    print*,'=============='
    IF(.NOT. ASSOCIATED(part_tmpl)) THEN
      print*,'Empty partition list'
      RETURN
    ENDIF
 
    p => part_tmpl
    DO
      print*,p%TMPL
      DO i=0,p%NP - 1
        print*,'   - ',i,trim(p%PART_TEXT(i))
      ENDDO
      IF(.NOT. ASSOCIATED(p%NEXT)) EXIT
      p => p%NEXT
    ENDDO
    print*,'=============='
  END SUBROUTINE print_part_tmpl
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE add_partno(META, IPART)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         02-Nov-2020 |
    !/                  +-----------------------------------+
    !/
    !/    09-Nov-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Adds partition number to meta-data.
    !
    !  2. Method :
    !
    !     Replaces all instances of "<IPART>" in the provided meta data with
    !     the partition number IPART.
    !
    !  3. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       META  META_T  I  Meta data type
    !       IPART   Int.  I  Partition number
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
 
    TYPE(meta_t), INTENT(INOUT) :: META
    INTEGER, INTENT(IN) :: IPART
    !/
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    CHARACTER(LEN=80) :: TMP
    INTEGER :: I, J
    TYPE(meta_pair_t), POINTER :: P
 
    CALL partno_string_sub(meta%ENAME, ipart)
    CALL partno_string_sub(meta%VARNM, ipart)
    CALL partno_string_sub(meta%VARNL, ipart)
    CALL partno_string_sub(meta%VARNS, ipart)
    CALL partno_string_sub(meta%VARNG, ipart)
    CALL partno_string_sub(meta%VARNC, ipart)
    CALL partno_string_sub(meta%VARND, ipart)
    IF(meta%EXTRA%N .GT. 0) THEN
      p => meta%EXTRA%HEAD
      DO
        CALL partno_string_sub(p%ATTNAME, ipart)
        IF(p%TYPE .EQ. "c") THEN
          CALL partno_string_sub(p%ATTVAL, ipart)
        ENDIF
        IF(.NOT. ASSOCIATED(p%NEXT)) EXIT
        p => p%NEXT
      ENDDO
    ENDIF
 
  END SUBROUTINE add_partno
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE partno_string_sub(INSTR, IPART)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         02-Nov-2020 |
    !/                  +-----------------------------------+
    !/
    !/    09-Nov-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Performs string substition of placeholder strings with partition
    !     number specfic values.
    !
    !     The placeholder <IPART> is automatically replaced with the
    !     partition number (0, 1, 2, etc).
    !
    !     Other template placeholders can be defined in the ounfmeta.inp
    !     file by the user.
    !
    !  2. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       INSTR  Char.  I/O  Input string
    !       IPART   Int.  I    Partition number
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
 
    CHARACTER(LEN=*), INTENT(INOUT) :: INSTR
    INTEGER, INTENT(IN) :: IPART
    !/
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    INTEGER :: I, J, ISN
    TYPE(part_tmpl_t), POINTER :: P
    CHARACTER(LEN=512) :: TMPL
 
    isn = ipart + 1
    IF(ptmeth .LE. 3) THEN
      IF (isn .GT. 5) isn = 5
    ELSE
      IF (isn .GT. 2) isn = 2
    ENDIF
 
    ! Set partition number (built-in IPART template)
    i = index(instr, ipart_token)
    j = i + len_trim(ipart_token)
    IF(i .GT. 0) THEN
      WRITE(tmpl, '(A,I1,A)') instr(1:i-1), ipart, instr(j:len(instr))
      instr = tmpl
    ENDIF
 
    ! Set standard name string (built-in SPART template)
    i = index(instr, spart_token)
    j = i + len_trim(spart_token)
 
    IF(i .GT. 0) THEN
      instr = instr(1:i-1) // trim(snamep(isn)) // instr(j:len(instr))
    ENDIF
 
    ! Also try underscore separated version: <SPART_>
    i = index(instr, spart_token_)
    j = i + len_trim(spart_token_)
 
    IF(i .GT. 0) THEN
      instr = instr(1:i-1) // trim(replace_char(snamep(isn), " ", "_")) &
           // instr(j:len(instr))
    ENDIF
 
    ! Merge in user defined partition templates (if any):
    IF(.NOT. ASSOCIATED(part_tmpl)) RETURN
 
    p => part_tmpl
    DO
      i = index(instr, trim(p%TMPL))
      j = i + len_trim(p%TMPL)
 
      IF(i .GT. 0) THEN
        IF(ipart .GE. p%NP) THEN
          WRITE(ndse, 1000) trim(p%TMPL), p%NP, ipart
          CALL extcde(10)
        ENDIF
        instr = instr(1:i-1) // trim(p%PART_TEXT(ipart)) // instr(j:len(instr))
      ENDIF
 
      ! Try "underscore" version <TMPL_>:
      i = len_trim(p%TMPL)
      tmpl = p%TMPL(1:i-1) // "_>"
      i = index(instr, trim(tmpl))
      j = i + len_trim(tmpl)
      IF(i .GT. 0) THEN
        instr = instr(1:i-1) // trim(replace_char(p%PART_TEXT(ipart), " ", "_")) &
             // instr(j:len(instr))
      ENDIF
 
      IF(.NOT. ASSOCIATED(p%NEXT)) EXIT
      p => p%NEXT
    ENDDO
 
    RETURN
 
1000 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '            / &
         '     NOT ENOUGH USER DEFINED ENTRIES FOR TEMPLATE'    / &
         '     TEMPLATE ID     : ',a                            / &
         '     NUM ENTRIES     : ',i2                           / &
         '     REQESTED IPART* : ',i2                           / &
         '     (*Note: IPART is zero-refernced)'                / &
         '     Please update your ounfmeta.inp file.'           /)
 
  END SUBROUTINE partno_string_sub
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE write_meta(NCID, VARID, META, ERR)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         02-Nov-2020 |
    !/                  +-----------------------------------+
    !/
    !/    09-Nov-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Writes the meta-data entries for a variable.
    !
    !  2. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       NCID    Int.  I/O  NetCDF file ID
    !       VARID   Int.  I/O  NetCDF variable ID
    !       META    Int.  I    Meta data type
    !       ERR     Int.  O    Error value
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
 
    INTEGER, INTENT(IN) :: NCID, VARID
    TYPE(meta_t), INTENT(IN) :: META
    INTEGER, INTENT(OUT) :: ERR
    !/
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    INTEGER :: IVAL, VAL
    REAL :: RVAL
    !/
    err = nf90_put_att(ncid, varid, 'long_name', meta%VARNL)
    IF(err /= nf90_noerr) RETURN
 
    IF(meta%VARNS .NE. '' .AND. meta%VARNS .NE. unsetc) THEN
      err = nf90_put_att(ncid, varid, 'standard_name', meta%VARNS)
      IF(err /= nf90_noerr) RETURN
    ENDIF
 
    IF(meta%VARNG .NE. '' .AND. meta%VARNG .NE. unsetc) THEN
      err = nf90_put_att(ncid, varid, 'globwave_name', meta%VARNG)
      IF(err /= nf90_noerr) RETURN
    ENDIF
 
    err = nf90_put_att(ncid, varid, 'units', meta%UNITS)
    IF(err /= nf90_noerr) RETURN
 
    ! Fill value dependent on variable type
    IF(ncvartype .EQ. 2) THEN
      err = nf90_put_att(ncid, varid, '_FillValue', nf90_fill_short)
    ELSE
      err = nf90_put_att(ncid, varid, '_FillValue', nf90_fill_float)
    END IF
    IF(err /= nf90_noerr) RETURN
 
    err = nf90_put_att(ncid, varid, 'scale_factor', meta%FSC)
    IF(err /= nf90_noerr) RETURN
 
    err = nf90_put_att(ncid, varid, 'add_offset', 0.)
    IF(err /= nf90_noerr) RETURN
 
    ! For variables with vartype SHORT, the valid min/max
    ! are scaled by scale_factor and converted to integers.
    ! If vartype is FLOAT, then no scaling is performed and
    ! valid min/max are written out directly as floats.
    IF(ncvartype .EQ. 2) THEN
      val = nint(meta%VMIN / meta%FSC)
      err = nf90_put_att(ncid, varid,'valid_min', val)
      IF(err /= nf90_noerr) RETURN
 
      val = nint(meta%VMAX / meta%FSC)
      err = nf90_put_att(ncid, varid,'valid_max', val)
      IF(err /= nf90_noerr) RETURN
    ELSE
      err = nf90_put_att(ncid, varid,'valid_min', meta%VMIN)
      IF(err /= nf90_noerr) RETURN
 
      err = nf90_put_att(ncid, varid,'valid_max', meta%VMAX)
      IF(err /= nf90_noerr) RETURN
    ENDIF
 
    IF(meta%VARNC .NE. '' .AND. meta%VARNC .NE. unsetc) THEN
      err = nf90_put_att(ncid, varid, 'comment', meta%VARNC)
      IF(err /= nf90_noerr) RETURN
    ENDIF
 
    IF(meta%VARND .NE. '' .AND. meta%VARND .NE. unsetc) THEN
      err = nf90_put_att(ncid, varid, 'direction_reference', meta%VARND)
      IF(err /= nf90_noerr) RETURN
    END IF
 
    IF(crs_name .NE. '' .AND. crs_name .NE. unsetc) THEN
      err = nf90_put_att(ncid, varid, 'grid_mapping', crs_name)
      IF(err /= nf90_noerr) RETURN
    ENDIF
 
    IF(coords_attr .NE. '' .AND. coords_attr .NE. unsetc) THEN
      err = nf90_put_att(ncid, varid, 'coordinates', adjustl(coords_attr))
      IF(err /= nf90_noerr) RETURN
    ENDIF
 
    IF (meta%EXTRA%N .GT. 0) THEN
      CALL write_freeform_meta_list(ncid, varid, meta%EXTRA, err)
      IF(err /= nf90_noerr) RETURN
    ENDIF
 
    RETURN
    !
  END SUBROUTINE write_meta
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE write_global_meta(NCID, ERR)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         09-Nov-2020 |
    !/                  +-----------------------------------+
    !/
    !/    09-Nov-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Writes the user meta-data entries for the global attributes
    !
    !  2. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       NCID    Int.  I/O  NetCDF file ID
    !       ERR     Int.  O    Error value
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
 
    INTEGER, INTENT(IN) :: NCID
    INTEGER, INTENT(OUT) :: ERR
    !/
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    CALL write_freeform_meta_list(ncid, nf90_global, global_meta, err)
  END SUBROUTINE write_global_meta
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE write_freeform_meta_list(NCID, VARID, METALIST, ERR)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         16-Dec-2020 |
    !/                  +-----------------------------------+
    !/
    !/    16-Dec-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Writes the freeform user meta-data entries for a NetCDF variable
    !
    !  2. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       NCID      Int.  I/O  NetCDF file ID
    !       VARID     Int.  I/O  NetCDF file ID
    !       METALIST  Type. I    META_LIST_T object to write
    !       ERR       Int.  O    Error value
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
 
    INTEGER, INTENT(IN) :: NCID, VARID
    TYPE(meta_list_t), INTENT(IN) :: METALIST
    INTEGER, INTENT(OUT) :: ERR
    !/
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    INTEGER :: I, IVAL
    REAL    :: RVAL
    TYPE(meta_pair_t), POINTER :: P
 
    IF(metalist%N .EQ. 0) RETURN
 
    p => metalist%HEAD
 
    ! Loop over global metadata pairs:
    DO
 
      IF (p%ATTNAME .EQ. '' .OR.                      &
           p%ATTNAME .EQ. unsetc) cycle
 
      SELECT CASE(p%TYPE)
 
      CASE('i')
        READ(p%ATTVAL, *) ival
        err = nf90_put_att(ncid, varid, p%ATTNAME, ival)
        IF(err /= nf90_noerr) RETURN
 
      CASE('r', 'f')
        READ(p%ATTVAL, *) rval
        err = nf90_put_att(ncid, varid, p%ATTNAME, rval)
        IF(err /= nf90_noerr) RETURN
 
      CASE('c')
        err = nf90_put_att(ncid, varid, p%ATTNAME,    &
             p%ATTVAL)
        IF(err /= nf90_noerr) RETURN
 
      CASE DEFAULT
        WRITE(ndse,1000) p%TYPE
        CALL extcde(10)
      END SELECT
 
      IF(.NOT. ASSOCIATED(p%NEXT)) EXIT
      p => p%NEXT
    ENDDO
    !
1000 FORMAT (/' *** WAVEWATCH III ERROR IN W3OUNFMETA : '            / &
         '     WRITE_FREEFORM_META: Unknown attribute'          / &
         '     data type: ', a1         / )
    !
  END SUBROUTINE write_freeform_meta_list
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE print_meta(META)
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         02-Nov-2020 |
    !/                  +-----------------------------------+
    !/
    !/    09-Nov-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Writes meta-data to the screen - for debugging purposes.
    !
    !  2. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       VARID   Int.  I/O  NetCDF variable ID
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
    TYPE(meta_t), INTENT(IN) :: META
    !/ ------------------------------------------------------------------- /
    !/ Local parameters
    !/
    TYPE(meta_pair_t), POINTER :: P
 
    WRITE(*,*) meta%VARNM
    WRITE(*,"(A20,':',1X,A)") "Standard name", trim(meta%VARNS)
    WRITE(*,"(A20,':',1X,A)") "Long name", trim(meta%VARNL)
    WRITE(*,"(A20,':',1X,A)") "Units", trim(meta%UNITS)
    WRITE(*,"(A20,':',1X,A)") "GLOBWAVE name", trim(meta%VARNG)
    WRITE(*,"(A20,':',1X,A)") "Direction conv", trim(meta%VARND)
    WRITE(*,"(A20,':',1X,A)") "Comment", trim(meta%VARNC)
    WRITE(*,"(A20,':',1X,2F12.3)") "Min/Max", meta%VMIN, meta%VMAX
    IF(meta%EXTRA%N .GT. 0) THEN
      p => meta%EXTRA%HEAD
      DO
        WRITE(*,"(A20,':',1X,A)") trim(p%ATTNAME), trim(p%ATTVAL)
        IF(.NOT. ASSOCIATED(p%NEXT)) EXIT
        p => p%NEXT
      ENDDO
    ENDIF
 
  END SUBROUTINE print_meta
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  FUNCTION meta_deep_copy(META) RESULT(COPY)
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         16-Dec-2020 |
    !/                  +-----------------------------------+
    !/
    !/    16-Dec-2020 : Creation                            ( version 7.12 )
    !/
    !
    !  1. Purpose :
    !
    !     Performs "Deep" copy of a META_T type. This ensures that the
    !     linked list data in the EXTRA field is copied, rather than just
    !     copying the pointer.
    !
    !  2. Parameters :
    !
    !     Parameter list
    !     ----------------------------------------------------------------
    !       META   META_T.  I   META data structure to copy
    !     ----------------------------------------------------------------
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
    TYPE(meta_t), INTENT(IN) :: meta
    TYPE(meta_t) :: copy
 
    ! Shallow copy first:
    copy = meta
 
    ! Now deep copy the EXTRA field (is pointer)
    copy%EXTRA = copy_meta_list(meta%EXTRA)
 
  END FUNCTION meta_deep_copy
 
  !/ ------------------------------------------------------------------- /
  !/ ------------------------------------------------------------------- /
  SUBROUTINE default_meta()
    !/
    !/                  +-----------------------------------+
    !/                  | WAVEWATCH III           NOAA/NCEP |
    !/                  |           C. Bunney               |
    !/                  |                                   |
    !/                  |                        FORTRAN 90 |
    !/                  | Last update :         22-Mar-2021 |
    !/                  +-----------------------------------+
    !/
    !/    02-Nov-2020 : Creation                            ( version 7.12 )
    !/    22-Mar-2021 : Adds extra coupling fields          ( version 7.13 )
    !/
    !
    !  1. Purpose :
    !
    !     Populates the default meta data for ww3_ounf.
    !
    !  2. Remarks :
    !
    !     VMIN and VMAX are now set in the units of the output field.
    !     Previously, they were set with scaled values based on the scaling
    !     factor FSC. The scaling is now performed (if necessary) in the
    !     WRITE_META subroutine.
    !
    !     FSC (scale factor) is only applied to data and valid_min/max if
    !     the netCDF variable type is NF90_SHORT.
    !
    !/ ------------------------------------------------------------------- /
    IMPLICIT NONE
    TYPE(meta_t), POINTER :: META(:)
    !
    !----------GROUP 1 ----------------
    !
    ! IFI=1, IFJ=1, DPT
    meta => group(1)%FIELD(1)%META
    meta(1)%FSC    = 0.5
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.dpt'
    meta(1)%VARNM = 'dpt'
    meta(1)%VARNL ='depth'
    meta(1)%VARNS ='depth'
    meta(1)%VARNG ='depth'
    meta(1)%VARNC =''
    meta(1)%VARND =''
    meta(1)%VMIN = -45000
    meta(1)%VMAX = 70000
    ! IFI=1, IFJ=2, CUR
    meta => group(1)%FIELD(2)%META
    meta(1)%ENAME  = '.cur'
    meta(1)%VARND = dircom
    IF(vector) THEN
      meta(1)%FSC    = 0.01
      meta(1)%UNITS  = 'm s-1'
      meta(1)%VMIN = -9.9
      meta(1)%VMAX =  9.9
      meta(1)%VARNM='ucur'
      meta(1)%VARNL='eastward current'
      meta(1)%VARNS='eastward_sea_water_velocity'
      meta(1)%VARNG='eastward_sea_water_velocity'
      meta(1)%VARNC='cur=sqrt(U**2+V**2)'
 
      ! Second component
      meta(2) = meta(1)
      meta(2)%VARNM='vcur'
      meta(2)%VARNL='northward current'
      meta(2)%VARNS='northward_sea_water_velocity'
      meta(2)%VARNG='northward_sea_water_velocity'
      meta(2)%VARNC='cur=sqrt(U**2+V**2)'
    ELSE
      meta(1)%FSC    = 0.01
      meta(1)%UNITS  = 'm s-1'
      meta(1)%VMIN = 0
      meta(1)%VMAX = 10.0
      meta(1)%VARNM='cspd'
      meta(1)%VARNL='current speed'
      meta(1)%VARNS='sea_water_speed'
      meta(1)%VARNG='sea_water_speed'
 
      ! Second component
      meta(2) = meta(1)
      meta(2)%FSC = 0.1
      meta(2)%VARNM='cdir'
      meta(2)%UNITS= 'degrees'
      meta(2)%VARNL='current direction (toward)'
      meta(2)%VARNS='direction_of_sea_water_velocity'
      meta(2)%VARNG='direction_of_sea_water_velocity'
      meta(2)%VMIN = 0
      meta(2)%VMAX = 360.0
    ENDIF ! VECTOR
    ! IFI=1, IFJ=3
    meta => group(1)%FIELD(3)%META
    meta(1)%ENAME  = '.wnd'
    meta(1)%VARND = dircom
    IF(vector) THEN
      meta(1)%FSC    = 0.1
      meta(1)%UNITS  = 'm s-1'
      meta(1)%VARNM='uwnd'
      meta(1)%VARNL='eastward_wind'
      meta(1)%VARNS='eastward_wind'
      meta(1)%VARNG='eastward_wind'
      meta(1)%VARNC='wind=sqrt(U10**2+V10**2)'
      meta(1)%VMIN = -99.0
      meta(1)%VMAX =  99.0
 
      ! Second component
      meta(2) = meta(1)
      meta(2)%VARNM='vwnd'
      meta(2)%VARNL='northward_wind'
      meta(2)%VARNS='northward_wind'
      meta(2)%VARNG='northward_wind'
      meta(2)%VARNC='wind=sqrt(U10**2+V10**2)'
    ELSE
      meta(1)%FSC = 0.01
      meta(1)%UNITS= 'm s-1'
      meta(1)%VARNM='wspd'
      meta(1)%VARNL='wind speed'
      meta(1)%VARNS='wind_speed'
      meta(1)%VARNG='wind_speed'
      meta(1)%VMIN = 0.0
      meta(1)%VMAX = 100.0
 
      ! Second component
      meta(2) = meta(1)
      meta(2)%FSC = 0.1
      meta(2)%VARNM='wdir'
      meta(2)%UNITS='degrees'
      meta(2)%VARNL='wind direction (from)'
      meta(2)%VARNS='wind_from_direction'
      meta(2)%VARNG='wind_from_direction'
      meta(2)%VMIN = 0.0
      meta(2)%VMAX = 360.0
    ENDIF ! VECTOR
    ! IFI=1, IFJ=4, AST
    meta => group(1)%FIELD(4)%META
    meta(1)%FSC    = 0.1
    meta(1)%ENAME = '.ast'
    meta(1)%UNITS = 'K'
    meta(1)%VARNM='ast'
    meta(1)%VARNL='air sea temperature difference'
    !META(1)%VARNS='air_sea_temperature_difference'
    meta(1)%VARNS=''
    meta(1)%VARNG='air_sea_temperature_difference'
    meta(1)%VMIN = -200.0
    meta(1)%VMAX = 200.0
    ! IFI=1, IFJ=5, WLV
    meta => group(1)%FIELD(5)%META
    meta(1)%FSC    = 0.01
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.wlv'
    meta(1)%VARNM='wlv'
    meta(1)%VARNL='sea surface height above sea level'
    meta(1)%VARNS='sea_surface_height_above_mean_sea_level'
    meta(1)%VARNG='sea_surface_height_above_sea_level'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 100
    ! IFI=1, IFJ=6, ICE
    meta => group(1)%FIELD(6)%META
    meta(1)%FSC    = 0.001
    meta(1)%UNITS  = '1'
    meta(1)%ENAME  = '.ice'
    meta(1)%VARNM='ice'
    meta(1)%VARNL='sea ice area fraction'
    meta(1)%VARNS='sea_ice_area_fraction'
    meta(1)%VARNG='sea_ice_area_fraction'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 1
    ! IFI=1, IFJ=7, IBG
    meta => group(1)%FIELD(7)%META
    meta(1)%FSC    = 0.0001
    meta(1)%UNITS  = 'km-1'
    meta(1)%ENAME  = '.ibg'
    meta(1)%VARNM='ibg'
    meta(1)%VARNL='icebergs_damping'
    !META(1)%VARNS='icebergs_induced_attenuation_scale_for_waves'
    meta(1)%VARNS=''
    meta(1)%VARNG='icebergs_damping'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 3.2
    ! IFI=1, IFJ=8, TAUA
    meta => group(1)%FIELD(8)%META
    meta(1)%FSC    = 0.01
    meta(1)%UNITS  = 'Pa'
    meta(1)%ENAME  = '.taua'
    meta(1)%VARNM='utaua'
    meta(1)%VARNL='surface_downward_eastward_stress'
    meta(1)%VARNS='surface_downward_eastward_stress'
    meta(1)%VARNG='surface_downward_eastward_stress'
    meta(2)%VARNC='taua=sqrt(utaua**2+vtaua**2)'
    meta(1)%VMIN = -320
    meta(1)%VMAX =  320
    meta(1)%VARND = dircom
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vtaua'
    meta(2)%VARNL='surface_downward_northward_stress'
    meta(2)%VARNS='surface_downward_northward_stress'
    meta(2)%VARNG='surface_downward_northward_stress'
    meta(2)%VARNC='taua=sqrt(utaua**2+vtaua**2)'
    ! IFI=1, IFJ=9, RHO
    meta => group(1)%FIELD(9)%META
    meta(1)%FSC    = 0.0001
    meta(1)%UNITS  = 'kg m-3'
    meta(1)%ENAME  = '.rhoa'
    meta(1)%VARNM='rhoa'
    meta(1)%VARNL='air_density'
    meta(1)%VARNS='air_density'
    meta(1)%VARNG='air_density'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 2
    ! IFI=1, IFJ=10, D50
#ifdef W3_BT4
    meta => group(1)%FIELD(10)%META
    meta(1)%FSC    = 0.001
    meta(1)%UNITS  = 'Krumbein phi scale'
    meta(1)%ENAME  = '.d50'
    meta(1)%VARNM='d50'
    meta(1)%VARNL='grain_size'
    !META(1)%VARNS='sediment_grain_size'
    meta(1)%VARNS=''
    meta(1)%VARNG='sediment_grain_size'
    meta(1)%VMIN = -10.0
    meta(1)%VMAX = 32.0
#endif
    ! IFI=1, IFJ=11, IC1
#ifdef W3_IS2
    meta => group(1)%FIELD(11)%META
    meta(1)%FSC = 0.001
    meta(1)%UNITS = 'm'
    meta(1)%ENAME = '.ic1'
    meta(1)%VARNM='ic1'
    meta(1)%VARNL='ice thickness'
    meta(1)%VARNS='sea_ice_thickness'
    meta(1)%VARNG='ice_thickness'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 30
    ! IFI=1, IFJ=12, IC5
    meta => group(1)%FIELD(12)%META
    meta(1)%FSC = 0.05
    meta(1)%UNITS = 'm'
    meta(1)%ENAME = '.ic5'
    meta(1)%VARNM='ic5'
    meta(1)%VARNL='maximum floe diameter'
    !META(1)%VARNS='maximum_ice_floe_diameter'
    meta(1)%VARNS=''
    meta(1)%VARNG='maximum_ice_floe_diameter'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 1500
#endif
    !
    !----------GROUP 2 ----------------
    !
    ! IFI=2, IFJ=1, HS
    meta => group(2)%FIELD(1)%META
    meta(1)%FSC    = 0.002
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.hs'
    meta(1)%VARNM='hs'
    meta(1)%VARNL='significant height of wind and swell waves'
    meta(1)%VARNS='sea_surface_wave_significant_height'
    meta(1)%VARNG='significant_wave_height'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 64
    ! IFI=2, IFJ=2, LM
    meta => group(2)%FIELD(2)%META
    meta(1)%FSC    = 1.
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.lm'
    meta(1)%VARNM='lm'
    meta(1)%VARNL='mean wave length'
    !META(1)%VARNS='mean_wave_length'
    meta(1)%VARNS=''
    meta(1)%VARNG='mean_wave_length'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 3200
    ! IFI=2, IFJ=3, T02
    meta => group(2)%FIELD(3)%META
    meta(1)%FSC    = 0.01
    meta(1)%UNITS  = 's'
    meta(1)%ENAME  = '.t02'
    meta(1)%VARNM='t02'
    meta(1)%VARNL='mean period T02'
    meta(1)%VARNS='sea_surface_wind_wave_mean_period' // &
         '_from_variance_spectral_density_second_frequency_moment'
    meta(1)%VARNG='mean_period_t02'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 50
    ! IFI=2, IFJ=4, T0M1
    meta => group(2)%FIELD(4)%META
    meta(1)%FSC    = 0.01
    meta(1)%UNITS  = 's'
    meta(1)%ENAME  = '.t0m1'
    meta(1)%VARNM='t0m1'
    meta(1)%VARNL='mean period T0m1'
    meta(1)%VARNS='sea_surface_wind_wave_mean_period_from_variance' // &
         '_spectral_density_inverse_frequency_moment'
    meta(1)%VARNG='mean_period_t0m1'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 50
    ! IFI=2, IFJ=5, T01
    meta => group(2)%FIELD(5)%META
    meta(1)%FSC    = 0.01
    meta(1)%UNITS  = 's'
    meta(1)%ENAME  = '.t01'
    meta(1)%VARNM='t01'
    meta(1)%VARNL='mean period T01'
    meta(1)%VARNS='sea_surface_wind_wave_mean_period_from_variance' // &
         '_spectral_density_first_frequency_moment'
    meta(1)%VARNG='mean_period_t01'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 50
    ! IFI=2, IFJ=6, FP
    meta => group(2)%FIELD(6)%META
    meta(1)%FSC    = 0.001
    meta(1)%UNITS  = 's-1'
    meta(1)%ENAME  = '.fp'
    meta(1)%VARNM='fp'
    meta(1)%VARNL='wave peak frequency'
    !META(1)%VARNS='sea_surface_wave_peak_frequency'
    meta(1)%VARNS=''
    meta(1)%VARNG='dominant_wave_frequency'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 10
    ! IFI=2, IFJ=7, DIR
    meta => group(2)%FIELD(7)%META
    meta(1)%FSC    = 0.1
    meta(1)%UNITS  = 'degree'
    meta(1)%ENAME  = '.dir'
    meta(1)%VARNM='dir'
    meta(1)%VARNL='wave mean direction'
    meta(1)%VARNS='sea_surface_wave_from_direction'
    meta(1)%VARNG='wave_from_direction'
    meta(1)%VARND=dircom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 360
    ! IFI=2, IFJ=8, SPR
    meta => group(2)%FIELD(8)%META
    meta(1)%FSC    = 0.1
    meta(1)%UNITS  = 'degree'
    meta(1)%ENAME  = '.spr'
    meta(1)%VARNM='spr'
    meta(1)%VARNL='directional spread'
    meta(1)%VARNS='sea_surface_wave_directional_spread'
    meta(1)%VARNG='directional_spread'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 90
    ! IFI=2, IFJ=9, DP
    meta => group(2)%FIELD(9)%META
    meta(1)%FSC    = 1.
    meta(1)%UNITS  = 'degree'
    meta(1)%ENAME  = '.dp'
    meta(1)%VARNM='dp'
    meta(1)%VARNL='peak direction'
    meta(1)%VARNS='sea_surface_wave_peak_direction'
    meta(1)%VARNG='dominant_wave_direction'
    meta(1)%VARND=dircom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 360
    ! IFI=2, IFJ=10, HIG
    meta => group(2)%FIELD(10)%META
    meta(1)%FSC    = 0.0002
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.hig'
    meta(1)%VARNM='hig'
    meta(1)%VARNL='infragravity_wave_height'
    !META(1)%VARNS='sea_surface_wave_infragravity_significant_height'
    meta(1)%VARNS=''
    meta(1)%VARNG='infragravity_significant_wave_height'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 1.0
    ! IFI=2, IFJ=11, MXE
    meta => group(2)%FIELD(11)%META
    meta(1)%FSC    = 0.002
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.mxe'
    meta(1)%VARNM='stmaxe'
    meta(1)%VARNL='expected maximum sea surface elevation (nonlinear,2nd order)'
    !META(1)%VARNS='expected maximum sea surface elevation (nonlinear,2nd order)'
    meta(1)%VARNS=''
    meta(1)%VARNG='expected maximum sea surface elevation (nonlinear,2nd order)'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 64
    ! IFI=2, IFJ=12, MXES
    meta => group(2)%FIELD(12)%META
    meta(1)%FSC    = 0.002
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.mxes'
    meta(1)%VARNM='stmaxd'
    meta(1)%VARNL='standard deviation of maximum sea surface elevation (nonlinear,2nd order)'
    !META(1)%VARNS='std of expected maximum sea surface elevation (nonlinear,2nd order)'
    meta(1)%VARNS=''
    meta(1)%VARNG='standard deviation of maximum sea surface elevation (nonlinear,2nd order)'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 64
    ! IFI=2, IFJ=13, MXH
    meta => group(2)%FIELD(13)%META
    meta(1)%FSC    = 0.002
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.mxh'
    meta(1)%VARNM='hmaxe'
    meta(1)%VARNL='expected maximum wave height (linear, 1st order)'
    !META(1)%VARNS='expected maximum wave height (linear, 1st order)'
    meta(1)%VARNS=''
    meta(1)%VARNG='expected maximum wave height (linear, 1st order)'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 64
    ! IFI=2, IFJ=14, MXHC
    meta => group(2)%FIELD(14)%META
    meta(1)%FSC    = 0.002
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.mxhc'
    meta(1)%VARNM='hcmaxe'
    meta(1)%VARNL='expected maximum wave height from crest (linear, 1st order)'
    !META(1)%VARNS='expected maximum wave height from crest (linear, 1st order)'
    meta(1)%VARNS=''
    meta(1)%VARNG='expected maximum wave height from crest (linear, 1st order)'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 64
    ! IFI=2, IFJ=15, SDMH
    meta => group(2)%FIELD(15)%META
    meta(1)%FSC    = 0.002
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.sdmh'
    meta(1)%VARNM='hmaxd'
    meta(1)%VARNL='STD of maximum wave height (linear, 1st order)'
    !META(1)%VARNS='STD of maximum wave height (linear, 1st order)'
    meta(1)%VARNS=''
    meta(1)%VARNG='STD of maximum wave height (linear, 1st order)'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 64
    ! IFI=2, IFJ=16, SDMHC
    meta => group(2)%FIELD(16)%META
    meta(1)%FSC    = 0.002
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.sdmhc'
    meta(1)%VARNM='hcmaxd'
    meta(1)%VARNL='STD of maximum wave height from crest (linear, 1st order)'
    !META(1)%VARNS='STD of maximum wave height from crest (linear, 1st order)'
    meta(1)%VARNS=''
    meta(1)%VARNG='STD of maximum wave height from crest (linear, 1st order)'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 64
    ! IFI=2, IFJ=17, WBT
    meta => group(2)%FIELD(17)%META
    meta(1)%FSC    = 0.001
    meta(1)%UNITS  = '1'
    meta(1)%ENAME  = '.wbt'
    meta(1)%VARNM='wbt'
    meta(1)%VARNL='dominant wave breaking probability'
    !META(1)%VARNS='dominant_wave_breaking_probability'
    meta(1)%VARNS=''
    meta(1)%VARNG='dominant_wave_breaking_probability'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 1
    ! IFI=2, IFJ=18, TP
    meta => group(2)%FIELD(18)%META
    meta(1)%FSC    = 0.01
    meta(1)%UNITS  = 's'
    meta(1)%ENAME  = '.tp'
    meta(1)%VARNM='tp'
    meta(1)%VARNL='wave peak period'
    meta(1)%VARNS='sea_surface_wave_peak_period'
    meta(1)%VARNG='dominant_wave_period'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 50
    ! IFI=2, IFJ=19
    meta => group(2)%FIELD(19)%META
    meta(1)%FSC    = 0.001
    meta(1)%UNITS  = 'm-1'
    meta(1)%ENAME  = '.wnm'
    meta(1)%VARNM='wnm'
    meta(1)%VARNL='mean wave number'
    meta(1)%VARNS=''
    meta(1)%VARNG=''
    meta(1)%VMIN = 0
    meta(1)%VMAX = 32
    !
    !---------- GROUP 3 ----------------
    !
    ! IFI=3, IFJ=1, EF
    meta => group(3)%FIELD(1)%META
    meta(1)%VARNM='ef'
    meta(1)%VARNL='wave_elevation_spectrum'
    meta(1)%VARNS='sea_surface_wave_variance_spectral_density'
    IF (ncvartype.LE.3) THEN
      meta(1)%UNITS  = 'log10(m2 s+1E-12)'
      !META(1)%VARNS='base_ten_logarithm_of_power_spectral_density_of_surface_elevation'
      meta(1)%VARNC='base_ten_logarithm'
      meta(1)%FSC = 0.0004
      meta(1)%VMIN = -12.
      meta(1)%VMAX = 12.
    ELSE
      meta(1)%UNITS  = 'm2 s'
      meta(1)%FSC = 1.
      meta(1)%VMIN = 0.
      meta(1)%VMAX = 1.e12
    END IF
    meta(1)%ENAME  = '.ef'
    meta(1)%VARNG = meta(1)%VARNS
    ! IFI=3, IFJ=2, TH1M
    meta => group(3)%FIELD(2)%META
    ! Information for spectral
    meta(1)%FSC     = 0.1
    meta(1)%VARNM='th1m'
    meta(1)%VARNL='mean wave direction frequency spectrum'
    !META(1)%VARNS='sea_surface_wave_from_direction_frequency_spectrum'
    meta(1)%VARNS=''
    meta(1)%VARNG = meta(1)%VARNS
    meta(1)%VARND=dircom
    meta(1)%UNITS  = 'degree'
    meta(1)%ENAME  = '.th1m'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 360
    ! IFI=3, IFJ=3, STH1M
    meta => group(3)%FIELD(3)%META
    ! Information for spectral
    meta(1)%FSC     = 0.01
    meta(1)%VARNM='sth1m'
    meta(1)%VARNL='spreading frequency spectrum'
    !META(1)%VARNS='sea_surface_wave_spreading_spectrum'
    meta(1)%VARNS=''
    meta(1)%VARNG = meta(1)%VARNS
    meta(1)%UNITS  = 'degree'
    meta(1)%ENAME  = '.sth1m'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 90
    ! IFI=3, IFJ=4, TH2M
    meta => group(3)%FIELD(4)%META
    meta(1)%FSC     = 0.1
    meta(1)%VARNM='th2m'
    meta(1)%VARNL='second mean wave direction frequency spectrum'
    !META(1)%VARNS='sea_surface_wave_from_direction_frequency_spectrum_from_second_moments'
    meta(1)%VARNS=''
    meta(1)%VARNG = meta(1)%VARNS
    meta(1)%VARND=dircom
    meta(1)%UNITS  = 'degree'
    meta(1)%ENAME  = '.th2m'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 360
    ! IFI=3, IFJ=5, STH2M
    meta => group(3)%FIELD(5)%META
    meta(1)%FSC     = 0.01
    meta(1)%VARNM='sth2m'
    meta(1)%VARNL='second spreading frequency spectrum'
    !META(1)%VARNS='sea_surface_wave_spreading_spectrum_from_second_moments'
    meta(1)%VARNS=''
    meta(1)%VARNG = meta(1)%VARNS
    meta(1)%UNITS  = 'degree'
    meta(1)%ENAME  = '.sth2m'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 90
    ! IFI=3, IFJ=6, WN
    meta => group(3)%FIELD(6)%META
    ! Information for spectral
    meta(1)%FSC    = 0.001
    meta(1)%UNITS  = 'm-1'
    meta(1)%ENAME  = '.wn'
    meta(1)%VARNM='wn'
    meta(1)%VARNL='wave numbers'
    !META(1)%VARNS='wave_numbers'
    meta(1)%VARNS=''
    meta(1)%VARNG='wave_numbers'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 32
    !
    !---------- GROUP 4 ----------------
    !
    ! IFI=4, IFJ=1, PHS
    meta => group(4)%FIELD(1)%META
    meta(1)%FSC    = 0.002
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME = '.phs'// ipart_token
    meta(1)%VARNM = 'phs'// ipart_token
    meta(1)%VARNL = 'wave significant height partition '// ipart_token
    meta(1)%VARNS = 'sea_surface_'// spart_token_ //'_wave_significant_height'
    meta(1)%VARNG = 'significant_wave_height_partition_'// ipart_token
    meta(1)%VARNC = partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 64
    ! IFI=4, IFJ=2, PTP
    meta => group(4)%FIELD(2)%META
    meta(1)%FSC = 0.01
    meta(1)%UNITS = 's'
    meta(1)%ENAME = '.ptp'// ipart_token
    meta(1)%VARNM = 'ptp'// ipart_token
    meta(1)%VARNL = 'peak period partition '// ipart_token
    meta(1)%VARNS = 'sea_surface_'// spart_token_ //'_wave_period_at_variance' // &
         '_spectral_density_maximum'
    meta(1)%VARNG = 'dominant_wave_period_partition_'// ipart_token
    meta(1)%VARNC = partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 100
    ! IFI=4, IFJ=3, PLP
    meta => group(4)%FIELD(3)%META
    meta(1)%FSC = 1.
    meta(1)%UNITS = 'm'
    meta(1)%ENAME = '.plp'// ipart_token
    meta(1)%VARNM = 'plp'// ipart_token
    meta(1)%VARNL = 'peak wave length partition '// ipart_token
    !META(1)%VARNS = 'peak_wave_length_partition_'// SPART_TOKEN_
    meta(1)%VARNS = ''
    meta(1)%VARNG = 'peak_wave_length_partition_'// ipart_token
    meta(1)%VARNC = partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 10000
    ! IFI=4, IFJ=4, PDIR
    meta => group(4)%FIELD(4)%META
    meta(1)%FSC = 0.1
    meta(1)%UNITS = 'degree'
    meta(1)%ENAME =  '.pdir'// ipart_token
    meta(1)%VARNM =  'pdir'// ipart_token
    meta(1)%VARNL = 'wave mean direction partition '// ipart_token
    meta(1)%VARNS = 'sea_surface_'// spart_token_ //'_wave_from_direction'
    meta(1)%VARNG = 'wave_from_direction_partition_'// ipart_token
    meta(1)%VARNC = partcom
    meta(1)%VARND = dircom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 360
    ! IFI=4, IFJ=5, PSPR
    meta => group(4)%FIELD(5)%META
    meta(1)%FSC = 0.1
    meta(1)%UNITS = 'degree'
    meta(1)%ENAME = '.pspr'// ipart_token
    meta(1)%VARNM = 'pspr'// ipart_token
    meta(1)%VARNL = 'directional spread partition '// ipart_token
    meta(1)%VARNS = 'sea_surface_'// spart_token_ //'_wave_diectional_spread'
    meta(1)%VARNG = 'directional_spread_partition_'// ipart_token
    meta(1)%VARNC = partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 90
    ! IFI=4, IFJ=6, PWS
    meta => group(4)%FIELD(6)%META
    meta(1)%FSC = 0.001
    meta(1)%UNITS = '1'
    meta(1)%ENAME = '.pws'// ipart_token
    meta(1)%VARNM = 'pws'// ipart_token
    meta(1)%VARNL = 'wind sea fraction in partition '// ipart_token
    !META(1)%VARNS = 'wind_sea_fraction_in_partition_'// IPART_TOKEN
    meta(1)%VARNS = ''
    meta(1)%VARNG = 'wind_sea_fraction_in_partition_'// ipart_token
    meta(1)%VARNC = partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 1
    ! IFI=4, IFJ=7, PDP
    meta => group(4)%FIELD(7)%META
    meta(1)%FSC = 0.1
    meta(1)%UNITS = 'degree'
    meta(1)%ENAME = '.pdp'// ipart_token
    meta(1)%VARNM = 'pdp'// ipart_token
    meta(1)%VARNL = 'peak direction partition '// ipart_token
    meta(1)%VARNS = 'sea_surface_'// spart_token_ //'_wave_from_direction_at_variance' // &
         '_spectral_density_maximum'
    meta(1)%VARNG = 'dominant_wave_from_direction_partition_'// ipart_token
    meta(1)%VARNC = partcom
    meta(1)%VARND = dircom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 360
    ! IFI=4, IFJ=8, PQP
    meta => group(4)%FIELD(8)%META
    meta(1)%FSC = 0.01
    meta(1)%UNITS = '1'
    meta(1)%ENAME = '.pqp'// ipart_token
    meta(1)%VARNM = 'pqp'// ipart_token
    meta(1)%VARNL = 'peakedness partition '// ipart_token
    !META(1)%VARNS = 'sea_surface_wave_peakedness_partition_'// IPART_TOKEN
    meta(1)%VARNS = ''
    meta(1)%VARNG = 'wave_peakedness_partition_'// ipart_token
    meta(1)%VARNC = partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 320
    ! IFI=4, IFJ=9, PPE
    meta => group(4)%FIELD(9)%META
    meta(1)%FSC = 0.01
    meta(1)%UNITS = '1'
    meta(1)%ENAME = '.ppe'// ipart_token
    meta(1)%VARNM = 'ppe'// ipart_token
    meta(1)%VARNL = 'peak enhancement factor partition '// ipart_token
    !META(1)%VARNS = 'wave_peak_enhancement_factor_partition_'// IPART_TOKEN
    meta(1)%VARNS = ''
    meta(1)%VARNG = 'wave_peak_enhancement_factor_partition_'// ipart_token
    meta(1)%VARNC = 'JONSWAP peak enhancement factor; ' // partcom
    meta(1)%VARND = ''
    meta(1)%VMIN = 0
    meta(1)%VMAX = 320
    ! IFI=4, IFJ=10, PGW
    meta => group(4)%FIELD(10)%META
    meta(1)%FSC = 0.0001
    meta(1)%UNITS = 's-1'
    meta(1)%ENAME =  '.pgw'// ipart_token
    meta(1)%VARNM = 'pgw'// ipart_token
    meta(1)%VARNL = 'frequency width partition '// ipart_token
    !META(1)%VARNS = 'Gaussian_frequency_spread_partition_'// IPART_TOKEN
    meta(1)%VARNS = ''
    meta(1)%VARNG = 'Gaussian_frequency_spread_partition_'// ipart_token
    meta(1)%VARNC = 'Gaussian least-square fit to ' // &
         'omni-directional spectral partition; ' // partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 3.2
    ! IFI=4, IFJ=11, PSW
    meta => group(4)%FIELD(11)%META
    meta(1)%FSC = 0.0001
    meta(1)%UNITS = '1'
    meta(1)%ENAME = '.psw'// ipart_token
    meta(1)%VARNM = 'psw'// ipart_token
    meta(1)%VARNL = 'spectral width partition '// ipart_token
    !META(1)%VARNS = 'sea_surface_wave_spectral_width_partition_'// IPART_TOKEN
    meta(1)%VARNS = ''
    meta(1)%VARNG = 'wave_spectral_width_partition_'// ipart_token
    meta(1)%VARNC = partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 3.2
    ! IFI=4, IFJ=12, PTM10
    meta => group(4)%FIELD(12)%META
    meta(1)%FSC = 0.01
    meta(1)%UNITS = 's'
    meta(1)%ENAME = '.ptm10c'// ipart_token
    meta(1)%VARNM = 'ptm10c'// ipart_token
    meta(1)%VARNL = 'mean period Tm10 partition '// ipart_token
    meta(1)%VARNS = 'sea_surface_'// spart_token_ //'_wave_mean_period_from_variance' // &
         '_spectral_density_inverse_frequency_moment'
    meta(1)%VARNG = 'mean_wave_period_Tm10_partition_'// ipart_token
    meta(1)%VARNC = partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 100
    ! IFI=4, IFJ=13, PT01
    meta => group(4)%FIELD(13)%META
    meta(1)%FSC = 0.01
    meta(1)%UNITS = 's'
    meta(1)%ENAME = '.pt01c'// ipart_token
    meta(1)%VARNM = 'pt01c'// ipart_token
    meta(1)%VARNL = 'mean period T01 partition '// ipart_token
    meta(1)%VARNS = 'sea_surface_'// spart_token_ //'_wave_mean_period_from_variance' // &
         '_spectral_density_first_frequency_moment'
    meta(1)%VARNG = 'mean_wave_period_T01_partition_'// ipart_token
    meta(1)%VARNC = partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 100
    ! IFI=4, IFJ=14, PT02
    meta => group(4)%FIELD(14)%META
    meta(1)%FSC = 0.01
    meta(1)%UNITS = 's'
    meta(1)%ENAME = '.pt02c'// ipart_token
    meta(1)%VARNM = 'pt02c'// ipart_token
    meta(1)%VARNL = 'mean period T02 partition '// ipart_token
    meta(1)%VARNS = 'sea_surface_'// spart_token_ //'_wave_mean_period_from_variance' // &
         '_spectral_density_second_frequency_moment'
    meta(1)%VARNG = 'mean_wave_period_T02_partition_'// ipart_token
    meta(1)%VARNC = partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 100
    ! IFI=4, IFJ=15, PEP
    meta => group(4)%FIELD(15)%META
    meta(1)%FSC = 0.02
    meta(1)%UNITS = 'm2 s rad-1'
    meta(1)%ENAME = '.pep'// ipart_token
    meta(1)%VARNM = 'pep'// ipart_token
    meta(1)%VARNL = 'energy at peak frequency partition '// ipart_token
    !META(1)%VARNS = 'sea_surface_wave_energy_at_variance_spectral_density_maximum_partition_'// IPART_TOKEN
    meta(1)%VARNS = ''
    meta(1)%VARNG = 'wave_energy_at_variance_spectral_density_maximum_partition_'// ipart_token
    meta(1)%VARNC = partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 200
    ! IFI=4, IFJ=16, TWS
    meta => group(4)%FIELD(16)%META
    meta(1)%FSC = 0.001
    meta(1)%UNITS = '1'
    meta(1)%ENAME = '.tws'
    meta(1)%VARNM = 'tws'
    meta(1)%VARNL = 'wind sea fraction'
    !META(1)%VARNS = 'wind_sea_fraction'
    meta(1)%VARNS = ''
    meta(1)%VARNG = 'wind_sea_fraction'
    meta(1)%VARNC = partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 1
    ! IFI=4, IFJ=17, PNR
    meta => group(4)%FIELD(17)%META
    meta(1)%FSC = 1.
    meta(1)%UNITS = '1'
    meta(1)%ENAME = '.pnr'
    meta(1)%VARNM = 'pnr'
    meta(1)%VARNL = 'number of wave partitions'
    !META(1)%VARNS = 'number_of_wave_partitions'
    meta(1)%VARNS = ''
    meta(1)%VARNG = 'number_of_wave_partitions'
    meta(1)%VARNC = partcom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 100
    !
    !---------- GROUP 5 ----------------
    !
    ! IFI=5, IFJ=1, UST
    meta => group(5)%FIELD(1)%META
    ! First component
    meta(1)%FSC    = 0.001
    meta(1)%ENAME  = '.ust'
    meta(1)%UNITS  = 'm s-1'
    meta(1)%VARNM='uust'
    meta(1)%VARNL='eastward friction velocity'
    !META(1)%VARNS='eastward_friction_velocity'
    meta(1)%VARNS=''
    meta(1)%VARNG='eastward_friction_velocity'
    meta(1)%VARNC='ust=sqrt(uust**2+vust**2)'
    meta(1)%VARND=dircom
    meta(1)%VMIN = -99.0
    meta(1)%VMAX =  99.0
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vust'
    meta(2)%VARNL='northward friction velocity'
    !META(2)%VARNS='northward_friction_velocity'
    meta(2)%VARNS=''
    meta(2)%VARNG='northward_friction_velocity'
    ! IFI=5, IFJ=2, CHA
    meta => group(5)%FIELD(2)%META
    meta(1)%FSC    = 1.e-5
    meta(1)%UNITS  = '1'
    meta(1)%ENAME  = '.cha'
    meta(1)%VARNM='cha'
    meta(1)%VARNL='charnock coefficient for surface roughness length for momentum in air'
    meta(1)%VARNS='charnock_coefficient_for_surface_roughness_length_for_momentum_in_air'
    meta(1)%VARNG='charnock_coefficient'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 0.327
    ! IFI=5, IFJ=3, CGE
    meta => group(5)%FIELD(3)%META
    meta(1)%FSC    = 0.1           !0.01
    meta(1)%UNITS  = 'kW m-1'
    meta(1)%ENAME  = '.cge'
    meta(1)%VARNM='cge'
    meta(1)%VARNL='wave energy flux'
    !META(1)%VARNS='sea_surface_wind_wave_energy_flux'
    meta(1)%VARNS=''
    meta(1)%VARNG='wave_energy_flux'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 999
    ! IFI=5, IFJ=4, FAW
    meta => group(5)%FIELD(4)%META
    meta(1)%FSC    = 0.1
    meta(1)%UNITS  = 'W m-2'
    meta(1)%ENAME  = '.faw'
    meta(1)%VARNM='faw'
    meta(1)%VARNL='wind to wave energy flux'
    meta(1)%VARNS='wind_mixing_energy_flux_into_sea_water'
    meta(1)%VARNG='wind_to_wave_energy_flux'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 999
    ! IFI=5, IFJ=5, TAW
    meta => group(5)%FIELD(5)%META
    ! First component
    meta(1)%FSC    = 0.000001
    meta(1)%UNITS  = 'm2 s-2'
    meta(1)%ENAME  = '.taw'
    meta(1)%VARNM='utaw'
    meta(1)%VARNL='eastward wave supported wind stress'
    !META(1)%VARNS='eastward_wave_supported_wind_stress'
    meta(1)%VARNS=''
    meta(1)%VARNC='taw=sqrt(utaw**2+vtaw**2)'
    meta(1)%VARNG='eastward_wave_supported_wind_stress'
    meta(1)%VARND=dircom
    meta(1)%VMIN = -0.032
    meta(1)%VMAX =  0.032
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vtaw'
    meta(2)%VARNL='northward wave supported wind stress'
    !META(2)%VARNS='northward_wave_supported_wind_stress'
    meta(2)%VARNS=''
    meta(2)%VARNG='northward_wave_supported_wind_stress'
    meta(2)%VARNC='taw=sqrt(utaw**2+vtaw**2)'
    ! IFI=5, IFJ=6, TWA
    meta => group(5)%FIELD(6)%META
    ! First component
    meta(1)%FSC    = 0.0001
    meta(1)%ENAME  = '.twa'
    meta(1)%UNITS  = 'm2 s-2'
    meta(1)%VARNM='utwa'
    meta(1)%VARNL='eastward wave to wind stress'
    !META(1)%VARNS='eastward_wave_to_wind_stress'
    meta(1)%VARNS=''
    meta(1)%VARNG='eastward_wave_to_wind_stress'
    meta(1)%VARNC='twa=sqrt(utwa**2+vtwa**2)'
    meta(1)%VARND=dircom
    meta(1)%VMIN = -3.2
    meta(1)%VMAX =  3.2
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vtwa'
    meta(2)%VARNL='northward wave to wind stress'
    !META(2)%VARNS='northward_wave_to_wind_stress'
    meta(2)%VARNS=''
    meta(2)%VARNG='northward_wave_to_wind_stress'
    meta(2)%VARNC='twa=sqrt(utwa**2+vtwa**2)'
    ! IFI=5, IFJ=7, WCC
    meta => group(5)%FIELD(7)%META
    meta(1)%FSC    = 0.0001
    meta(1)%UNITS  = '1'
    meta(1)%ENAME  = '.wcc'
    meta(1)%VARNM='wcc'
    meta(1)%VARNL='whitecap coverage'
    !META(1)%VARNS='whitecap_coverage'
    meta(1)%VARNS=''
    meta(1)%VARNG='whitecap_coverage'
    meta(1)%VARNC=''
    meta(1)%VARND=''
    meta(1)%VMIN = 0
    meta(1)%VMAX = 1
    ! IFI=5, IFJ=8, WCF
    meta => group(5)%FIELD(8)%META
    meta(1)%FSC    = 0.001
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.wcf'
    meta(1)%VARNM='wcf'
    meta(1)%VARNL='whitecap foam thickness'
    !META(1)%VARNS='whitecap_foam_thickness'
    meta(1)%VARNS=''
    meta(1)%VARNG='whitecap_foam_thickness'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 10
    ! IFI=5, IFJ=9, WCH
    meta => group(5)%FIELD(9)%META
    meta(1)%FSC    = 0.002
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.wch'
    meta(1)%VARNM='wch'
    meta(1)%VARNL='significant breaking wave height'
    !META(1)%VARNS='significant_breaking_wave_height'
    meta(1)%VARNS=''
    meta(1)%VARNG='significant_breaking_wave_height'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 64
    ! IFI=5, IFJ=10, WCM
    meta => group(5)%FIELD(10)%META
    meta(1)%FSC    = 0.0001
    meta(1)%UNITS  = '1'
    meta(1)%ENAME  = '.wcm'
    meta(1)%VARNM='wcm'
    meta(1)%VARNL='whitecap moment'
    !META(1)%VARNS='whitecap_moment'
    meta(1)%VARNS=''
    meta(1)%VARNG='whitecap_moment'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 1
    ! IFI=5, IFJ=11, FWS
    meta => group(5)%FIELD(11)%META
    meta(1)%FSC    = 0.002
    meta(1)%UNITS  = 's'
    meta(1)%ENAME  = '.fws'
    meta(1)%VARNM='fws'
    meta(1)%VARNL='Wind_sea_mean_period_T0M1'
    meta(1)%VARNS='sea_surface_wind_wave_mean_period_from_variance' // &
         '_spectral_density_inverse_frequency_moment'
    meta(1)%VARNG='Wind_sea_mean_period_T0M1'
    meta(1)%VARNC=''
    meta(1)%VARND=''
    meta(1)%VMIN = 0
    meta(1)%VMAX = 64
    !
    !---------- GROUP 6 ----------------
    !
    ! IFI=6, IFJ=1, SXY
    meta => group(6)%FIELD(1)%META
    meta(1)%FSC  = 10.
    meta(1)%UNITS = 'N m-1'
    meta(1)%ENAME = '.sxy'
    meta(1)%VARND = dircom
    meta(1)%VMIN = -30000
    meta(1)%VMAX = 30000
 
    ! First component
    meta(1)%VARNM='sxx'
    meta(1)%VARNL='Radiation stress component Sxx'
    !META(1)%VARNS='radiation_stress_component_sxx'
    meta(1)%VARNS=''
 
    ! S6cond component
    meta(2) = meta(1)
    meta(2)%VARNM='syy'
    meta(2)%VARNL='Radiation stress component Syy'
    !META(2)%VARNS='radiation_stress_component_syy'
    meta(2)%VARNS=''
 
    ! Third component
    meta(3) = meta(1)
    meta(3)%FSC = 1.
    meta(3)%VARNM='sxy'
    meta(3)%VARNL='Radiation stress component Sxy'
    !META(3)%VARNS='radiation_stress_component_sxy'
    meta(3)%VARNS=''
    ! IFI=6, IFJ=2, TWO
    meta => group(6)%FIELD(2)%META
    meta(1)%FSC    = 0.000001
    meta(1)%UNITS  = 'm2 s-2'
    meta(1)%ENAME  = '.two'
    meta(1)%VMIN = -0.032
    meta(1)%VMAX =  0.032
    meta(1)%VARND = dircom
 
    ! First component
    meta(1)%VARNM='utwo'
    meta(1)%VARNL='eastward wave to ocean stress'
    !META(1)%VARNS='eastward_wave_to_ocean_stress'
    meta(1)%VARNS=''
    meta(1)%VARNG='eastward_wave_to_ocean_stress'
    meta(1)%VARNC='two=sqrt(utwo**2+vtwo**2)'
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vtwo'
    meta(2)%VARNL='northward wave to ocean stress'
    !META(2)%VARNS='northward_wave_to_ocean_stress'
    meta(2)%VARNS=''
    meta(2)%VARNG='northward_wave_to_ocean_stress'
    meta(2)%VARNC='two=sqrt(utwo**2+vtwo**2)'
    ! IFI=6, IFJ=3, BHD
    meta => group(6)%FIELD(3)%META
    meta(1)%FSC    = 0.1
    meta(1)%UNITS  = 'm2 s-2'
    meta(1)%ENAME  = '.bhd'
    meta(1)%VARNM='bhd'
    meta(1)%VARNL='radiation pressure (Bernouilli Head)'
    !META(1)%VARNS='radiation_pressure'
    meta(1)%VARNS=''
    meta(1)%VARNG='radiation_pressure'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 100
    ! IFI=6, IFJ=4, FOC
    meta => group(6)%FIELD(4)%META
    meta(1)%FSC    = 0.1
    meta(1)%UNITS  = 'W m-2'
    meta(1)%ENAME  = '.foc'
    meta(1)%VARNM='foc'
    meta(1)%VARNL='wave to ocean energy flux'
    !META(1)%VARNS='wave_to_ocean_energy_flux'
    meta(1)%VARNS=''
    meta(1)%VARNG='wave_to_ocean_energy_flux'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 999
    ! IFI=6, IFJ=5, TUS
    meta => group(6)%FIELD(5)%META
    meta(1)%FSC    = 0.001
    meta(1)%UNITS  = 'm2 s-1'
    meta(1)%ENAME  = '.tus'
    meta(1)%VARND = dircom
    meta(1)%VMIN = -32.0 ! C Hansen: The former values of +-9.9 might be
    meta(1)%VMAX = 32.0  ! exceeded more frequently in real storms
 
    ! First component
    meta(1)%VARNM='utus'
    meta(1)%VARNL='eastward stokes transport'
    !META(1)%VARNS='eastward_stokes_transport'
    meta(1)%VARNS=''
    meta(1)%VARNG='eastward_stokes_transport'
    meta(1)%VARNC='tus=sqrt(utus**2+vtus**2)'
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vtus'
    meta(2)%VARNL='northward stokes transport'
    !META(2)%VARNS='northward_stokes_transport'
    meta(2)%VARNS=''
    meta(2)%VARNG='northward_stokes_transport'
    meta(2)%VARNC='tus=sqrt(utus**2+vtus**2)'
 
    ! IFI=6, IFJ=6, USS
    meta => group(6)%FIELD(6)%META
    meta(1)%FSC    = 0.0005
    meta(1)%UNITS  = 'm s-1'
    meta(1)%ENAME  = '.uss'
 
    ! First component
    meta(1)%VARNM='uuss'
    meta(1)%VARNL='eastward surface stokes drift'
    meta(1)%VARNS='sea_surface_wave_stokes_drift_eastward_velocity'
    meta(1)%VARNG='eastward_surface_stokes_drift'
    meta(1)%VARNC='uss=sqrt(uuss**2+vuss**2)'
    meta(1)%VARND=dircom
    meta(1)%VMIN = -4.95
    meta(1)%VMAX =  4.95
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vuss'
    meta(2)%VARNL='northward surface stokes drift'
    meta(2)%VARNS='sea_surface_wave_stokes_drift_northward_velocity'
    meta(2)%VARNG='northward_surface_stokes_drift'
    WRITE(meta(2)%VARNC,'(A,F8.4,A,F8.4,A)') 'Frequency range ',sig(1)*tpiinv,' to ',sig(nk)*tpiinv,' Hz'
    ! IFI=6, IFJ=7, P2S
    meta => group(6)%FIELD(7)%META
    meta(1)%FSC    = 0.01
    meta(1)%ENAME  = '.p2s'
    meta(1)%UNITS  = 'm4'
    meta(1)%VMIN = -150
    meta(1)%VMAX = 320
 
    ! First component
    meta(1)%VARNL='power spectral density of equivalent surface pressure'
    !META(1)%VARNS='power_spectral_density_of_equivalent_surface_pressure'
    meta(1)%VARNS=''
    meta(1)%VARNG='power_spectral_density_of_equivalent_surface_pressure'
    meta(1)%VARNM='fp2s'
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='pp2s'
    meta(2)%UNITS= 's-1'
    meta(2)%VARNL='peak period of power spectral density of equivalent surface pressure'
    !META(2)%VARNS='peak_period_of_power_spectral_density_of_equivalent_surface_pressure'
    meta(2)%VARNS=''
    meta(2)%VARNG='peak_period_of_power_spectral_density_of_equivalent_surface_pressure'
 
    ! IFI=6, IFJ=8, USF
    meta => group(6)%FIELD(8)%META
    meta(1)%UNITS = 'm s-1 Hz-1'
    meta(1)%FSC = 0.0005
    meta(1)%ENAME = '.usf'
    meta(1)%VMIN = -4.95
    meta(1)%VMAX =  4.95
    meta(1)%VARND = dircom
 
    ! First component
    meta(1)%VARNM='uusf'
    meta(1)%VARNL='eastward spectral variance of surface stokes drift'
    !META(1)%VARNS='eastward_spectral_variance_of_surface_stokes_drift'
    meta(1)%VARNS=''
    meta(1)%VARNC='usf=sqrt(uusf**2+vusf**2)'
    meta(1)%VARNG='eastward_spectral_variance_of_surface_stokes_drift'
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vusf'
    meta(2)%VARNL='northward spectral variance of surface stokes drift'
    !META(2)%VARNS='northward_spectral_variance_of_surface_stokes_drift'
    meta(2)%VARNS=''
    meta(2)%VARNG='northward_spectral_variance_of_surface_stokes_drift'
    meta(2)%VARNC='usf=sqrt(uusf**2+vusf**2)'
    ! IFI=6, IFJ=9, P2L
    meta => group(6)%FIELD(9)%META
    ! Information for spectral microseismic generation data (2nd file)
    meta(1)%FSC = 0.0004
    meta(1)%VARNM='p2l'
    meta(1)%VARNL='base ten logarithm of power spectral density of equivalent surface pressure'
    !META(1)%VARNS='base_ten_logarithm_of_power_spectral_density_of_equivalent_surface_pressure'
    meta(1)%VARNS=''
    meta(1)%VARNG='base_ten_logarithm_of_power_spectral_density_of_equivalent_surface_pressure'
    IF (ncvartype.EQ.2) THEN
      meta(1)%UNITS='log10(Pa2 m2 s+1E-12)'
      meta(1)%VMIN = -12.
      meta(1)%VMAX = 12.
    ELSE
      meta(1)%UNITS='Pa2 m2 s'
      meta(1)%VARNL='power spectral density of equivalent surface pressure'
      !META(1)%VARNS='power_spectral_density_of_equivalent_surface_pressure'
      meta(1)%VARNG='power_spectral_density_of_equivalent_surface_pressure'
      meta(1)%VMIN = 0.
      meta(1)%VMAX = 1.e12
    ENDIF
    meta(1)%VARNC=''
    meta(1)%VARND=''
    meta(1)%ENAME='.p2l'
    ! IFI=6, IFJ=10, TWI
    meta => group(6)%FIELD(10)%META
    meta(1)%FSC = 0.000001
    meta(1)%UNITS = 'm2 s-2'
    meta(1)%ENAME = '.tic'
    meta(1)%VMIN = -0.032
    meta(1)%VMAX =  0.032
    meta(1)%VARND = dircom
 
    ! First component
    meta(1)%VARNL='eastward wave to sea ice stress'
    meta(1)%VARNM='utic'
    !META(1)%VARNS='eastward_wave_to_sea_ice_stress'
    meta(1)%VARNS=''
    meta(1)%VARNG='eastward_wave_to_sea_ice_stress'
    meta(1)%VARNC='two=sqrt(utwo**2+vtwo**2)'
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vtic'
    meta(2)%VARNL='northward wave to sea ice stress'
    !META(2)%VARNS='northward_wave_to_sea_ice_stress'
    meta(2)%VARNS=''
    meta(2)%VARNG='northward_wave_to_sea_ice_stress'
    meta(2)%VARNC='two=sqrt(utwo**2+vtwo**2)'
    ! IFI=6, IFJ=11, FIC
    meta => group(6)%FIELD(11)%META
    meta(1)%FSC    = 0.1
    meta(1)%UNITS  = 'W m-2'
    meta(1)%ENAME  = '.fic'
    meta(1)%VARNM='fic'
    meta(1)%VARNL='wave to sea ice energy flux'
    !META(1)%VARNS='wave_to_sea_ice_energy_flux'
    meta(1)%VARNS=''
    meta(1)%VARNG='wave_to_sea_ice_energy_flux'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 999
    ! IFI=6, IFJ=12, USP
    meta => group(6)%FIELD(12)%META
    meta(1)%UNITS   = 'm s-1'
    meta(1)%FSC    = 0.0005
    meta(1)%ENAME  = '.usp'
    meta(1)%VARND  = dircom
    meta(1)%VMIN = -9.99
    meta(1)%VMAX = 9.98
 
    ! First component
    meta(1)%VARNM='ussp'
    meta(1)%VARNL='eastward partitioned surface stokes drift'
    !META(1)%VARNS='eastward_partitioned_surface_stokes_drift'
    meta(1)%VARNS=''
    meta(1)%VARNG='eastward_partitioned_surface_stokes_drift'
    meta(1)%VARNC='usp=sqrt(ussp**2+vssp**2)'
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vssp'
    meta(2)%VARNL='northward partitioned surface stokes drift'
    !META(2)%VARNS='northward_partitioned_surface_stokes_drift'
    meta(2)%VARNS=''
    meta(2)%VARNG='northward_partitioned_surface_stokes_drift'
    meta(2)%VARNC='usp=sqrt(ussp**2+vssp**2)'
    ! IFI=6, IFJ=13
    meta => group(6)%FIELD(13)%META
    meta(1)%UNITS  = 'Pa'
    meta(1)%FSC    = 0.01
    meta(1)%ENAME  = '.toc'
    meta(1)%VMIN   = -320
    meta(1)%VMAX   =  320
    meta(1)%VARND  = dircom
 
    ! First component
    meta(1)%VARNM='utoc'
    meta(1)%VARNL='eastward total wave to ocean stres'
    meta(1)%VARNS=''
    meta(1)%VARNG=''
    meta(1)%VARNC='toc=sqrt(utoc**2+vtoc**2)'
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vtoc'
    meta(2)%VARNL='northward total wave to ocean stres'
    meta(2)%VARNS=''
    meta(2)%VARNG=''
    meta(2)%VARNC='toc=sqrt(utoc**2+vtoc**2)'
    !
    !---------- GROUP 7 ----------------
    !
    ! IFI=7, IFJ=1, ABR
    meta => group(7)%FIELD(1)%META
    meta(1)%FSC    = 0.01
    meta(1)%ENAME  = '.abr'
    meta(1)%UNITS  = 'm'
    meta(1)%VMIN = -180
    meta(1)%VMAX = 180
    meta(1)%VARND = dircom
 
    ! First component
    meta(1)%VARNM='uabr'
    meta(1)%VARNL='rms of bottom displacement amplitude zonal'
    !META(1)%VARNS='rms_of_bottom_displacement_amplitude_zonal'
    meta(1)%VARNS=''
    meta(1)%VARNG='rms_of_bottom_displacement_amplitude_zonal'
    meta(1)%VARNC='abr=sqrt(uabr**2+vabr**2)'
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vabr'
    meta(2)%VARNL='rms of bottom displacement amplitude meridional'
    !META(2)%VARNS='rms_of_bottom_displacement_amplitude_meridional'
    meta(2)%VARNS=''
    meta(2)%VARNG='rms_of_bottom_displacement_amplitude_meridional'
    meta(2)%VARNC='abr=sqrt(uabr**2+vabr**2)'
    ! IFI=7, IFJ=2, UBR
    meta => group(7)%FIELD(2)%META
    meta(1)%FSC    = 0.01
    meta(1)%ENAME  = '.ubr'
    meta(1)%UNITS  = 'm s-1'
    meta(1)%VMIN = -180
    meta(1)%VMAX = 180
    meta(1)%VARND = dircom
 
    ! First component
    meta(1)%VARNM='uubr'
    meta(1)%VARNL='rms of bottom velocity amplitude zonal'
    !META(1)%VARNS='rms_of_bottom_velocity_amplitude_zonal'
    meta(1)%VARNS=''
    meta(1)%VARNG='rms_of_bottom_velocity_amplitude_zonal'
    meta(1)%VARNC='ubr=sqrt(uubr**2+vubr**2)'
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vubr'
    meta(2)%VARNL='rms of bottom velocity amplitude meridional'
    !META(2)%VARNS='rms_of_bottom_velocity_amplitude_meridional'
    meta(2)%VARNS=''
    meta(2)%VARNG='rms_of_bottom_velocity_amplitude_meridional'
    ! IFI=7, IFJ=3, BED
    meta => group(7)%FIELD(3)%META
    meta(1)%FSC    = 0.001
    meta(1)%UNITS  = 'm'
    meta(1)%ENAME  = '.bed'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 30
    meta(1)%VARND = dircom
 
    ! First component
    meta(1)%VARNM='bed'
    meta(1)%VARNL='bottom roughness'
    !META(1)%VARNS='sea bottom roughness length'
    meta(1)%VARNS=''
    meta(1)%VARNG='ripple_wavelength'
    meta(1)%VARNC='ripple_length=sqrt(ripplex**2+rippley**2)'
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='ripplex'
    meta(2)%VARNL='eastward sea bottom ripple wavelength'
    !META(2)%VARNS='eastward_ripple_wavelength'
    meta(2)%VARNS=''
    meta(2)%VARNG='eastward_ripple_wavelength'
    meta(2)%VARNC='ripple_length=sqrt(ripplex**2+rippley**2)'
 
    ! Third component
    meta(3) = meta(1)
    meta(3)%VARNM='rippley'
    meta(3)%VARNL='northward sea bottom ripple wavelength'
    !META(3)%VARNS='northward_ripple_wavelength'
    meta(3)%VARNS=''
    meta(3)%VARNG='northward_ripple_wavelength'
    meta(3)%VARNC='ripple_length=sqrt(ripplex**2+rippley**2)'
    ! IFI=7, IFJ=4, FBB
    meta => group(7)%FIELD(4)%META
    meta(1)%FSC    = 0.1
    meta(1)%UNITS  = 'W m-2'
    meta(1)%ENAME  = '.fbb'
    meta(1)%VARNM='fbb'
    meta(1)%VARNL='wave dissipation in bbl'
    !META(1)%VARNS='wave_energy_dissipation_in_bottom_boundary_layer'
    meta(1)%VARNS=''
    meta(1)%VARNG='wave_dissipation_in_bbl'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 999
    ! IFI=7, IFJ=5, TBB
    meta => group(7)%FIELD(5)%META
    meta(1)%FSC    = 0.000001
    meta(1)%UNITS  = 'm2 s-2'
    meta(1)%ENAME  = '.tbb'
    meta(1)%VMIN = -0.032
    meta(1)%VMAX =  0.032
    meta(1)%VARND = dircom
 
    ! First component
    meta(1)%VARNM='utbb'
    meta(1)%VARNL='eastward wave to bbl stress'
    !META(1)%VARNS='eastward_wave_to_bottom_boundary_layer_stress'
    meta(1)%VARNS=''
    meta(1)%VARNG='eastward_wave_to_bbl_stress'
    meta(1)%VARNC='tbb=sqrt(utbb**2+vtbb**2)'
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='vtbb'
    meta(2)%VARNL='northward wave to bbl stress'
    !META(2)%VARNS='northward_wave_to_bottom_boundary_layer_stress'
    meta(2)%VARNS=''
    meta(2)%VARNG='northward_wave_to_bbl_stress'
    meta(2)%VARNC='tbb=sqrt(utbb**2+vtbb**2)'
    !
    !---------- GROUP 8 ----------------
    ! IFI=8, IFJ=1, MSS
    meta => group(8)%FIELD(1)%META
    meta(1)%FSC    = 0.00001
    meta(1)%ENAME  = '.mss'
    meta(1)%UNITS  = '1'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 0.3
    meta(1)%VARND = dircom
    WRITE(meta(1)%VARNC,'(A,F8.4,A,F8.4,A)') 'Frequency range ',sig(1)*tpiinv,' to ',sig(nk)*tpiinv,' Hz'
 
    ! First component
    meta(1)%VARNM='mssu'
    meta(1)%VARNL='downwave mean square slope'
    meta(1)%VARNS='sea_surface_wave_mean_square_upwave_slope'
    meta(1)%VARNG='x_mean_square_slope'
    meta(1)%VARNC='mss=mssu+mssc'
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='mssc'
    meta(2)%VARNL='crosswave mean square slope'
    meta(2)%VARNS='sea_surface_wave_mean_square_crosswave_slope'
    meta(2)%VARNG='y_mean_square_slope'
    ! IFI=8, IFJ=2, MSC
    meta => group(8)%FIELD(2)%META
    meta(1)%FSC    = 1e-7
    meta(1)%ENAME  = '.msc'
    meta(1)%UNITS  = '1'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 0.003
    meta(1)%VARND = dircom
 
    ! First component
    meta(1)%VARNM='mscx'
    meta(1)%VARNL='eastward phillips constant'
    !META(1)%VARNS='eastward_phillips_constant'
    meta(1)%VARNS=''
    meta(1)%VARNG='eastward_phillips_constant'
    meta(1)%VARNC='msc=mscx+mscy'
 
    ! Second component
    meta(2) = meta(1)
    meta(2)%VARNM='mscy'
    meta(2)%VARNL='northward phillips constant'
    !META(2)%VARNS='northward_phillips_constant'
    meta(2)%VARNS=''
    meta(2)%VARNG='northward_phillips_constant'
    meta(2)%VARNC='msc=mscx+mscy'
    ! IFI=8, IFJ=3, MSD
    meta => group(8)%FIELD(3)%META
    meta(1)%FSC    = 0.1
    meta(1)%UNITS  = 'degree'
    meta(1)%ENAME  = '.msd'
    meta(1)%VARNM='mssd'
    meta(1)%VARNL='u direction for mss'
    meta(1)%VARNS='sea_surface_mean_square_upwave_slope_direction'
    meta(1)%VARNG='sea_surface_wave_dominant_mean_square_slope_direction'
    WRITE(meta(1)%VARNC,'(A,F8.4,A,F8.4,A)') 'Frequency range ',sig(1)*tpiinv,' to ',sig(nk)*tpiinv,' Hz'
    meta(1)%VARND = dircom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 360
    ! IFI=8, IFJ=4, MCD
    meta => group(8)%FIELD(4)%META
    meta(1)%FSC    = 0.1
    meta(1)%UNITS  = 'degree'
    meta(1)%ENAME  = '.mcd'
    meta(1)%VARNM='mscd'
    meta(1)%VARNL='x direction for msc'
    !META(1)%VARNS='sea_surface_wave_dominant_mean_square_slope_direction_in_highest_frequency'
    meta(1)%VARNS=''
    meta(1)%VARNG='sea_surface_wave_dominant_mean_square_slope_direction_in_highest_frequency'
    meta(1)%VARND = dircom
    meta(1)%VMIN = 0
    meta(1)%VMAX = 360
    ! IFI=8, IFJ=5, QP
    meta => group(8)%FIELD(5)%META
    meta(1)%FSC    = 0.001
    meta(1)%UNITS  = '1'
    meta(1)%ENAME  = '.qp'
    meta(1)%VARNM='qp'
    meta(1)%VARNL='peakedness'
    !META(1)%VARNS='sea_surface_wave_peakedness'
    meta(1)%VARNS=''
    meta(1)%VARNG='wave_peakedness'
    meta(1)%VARNC='Goda wave peakedness parameter'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 32
    ! IFI=8, IFJ=6, QKK
    meta => group(8)%FIELD(6)%META
    meta(1)%FSC    = 0.05
    meta(1)%UNITS  = 'm/rad'
    meta(1)%ENAME  = '.qkk'
    meta(1)%VARNM='qkk'
    meta(1)%VARNL='k-peakedness'
    !META(1)%VARNS='sea_surface_wave_peakedness'
    meta(1)%VARNS=''
    meta(1)%VARNG='wavenumber_peakedness'
    meta(1)%VARNC='2D wavenumber peakedness'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 1600
    ! IFI=8, IFJ=7, SKW
    meta => group(8)%FIELD(7)%META
    meta(1)%FSC    = 0.00001
    meta(1)%UNITS  = '1'
    meta(1)%ENAME  = '.skw'
    meta(1)%VARNM='skw'
    meta(1)%VARNL='skewness'
    !META(1)%VARNS='sea_surface_wave_peakedness'
    meta(1)%VARNS=''
    meta(1)%VARNG='skewness of P(z,sx,sy=0)'
    meta(1)%VARNC='skewness of P(z,sx,sy=0)'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 1
    ! IFI=8, IFJ=8, EMB
    meta => group(8)%FIELD(8)%META
    meta(1)%FSC    = 0.00001
    meta(1)%UNITS  = '1'
    meta(1)%ENAME  = '.emb'
    meta(1)%VARNM='emb'
    meta(1)%VARNL='EM-bias'
    !META(1)%VARNS='sea_surface_wave_peakedness'
    meta(1)%VARNS=''
    meta(1)%VARNG='EM bias coefficient'
    meta(1)%VARNC='EM bias coefficient'
    meta(1)%VMIN = -1
    meta(1)%VMAX = 1    
    ! IFI=8, IFJ=7, SKW
    meta => group(8)%FIELD(9)%META
    meta(1)%FSC    = 0.00001
    meta(1)%UNITS  = '1'
    meta(1)%ENAME  = '.emc'
    meta(1)%VARNM='emc'
    meta(1)%VARNL='trackerbias'
    !META(1)%VARNS='sea_surface_wave_peakedness'
    meta(1)%VARNS=''
    meta(1)%VARNG='tracker bias coefficient'
    meta(1)%VARNC='tracker bias coefficient'
    meta(1)%VMIN = -1
    meta(1)%VMAX = 1    !
    !
    !---------- GROUP 9 ----------------
    !
    ! IFI=9, IFJ=1, DTD
    meta => group(9)%FIELD(1)%META
    meta(1)%FSC    = 0.1
    meta(1)%UNITS  = 'min.'
    meta(1)%ENAME  = '.dtd'
    meta(1)%VARNM='dtd'
    meta(1)%VARNL='dynamic time step'
    !META(1)%VARNS='dynamic_time_step'
    meta(1)%VARNS=''
    meta(1)%VARNG='dynamic_time_step'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 3200
    ! IFI=9, IFJ=2, FC
    meta => group(9)%FIELD(2)%META
    meta(1)%FSC    = 0.001
    meta(1)%UNITS  = 's-1'
    meta(1)%ENAME  = '.fc'
    meta(1)%VARNM='fc'
    meta(1)%VARNL='cut off frequency'
    !META(1)%VARNS='cut_off_frequency'
    meta(1)%VARNS=''
    meta(1)%VARNG='cut_off_frequency'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 8
    ! IFI=9, IFJ=3, CFX
    meta => group(9)%FIELD(3)%META
    meta(1)%FSC    = 0.01
    meta(1)%UNITS  = '1'
    meta(1)%ENAME  = '.cfx'
    meta(1)%VARNM='cfx'
    meta(1)%VARNL='maximum cfl for spatial advection'
    !META(1)%VARNS='maximum_cfl_for_spatial_advection'
    meta(1)%VARNS=''
    meta(1)%VARNG='maximum_cfl_for_spatial_advection'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 320
    ! IFI=9, IFJ=4, CFD
    meta => group(9)%FIELD(4)%META
    meta(1)%FSC    = 0.01
    meta(1)%UNITS  = '1'
    meta(1)%ENAME  = '.cfd'
    meta(1)%VARNM='cfd'
    meta(1)%VARNL='maximum cfl for direction advection'
    !META(1)%VARNS='maximum_cfl_for_direction_advection'
    meta(1)%VARNS=''
    meta(1)%VARNG='maximum_cfl_for_direction_advection'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 320
    ! IFI=9, IFJ=5, CFK
    meta => group(9)%FIELD(5)%META
    meta(1)%FSC    = 0.01
    meta(1)%UNITS  = '1'
    meta(1)%ENAME  = '.cfk'
    meta(1)%VARNM='cfk'
    meta(1)%VARNL='maximum cfl for frequency advection'
    !META(1)%VARNS='maximum_cfl_for_frequency_advection'
    meta(1)%VARNS=''
    meta(1)%VARNG='maximum_cfl_for_frequency_advection'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 320
    !
    ! ------ Group 10 (User defined) -------
    !
    ! IFI=10, IFJ=1
    meta => group(10)%FIELD(1)%META
    meta(1)%FSC    = 0.1
    meta(1)%UNITS  = 'm'
    meta(1)%VMIN = 0
    meta(1)%VMAX = 0
    WRITE (meta(1)%ENAME,'(A2,I2.2)') '.u'
    WRITE (meta(1)%VARNM,'(A1,I2.2)') 'u'
    WRITE (meta(1)%VARNL,'(A12,I2.2)') 'User_defined'
    WRITE (meta(1)%VARNS,'(A12,I2.2)') 'User_defined'
    WRITE (meta(1)%VARNG,'(A12,I2.2)') 'user_defined'
    !
  END SUBROUTINE default_meta
  !/ ------------------------------------------------------------------- /
 
END MODULE w3ounfmetamd