AOvar_Define.f90

Go to the documentation of this file.

! 
! AOvar_Define
!
! Module defining the CRTM AtmOptics module internal
! variable object.
! 
!
! CREATION HISTORY:
!       Written by:     Paul van Delst, 01-Jul-2013
!                       paul.vandelst@noaa.gov
!                       

MODULE aovar_define

  ! -----------------
  ! Environment setup
  ! -----------------
  ! Module use
  USE type_kinds           , ONLY: fp
  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
  ! Disable implicit typing
  IMPLICIT NONE


  ! ------------
  ! Visibilities
  ! ------------
  ! Everything private by default
  PRIVATE
  ! Datatypes
  PUBLIC :: aovar_type
  ! Operators
  PUBLIC :: OPERATOR(==)
  ! Procedures
  PUBLIC :: aovar_associated
  PUBLIC :: aovar_destroy
  PUBLIC :: aovar_create
  PUBLIC :: aovar_inspect
  PUBLIC :: aovar_validrelease
  PUBLIC :: aovar_info
  PUBLIC :: aovar_defineversion
  PUBLIC :: aovar_inquirefile
  PUBLIC :: aovar_readfile
  PUBLIC :: aovar_writefile


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


  ! -----------------
  ! Module parameters
  ! -----------------
  CHARACTER(*), PARAMETER :: module_version_id = &
    '$Id: AOvar_Define.f90 60152 2015-08-13 19:19:13Z paul.vandelst@noaa.gov $'
  ! Release and version
  INTEGER, PARAMETER :: aovar_release = 1  ! This determines structure and file formats.
  INTEGER, PARAMETER :: aovar_version = 1  ! This is just the default data version.
  ! Close status for write errors
  CHARACTER(*), PARAMETER :: write_error_status = 'DELETE'
  ! Literal constants
  REAL(fp), PARAMETER :: zero = 0.0_fp
  REAL(fp), PARAMETER :: one  = 1.0_fp
  ! String lengths
  INTEGER,  PARAMETER :: ml = 256 ! Message length
  INTEGER,  PARAMETER :: sl =  80 ! String length

  
  ! ---------------------
  ! Structure definitions
  ! ---------------------
  ! The internal variable definition to hold information
  ! between FWD, TL, AD, and K-matrix calls
  TYPE :: aovar_type
    ! Allocation indicator
    LOGICAL :: is_allocated = .false.
    ! Release and version information
    INTEGER :: release = aovar_release
    INTEGER :: version = aovar_version
    ! Dimensions
    INTEGER :: n_layers = 0
    ! The total atmospheric transmittance
    REAL(fp) :: transmittance = zero
    ! The profile data
    REAL(fp), ALLOCATABLE :: optical_depth(:)
    REAL(fp), ALLOCATABLE :: bs(:)
    REAL(fp), ALLOCATABLE :: w(:)
  END TYPE aovar_type


CONTAINS


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

  ELEMENTAL FUNCTION aovar_associated( self ) RESULT( Status )
    TYPE(aovar_type), INTENT(IN) :: self
    LOGICAL :: status
    status = self%Is_Allocated
  END FUNCTION aovar_associated

 
  ELEMENTAL SUBROUTINE aovar_destroy( self )
    TYPE(aovar_type), INTENT(OUT) :: self
    self%Is_Allocated = .false.
    self%n_Layers         = 0
  END SUBROUTINE aovar_destroy
  
  
  ELEMENTAL SUBROUTINE aovar_create( &
    self    , &  ! Output
    n_Layers  )  ! Input
    ! Arguments
    TYPE(aovar_type), INTENT(OUT) :: self
    INTEGER         , INTENT(IN)  :: n_layers                
    ! Local variables
    INTEGER :: alloc_stat

    ! Check input
    IF ( n_layers < 1 ) RETURN

    ! Perform the allocation
    ALLOCATE( self%optical_depth(n_layers), &
              self%bs(n_layers), &
              self%w(n_layers), &
              stat = alloc_stat )
    IF ( alloc_stat /= 0 ) RETURN

    ! Initialise dimensions only!
    self%n_Layers = n_layers        

    ! Set allocation indicator
    self%Is_Allocated = .true.
  END SUBROUTINE aovar_create
  
  
  SUBROUTINE aovar_inspect( self)
    TYPE(aovar_type), INTENT(IN) :: self
    WRITE(*,'(1x,"AOvar OBJECT")')

    ! Release/version info
    WRITE(*,'(3x,"Release.Version     :",1x,i0,".",i0)') self%Release, self%Version

    ! Dimensions
    WRITE(*,'(3x,"n_Layers            :",1x,i0)') self%n_Layers        
    IF ( .NOT. aovar_associated(self) ) RETURN

    ! Data
    WRITE(*,'(3x,"Total transmittance :",1x,es13.6)') self%transmittance
    WRITE(*,'(3x,"Optical depth (sigma) :")')
    WRITE(*,'(5(1x,es13.6,:))') self%optical_depth
    WRITE(*,'(3x,"Volume scattering coefficient (bs) :")')
    WRITE(*,'(5(1x,es13.6,:))') self%bs
    WRITE(*,'(3x,"Single scatter albedo (w) :")')
    WRITE(*,'(5(1x,es13.6,:))') self%w
  END SUBROUTINE aovar_inspect


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

    ! Set up
    isvalid = .true.

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

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


  SUBROUTINE aovar_info( self, Info )
    ! Arguments
    TYPE(aovar_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,"AOvar RELEASE.VERSION: ",i2,".",i2.2,a,3x, &
           &"N_LAYERS=",i0 )' ) &
           achar(carriage_return)//achar(linefeed), &
           self%Release, self%Version, &
           achar(carriage_return)//achar(linefeed), &
           self%n_Layers      
                       
    ! Trim the output based on the
    ! dummy argument string length
    info = long_string(1:min(len(info), len_trim(long_string)))
  END SUBROUTINE aovar_info


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


  FUNCTION aovar_inquirefile( &
    Filename, &  ! Input
    n_Layers, &  ! Optional output  
    Release , &  ! Optional output
    Version , &  ! Optional output
    Title   , &  ! Optional output
    History , &  ! Optional output
    Comment ) &  ! Optional output
  result( err_stat )
    ! Arguments
    CHARACTER(*),           INTENT(IN)  :: filename
    INTEGER     , OPTIONAL, INTENT(OUT) :: n_layers        
    INTEGER     , OPTIONAL, INTENT(OUT) :: release        
    INTEGER     , OPTIONAL, INTENT(OUT) :: version        
    CHARACTER(*), OPTIONAL, INTENT(OUT) :: title           
    CHARACTER(*), OPTIONAL, INTENT(OUT) :: history         
    CHARACTER(*), OPTIONAL, INTENT(OUT) :: comment         
    ! Function result
    INTEGER :: err_stat
    ! Function parameters
    CHARACTER(*), PARAMETER :: routine_name = 'AOvar_InquireFile'
    ! Function variables
    CHARACTER(ML) :: msg
    CHARACTER(ML) :: io_msg
    INTEGER :: io_stat
    INTEGER :: fid
    TYPE(aovar_type) :: aovar

 
    ! Setup
    err_stat = success
    ! ...Check that the file exists
    IF ( .NOT. file_exists( filename ) ) THEN
      msg = 'File '//trim(filename)//' not found.'
      CALL inquire_cleanup(); RETURN
    END IF


    ! Open the file
    err_stat = open_binary_file( filename, fid )
    IF ( err_stat /= success ) THEN
      msg = 'Error opening '//trim(filename)
      CALL inquire_cleanup(); RETURN
    END IF


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


    ! Read the dimensions
    READ( fid, iostat=io_stat, iomsg=io_msg ) aovar%n_Layers      
    IF ( io_stat /= 0 ) THEN
      msg = 'Error reading dimension values from '//trim(filename)//' - '//trim(io_msg)
      CALL inquire_cleanup(); RETURN
    END IF


    ! 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 inquire_cleanup(); RETURN
    END IF


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


    ! Assign the return arguments
    IF ( PRESENT(n_layers) ) n_layers = aovar%n_Layers        
    IF ( PRESENT(release ) ) release  = aovar%Release        
    IF ( PRESENT(version ) ) version  = aovar%Version        
    
  CONTAINS
  
    SUBROUTINE inquire_cleanup()
      ! Close file if necessary
      IF ( file_open(fid) ) 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
      ! Set error status and print error message
      err_stat = failure
      CALL display_message( routine_name, msg, err_stat )
    END SUBROUTINE inquire_cleanup
    
  END FUNCTION aovar_inquirefile


  FUNCTION aovar_readfile( &
    AOvar     , &  ! 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(aovar_type),       INTENT(OUT) :: aovar
    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 = 'AOvar_ReadFile'
    ! Function variables
    CHARACTER(ML) :: msg
    CHARACTER(ML) :: io_msg
    LOGICAL :: close_file
    LOGICAL :: noisy
    INTEGER :: io_stat
    INTEGER :: fid
    TYPE(aovar_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. aovar_validrelease( dummy ) ) THEN
      msg = 'AOvar Release check failed.'
      CALL read_cleanup(); RETURN
    END IF


    ! Read the dimensions
    READ( fid, iostat=io_stat, iomsg=io_msg ) dummy%n_Layers
    IF ( io_stat /= 0 ) THEN
      msg = 'Error reading data dimensions - '//trim(io_msg)
      CALL read_cleanup(); RETURN
    END IF
    ! ...Allocate the object
    CALL aovar_create( aovar, dummy%n_Layers )
    IF ( .NOT. aovar_associated( aovar ) ) THEN
      msg = 'AOvar object allocation failed.'
      CALL read_cleanup(); RETURN
    END IF
    ! ...Explicitly assign the version number
    aovar%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 data
    ! ...total transmittance
    READ( fid, iostat=io_stat, iomsg=io_msg ) aovar%transmittance
    IF ( io_stat /= 0 ) THEN
      msg = 'Error reading total transmittance - '//trim(io_msg)
      CALL read_cleanup(); RETURN
    END IF
    ! ...optical depth
    READ( fid, iostat=io_stat, iomsg=io_msg ) aovar%optical_depth
    IF ( io_stat /= 0 ) THEN
      msg = 'Error reading optical depth - '//trim(io_msg)
      CALL read_cleanup(); RETURN
    END IF
    ! ...volume scattering coefficient
    READ( fid, iostat=io_stat, iomsg=io_msg ) aovar%bs
    IF ( io_stat /= 0 ) THEN
      msg = 'Error reading volume scattering coefficient - '//trim(io_msg)
      CALL read_cleanup(); RETURN
    END IF
    ! ...single scatter albedo
    READ( fid, iostat=io_stat, iomsg=io_msg ) aovar%w
    IF ( io_stat /= 0 ) THEN
      msg = 'Error reading single scatter albedo - '//trim(io_msg)
      CALL read_cleanup(); RETURN
    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 aovar_info( aovar, 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 aovar_destroy( aovar )
       err_stat = failure
       CALL display_message( routine_name, msg, err_stat )
     END SUBROUTINE read_cleanup

  END FUNCTION aovar_readfile


  FUNCTION aovar_writefile( &
    AOvar     , &  ! 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(aovar_type),  INTENT(IN) :: aovar
    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 = 'AOvar_WriteFile'
    ! Function variables
    CHARACTER(ML) :: msg
    CHARACTER(ML) :: io_msg
    LOGICAL :: close_file
    LOGICAL :: noisy
    INTEGER :: io_stat
    INTEGER :: fid
    

    ! 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. aovar_associated( aovar ) ) THEN
      msg = 'AOvar 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 ) &
      aovar%Release, &
      aovar%Version
    IF ( io_stat /= 0 ) THEN
      msg = 'Error writing Release/Version - '//trim(io_msg)
      CALL write_cleanup(); RETURN
    END IF


    ! Write the dimensions
    WRITE( fid, iostat=io_stat, iomsg=io_msg ) aovar%n_Layers
    IF ( io_stat /= 0 ) THEN
      msg = 'Error writing data dimensions - '//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 data
    ! ...total transmittance
    WRITE( fid, iostat=io_stat, iomsg=io_msg ) aovar%transmittance
    IF ( io_stat /= 0 ) THEN
      msg = 'Error writing total transmittance - '//trim(io_msg)
      CALL write_cleanup(); RETURN
    END IF
    ! ...optical depth
    WRITE( fid, iostat=io_stat, iomsg=io_msg ) aovar%optical_depth
    IF ( io_stat /= 0 ) THEN
      msg = 'Error writing optical depth - '//trim(io_msg)
      CALL write_cleanup(); RETURN
    END IF
    ! ...volume scattering coefficient
    WRITE( fid, iostat=io_stat, iomsg=io_msg ) aovar%bs
    IF ( io_stat /= 0 ) THEN
      msg = 'Error writing volume scattering coefficient - '//trim(io_msg)
      CALL write_cleanup(); RETURN
    END IF
    ! ...single scatter albedo
    WRITE( fid, iostat=io_stat, iomsg=io_msg ) aovar%w
    IF ( io_stat /= 0 ) THEN
      msg = 'Error writing single scatter albedo - '//trim(io_msg)
      CALL write_cleanup(); RETURN
    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 aovar_info( aovar, 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, 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 aovar_writefile

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

  ELEMENTAL FUNCTION aovar_equal( x, y ) RESULT( is_equal )
    TYPE(aovar_type), INTENT(IN) :: x, y
    LOGICAL :: is_equal

    ! Set up
    is_equal = .false.
   
    ! Check the object association status
    IF ( (.NOT. aovar_associated(x)) .OR. &
         (.NOT. aovar_associated(y))      ) RETURN

   ! Check contents
    ! ...Release/version info
    IF ( (x%Release /= y%Release) .OR. &
         (x%Version /= y%Version) ) RETURN
    ! ...Dimensions
    IF ( (x%n_Layers /= y%n_Layers ) ) RETURN
    ! ...Data
    IF (    (x%transmittance .equalto. y%transmittance ) .AND. &
         all(x%optical_depth .equalto. y%optical_depth ) .AND. &
         all(x%bs            .equalto. y%bs            ) .AND. &
         all(x%w             .equalto. y%w             ) ) &
      is_equal = .true.
  END FUNCTION aovar_equal

END MODULE aovar_define