MWwaterCoeff_Define.f90

Go to the documentation of this file.

!
! MWwaterCoeff_Define
!
! Module defining the MWwaterCoeff object.
!
!
! CREATION HISTORY:
!       Written by:     Paul van Delst, 10-Nov-2011
!                       paul.vandelst@noaa.gov

MODULE mwwatercoeff_define

  ! -----------------
  ! Environment setup
  ! -----------------
  ! Module use
  USE type_kinds           , ONLY: fp, long, double
  USE message_handler      , ONLY: success, failure, information, display_message
  USE compare_float_numbers, ONLY: OPERATOR(.equalto.)
  USE file_utility         , ONLY: file_open, file_exists
  USE binary_file_utility  , ONLY: open_binary_file      , &
                                   writegatts_binary_file, &
                                   readgatts_binary_file
  USE fitcoeff_define
  USE mwwaterlut_define
  ! Disable implicit typing
  IMPLICIT NONE


  ! ------------
  ! Visibilities
  ! ------------
  ! Everything private by default
  PRIVATE
  ! Datatypes
  PUBLIC :: mwwatercoeff_type
  ! Operators
  PUBLIC :: OPERATOR(==)
  ! Procedures
  PUBLIC :: mwwatercoeff_associated
  PUBLIC :: mwwatercoeff_destroy
  PUBLIC :: mwwatercoeff_create
  PUBLIC :: mwwatercoeff_inspect
  PUBLIC :: mwwatercoeff_validrelease
  PUBLIC :: mwwatercoeff_info
  PUBLIC :: mwwatercoeff_defineversion
  PUBLIC :: mwwatercoeff_readfile
  PUBLIC :: mwwatercoeff_writefile


  ! ---------------------
  ! Procedure overloading
  ! ---------------------
  INTERFACE OPERATOR(==)
    MODULE PROCEDURE mwwatercoeff_equal
  END INTERFACE OPERATOR(==)


  ! -----------------
  ! Module parameters
  ! -----------------
  CHARACTER(*), PARAMETER :: module_version_id = &
    '$Id: MWwaterCoeff_Define.f90 60152 2015-08-13 19:19:13Z paul.vandelst@noaa.gov $'
  ! Release and version
  INTEGER, PARAMETER :: mwwatercoeff_release = 1  ! This determines structure and file formats.
  INTEGER, PARAMETER :: mwwatercoeff_version = 1  ! This is just the default data version.
  ! Close status for write errors
  CHARACTER(*), PARAMETER :: write_error_status = 'DELETE'
  ! Data indicators
  INTEGER, PARAMETER :: data_missing = 0
  INTEGER, PARAMETER :: data_present = 1
  ! String lengths
  INTEGER,  PARAMETER :: ml = 256 ! Message length
  INTEGER,  PARAMETER :: sl =  80 ! String length


  ! ---------------------------------
  ! MWwaterCoeff data type definition
  ! ---------------------------------
  !:tdoc+:
  TYPE :: mwwatercoeff_type
    ! Allocation indicator
    LOGICAL :: is_allocated = .false.
    ! Release and version information
    INTEGER(Long) :: release = mwwatercoeff_release
    INTEGER(Long) :: version = mwwatercoeff_version
    ! Derived type components
    TYPE(fitcoeff_1d_type) :: fccoeff     ! Foam coverage          fitting coefficients
    TYPE(fitcoeff_1d_type) :: frcoeff     ! Foam reflectivity      fitting coefficients
    TYPE(fitcoeff_3d_type) :: rccoeff     ! Reflection correction  fitting coefficients
    TYPE(fitcoeff_3d_type) :: azcoeff     ! Azimuth emissivity     fitting coefficients
    TYPE(fitcoeff_1d_type) :: ssccoeff    ! Small-scale correction fitting coefficients
    TYPE(fitcoeff_3d_type) :: lsccoeff    ! Large-scale correction fitting coefficients
    TYPE(mwwaterlut_type)  :: lut         ! Emissivity look-up table
  END TYPE mwwatercoeff_type
  !:tdoc-:


CONTAINS


!################################################################################
!################################################################################
!##                                                                            ##
!##                           ## PUBLIC PROCEDURES ##                          ##
!##                                                                            ##
!################################################################################
!################################################################################

!--------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       MWwaterCoeff_Associated
!
! PURPOSE:
!       Pure function to test the status of the allocatable components
!       of the MWwaterCoeff structure.
!
! CALLING SEQUENCE:
!       Status = MWwaterCoeff_Associated( MWwaterCoeff )
!
! OBJECTS:
!       MWwaterCoeff:  Structure which is to have its member's
!                      status tested.
!                      UNITS:      N/A
!                      TYPE:       MWwaterCoeff_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN)
!
! FUNCTION RESULT:
!       Status:        The return value is a logical value indicating the
!                      status of the components.
!                       .TRUE.  - if ALL of the MWwaterCoeff allocatable members
!                                 are in use.
!                       .FALSE. - if ANY of the MWwaterCoeff allocatable members
!                                 are not in use.
!                      UNITS:      N/A
!                      TYPE:       LOGICAL
!                      DIMENSION:  Scalar
!
!:sdoc-:
!--------------------------------------------------------------------------------

  PURE FUNCTION mwwatercoeff_associated( self ) RESULT( Status )
    TYPE(mwwatercoeff_type), INTENT(IN) :: self
    LOGICAL :: status
    status = self%Is_Allocated
  END FUNCTION mwwatercoeff_associated


!--------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       MWwaterCoeff_Destroy
!
! PURPOSE:
!       Pure subroutine to re-initialize MWwaterCoeff objects.
!
! CALLING SEQUENCE:
!       CALL MWwaterCoeff_Destroy( MWwaterCoeff )
!
! OBJECTS:
!       MWwaterCoeff: Re-initialized MWwaterCoeff structure.
!                     UNITS:      N/A
!                     TYPE:       MWwaterCoeff_type
!                     DIMENSION:  Scalar
!                     ATTRIBUTES: INTENT(OUT)
!
!:sdoc-:
!--------------------------------------------------------------------------------

  PURE SUBROUTINE mwwatercoeff_destroy( self )
    TYPE(mwwatercoeff_type), INTENT(OUT) :: self
    self%Is_Allocated = .false.
  END SUBROUTINE mwwatercoeff_destroy


!--------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       MWwaterCoeff_Create
!
! PURPOSE:
!       Pure subroutine to create a valid instance of an MWwaterCoeff object.
!
! CALLING SEQUENCE:
!       CALL MWwaterCoeff_Create( MWwaterCoeff )
!
! OBJECTS:
!       MWwaterCoeff:       MWwaterCoeff object structure.
!                           UNITS:      N/A
!                           TYPE:       MWwaterCoeff_type
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN OUT)
!
!:sdoc-:
!--------------------------------------------------------------------------------

  PURE SUBROUTINE mwwatercoeff_create( &
    self )  ! Output
    ! Arguments
    TYPE(mwwatercoeff_type), INTENT(IN OUT) :: self
    ! Set allocation indicator
    self%Is_Allocated = .true.
  END SUBROUTINE mwwatercoeff_create


!--------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       MWwaterCoeff_Inspect
!
! PURPOSE:
!       Subroutine to print the contents of a MWwaterCoeff object to stdout.
!
! CALLING SEQUENCE:
!       CALL MWwaterCoeff_Inspect( MWwaterCoeff )
!
! OBJECTS:
!       MWwaterCoeff:  MWwaterCoeff object to display.
!                      UNITS:      N/A
!                      TYPE:       MWwaterCoeff_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN)
!
!:sdoc-:
!--------------------------------------------------------------------------------

  SUBROUTINE mwwatercoeff_inspect( self, pause )
    TYPE(mwwatercoeff_type), INTENT(IN) :: self
    LOGICAL,       OPTIONAL, INTENT(IN) :: pause
    LOGICAL :: wait

    wait = .false.
    IF ( PRESENT(pause) ) wait = pause

    WRITE(*,'(1x,"MWwaterCoeff OBJECT")')
    ! Release/version info
    WRITE(*,'(3x,"Release.Version : ",i0,".",i0)') self%Release, self%Version
    IF ( .NOT. mwwatercoeff_associated(self) ) RETURN
    ! Derived types
    IF ( fitcoeff_associated(self%FCCoeff ) ) THEN
      WRITE(*,'(1x,"Foam coverage cofficients, ")',advance='NO')
      CALL fitcoeff_inspect(self%FCCoeff )
    END IF
    IF ( fitcoeff_associated(self%FRCoeff ) ) THEN
      WRITE(*,'(1x,"Foam reflectivity cofficients, ")',advance='NO')
      CALL fitcoeff_inspect(self%FRCoeff )
    END IF
    IF ( fitcoeff_associated(self%RCCoeff) ) THEN
      WRITE(*,'(1x,"Reflection correction cofficients, ")',advance='NO')
      CALL fitcoeff_inspect(self%RCCoeff)
    END IF
    IF ( fitcoeff_associated(self%AZCoeff) ) THEN
      WRITE(*,'(1x,"Azimuth emissivity coefficients, ")',advance='NO')
      CALL fitcoeff_inspect(self%AZCoeff)
    END IF
    IF ( fitcoeff_associated(self%SSCCoeff) ) THEN
      WRITE(*,'(1x,"Small-scale correction coefficients, ")',advance='NO')
      CALL fitcoeff_inspect(self%SSCCoeff)
    END IF
    IF ( fitcoeff_associated(self%LSCCoeff) ) THEN
      WRITE(*,'(1x,"Large-scale correction coefficients, ")',advance='NO')
      CALL fitcoeff_inspect(self%LSCCoeff)
    END IF
    IF ( mwwaterlut_associated(self%LUT) ) THEN
      WRITE(*,'(1x,"Emissivity look-up table, ")',advance='NO')
      CALL mwwaterlut_inspect(self%LUT,pause=pause)
    END IF
  END SUBROUTINE mwwatercoeff_inspect



!----------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       MWwaterCoeff_ValidRelease
!
! PURPOSE:
!       Function to check the MWwaterCoeff Release value.
!
! CALLING SEQUENCE:
!       IsValid = MWwaterCoeff_ValidRelease( MWwaterCoeff )
!
! INPUTS:
!       MWwaterCoeff:  MWwaterCoeff object for which the Release component
!                      is to be checked.
!                      UNITS:      N/A
!                      TYPE:       MWwaterCoeff_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN)
!
! FUNCTION RESULT:
!       IsValid:       Logical value defining the release validity.
!                      UNITS:      N/A
!                      TYPE:       LOGICAL
!                      DIMENSION:  Scalar
!
!:sdoc-:
!----------------------------------------------------------------------------------

  FUNCTION mwwatercoeff_validrelease( self ) RESULT( IsValid )
    ! Arguments
    TYPE(mwwatercoeff_type), INTENT(IN) :: self
    ! Function result
    LOGICAL :: isvalid
    ! Local parameters
    CHARACTER(*), PARAMETER :: routine_name = 'MWwaterCoeff_ValidRelease'
    ! Local variables
    CHARACTER(ML) :: msg

    ! Set up
    isvalid = .true.


    ! Check release is not too old
    IF ( self%Release < mwwatercoeff_release ) THEN
      isvalid = .false.
      WRITE( msg,'("An MWwaterCoeff data update is needed. ", &
                  &"MWwaterCoeff release is ",i0,". Valid release is ",i0,"." )' ) &
                  self%Release, mwwatercoeff_release
      CALL display_message( routine_name, msg, information ); RETURN
    END IF


    ! Check release is not too new
    IF ( self%Release > mwwatercoeff_release ) THEN
      isvalid = .false.
      WRITE( msg,'("An MWwaterCoeff software update is needed. ", &
                  &"MWwaterCoeff release is ",i0,". Valid release is ",i0,"." )' ) &
                  self%Release, mwwatercoeff_release
      CALL display_message( routine_name, msg, information ); RETURN
    END IF

  END FUNCTION mwwatercoeff_validrelease


!--------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       MWwaterCoeff_Info
!
! PURPOSE:
!       Subroutine to return a string containing version and dimension
!       information about a MWwaterCoeff object.
!
! CALLING SEQUENCE:
!       CALL MWwaterCoeff_Info( MWwaterCoeff, Info )
!
! OBJECTS:
!       MWwaterCoeff:  MWwaterCoeff object about which info is required.
!                      UNITS:      N/A
!                      TYPE:       MWwaterCoeff_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN)
!
! OUTPUTS:
!       Info:          String containing version and dimension information
!                      about the MWwaterCoeff object.
!                      UNITS:      N/A
!                      TYPE:       CHARACTER(*)
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(OUT)
!
!:sdoc-:
!--------------------------------------------------------------------------------

  SUBROUTINE mwwatercoeff_info( self, Info )
    ! Arguments
    TYPE(mwwatercoeff_type), INTENT(IN)  :: self
    CHARACTER(*),            INTENT(OUT) :: info
    ! Parameters
    INTEGER, PARAMETER :: carriage_return = 13
    INTEGER, PARAMETER :: linefeed = 10
    ! Local variables
    CHARACTER(2000) :: long_string

    ! Write the required data to the local string
    WRITE( long_string, &
           '(a,1x,"MWwaterCoeff RELEASE.VERSION: ",i0,".",i0 )' ) &
           achar(carriage_return)//achar(linefeed), &
           self%Release, self%Version

    ! Trim the output based on the
    ! dummy argument string length
    info = long_string(1:min(len(info), len_trim(long_string)))

  END SUBROUTINE mwwatercoeff_info


!--------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       MWwaterCoeff_DefineVersion
!
! PURPOSE:
!       Subroutine to return the module version information.
!
! CALLING SEQUENCE:
!       CALL MWwaterCoeff_DefineVersion( Id )
!
! OUTPUTS:
!       Id:    Character string containing the version Id information
!              for the module.
!              UNITS:      N/A
!              TYPE:       CHARACTER(*)
!              DIMENSION:  Scalar
!              ATTRIBUTES: INTENT(OUT)
!
!:sdoc-:
!--------------------------------------------------------------------------------

  SUBROUTINE mwwatercoeff_defineversion( Id )
    CHARACTER(*), INTENT(OUT) :: id
    id = module_version_id
  END SUBROUTINE mwwatercoeff_defineversion


!--------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       MWwaterCoeff_SetValue
!
! PURPOSE:
!       Subroutine to set the contents of a valid MWwaterCoeff object.
!
! CALLING SEQUENCE:
!       CALL MWwaterCoeff_SetValue( MWwaterCoeff, &
!                                   FCCoeff    = FCCoeff   , &
!                                   FRCoeff    = FRCoeff   , &
!                                   RCCoeff    = RCCoeff   , &
!                                   AZCoeff    = AZCoeff   , &
!                                   SSCCoeff   = SSCCoeff  , &
!                                   LSCCoeff   = LSCCoeff  , &
!                                   MWwaterLUT = MWwaterLUT  )
! OBJECTS:
!       MWwaterCoeff:  Valid, allocated MWwaterCoeff object for which
!                      values are to be set.
!                      UNITS:      N/A
!                      TYPE:       MWwaterCoeff_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN OUT)
!
! OPTIONAL INPUTS:
!       FCCoeff:       Object containing the foam coverage fitting coefficients.
!                      UNITS:      N/A
!                      TYPE:       FitCoeff_1D_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       FRCoeff:       Object containing the foam reflectivity fitting
!                      coefficients.
!                      UNITS:      N/A
!                      TYPE:       FitCoeff_1D_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       RCCoeff:       Object containing the reflection correction fitting
!                      coefficients.
!                      UNITS:      N/A
!                      TYPE:       FitCoeff_3D_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       AZCoeff:       Object containing the azimuth emissivity fitting
!                      coefficients.
!                      UNITS:      N/A
!                      TYPE:       FitCoeff_3D_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       SSCCoeff:      Object containing the small-scale correction fitting
!                      coefficients.
!                      UNITS:      N/A
!                      TYPE:       FitCoeff_1D_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       LSCCoeff:      Object containing the large-scale correction fitting
!                      coefficients.
!                      UNITS:      N/A
!                      TYPE:       FitCoeff_3D_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       MWwaterLUT:    Object containing the emissivity look-up table.
!                      UNITS:      N/A
!                      TYPE:       MWwaterLUT_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN), OPTIONAL
!
!:sdoc-:
!--------------------------------------------------------------------------------

  SUBROUTINE mwwatercoeff_setvalue( &
    self      , &  ! In/Output
    FCCoeff   , &  ! Optional input
    FRCoeff   , &  ! Optional input
    RCCoeff   , &  ! Optional input
    AZCoeff   , &  ! Optional input
    SSCCoeff  , &  ! Optional input
    LSCCoeff  , &  ! Optional input
    MWwaterLUT  )  ! Optional input
    ! Arguments
    TYPE(MWwaterCoeff_type)         , INTENT(IN OUT) :: self
    TYPE(FitCoeff_1D_type), OPTIONAL, INTENT(IN)     :: FCCoeff
    TYPE(FitCoeff_1D_type), OPTIONAL, INTENT(IN)     :: FRCoeff
    TYPE(FitCoeff_3D_type), OPTIONAL, INTENT(IN)     :: RCCoeff
    TYPE(FitCoeff_3D_type), OPTIONAL, INTENT(IN)     :: AZCoeff
    TYPE(FitCoeff_1D_type), OPTIONAL, INTENT(IN)     :: SSCCoeff
    TYPE(FitCoeff_3D_type), OPTIONAL, INTENT(IN)     :: LSCCoeff
    TYPE(MWwaterLUT_type) , OPTIONAL, INTENT(IN)     :: MWwaterLUT

    IF ( PRESENT(fccoeff   ) ) self%FCCoeff  = fccoeff
    IF ( PRESENT(frcoeff   ) ) self%FRCoeff  = frcoeff
    IF ( PRESENT(rccoeff   ) ) self%RCCoeff  = rccoeff
    IF ( PRESENT(azcoeff   ) ) self%AZCoeff  = azcoeff
    IF ( PRESENT(ssccoeff  ) ) self%SSCCoeff = ssccoeff
    IF ( PRESENT(lsccoeff  ) ) self%LSCCoeff = lsccoeff
    IF ( PRESENT(mwwaterlut) ) self%LUT      = mwwaterlut

  END SUBROUTINE mwwatercoeff_setvalue


!--------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       MWwaterCoeff_GetValue
!
! PURPOSE:
!       Subroutine to get the contents of a valid MWwaterCoeff object.
!
! CALLING SEQUENCE:
!       CALL MWwaterCoeff_GetValue( MWwaterCoeff, &
!                                   FCCoeff    = FCCoeff   , &
!                                   FRCoeff    = FRCoeff   , &
!                                   RCCoeff    = RCCoeff   , &
!                                   AZCoeff    = AZCoeff   , &
!                                   SSCCoeff   = SSCCoeff  , &
!                                   LSCCoeff   = LSCCoeff  , &
!                                   MWwaterLUT = MWwaterLUT  )
!
! OBJECTS:
!       MWwaterCoeff:  Valid, allocated MWwaterCoeff object from which
!                      values are to be retrieved.
!                      UNITS:      N/A
!                      TYPE:       MWwaterCoeff_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN OUT)
!
! OPTIONAL OUTPUTS:
!       FCCoeff:       Object containing the foam coverage fitting coefficients.
!                      UNITS:      N/A
!                      TYPE:       FitCoeff_1D_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(OUT), OPTIONAL
!
!       FRCoeff:       Object containing the foam reflectivity fitting
!                      coefficients.
!                      UNITS:      N/A
!                      TYPE:       FitCoeff_1D_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(OUT), OPTIONAL
!
!       RCCoeff:       Object containing the reflection correction fitting
!                      coefficients.
!                      UNITS:      N/A
!                      TYPE:       FitCoeff_3D_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(OUT), OPTIONAL
!
!       AZCoeff:       Object containing the azimuth emissivity fitting
!                      coefficients.
!                      UNITS:      N/A
!                      TYPE:       FitCoeff_3D_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(OUT), OPTIONAL
!
!       SSCCoeff:      Object containing the small-scale correction fitting
!                      coefficients.
!                      UNITS:      N/A
!                      TYPE:       FitCoeff_1D_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(OUT), OPTIONAL
!
!       LSCCoeff:      Object containing the large-scale correction fitting
!                      coefficients.
!                      UNITS:      N/A
!                      TYPE:       FitCoeff_3D_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(OUT), OPTIONAL
!
!       MWwaterLUT:    Object containing the emissivity look-up table.
!                      UNITS:      N/A
!                      TYPE:       MWwaterLUT_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(OUT), OPTIONAL
!
!:sdoc-:
!--------------------------------------------------------------------------------

  SUBROUTINE mwwatercoeff_getvalue( &
    self      , &  ! Input
    FCCoeff   , &  ! Optional output
    FRCoeff   , &  ! Optional output
    RCCoeff   , &  ! Optional output
    AZCoeff   , &  ! Optional output
    SSCCoeff  , &  ! Optional output
    LSCCoeff  , &  ! Optional output
    MWwaterLUT  )  ! Optional output
    ! Arguments
    TYPE(MWwaterCoeff_type)         , INTENT(IN)  :: self
    TYPE(FitCoeff_1D_type), OPTIONAL, INTENT(OUT) :: FCCoeff
    TYPE(FitCoeff_1D_type), OPTIONAL, INTENT(OUT) :: FRCoeff
    TYPE(FitCoeff_3D_type), OPTIONAL, INTENT(OUT) :: RCCoeff
    TYPE(FitCoeff_3D_type), OPTIONAL, INTENT(OUT) :: AZCoeff
    TYPE(FitCoeff_1D_type), OPTIONAL, INTENT(OUT) :: SSCCoeff
    TYPE(FitCoeff_3D_type), OPTIONAL, INTENT(OUT) :: LSCCoeff
    TYPE(MWwaterLUT_type) , OPTIONAL, INTENT(OUT) :: MWwaterLUT

    IF ( PRESENT(fccoeff   ) ) fccoeff    = self%FCCoeff
    IF ( PRESENT(frcoeff   ) ) frcoeff    = self%FRCoeff
    IF ( PRESENT(rccoeff   ) ) rccoeff    = self%RCCoeff
    IF ( PRESENT(azcoeff   ) ) azcoeff    = self%AZCoeff
    IF ( PRESENT(ssccoeff  ) ) ssccoeff   = self%SSCCoeff
    IF ( PRESENT(lsccoeff  ) ) lsccoeff   = self%LSCCoeff
    IF ( PRESENT(mwwaterlut) ) mwwaterlut = self%LUT

  END SUBROUTINE mwwatercoeff_getvalue


!--------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       MWwaterCoeff_ReadFile
!
! PURPOSE:
!       Function to read MWwaterCoeff object files.
!
! CALLING SEQUENCE:
!       Error_Status = MWwaterCoeff_ReadFile( &
!                        MWwaterCoeff       , &
!                        Filename           , &
!                        No_Close = No_Close, &
!                        Quiet    = Quiet   , &
!                        Title    = Title   , &
!                        History  = History , &
!                        Comment  = Comment   )
!
! OBJECTS:
!       MWwaterCoeff:   MWwaterCoeff object containing the data read from file.
!                       UNITS:      N/A
!                       TYPE:       MWwaterCoeff_type
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(OUT)
!
! INPUTS:
!       Filename:       Character string specifying the name of a
!                       MWwaterCoeff data file to read.
!                       UNITS:      N/A
!                       TYPE:       CHARACTER(*)
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(IN)
!
! OPTIONAL INPUTS:
!       No_Close:       Set this logical argument to *NOT* close the datafile
!                       upon exiting this routine. This option is required if
!                       the MWwaterCoeff data is embedded within another file.
!                       If == .FALSE., File is closed upon function exit [DEFAULT].
!                          == .TRUE.,  File is NOT closed upon function exit
!                       If not specified, default is .FALSE.
!                       UNITS:      N/A
!                       TYPE:       LOGICAL
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       Quiet:          Set this logical argument to suppress INFORMATION
!                       messages being printed to stdout
!                       If == .FALSE., INFORMATION messages are OUTPUT [DEFAULT].
!                          == .TRUE.,  INFORMATION messages are SUPPRESSED.
!                       If not specified, default is .FALSE.
!                       UNITS:      N/A
!                       TYPE:       LOGICAL
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(IN), OPTIONAL
!
! OPTIONAL OUTPUTS:
!       Title:          Character string containing a succinct description
!                       of what is in the dataset.
!                       UNITS:      N/A
!                       TYPE:       CHARACTER(*)
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(OUT), OPTIONAL
!
!       History:        Character string containing dataset creation
!                       history.
!                       UNITS:      N/A
!                       TYPE:       CHARACTER(*)
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(OUT), OPTIONAL
!
!       Comment:        Character string containing any comments about
!                       the dataset.
!                       UNITS:      N/A
!                       TYPE:       CHARACTER(*)
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(OUT), OPTIONAL
!
! FUNCTION RESULT:
!       Error_Status:   The return value is an integer defining the error status.
!                       The error codes are defined in the Message_Handler module.
!                       If == SUCCESS, the file read was successful
!                          == FAILURE, an unrecoverable error occurred.
!                       UNITS:      N/A
!                       TYPE:       INTEGER
!                       DIMENSION:  Scalar
!
!:sdoc-:
!------------------------------------------------------------------------------

  FUNCTION mwwatercoeff_readfile( &
    MWwaterCoeff, &  ! Output
    Filename    , &  ! Input
    No_Close    , &  ! Optional input
    Quiet       , &  ! Optional input
    Title       , &  ! Optional output
    History     , &  ! Optional output
    Comment     , &  ! Optional output
    Debug       ) &  ! Optional input (Debug output control)
  result( err_stat )
    ! Arguments
    TYPE(mwwatercoeff_type), INTENT(OUT) :: mwwatercoeff
    CHARACTER(*),            INTENT(IN)  :: filename
    LOGICAL     ,  OPTIONAL, INTENT(IN)  :: no_close
    LOGICAL     ,  OPTIONAL, INTENT(IN)  :: quiet
    CHARACTER(*),  OPTIONAL, INTENT(OUT) :: title
    CHARACTER(*),  OPTIONAL, INTENT(OUT) :: history
    CHARACTER(*),  OPTIONAL, INTENT(OUT) :: comment
    LOGICAL     ,  OPTIONAL, INTENT(IN)  :: debug
    ! Function result
    INTEGER :: err_stat
    ! Function parameters
    CHARACTER(*), PARAMETER :: routine_name = 'MWwaterCoeff_ReadFile'
    ! Function variables
    CHARACTER(ML) :: msg
    CHARACTER(ML) :: io_msg
    LOGICAL :: close_file
    LOGICAL :: noisy
    INTEGER :: io_stat
    INTEGER :: fid
    INTEGER(Long) :: lut_present
    TYPE(mwwatercoeff_type) :: dummy

    ! Setup
    err_stat = success
    ! ...Check No_Close argument
    close_file = .true.
    IF ( PRESENT(no_close) ) close_file = .NOT. no_close
    ! ...Check Quiet argument
    noisy = .true.
    IF ( PRESENT(quiet) ) noisy = .NOT. quiet
    ! ...Override Quiet settings if debug set.
    IF ( PRESENT(debug) ) THEN
      IF ( debug ) noisy = .true.
    END IF


    ! Check if the file is open.
    IF ( file_open( filename ) ) THEN
      ! ...Inquire for the logical unit number
      INQUIRE( file=filename, number=fid )
      ! ...Ensure it's valid
      IF ( fid < 0 ) THEN
        msg = 'Error inquiring '//trim(filename)//' for its FileID'
        CALL read_cleanup(); RETURN
      END IF
    ELSE
      ! ...Open the file if it exists
      IF ( file_exists( filename ) ) THEN
        err_stat = open_binary_file( filename, fid )
        IF ( err_stat /= success ) THEN
          msg = 'Error opening '//trim(filename)
          CALL read_cleanup(); RETURN
        END IF
      ELSE
        msg = 'File '//trim(filename)//' not found.'
        CALL read_cleanup(); RETURN
      END IF
    END IF


    ! Read and check the release and version
    READ( fid, iostat=io_stat, iomsg=io_msg ) &
      dummy%Release, &
      dummy%Version
    IF ( io_stat /= 0 ) THEN
      msg = 'Error reading Release/Version - '//trim(io_msg)
      CALL read_cleanup(); RETURN
    END IF
    IF ( .NOT. mwwatercoeff_validrelease( dummy ) ) THEN
      msg = 'MWwaterCoeff Release check failed.'
      CALL read_cleanup(); RETURN
    END IF


    ! Create instance of the structure
    CALL mwwatercoeff_create(mwwatercoeff)
    ! ...Explicitly set the version
    mwwatercoeff%Version = dummy%Version


    ! Read the global attributes
    err_stat = readgatts_binary_file( &
                 fid, &
                 title   = title  , &
                 history = history, &
                 comment = comment  )
    IF ( err_stat /= success ) THEN
      msg = 'Error reading global attributes'
      CALL read_cleanup(); RETURN
    END IF


    ! Read the coefficient data
    ! ...Read the foam coverage coefficients
    err_stat = fitcoeff_readfile( &
                 mwwatercoeff%FCCoeff, &
                 filename, &
                 no_close = .true., &
                 quiet    = quiet , &
                 debug    = debug   )
    IF ( err_stat /= success ) THEN
      msg = 'Error reading foam coverage coefficient data.'
      CALL read_cleanup(); RETURN
    END IF
    ! ...Read the foam reflectivity coefficients
    err_stat = fitcoeff_readfile( &
                 mwwatercoeff%FRCoeff, &
                 filename, &
                 no_close = .true., &
                 quiet    = quiet , &
                 debug    = debug   )
    IF ( err_stat /= success ) THEN
      msg = 'Error reading foam reflectivity coefficient data.'
      CALL read_cleanup(); RETURN
    END IF
    ! ...Read the reflection correction coefficients
    err_stat = fitcoeff_readfile( &
                 mwwatercoeff%RCCoeff, &
                 filename, &
                 no_close = .true., &
                 quiet    = quiet , &
                 debug    = debug   )
    IF ( err_stat /= success ) THEN
      msg = 'Error reading reflection correction coefficient data'
      CALL read_cleanup(); RETURN
    END IF
    ! ...Read the azimuth emissivity coefficients
    err_stat = fitcoeff_readfile( &
                 mwwatercoeff%AZCoeff, &
                 filename, &
                 no_close = .true., &
                 quiet    = quiet , &
                 debug    = debug   )
    IF ( err_stat /= success ) THEN
      msg = 'Error reading azimuth emissivity coefficient data.'
      CALL read_cleanup(); RETURN
    END IF
    ! ...Read the small-scale correction coefficients
    err_stat = fitcoeff_readfile( &
                 mwwatercoeff%SSCCoeff, &
                 filename, &
                 no_close = .true., &
                 quiet    = quiet , &
                 debug    = debug   )
    IF ( err_stat /= success ) THEN
      msg = 'Error reading small-scale correction coefficient data.'
      CALL read_cleanup(); RETURN
    END IF
    ! ...Read the large-scale correction coefficients
    err_stat = fitcoeff_readfile( &
                 mwwatercoeff%LSCCoeff, &
                 filename, &
                 no_close = .true., &
                 quiet    = quiet , &
                 debug    = debug   )
    IF ( err_stat /= success ) THEN
      msg = 'Error reading large-scale correction coefficient data.'
      CALL read_cleanup(); RETURN
    END IF


    ! Read the emissivity look-up table if it's present
    ! ...Read the data indicator
    READ( fid, iostat=io_stat, iomsg=io_msg ) lut_present
    IF ( io_stat /= 0 ) THEN
      msg = 'Error reading emissivity look-up table data indicator - '//trim(io_msg)
      CALL read_cleanup(); RETURN
    END IF
    ! ...Read the LUT data
    IF ( lut_present == data_present ) THEN
      err_stat = mwwaterlut_readfile( &
                   mwwatercoeff%LUT, &
                   filename, &
                   no_close = .true., &
                   quiet    = quiet , &
                   debug    = debug   )
      IF ( err_stat /= success ) THEN
        msg = 'Error reading emissivity look-up table data.'
        CALL read_cleanup(); RETURN
      END IF
    END IF


    ! Close the file
    IF ( close_file ) THEN
      CLOSE( fid, iostat=io_stat, iomsg=io_msg )
      IF ( io_stat /= 0 ) THEN
        msg = 'Error closing '//trim(filename)//' - '//trim(io_msg)
        CALL read_cleanup(); RETURN
      END IF
    END IF


    ! Output an info message
    IF ( noisy ) THEN
      CALL mwwatercoeff_info( mwwatercoeff, msg )
      CALL display_message( routine_name, 'FILE: '//trim(filename)//'; '//trim(msg), information )
    END IF

  CONTAINS

    SUBROUTINE read_cleanup()
      IF ( file_open(filename) ) THEN
        CLOSE( fid, iostat=io_stat, iomsg=io_msg )
        IF ( io_stat /= 0 ) &
          msg = trim(msg)//'; Error closing input file during error cleanup - '//trim(io_msg)
      END IF
      CALL mwwatercoeff_destroy( mwwatercoeff )
      err_stat = failure
      CALL display_message( routine_name, msg, err_stat )
    END SUBROUTINE read_cleanup

  END FUNCTION mwwatercoeff_readfile


!--------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       MWwaterCoeff_WriteFile
!
! PURPOSE:
!       Function to write MWwaterCoeff object files.
!
! CALLING SEQUENCE:
!       Error_Status = MWwaterCoeff_WriteFile( &
!                        MWwaterCoeff       , &
!                        Filename           , &
!                        No_Close = No_Close, &
!                        Quiet    = Quiet   , &
!                        Title    = Title   , &
!                        History  = History , &
!                        Comment  = Comment   )
!
! OBJECTS:
!       MWwaterCoeff:   MWwaterCoeff object containing the data to write to file.
!                       UNITS:      N/A
!                       TYPE:       MWwaterCoeff_type
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(IN)
!
! INPUTS:
!       Filename:       Character string specifying the name of a
!                       MWwaterCoeff format data file to write.
!                       UNITS:      N/A
!                       TYPE:       CHARACTER(*)
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(IN)
!
! OPTIONAL INPUTS:
!       No_Close:       Set this logical argument to *NOT* close the datafile
!                       upon exiting this routine. This option is required if
!                       the MWwaterCoeff data is to be embedded within another file.
!                       If == .FALSE., File is closed upon function exit [DEFAULT].
!                          == .TRUE.,  File is NOT closed upon function exit
!                       If not specified, default is .FALSE.
!                       UNITS:      N/A
!                       TYPE:       LOGICAL
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       Quiet:          Set this logical argument to suppress INFORMATION
!                       messages being printed to stdout
!                       If == .FALSE., INFORMATION messages are OUTPUT [DEFAULT].
!                          == .TRUE.,  INFORMATION messages are SUPPRESSED.
!                       If not specified, default is .FALSE.
!                       UNITS:      N/A
!                       TYPE:       LOGICAL
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       Title:          Character string containing a succinct description
!                       of what is in the dataset.
!                       UNITS:      N/A
!                       TYPE:       CHARACTER(*)
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       History:        Character string containing dataset creation
!                       history.
!                       UNITS:      N/A
!                       TYPE:       CHARACTER(*)
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       Comment:        Character string containing any comments about
!                       the dataset.
!                       UNITS:      N/A
!                       TYPE:       CHARACTER(*)
!                       DIMENSION:  Scalar
!                       ATTRIBUTES: INTENT(IN), OPTIONAL
!
! FUNCTION RESULT:
!       Error_Status:   The return value is an integer defining the error status.
!                       The error codes are defined in the Message_Handler module.
!                       If == SUCCESS, the file write was successful
!                          == FAILURE, an unrecoverable error occurred.
!                       UNITS:      N/A
!                       TYPE:       INTEGER
!                       DIMENSION:  Scalar
!
!:sdoc-:
!------------------------------------------------------------------------------

  FUNCTION mwwatercoeff_writefile( &
    MWwaterCoeff, &  ! Input
    Filename    , &  ! Input
    No_Close    , &  ! Optional input
    Quiet       , &  ! Optional input
    Title       , &  ! Optional input
    History     , &  ! Optional input
    Comment     , &  ! Optional input
    Debug       ) &  ! Optional input (Debug output control)
  result( err_stat )
    ! Arguments
    TYPE(mwwatercoeff_type), INTENT(IN) :: mwwatercoeff
    CHARACTER(*),            INTENT(IN) :: filename
    LOGICAL     ,  OPTIONAL, INTENT(IN) :: no_close
    LOGICAL     ,  OPTIONAL, INTENT(IN) :: quiet
    CHARACTER(*),  OPTIONAL, INTENT(IN) :: title
    CHARACTER(*),  OPTIONAL, INTENT(IN) :: history
    CHARACTER(*),  OPTIONAL, INTENT(IN) :: comment
    LOGICAL     ,  OPTIONAL, INTENT(IN) :: debug
    ! Function result
    INTEGER :: err_stat
    ! Function parameters
    CHARACTER(*), PARAMETER :: routine_name = 'MWwaterCoeff_WriteFile'
    ! Function variables
    CHARACTER(ML) :: msg
    CHARACTER(ML) :: io_msg
    LOGICAL :: close_file
    LOGICAL :: noisy
    INTEGER :: io_stat
    INTEGER :: fid
    INTEGER(Long) :: lut_present


    ! Setup
    err_stat = success
    ! ...Check No_Close argument
    close_file = .true.
    IF ( PRESENT(no_close) ) close_file = .NOT. no_close
    ! ...Check Quiet argument
    noisy = .true.
    IF ( PRESENT(quiet) ) noisy = .NOT. quiet
    ! ...Override Quiet settings if debug set.
    IF ( PRESENT(debug) ) THEN
      IF ( debug ) noisy = .true.
    END IF
    ! ...Check there is data to write
    IF ( .NOT. mwwatercoeff_associated( mwwatercoeff ) ) THEN
      msg = 'MWwaterCoeff object is empty.'
      CALL write_cleanup(); RETURN
    END IF


    ! Check if the file is open.
    IF ( file_open( filename ) ) THEN
      ! ...Inquire for the logical unit number
      INQUIRE( file=filename, number=fid )
      ! ...Ensure it's valid
      IF ( fid < 0 ) THEN
        msg = 'Error inquiring '//trim(filename)//' for its FileID'
        CALL write_cleanup(); RETURN
      END IF
    ELSE
      ! ...Open the file for output
      err_stat = open_binary_file( filename, fid, for_output=.true. )
      IF ( err_stat /= success ) THEN
        msg = 'Error opening '//trim(filename)
        CALL write_cleanup(); RETURN
      END IF
    END IF


    ! Write the release and version
    WRITE( fid, iostat=io_stat, iomsg=io_msg ) &
      mwwatercoeff%Release, &
      mwwatercoeff%Version
    IF ( io_stat /= 0 ) THEN
      msg = 'Error writing Release/Version - '//trim(io_msg)
      CALL write_cleanup(); RETURN
    END IF


    ! Write the global attributes
    err_stat = writegatts_binary_file( &
                 fid, &
                 write_module = module_version_id, &
                 title        = title  , &
                 history      = history, &
                 comment      = comment  )
    IF ( err_stat /= success ) THEN
      msg = 'Error writing global attributes'
      CALL write_cleanup(); RETURN
    END IF


    ! Write the coefficient data
    ! ...Write the foam coverage coefficients
    err_stat = fitcoeff_writefile( &
                 mwwatercoeff%FCCoeff, &
                 filename, &
                 no_close = .true., &
                 quiet    = quiet , &
                 debug    = debug   )
    IF ( err_stat /= success ) THEN
      msg = 'Error writing foam coverage coefficient data.'
      CALL write_cleanup(); RETURN
    END IF
    ! ...Write the foam reflectivity coefficients
    err_stat = fitcoeff_writefile( &
                 mwwatercoeff%FRCoeff, &
                 filename, &
                 no_close = .true., &
                 quiet    = quiet , &
                 debug    = debug   )
    IF ( err_stat /= success ) THEN
      msg = 'Error writing foam reflectivity coefficient data.'
      CALL write_cleanup(); RETURN
    END IF
    ! ...Write the reflection correction coefficients
    err_stat = fitcoeff_writefile( &
                 mwwatercoeff%RCCoeff, &
                 filename, &
                 no_close = .true., &
                 quiet    = quiet , &
                 debug    = debug   )
    IF ( err_stat /= success ) THEN
      msg = 'Error writing reflection correction coefficient data'
      CALL write_cleanup(); RETURN
    END IF
    ! ...Write the azimuth emissivity coefficients
    err_stat = fitcoeff_writefile( &
                 mwwatercoeff%AZCoeff, &
                 filename, &
                 no_close = .true., &
                 quiet    = quiet , &
                 debug    = debug   )
    IF ( err_stat /= success ) THEN
      msg = 'Error writing azimuth emissivity coefficient data.'
      CALL write_cleanup(); RETURN
    END IF
    ! ...Write the small-scale correction coefficients
    err_stat = fitcoeff_writefile( &
                 mwwatercoeff%SSCCoeff, &
                 filename, &
                 no_close = .true., &
                 quiet    = quiet , &
                 debug    = debug   )
    IF ( err_stat /= success ) THEN
      msg = 'Error writing small-scale correction coefficient data.'
      CALL write_cleanup(); RETURN
    END IF
    ! ...Write the large-scale correction coefficients
    err_stat = fitcoeff_writefile( &
                 mwwatercoeff%LSCCoeff, &
                 filename, &
                 no_close = .true., &
                 quiet    = quiet , &
                 debug    = debug   )
    IF ( err_stat /= success ) THEN
      msg = 'Error writing large-scale correction coefficient data.'
      CALL write_cleanup(); RETURN
    END IF


    ! Write the emissivity look-up table if it's present
    IF ( mwwaterlut_associated( mwwatercoeff%LUT ) ) THEN
      lut_present = data_present
    ELSE
      lut_present = data_missing
    END IF
    ! ...Write the data indicator
    WRITE( fid, iostat=io_stat, iomsg=io_msg ) lut_present
    IF ( io_stat /= 0 ) THEN
      msg = 'Error writing emissivity look-up table data indicator - '//trim(io_msg)
      CALL write_cleanup(); RETURN
    END IF
    ! ...Write the LUT data
    IF ( lut_present == data_present ) THEN
      err_stat = mwwaterlut_writefile( &
                   mwwatercoeff%LUT, &
                   filename, &
                   no_close = .true., &
                   quiet    = quiet , &
                   debug    = debug   )
      IF ( err_stat /= success ) THEN
        msg = 'Error writing emissivity look-up table data.'
        CALL write_cleanup(); RETURN
      END IF
    END IF


    ! Close the file
    IF ( close_file ) THEN
      CLOSE( fid, iostat=io_stat, iomsg=io_msg )
      IF ( io_stat /= 0 ) THEN
        msg = 'Error closing '//trim(filename)//' - '//trim(io_msg)
        CALL write_cleanup(); RETURN
      END IF
    END IF


    ! Output an info message
    IF ( noisy ) THEN
      CALL mwwatercoeff_info( mwwatercoeff, msg )
      CALL display_message( routine_name, 'FILE: '//trim(filename)//'; '//trim(msg), information )
    END IF

  CONTAINS

    SUBROUTINE write_cleanup()
      IF ( file_open(filename) ) THEN
        CLOSE( fid, status=write_error_status, iostat=io_stat, iomsg=io_msg )
        IF ( io_stat /= 0 ) &
          msg = trim(msg)//'; Error closing output file during error cleanup - '//trim(io_msg)
      END IF
      err_stat = failure
      CALL display_message( routine_name, msg, err_stat )
    END SUBROUTINE write_cleanup

  END FUNCTION mwwatercoeff_writefile


!################################################################################
!################################################################################
!##                                                                            ##
!##                          ## PRIVATE PROCEDURES ##                          ##
!##                                                                            ##
!################################################################################
!################################################################################

!--------------------------------------------------------------------------------
!
! NAME:
!       MWwaterCoeff_Equal
!
! PURPOSE:
!       Pure function to test the equality of two MWwaterCoeff objects.
!       Used in OPERATOR(==) interface block.
!
! CALLING SEQUENCE:
!       is_equal = MWwaterCoeff_Equal( x, y )
!
!         or
!
!       IF ( x == y ) THEN
!         ...
!       END IF
!
! OBJECTS:
!       x, y:          Two MWwaterCoeff objects to be compared.
!                      UNITS:      N/A
!                      TYPE:       MWwaterCoeff_type
!                      DIMENSION:  Scalar
!                      ATTRIBUTES: INTENT(IN)
!
! FUNCTION RESULT:
!       is_equal:      Logical value indicating whether the inputs are equal.
!                      UNITS:      N/A
!                      TYPE:       LOGICAL
!                      DIMENSION:  Scalar
!
!--------------------------------------------------------------------------------

  PURE FUNCTION mwwatercoeff_equal( x, y ) RESULT( is_equal )
    TYPE(mwwatercoeff_type), INTENT(IN) :: x, y
    LOGICAL :: is_equal

    ! Set up
    is_equal = .false.

    ! Check the object association status
    IF ( (.NOT. mwwatercoeff_associated(x)) .OR. &
         (.NOT. mwwatercoeff_associated(y))      ) RETURN

    ! Check contents
    ! ...Release/version info
    IF ( (x%Release /= y%Release) .OR. &
         (x%Version /= y%Version) ) RETURN
    ! ...Structures
    IF ( fitcoeff_associated( x%FCCoeff ) .NEQV. fitcoeff_associated( x%FCCoeff ) ) RETURN
    IF ( fitcoeff_associated( x%FCCoeff ) .AND.  fitcoeff_associated( x%FCCoeff ) ) THEN
      IF ( .NOT. (x%FCCoeff == y%FCCoeff) ) RETURN
    END IF
    IF ( fitcoeff_associated( x%FRCoeff ) .NEQV. fitcoeff_associated( x%FRCoeff ) ) RETURN
    IF ( fitcoeff_associated( x%FRCoeff ) .AND.  fitcoeff_associated( x%FRCoeff ) ) THEN
      IF ( .NOT. (x%FRCoeff == y%FRCoeff) ) RETURN
    END IF
    IF ( fitcoeff_associated( x%RCCoeff ) .NEQV. fitcoeff_associated( x%RCCoeff ) ) RETURN
    IF ( fitcoeff_associated( x%RCCoeff ) .AND.  fitcoeff_associated( x%RCCoeff ) ) THEN
      IF ( .NOT. (x%RCCoeff == y%RCCoeff) ) RETURN
    END IF
    IF ( fitcoeff_associated( x%AZCoeff ) .NEQV. fitcoeff_associated( x%AZCoeff ) ) RETURN
    IF ( fitcoeff_associated( x%AZCoeff ) .AND.  fitcoeff_associated( x%AZCoeff ) ) THEN
      IF ( .NOT. (x%AZCoeff == y%AZCoeff) ) RETURN
    END IF
    IF ( fitcoeff_associated( x%SSCCoeff ) .NEQV. fitcoeff_associated( x%SSCCoeff ) ) RETURN
    IF ( fitcoeff_associated( x%SSCCoeff ) .AND.  fitcoeff_associated( x%SSCCoeff ) ) THEN
      IF ( .NOT. (x%SSCCoeff == y%SSCCoeff) ) RETURN
    END IF
    IF ( fitcoeff_associated( x%LSCCoeff ) .NEQV. fitcoeff_associated( x%LSCCoeff ) ) RETURN
    IF ( fitcoeff_associated( x%LSCCoeff ) .AND.  fitcoeff_associated( x%LSCCoeff ) ) THEN
      IF ( .NOT. (x%LSCCoeff == y%LSCCoeff) ) RETURN
    END IF
    IF ( mwwaterlut_associated( x%LUT ) .NEQV. mwwaterlut_associated( x%LUT ) ) RETURN
    IF ( mwwaterlut_associated( x%LUT ) .AND.  mwwaterlut_associated( x%LUT ) ) THEN
      IF ( .NOT. (x%LUT == y%LUT) ) RETURN
    END IF

    ! If we get here everything is equal
    is_equal = .true.

  END FUNCTION mwwatercoeff_equal

END MODULE mwwatercoeff_define