w3ounfmetamd Module Reference

Manages user configurable netCDF meta-data for ww3_ounf program. More...

Data Types
type	field_t
	Type for storage of meta data aggregated by component (NFIELD) More...
type	group_t
	Type for storage of meta data aggregated by field (IFI) More...
type	meta_t
	Type for storing WW3 netCDF metadata for a variable. More...
type	part_tmpl_t
	Type for storing partitioned parameter template strings. More...

Functions/Subroutines
subroutine	init_meta (VEC)
	Allocates space for the META_T arrays and sets some defaults. More...
subroutine	teardown_meta ()
	De-allocates memory used for the META_T arrays. More...
subroutine	next_line (NDMI, BUF, ILINE, EOF, NEW_SECTION)
	Reads the next valid line from the user meta input file. More...
subroutine	notabs (STR)
	Replaces tab characters in a string with a space. More...
character(len(str)) function	replace_char (STR, C, REP)
	Replaces single characters in a string. More...
subroutine	read_meta ()
	Reads meta data entries from the ountmeta.inp file. More...
subroutine	decode_header (BUF, ILINE, IFI, IFJ, IFC)
	Decode the META header line. More...
subroutine	read_meta_pairs (NDMI, META, ILINE)
	Reads in attribute name/value pairs and updates the relevant values in the META type. More...
subroutine	get_attval_type (BUF, ILINE, ATTV, ATT_TYPE)
	Gets the attribute value and optional variable type from the passed in string. More...
subroutine	read_freeform_meta_list (NDMI, ILINE, METALIST)
	Reads in freeform attribute name/value pairs. More...
subroutine	read_crs_meta (NDMI, ILINE)
	Reads in metadata for the coordinate reference system (CRS). More...
subroutine	default_crs_meta ()
	Set up a default coordinate reference system for the grid. More...
type(meta_t) function	getmeta (IFI, IFJ, ICOMP, IPART)
	Get the meta data for a particular field. More...
subroutine	read_part_tmpl (NDMI, ILINE)
	Reads in a TEMPLATE section from file. More...
subroutine	print_part_tmpl ()
	Prints the patition templates to screen (for debug use). More...
subroutine	add_partno (META, IPART)
	Adds partition number to meta-data. More...
subroutine	partno_string_sub (INSTR, IPART)
	Performs string substition of placeholder strings with partition number specfic values. More...
subroutine	write_meta (NCID, VARID, META, ERR)
	Writes the meta-data entries for a variable. More...
subroutine	write_global_meta (NCID, ERR)
	Writes the user meta-data entries for the global attributes. More...
subroutine	write_freeform_meta_list (NCID, VARID, METALIST, ERR)
	Writes the freeform user meta-data entries for a NetCDF variable. More...
subroutine	print_meta (META)
	Writes meta-data to the screen - for debugging purposes. More...
type(meta_t) function	meta_deep_copy (META)
	Performs "deep" copy of a META_T type. More...
subroutine	default_meta ()
	Populates the default meta data for ww3_ounf output fields. More...

Variables
character(len= *), parameter	fn_meta = "ounfmeta.inp"
	Meta-data input filename. More...
character(len= *), parameter	ipart_token = "<IPART>"
	String token for integer partition number. More...
character(len= *), parameter	spart_token = "<SPART>"
	String token for partition descriptive string (space separated) More...
character(len= *), parameter	spart_token_ = "<SPART_>"
	String token for partition descriptive string (underscore separated) More...
type(group_t), dimension(:), allocatable	group
	Storage for meta data aggregated by group (IFJ) More...
type(meta_list_t)	global_meta
	Storage for the Global meta data (free form) More...
logical	fl_default_gbl_meta = .TRUE.
	Flag for using default (true) or user-defined (false) global meta data. More...
character(len=128)	crs_name = ''
	Coordinate reference system (CRS) name. More...
type(meta_list_t)	crs_meta
	Meta data list for CRS. More...
logical	crs_is_default = .FALSE.
	True if CRS set by this module. More...
character(len=256)	coords_attr = ''
	"coordinates" attribute - for defining auxiliary coordinates (for all variables) More...
type(part_tmpl_t), pointer	part_tmpl
	User-defined partitionted paratmeters template strings. More...
integer	ncvartype
	NetCDF variable type (2=int, 3=real, 4=depends) More...
character(len=30)	dircom
	Directional convention comment. More...
character(len=128)	partcom
	Partitioning method comment. More...
character(len=15), dimension(5)	snamep
	Part. More...
logical	flrtd = .FALSE.
	Flag for rototed pole grid. More...

Detailed Description

Manages user configurable netCDF meta-data for ww3_ounf program.

Default netCDF meta data is provided for each WW3 output variable and is stored intentally via the META_T type. The meta values are grouped by component (max 3), field (IFI) and group (IFJ).

The user can override this meta data via an input text file with the filename ounfmeta.inp.

Entries in the file are formatted as follows:

    META [ IFI [ IFJ ]  |  FLDID ]   [ IFC ]
      attr_name = attr_value
      attr_name = attr_value
      extra_attr = extra_value [type]
    ... repeated as many times as required.

An output field is selected using the META keyword followed by either an [IFI, IFJ] integer pair or a FieldID string. Optionally, either form may be followed by an integer value to select the component in multi-component fields (such as wind).

Blank lines and comments lines (starting with $) are ignored.

attr_name is the netCDF attribute name that you wish to override. This can be one of the following:

• "varnm"
• "ename"
• "standard_name", or "varns"
• "long_name" or "varnl"
• "globwave_name" or "varng"
• "direction_reference", "dir_ref" or "varnd"
• "comment" or "varnc"
• "units"
• "valid_min" or "vmin"
• "valid_max" or "vmax"
• "scale_factor" or "fsc"

Any other attribute name is assumed to be an optional "extra" attribute. This extra attribute can take an optional "type" keyworkd to specify the variable tpye of the metadata. If no type is supplied, it defaults to a characer type. Valid types are one of ["c", "r", "i"] for character/string, real/float or integer values respectively.

Global meta data can be specified with a special "META global" line:

    META global
      extra_attr = extra_value [type]
      extra_attr = extra_value [type]

A "coordinate reference system" (CRS) can be specified for all output fields using the "CRS" keyword. As a minimum, the "grid_mapping_name" attribute must be specified. If the CRS section is defined, all output fields will have a "grid_mapping" attribute added referencing the CRS variable. "crs_vaname" will be created as a scalar NF90_CHAR variable in the output file.

    CRS <crs_varname>
      grid_mapping_name = <mapping name>
      attr = value
      attr = value

Note: ALL keywords and "Field Name ID" strings (e.g. HS) are case insensitive. All netCDF attribute names are case sensitive.

Partitioned outputs are handles slightly differently; one meta data entry is used for all partitions of a field. The metadata is made specific to a particular partition via template strings. There are two built-in template strings: SPART and IPART. These provide a "string" description (e.g. "wind sea", "primary swell", etc) or an integer partition number. These can be references in the meta data using the template name surrounded by < .. >, e.g. <SPART>

It is also possible to supply user defined partitioned parameter template strings in the ounfmeta.inp file using the TEMPLATE keyword, as below:

    TEMPLATE <template-name>
      String for partition 0
      String for partition 1
      String for partition 2
      String for partition 3
      ... etc

Specifying the <template-name> with a trailing underscore will provide an underscore seperated (_) string, rather than space seperated.

Example ounfmeta.inp file:

    $ Lines starting with dollars are comments.
    $ The line starts a meta-data section for the depth field
    META DPT
      standard_name = depth
      long_name = "can be quoted string"
      comment = or an unquoted string
      vmax = 999.9

    $ Next one is HSig (group 2, field 1)
    META 2 1
      varns = "sig. wave height"
      varnl = "this is long name"

    $ Next one is second component of wind. It also sets an
    $ "extra" meta data value (height - a float)
    META WND 2
      standard_name = "v-wind"
      height = 10.0 "r"

    $ User defined partitioned parameters template strings:
    TEMPLATE PARTSTR
      wind wave
      primary swell
      secondary swell

    $ Use partition templates in partitioned Hs field:
    $ (SPART and IPART are built-in)
    META PHS
      standard_name = "<SPART_>_sigificant_wave_height"
      long_name = "<PARTSTR>"
      partition_number = "<IPART>"

    $ Coordinate reference system:
    CRS crs
      grid_mapping_name = "latitude_longitude"
      semi_major_axis = 6371000.0 f
      inverse_flattening = 0 f

    $ Global metadata:
    META global
      institution = UKMO
      comment "space seperated strings should be quoted" c
      version = 1.0 r

Author

Chris Bunney

Date

02-Nov-2020

Change log

Date	Ver	Comments
02-Nov-2020	7.12	Creation
26-Jan-2021	7.12	Added Tp and alternative dir/mag metadata for directional fields.
16-Dec-2020	7.12	Added user partition templates and coordinate reference system.
02-Feb-2021	7.12	Improved partitioned parameter template string implementation.
22-Mar-2021	7.12	Add extra coupling fields
02-Sep-2021	7.12	Add coordinates attribute

Function/Subroutine Documentation

add_partno()

subroutine w3ounfmetamd::add_partno	(	type(meta_t), intent(inout)	META,
		integer, intent(in)	IPART
	)

Adds partition number to meta-data.

Replaces all instances of "<IPART>" in the provided meta data with the partition number IPART.

Parameters

[in]	META	Meta data type
[in]	IPART	Partition number

Author

Chris Bunney

Date

02-Nov-2020

Definition at line 1898 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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
 

References partno_string_sub().

Referenced by getmeta().

decode_header()

subroutine w3ounfmetamd::decode_header	(	character(*), intent(in)	BUF,
		integer, intent(in)	ILINE,
		integer, intent(out)	IFI,
		integer, intent(out)	IFJ,
		integer, intent(out)	IFC
	)

Decode the META header line.

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.

Parameters

[in]	BUF	Input header string (without leading META tag)
[in]	ILINE	Line number (for error reporting)
[out]	IFI	Output group number
[out]	IFJ	Output field number
[out]	IFC	Component number (defaults to 1)

Author

Chris Bunney

Date

02-Feb-2021

Definition at line 906 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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 /)
    !

References fl_default_gbl_meta, fn_meta, and w3iogomd::w3fldtoij().

Referenced by read_meta().

default_crs_meta()

subroutine w3ounfmetamd::default_crs_meta

Set up a default coordinate reference system for the grid.

The default coordinate reference system (CRS) will be defined based on the type of grid the model is formulated on, e.g. regular lat-lon, rotated pole, etc.

Remarks

See "Grid Mappings" section of CF conventions:

• https://cfconventions.org/Data/cf-conventions/cf-conventions-1.7/build/ch05s06.html
• https://cfconventions.org/Data/cf-conventions/cf-conventions-1.7/build/apf.html

Author

Chris Bunney

Date

25-May-2021

Definition at line 1558 of file w3ounfmetamd.F90.

    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
 

References crs_is_default, crs_meta, crs_name, and flrtd.

Referenced by init_meta().

default_meta()

subroutine w3ounfmetamd::default_meta

Populates the default meta data for ww3_ounf output fields.

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.

Author

Chris Bunney

Date

22-Mar-2021

Definition at line 2499 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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'
    !

References dircom, group, ipart_token, ncvartype, partcom, and spart_token_.

Referenced by init_meta().

get_attval_type()

subroutine w3ounfmetamd::get_attval_type	(	character(*), intent(in)	BUF,
		integer, intent(in)	ILINE,
		character(*), intent(out)	ATTV,
		character(*), intent(out)	ATT_TYPE
	)

Gets the attribute value and optional variable type from the passed in string.

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.

Parameters

[in]	BUF	Input string to process
[in]	ILINE	Line number (for error reporting)
[out]	ATTV	Attribute value
[out]	ATT_TYPE	Attribute type

Author

Chris Bunney

Date

09-Nov-2020

Definition at line 1225 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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 /)
    !

References fn_meta.

Referenced by read_freeform_meta_list(), and read_meta_pairs().

getmeta()

type(meta_t) function w3ounfmetamd::getmeta	(	integer, intent(in)	IFI,
		integer, intent(in)	IFJ,
		integer, intent(in), optional	ICOMP,
		integer, intent(in), optional	IPART
	)

Get the meta data for a particular field.

The required field is specified using the group (IFI) and field (IFJ) index. Optionally, the component (ICOMP) and partition (IPART) numbers can be specified for vector/tensor or partitioned parameter fields. If not specified, these default to 1. 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.

Parameters

[in]	IFI	Output group number
[in]	IFJ	Output field number
[in]	ICOMP	Component number (defaults to 1)
[in]	IPART	Partition number (defaults to 1)

Author

Chris Bunney

Date

02-Nov-2020

Definition at line 1610 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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'      / )
    !

References add_partno(), group, and meta_deep_copy().

Referenced by w3ounf().

init_meta()

subroutine w3ounfmetamd::init_meta

(

logical, intent(in), optional

VEC

)

Allocates space for the META_T arrays and sets some defaults.

By default, directional fields will be set up to output a magnitude and direction field. Alternatively, if VEC is set to True, then u/v vectors will be generated instead.

Note

- vector output is currently only implemented for the "current" and "wind" fields.

Parameters

VEC	Output vectors for directional fields rather than direction/magnitude.

Author

Chris Bunney

Date

09-Mar-2020

Definition at line 299 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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
 

References default_crs_meta(), default_meta(), dircom, file(), w3gdatmd::flagunr, flrtd, w3odatmd::fnmpre, group, w3odatmd::noge, w3odatmd::nogrp, part_tmpl, partcom, w3gdatmd::polat, w3odatmd::ptfcut, w3odatmd::ptmeth, read_meta(), w3smcomd::smcotype, and snamep.

Referenced by w3ounf().

meta_deep_copy()

type(meta_t) function w3ounfmetamd::meta_deep_copy

(

type(meta_t), intent(in)

META

)

Performs "deep" copy of a META_T type.

A "deep" copy ensures that the linked list data in the EXTRA field is copied, rather than just copying the pointer.

Calls copy_meta_list() internally to copy the EXTRA linked list.

Returns

A new META_T variable

Parameters

[in]

META

META data structure to copy

Author

Chris Bunney

Date

16-Dec-2020

Definition at line 2448 of file w3ounfmetamd.F90.

    !/                  +-----------------------------------+
    !/                  | 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)
 

Referenced by getmeta().

next_line()

subroutine w3ounfmetamd::next_line	(	integer, intent(in)	NDMI,
		character(*), intent(out)	BUF,
		integer, intent(inout)	ILINE,
		logical, intent(out)	EOF,
		logical, intent(out), optional	NEW_SECTION
	)

Reads the next valid line from the user meta input file.

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.

Parameters

[in]	NDMI	Unit number of input file
[out]	BUF	Next input line read from file
[in,out]	ILINE	Line number of file
[out]	EOF	True if end-of-file is reached.
[out]	NEW_SECTION	True if new section marker found

Author

Chris Bunney

Date

09-Nov-2020

Definition at line 475 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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 /)
    !

References w3servmd::extcde(), fn_meta, w3odatmd::ndse, notabs(), and w3servmd::str_to_upper().

Referenced by read_freeform_meta_list(), read_meta(), read_meta_pairs(), and read_part_tmpl().

notabs()

subroutine w3ounfmetamd::notabs

(

character(*), intent(inout)

STR

)

Replaces tab characters in a string with a space.

Remarks

Assumes ASCII encoding (Tab character is ASCII value 9)

Parameters

[in,out]

STR

Character string to process

Author

Chris Bunney

Date

02-Nov-2020

Definition at line 615 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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
 

Referenced by next_line().

partno_string_sub()

subroutine w3ounfmetamd::partno_string_sub	(	character(len=*), intent(inout)	INSTR,
		integer, intent(in)	IPART
	)

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.

Parameters

[in,out]	INSTR	Input string
[in]	IPART	Partition number

Author

Chris Bunney

Date

02-Nov-2020

Definition at line 1977 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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.'           /)
 

References ipart_token, part_tmpl, replace_char(), snamep, spart_token, and spart_token_.

Referenced by add_partno().

print_meta()

subroutine w3ounfmetamd::print_meta

(

type(meta_t), intent(in)

META

)

Writes meta-data to the screen - for debugging purposes.

Parameters

[in]

META

Meta data type

Author

Chris Bunney

Date

09-Nov-2020

Definition at line 2383 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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
 

print_part_tmpl()

subroutine w3ounfmetamd::print_part_tmpl

Prints the patition templates to screen (for debug use).

Author

Chris Bunney

Date

04-Dec-2020

Definition at line 1844 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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*,'=============='

References part_tmpl.

read_crs_meta()

subroutine w3ounfmetamd::read_crs_meta	(	integer, intent(in)	NDMI,
		integer, intent(inout)	ILINE
	)

Reads in metadata for the coordinate reference system (CRS).

The "grid_mapping_name" must be supplied as an attribute.

Parameters

[in]	NDMI	Unit number of metadata input file
[in,out]	ILINE	Current line number in file

Author

Chris Bunney

Date

07-Dec-2020

Definition at line 1451 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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 /)
    !
 

References crs_is_default, crs_meta, crs_name, fn_meta, and read_freeform_meta_list().

Referenced by read_meta().

read_freeform_meta_list()

subroutine w3ounfmetamd::read_freeform_meta_list	(	integer, intent(in)	NDMI,
		integer, intent(inout)	ILINE,
		type(meta_list_t), intent(inout)	METALIST
	)

Reads in freeform attribute name/value pairs.

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.

Parameters

[in]	NDMI	Unit number of metadata input file
[in,out]	ILINE	Current line number in file
[in,out]	METALIST	A META_LIST_T object to append to

Author

Chris Bunney

Date

16-Dec-2020

Definition at line 1338 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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 /)
    !

References fn_meta, get_attval_type(), and next_line().

Referenced by read_crs_meta(), and read_meta().

read_meta()

subroutine w3ounfmetamd::read_meta

Reads meta data entries from the ountmeta.inp file.

This is the main entry routine for parsing the ounfmeta.inp file. Values read from the file will be used to update or add to the default meta data values set in the default_meta() subroutine.

Author

Chris Bunney

Date

26-Jan-2021

Definition at line 729 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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 /)
    !

References decode_header(), file(), fn_meta, global_meta, group, next_line(), read_crs_meta(), read_freeform_meta_list(), read_meta_pairs(), and read_part_tmpl().

Referenced by init_meta().

read_meta_pairs()

subroutine w3ounfmetamd::read_meta_pairs	(	integer, intent(in)	NDMI,
		type(meta_t), intent(inout), pointer	META,
		integer, intent(inout)	ILINE
	)

Reads in attribute name/value pairs and updates the relevant values in the META type.

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)

Parameters

[in]	NDMI	Unit number of metadata input file
[out]	META	Pointer to META type
[in,out]	ILINE	Current line number in file

Author

Chris Bunney

Date

09-11-2020

Definition at line 1043 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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 /)
    !

References fn_meta, get_attval_type(), and next_line().

Referenced by read_meta().

read_part_tmpl()

subroutine w3ounfmetamd::read_part_tmpl	(	integer, intent(in)	NDMI,
		integer, intent(inout)	ILINE
	)

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

Parameters

[in,out]	NDMI	Unit number
[in,out]	ILINE	Line number

Author

Chris Bunney

Date

04-Dec-2020

Definition at line 1717 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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 /)
    !

References fn_meta, next_line(), and part_tmpl.

Referenced by read_meta().

replace_char()

character(len(str)) function w3ounfmetamd::replace_char	(	character(*)	STR,
		character	C,
		character	REP
	)

Replaces single characters in a string.

Returns

A new string

Parameters

[in]	STR	Character string to process
[in]	C	Character to search for
[in]	REP	Character to substitute

Author

Chris Bunney

Date

02-Feb-2021

Definition at line 674 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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
 

Referenced by partno_string_sub().

teardown_meta()

subroutine w3ounfmetamd::teardown_meta

De-allocates memory used for the META_T arrays.

Author

Chris Bunney

Date

09-Nov-2020

Definition at line 415 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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)
 

References crs_meta, w3metamd::del_meta_list(), global_meta, group, w3odatmd::noge, and w3odatmd::nogrp.

Referenced by w3ounf().

write_freeform_meta_list()

subroutine w3ounfmetamd::write_freeform_meta_list	(	integer, intent(in)	NCID,
		integer, intent(in)	VARID,
		type(meta_list_t), intent(in)	METALIST,
		integer, intent(out)	ERR
	)

Writes the freeform user meta-data entries for a NetCDF variable.

Parameters

[in,out]	NCID	NetCDF file ID
[in,out]	VARID	NetCDF variable ID
[in]	METALIST	META_LIST_T object to write
[out]	ERR	Error value

Author

Chris Bunney

Date

16-Dec-2020

Definition at line 2293 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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         / )
    !

Referenced by w3ounf(), write_global_meta(), and write_meta().

write_global_meta()

subroutine w3ounfmetamd::write_global_meta	(	integer, intent(in)	NCID,
		integer, intent(out)	ERR
	)

Writes the user meta-data entries for the global attributes.

Global meta-data is stored as a meta-data list, so this is essentially a convenience/legacy function that calls the write_freeform_meta_list() subroutine.

Parameters

[in]	NCID	NetCDF file ID
[out]	ERR	Error value

Author

Chris Bunney

Date

09-Nov-2020

Definition at line 2246 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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)

References global_meta, and write_freeform_meta_list().

Referenced by w3ounf().

write_meta()

subroutine w3ounfmetamd::write_meta	(	integer, intent(in)	NCID,
		integer, intent(in)	VARID,
		type(meta_t), intent(in)	META,
		integer, intent(out)	ERR
	)

Writes the meta-data entries for a variable.

Attribute pairs defined in META are written to the netCDF variable specificed in the VARID handle.

There are two stages to the write - first all "mandatory" or "pre-defined" attributes are written out (those defined in the META_T type). Secondly, if there is any user-defined "extra" freeform meta data defined, this is written out via a separate call to write_freeform_meta_list().

Parameters

[in,out]	NCID	NetCDF file ID
[in,out]	VARID	NetCDF variable ID
[in]	META	Meta data type
[out]	ERR	Error value

Author

Chris Bunney

Date

02-Nov-2020

Definition at line 2115 of file w3ounfmetamd.F90.

    !/
    !/                  +-----------------------------------+
    !/                  | 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
    !

References coords_attr, crs_name, ncvartype, and write_freeform_meta_list().

Referenced by w3ounf().

Variable Documentation

coords_attr

character(len=256) w3ounfmetamd::coords_attr = ''

"coordinates" attribute - for defining auxiliary coordinates (for all variables)

Definition at line 257 of file w3ounfmetamd.F90.

  CHARACTER(LEN=256) :: COORDS_ATTR = ''

Referenced by w3ounf(), and write_meta().

crs_is_default

logical w3ounfmetamd::crs_is_default = .FALSE.

True if CRS set by this module.

Definition at line 254 of file w3ounfmetamd.F90.

  LOGICAL            :: CRS_IS_DEFAULT = .false. 

Referenced by default_crs_meta(), and read_crs_meta().

crs_meta

type(meta_list_t) w3ounfmetamd::crs_meta

Meta data list for CRS.

Definition at line 253 of file w3ounfmetamd.F90.

  TYPE(META_LIST_T)  :: CRS_META

Referenced by default_crs_meta(), read_crs_meta(), teardown_meta(), and w3ounf().

crs_name

character(len=128) w3ounfmetamd::crs_name = ''

Coordinate reference system (CRS) name.

Definition at line 252 of file w3ounfmetamd.F90.

  CHARACTER(LEN=128) :: CRS_NAME = ''

Referenced by default_crs_meta(), read_crs_meta(), w3ounf(), and write_meta().

dircom

character(len=30) w3ounfmetamd::dircom

Directional convention comment.

Definition at line 272 of file w3ounfmetamd.F90.

  CHARACTER(LEN=30)  :: DIRCOM

Referenced by default_meta(), and init_meta().

fl_default_gbl_meta

logical w3ounfmetamd::fl_default_gbl_meta = .TRUE.

Flag for using default (true) or user-defined (false) global meta data.

Definition at line 249 of file w3ounfmetamd.F90.

  LOGICAL           :: FL_DEFAULT_GBL_META = .true.

Referenced by decode_header(), and w3ounf().

flrtd

logical w3ounfmetamd::flrtd = .FALSE.

Flag for rototed pole grid.

Definition at line 279 of file w3ounfmetamd.F90.

  LOGICAL            :: FLRTD = .false. 

Referenced by default_crs_meta(), and init_meta().

fn_meta

character(len=*), parameter w3ounfmetamd::fn_meta = "ounfmeta.inp"

Meta-data input filename.

Definition at line 199 of file w3ounfmetamd.F90.

  CHARACTER(LEN=*), PARAMETER :: FN_META = "ounfmeta.inp"

Referenced by decode_header(), get_attval_type(), next_line(), read_crs_meta(), read_freeform_meta_list(), read_meta(), read_meta_pairs(), and read_part_tmpl().

global_meta

type(meta_list_t) w3ounfmetamd::global_meta

Storage for the Global meta data (free form)

Definition at line 246 of file w3ounfmetamd.F90.

  TYPE(META_LIST_T) :: GLOBAL_META

Referenced by read_meta(), teardown_meta(), and write_global_meta().

group

type(group_t), dimension(:), allocatable w3ounfmetamd::group

Storage for meta data aggregated by group (IFJ)

Definition at line 243 of file w3ounfmetamd.F90.

  TYPE(GROUP_T), ALLOCATABLE :: GROUP(:)

Referenced by default_meta(), getmeta(), init_meta(), read_meta(), and teardown_meta().

ipart_token

character(len=*), parameter w3ounfmetamd::ipart_token = "<IPART>"

String token for integer partition number.

Definition at line 202 of file w3ounfmetamd.F90.

  CHARACTER(LEN=*), PARAMETER :: IPART_TOKEN  = "<IPART>"

Referenced by default_meta(), and partno_string_sub().

ncvartype

integer w3ounfmetamd::ncvartype

NetCDF variable type (2=int, 3=real, 4=depends)

Definition at line 271 of file w3ounfmetamd.F90.

  INTEGER            :: NCVARTYPE

Referenced by default_meta(), w3ounf(), and write_meta().

part_tmpl

type(part_tmpl_t), pointer w3ounfmetamd::part_tmpl

User-defined partitionted paratmeters template strings.

Definition at line 269 of file w3ounfmetamd.F90.

  TYPE(PART_TMPL_T), POINTER :: PART_TMPL

Referenced by init_meta(), partno_string_sub(), print_part_tmpl(), and read_part_tmpl().

partcom

character(len=128) w3ounfmetamd::partcom

Partitioning method comment.

Definition at line 273 of file w3ounfmetamd.F90.

  CHARACTER(LEN=128) :: PARTCOM

Referenced by default_meta(), and init_meta().

snamep

character(len=15), dimension(5) w3ounfmetamd::snamep

Part.

standard name templates

Definition at line 274 of file w3ounfmetamd.F90.

  CHARACTER(LEN=15)  :: SNAMEP(5)

Referenced by init_meta(), and partno_string_sub().

spart_token

character(len=*), parameter w3ounfmetamd::spart_token = "<SPART>"

String token for partition descriptive string (space separated)

Definition at line 205 of file w3ounfmetamd.F90.

  CHARACTER(LEN=*), PARAMETER :: SPART_TOKEN  = "<SPART>"

Referenced by partno_string_sub().

spart_token_

character(len=*), parameter w3ounfmetamd::spart_token_ = "<SPART_>"

String token for partition descriptive string (underscore separated)

Definition at line 208 of file w3ounfmetamd.F90.

  CHARACTER(LEN=*), PARAMETER :: SPART_TOKEN_ = "<SPART_>"

Referenced by default_meta(), and partno_string_sub().