CRTM_LifeCycle.f90

Go to the documentation of this file.

!
! CRTM_LifeCycle
!
! Module containing CRTM life cycle functions to initialize and destroy
! the CRTM space.
!
! Written by:     Paul van Delst, 21-May-2004
!                 paul.vandelst@noaa.gov
!

MODULE crtm_lifecycle


  ! -----------------
  ! Environment setup
  ! -----------------
  ! Module usage
  USE message_handler
  USE crtm_channelinfo_define, ONLY: crtm_channelinfo_type, &
                                     crtm_channelinfo_associated, &
                                     crtm_channelinfo_destroy, &
                                     crtm_channelinfo_create
  ! ...Spectral coefficients
  USE crtm_spccoeff          , ONLY: sc, &
                                     crtm_spccoeff_load   , &
                                     crtm_spccoeff_destroy, &
                                     spccoeff_ismicrowavesensor  , &
                                     spccoeff_isinfraredsensor   , &
                                     spccoeff_isvisiblesensor    , &
                                     spccoeff_isultravioletsensor
  ! ...Transmittance model coefficients
  USE crtm_taucoeff
  ! ...Aerosol optical properties
  USE crtm_aerosolcoeff      , ONLY: crtm_aerosolcoeff_load, &
                                     crtm_aerosolcoeff_destroy
  ! ...Cloud optical properties
  USE crtm_cloudcoeff        , ONLY: crtm_cloudcoeff_load, &
                                     crtm_cloudcoeff_destroy
  ! ...Infrared surface emissivities
  USE crtm_irwatercoeff      , ONLY: crtm_irwatercoeff_load, &
                                     crtm_irwatercoeff_destroy
  USE crtm_irlandcoeff       , ONLY: crtm_irlandcoeff_load, &
                                     crtm_irlandcoeff_destroy
  USE crtm_irsnowcoeff       , ONLY: crtm_irsnowcoeff_load, &
                                     crtm_irsnowcoeff_destroy
  USE crtm_iricecoeff        , ONLY: crtm_iricecoeff_load, &
                                     crtm_iricecoeff_destroy
  ! ...Visible surface emissivities
  USE crtm_viswatercoeff     , ONLY: crtm_viswatercoeff_load, &
                                     crtm_viswatercoeff_destroy
  USE crtm_vislandcoeff      , ONLY: crtm_vislandcoeff_load, &
                                     crtm_vislandcoeff_destroy
  USE crtm_vissnowcoeff      , ONLY: crtm_vissnowcoeff_load, &
                                     crtm_vissnowcoeff_destroy
  USE crtm_visicecoeff       , ONLY: crtm_visicecoeff_load, &
                                     crtm_visicecoeff_destroy
  ! ...Microwave surface emissivities
  USE crtm_mwwatercoeff      , ONLY: crtm_mwwatercoeff_load, &
                                     crtm_mwwatercoeff_destroy
  ! Disable all implicit typing
  IMPLICIT NONE


  ! ------------
  ! Visibilities
  ! ------------
  ! Everything private by default
  PRIVATE
  ! Public procedures
  PUBLIC :: crtm_init
  PUBLIC :: crtm_destroy
  PUBLIC :: crtm_isinitialized
  PUBLIC :: crtm_lifecycleversion

  ! -----------------
  ! Module parameters
  ! -----------------
  CHARACTER(*), PARAMETER :: module_version_id = &
  '$Id: CRTM_LifeCycle.f90 60152 2015-08-13 19:19:13Z paul.vandelst@noaa.gov $'
  ! String lengths
  INTEGER, PARAMETER :: ml = 256   ! Error message length
  INTEGER, PARAMETER :: sl = 5000  ! Maximum length for path+filenames


CONTAINS


!------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       CRTM_Init
!
! PURPOSE:
!       Function to initialise the CRTM.
!
! CALLING SEQUENCE:
!       Error_Status = CRTM_Init( Sensor_ID  , &
!                                 ChannelInfo, &
!                                 CloudCoeff_File    = CloudCoeff_File   , &
!                                 AerosolCoeff_File  = AerosolCoeff_File , &
!                                 Load_CloudCoeff    = Load_CloudCoeff   , &
!                                 Load_AerosolCoeff  = Load_AerosolCoeff , &
!                                 IRwaterCoeff_File  = IRwaterCoeff_File , &
!                                 IRlandCoeff_File   = IRlandCoeff_File  , &
!                                 IRsnowCoeff_File   = IRsnowCoeff_File  , &
!                                 IRiceCoeff_File    = IRiceCoeff_File   , &
!                                 VISwaterCoeff_File = VISwaterCoeff_File, &
!                                 VISlandCoeff_File  = VISlandCoeff_File , &
!                                 VISsnowCoeff_File  = VISsnowCoeff_File , &
!                                 VISiceCoeff_File   = VISiceCoeff_File  , &
!                                 MWwaterCoeff_File  = MWwaterCoeff_File , &
!                                 File_Path          = File_Path         , &
!                                 Quiet              = Quiet             , &
!                                 Process_ID         = Process_ID        , &
!                                 Output_Process_ID  = Output_Process_ID   )
!
! INPUTS:
!       Sensor_ID:          List of the sensor IDs (e.g. hirs3_n17, amsua_n18,
!                           ssmis_f16, etc) with which the CRTM is to be
!                           initialised. These sensor ids are used to construct
!                           the sensor specific SpcCoeff and TauCoeff filenames
!                           containing the necessary coefficient data, i.e.
!                             <Sensor_ID>.SpcCoeff.bin
!                           and
!                             <Sensor_ID>.TauCoeff.bin
!                           for each sensor Id in the list.
!                           UNITS:      N/A
!                           TYPE:       CHARACTER(*)
!                           DIMENSION:  Rank-1 (n_Sensors)
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
! OUTPUTS:
!       ChannelInfo:        ChannelInfo structure array populated based on
!                           the contents of the coefficient files and the
!                           user inputs.
!                           UNITS:      N/A
!                           TYPE:       CRTM_ChannelInfo_type
!                           DIMENSION:  Same as input Sensor_Id argument
!                           ATTRIBUTES: INTENT(OUT)
!
! OPTIONAL INPUTS:
!       CloudCoeff_File:    Name of the data file containing the cloud optical
!                           properties data for scattering calculations.
!                           Available datafiles:
!                           - CloudCoeff.bin  [DEFAULT]
!                           UNITS:      N/A
!                           TYPE:       CHARACTER(*)
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       AerosolCoeff_File:  Name of the data file containing the aerosol optical
!                           properties data for scattering calculations.
!                           Available datafiles:
!                           - AerosolCoeff.bin  [DEFAULT]
!                           UNITS:      N/A
!                           TYPE:       CHARACTER(*)
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       Load_CloudCoeff:    Set this logical argument for not loading the CloudCoeff data
!                           to save memory space under the clear conditions
!                           If == .FALSE., the CloudCoeff data will not be loaded;
!                              == .TRUE.,  the CloudCoeff data will be loaded.
!                           If not specified, default is .TRUE. (will be loaded)
!                           UNITS:      N/A
!                           TYPE:       LOGICAL
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       Load_AerosolCoeff:  Set this logical argument for not loading the AerosolCoeff data
!                           to save memory space under the clear conditions
!                           If == .FALSE., the AerosolCoeff data will not be loaded;
!                              == .TRUE.,  the AerosolCoeff data will be loaded.
!                           If not specified, default is .TRUE. (will be loaded)
!                           UNITS:      N/A
!                           TYPE:       LOGICAL
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       MWwaterCoeff_File:  Name of the data file containing the coefficient
!                           data for the microwave water emissivity model.
!                           Available datafiles:
!                           - FASTEM6.MWwater.EmisCoeff.bin  [DEFAULT]
!                           - FASTEM5.MWwater.EmisCoeff.bin
!                           - FASTEM4.MWwater.EmisCoeff.bin
!                           UNITS:      N/A
!                           TYPE:       CHARACTER(*)
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       IRwaterCoeff_File:  Name of the data file containing the coefficient
!                           data for the infrared water emissivity model.
!                           Available datafiles:
!                           - Nalli.IRwater.EmisCoeff.bin  [DEFAULT]
!                           - WuSmith.IRwater.EmisCoeff.bin
!                           If not specified the Nalli datafile is read.
!                           UNITS:      N/A
!                           TYPE:       CHARACTER(*)
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       IRlandCoeff_File:   Name of the data file containing the coefficient
!                           data for the infrared land emissivity model.
!                           Available datafiles:
!                           - NPOESS.IRland.EmisCoeff.bin  [DEFAULT]
!                           - IGBP.IRland.EmisCoeff.bin
!                           - USGS.IRland.EmisCoeff.bin
!                           UNITS:      N/A
!                           TYPE:       CHARACTER(*)
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       IRsnowCoeff_File:   Name of the data file containing the coefficient
!                           data for the infrared snow emissivity model.
!                           Available datafiles:
!                           - NPOESS.IRsnow.EmisCoeff.bin  [DEFAULT]
!                           - IGBP.IRsnow.EmisCoeff.bin
!                           - USGS.IRsnow.EmisCoeff.bin
!                           UNITS:      N/A
!                           TYPE:       CHARACTER(*)
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       IRiceCoeff_File:    Name of the data file containing the coefficient
!                           data for the infrared ice emissivity model.
!                           Available datafiles:
!                           - NPOESS.IRice.EmisCoeff.bin  [DEFAULT]
!                           - IGBP.IRice.EmisCoeff.bin
!                           - USGS.IRice.EmisCoeff.bin
!                           UNITS:      N/A
!                           TYPE:       CHARACTER(*)
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       VISwaterCoeff_File: Name of the data file containing the coefficient
!                           data for the visible water emissivity model.
!                           Available datafiles:
!                           - NPOESS.VISwater.EmisCoeff.bin  [DEFAULT]
!                           - IGBP.VISwater.EmisCoeff.bin
!                           - USGS.VISwater.EmisCoeff.bin
!                           UNITS:      N/A
!                           TYPE:       CHARACTER(*)
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       VISlandCoeff_File:  Name of the data file containing the coefficient
!                           data for the visible land emissivity model.
!                           Available datafiles:
!                           - NPOESS.VISland.EmisCoeff.bin  [DEFAULT]
!                           - IGBP.VISland.EmisCoeff.bin
!                           - USGS.VISland.EmisCoeff.bin
!                           UNITS:      N/A
!                           TYPE:       CHARACTER(*)
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       VISsnowCoeff_File:  Name of the data file containing the coefficient
!                           data for the visible snow emissivity model.
!                           Available datafiles:
!                           - NPOESS.VISsnow.EmisCoeff.bin  [DEFAULT]
!                           - IGBP.VISsnow.EmisCoeff.bin
!                           - USGS.VISsnow.EmisCoeff.bin
!                           UNITS:      N/A
!                           TYPE:       CHARACTER(*)
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       VISiceCoeff_File:   Name of the data file containing the coefficient
!                           data for the visible ice emissivity model.
!                           Available datafiles:
!                           - NPOESS.VISice.EmisCoeff.bin  [DEFAULT]
!                           - IGBP.VISice.EmisCoeff.bin
!                           - USGS.VISice.EmisCoeff.bin
!                           UNITS:      N/A
!                           TYPE:       CHARACTER(*)
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       File_Path:          Character string specifying a file path for the
!                           input data files. If not specified, the current
!                           directory is the default.
!                           UNITS:      N/A
!                           TYPE:       CHARACTER(*)
!                           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
!
!       Process_ID:         Set this argument to the MPI process ID that this
!                           function call is running under. This value is used
!                           solely for controlling INFORMATION message output.
!                           If MPI is not being used, ignore this argument.
!                           This argument is ignored if the Quiet argument is set.
!                           UNITS:      N/A
!                           TYPE:       INTEGER
!                           DIMENSION:  Scalar
!                           ATTRIBUTES: INTENT(IN), OPTIONAL
!
!       Output_Process_ID:  Set this argument to the MPI process ID in which
!                           all INFORMATION messages are to be output. If
!                           the passed Process_ID value agrees with this value
!                           the INFORMATION messages are output.
!                           This argument is ignored if the Quiet argument
!                           is set.
!                           UNITS:      N/A
!                           TYPE:       INTEGER
!                           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 CRTM initialisation was successful
!                              == FAILURE an unrecoverable error occurred.
!                           UNITS:      N/A
!                           TYPE:       INTEGER
!                           DIMENSION:  Scalar
!
! SIDE EFFECTS:
!       All public data arrays accessed by this module and its dependencies
!       are overwritten.
!
!:sdoc-:
!------------------------------------------------------------------------------

  FUNCTION crtm_init( &
    Sensor_ID         , &  ! Input
    ChannelInfo       , &  ! Output
    CloudCoeff_File   , &  ! Optional input
    AerosolCoeff_File , &  ! Optional input
    EmisCoeff_File    , &  ! Optional input  ! *** DEPRECATED. Replaced by IRwaterCoeff_File
    IRwaterCoeff_File , &  ! Optional input
    IRlandCoeff_File  , &  ! Optional input
    IRsnowCoeff_File  , &  ! Optional input
    IRiceCoeff_File   , &  ! Optional input
    VISwaterCoeff_File, &  ! Optional input
    VISlandCoeff_File , &  ! Optional input
    VISsnowCoeff_File , &  ! Optional input
    VISiceCoeff_File  , &  ! Optional input
    MWwaterCoeff_File , &  ! Optional input
    File_Path         , &  ! Optional input
    Load_CloudCoeff   , &  ! Optional input
    Load_AerosolCoeff , &  ! Optional input
    Quiet             , &  ! Optional input
    Process_ID        , &  ! Optional input
    Output_Process_ID ) &  ! Optional input
  result( err_stat )
    ! Arguments
    CHARACTER(*)               , INTENT(IN)  :: sensor_id(:)
    TYPE(crtm_channelinfo_type), INTENT(OUT) :: channelinfo(:)
    CHARACTER(*),      OPTIONAL, INTENT(IN)  :: cloudcoeff_file
    CHARACTER(*),      OPTIONAL, INTENT(IN)  :: aerosolcoeff_file
    CHARACTER(*),      OPTIONAL, INTENT(IN)  :: emiscoeff_file  ! *** DEPRECATED. Replaced by IRwaterCoeff_File
    CHARACTER(*),      OPTIONAL, INTENT(IN)  :: irwatercoeff_file
    CHARACTER(*),      OPTIONAL, INTENT(IN)  :: irlandcoeff_file
    CHARACTER(*),      OPTIONAL, INTENT(IN)  :: irsnowcoeff_file
    CHARACTER(*),      OPTIONAL, INTENT(IN)  :: iricecoeff_file
    CHARACTER(*),      OPTIONAL, INTENT(IN)  :: viswatercoeff_file
    CHARACTER(*),      OPTIONAL, INTENT(IN)  :: vislandcoeff_file
    CHARACTER(*),      OPTIONAL, INTENT(IN)  :: vissnowcoeff_file
    CHARACTER(*),      OPTIONAL, INTENT(IN)  :: visicecoeff_file
    CHARACTER(*),      OPTIONAL, INTENT(IN)  :: mwwatercoeff_file
    CHARACTER(*),      OPTIONAL, INTENT(IN)  :: file_path
    LOGICAL     ,      OPTIONAL, INTENT(IN)  :: load_cloudcoeff
    LOGICAL     ,      OPTIONAL, INTENT(IN)  :: load_aerosolcoeff
    LOGICAL     ,      OPTIONAL, INTENT(IN)  :: quiet
    INTEGER     ,      OPTIONAL, INTENT(IN)  :: process_id
    INTEGER     ,      OPTIONAL, INTENT(IN)  :: output_process_id
    ! Function result
    INTEGER :: err_stat
    ! Local parameters
    CHARACTER(*), PARAMETER :: routine_name = 'CRTM_Init'
    ! Local variables
    CHARACTER(ML) :: msg, pid_msg
    CHARACTER(SL) :: default_cloudcoeff_file
    CHARACTER(SL) :: default_aerosolcoeff_file
    CHARACTER(SL) :: default_irwatercoeff_file
    CHARACTER(SL) :: default_irlandcoeff_file
    CHARACTER(SL) :: default_irsnowcoeff_file
    CHARACTER(SL) :: default_iricecoeff_file
    CHARACTER(SL) :: default_viswatercoeff_file
    CHARACTER(SL) :: default_vislandcoeff_file
    CHARACTER(SL) :: default_vissnowcoeff_file
    CHARACTER(SL) :: default_visicecoeff_file
    CHARACTER(SL) :: default_mwwatercoeff_file

    INTEGER :: l, n, n_sensors
    LOGICAL :: local_load_cloudcoeff
    LOGICAL :: local_load_aerosolcoeff

    ! ******
    ! TEMPORARY UNTIL LOAD ROUTINE INTERFACES HAVE BEEN MODIFIED
    INTEGER :: iquiet
    iquiet = 0
    IF ( PRESENT(quiet) ) THEN
      IF ( quiet ) THEN
        iquiet = 1  ! Set
      END IF
    END IF
    ! ******

    ! Set up
    err_stat = success
    ! ...Create a process ID message tag for error messages
    IF ( PRESENT(process_id) ) THEN
      WRITE( pid_msg,'("; Process ID: ",i0)' ) process_id
    ELSE
      pid_msg = ''
    END IF
    ! ...Check coefficient loading flags
    local_load_cloudcoeff = .true.
    IF( PRESENT(load_cloudcoeff) ) local_load_cloudcoeff = load_cloudcoeff
    local_load_aerosolcoeff = .true.
    IF( PRESENT(load_aerosolcoeff) ) local_load_aerosolcoeff = load_aerosolcoeff
    ! ...Check dimensionality
    n_sensors = SIZE(sensor_id)
    IF ( SIZE(channelinfo) /= n_sensors ) THEN
      err_stat = failure
      msg = 'Inconsistent Sensor_ID and ChannelInfo dimensions'
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
      RETURN
    END IF
    ! ...Check for deprecated EmisCoeff_File argument
    IF ( PRESENT(emiscoeff_file) ) THEN
      err_stat = failure
      msg = 'The EmisCoeff_File optional argument is deprecated. Use IRwaterCoeff_File instead'
      CALL display_message( routine_name, trim(msg)//trim(pid_msg), err_stat )
      RETURN
    END IF
      

    ! Specify sensor-independent coefficient filenames
    ! ...Default filenames
    default_cloudcoeff_file    = 'CloudCoeff.bin'
    default_aerosolcoeff_file  = 'AerosolCoeff.bin'
    default_irwatercoeff_file  = 'Nalli.IRwater.EmisCoeff.bin'
    default_irlandcoeff_file   = 'NPOESS.IRland.EmisCoeff.bin'
    default_irsnowcoeff_file   = 'NPOESS.IRsnow.EmisCoeff.bin'
    default_iricecoeff_file    = 'NPOESS.IRice.EmisCoeff.bin'
    default_viswatercoeff_file = 'NPOESS.VISwater.EmisCoeff.bin'
    default_vislandcoeff_file  = 'NPOESS.VISland.EmisCoeff.bin'
    default_vissnowcoeff_file  = 'NPOESS.VISsnow.EmisCoeff.bin'
    default_visicecoeff_file   = 'NPOESS.VISice.EmisCoeff.bin'
    default_mwwatercoeff_file  = 'FASTEM6.MWwater.EmisCoeff.bin'
    ! ...Were other filenames specified?
    IF ( PRESENT(cloudcoeff_file   ) ) default_cloudcoeff_file    = trim(adjustl(cloudcoeff_file))
    IF ( PRESENT(aerosolcoeff_file ) ) default_aerosolcoeff_file  = trim(adjustl(aerosolcoeff_file))
    IF ( PRESENT(irwatercoeff_file ) ) default_irwatercoeff_file  = trim(adjustl(irwatercoeff_file))
    IF ( PRESENT(irlandcoeff_file  ) ) default_irlandcoeff_file   = trim(adjustl(irlandcoeff_file))
    IF ( PRESENT(irsnowcoeff_file  ) ) default_irsnowcoeff_file   = trim(adjustl(irsnowcoeff_file))
    IF ( PRESENT(iricecoeff_file   ) ) default_iricecoeff_file    = trim(adjustl(iricecoeff_file))
    IF ( PRESENT(viswatercoeff_file) ) default_viswatercoeff_file = trim(adjustl(viswatercoeff_file))
    IF ( PRESENT(vislandcoeff_file ) ) default_vislandcoeff_file  = trim(adjustl(vislandcoeff_file))
    IF ( PRESENT(vissnowcoeff_file ) ) default_vissnowcoeff_file  = trim(adjustl(vissnowcoeff_file))
    IF ( PRESENT(visicecoeff_file  ) ) default_visicecoeff_file   = trim(adjustl(visicecoeff_file))
    IF ( PRESENT(mwwatercoeff_file ) ) default_mwwatercoeff_file  = trim(adjustl(mwwatercoeff_file))
    ! ...Was a path specified?
    IF ( PRESENT(file_path) ) THEN
      default_cloudcoeff_file    = trim(adjustl(file_path)) // trim(default_cloudcoeff_file)
      default_aerosolcoeff_file  = trim(adjustl(file_path)) // trim(default_aerosolcoeff_file)
      default_irwatercoeff_file  = trim(adjustl(file_path)) // trim(default_irwatercoeff_file)
      default_irlandcoeff_file   = trim(adjustl(file_path)) // trim(default_irlandcoeff_file)
      default_irsnowcoeff_file   = trim(adjustl(file_path)) // trim(default_irsnowcoeff_file)
      default_iricecoeff_file    = trim(adjustl(file_path)) // trim(default_iricecoeff_file)
      default_viswatercoeff_file = trim(adjustl(file_path)) // trim(default_viswatercoeff_file)
      default_vislandcoeff_file  = trim(adjustl(file_path)) // trim(default_vislandcoeff_file)
      default_vissnowcoeff_file  = trim(adjustl(file_path)) // trim(default_vissnowcoeff_file)
      default_visicecoeff_file   = trim(adjustl(file_path)) // trim(default_visicecoeff_file)
      default_mwwatercoeff_file  = trim(adjustl(file_path)) // trim(default_mwwatercoeff_file)
    END IF


    ! Load the spectral coefficients
    err_stat = crtm_spccoeff_load( &
                 sensor_id                            , &
                 file_path         = file_path        , &
                 quiet             = quiet            , &
                 process_id        = process_id       , &
                 output_process_id = output_process_id  )
    IF ( err_stat /= success ) THEN
      CALL display_message( routine_name,'Error loading SpcCoeff data'//trim(pid_msg),err_stat )
      RETURN
    END IF


    ! Load the transmittance model coefficients
    err_stat = crtm_load_taucoeff( &
                 sensor_id         = sensor_id        , &
                 file_path         = file_path        , &
                 quiet             = iquiet           , &  ! *** Use of iQuiet temporary
                 process_id        = process_id       , &
                 output_process_id = output_process_id  )
    IF ( err_stat /= success ) THEN
      CALL display_message( routine_name,'Error loading TauCoeff data'//trim(pid_msg),err_stat )
      RETURN
    END IF


    ! Load the cloud coefficients
    IF ( local_load_cloudcoeff ) THEN
      err_stat = crtm_cloudcoeff_load( &
                   default_cloudcoeff_file, &
                   quiet             = quiet            , &
                   process_id        = process_id       , &
                   output_process_id = output_process_id  )
      IF ( err_stat /= success ) THEN
        msg = 'Error loading CloudCoeff data from '//trim(default_cloudcoeff_file)
        CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
        RETURN
      END IF
    END IF


    ! Load the aerosol coefficients
    IF ( local_load_aerosolcoeff ) THEN
      err_stat = crtm_aerosolcoeff_load( &
                   default_aerosolcoeff_file, &
                   quiet             = quiet            , &
                   process_id        = process_id       , &
                   output_process_id = output_process_id  )
      IF ( err_stat /= success ) THEN
        msg = 'Error loading AerosolCoeff data from '//trim(default_aerosolcoeff_file)
        CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
        RETURN
      END IF
    END IF


    ! Load the emissivity model coefficients
    ! ...Infrared
    infrared_sensor: IF ( any(spccoeff_isinfraredsensor(sc)) ) THEN
      ! ...IR land
      err_stat = crtm_irlandcoeff_load( &
                   default_irlandcoeff_file, &
                   quiet             = quiet            , &
                   process_id        = process_id       , &
                   output_process_id = output_process_id  )
      IF ( err_stat /= success ) THEN
        msg = 'Error loading IRlandCoeff data from '//trim(default_irlandcoeff_file)
        CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
        RETURN
      END IF
      ! ...IR Water
      err_stat = crtm_irwatercoeff_load( &
                   default_irwatercoeff_file, &
                   quiet             = quiet            , &
                   process_id        = process_id       , &
                   output_process_id = output_process_id  )
      IF ( err_stat /= success ) THEN
        msg = 'Error loading IRwaterCoeff data from '//trim(default_irwatercoeff_file)
        CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
        RETURN
      END IF
      ! ...IR snow
      err_stat = crtm_irsnowcoeff_load( &
                   default_irsnowcoeff_file, &
                   quiet             = quiet            , &
                   process_id        = process_id       , &
                   output_process_id = output_process_id  )
      IF ( err_stat /= success ) THEN
        msg = 'Error loading IRsnowCoeff data from '//trim(default_irsnowcoeff_file)
        CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
        RETURN
      END IF
      ! ...IR ice
      err_stat = crtm_iricecoeff_load( &
                   default_iricecoeff_file, &
                   quiet             = quiet            , &
                   process_id        = process_id       , &
                   output_process_id = output_process_id  )
      IF ( err_stat /= success ) THEN
        msg = 'Error loading IRiceCoeff data from '//trim(default_iricecoeff_file)
        CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
        RETURN
      END IF
    END IF infrared_sensor
    
    ! ...Visible
    visible_sensor: IF ( any(spccoeff_isvisiblesensor(sc)) ) THEN
      ! ...VIS land
      err_stat = crtm_vislandcoeff_load( &
                   default_vislandcoeff_file, &
                   quiet             = quiet            , &
                   process_id        = process_id       , &
                   output_process_id = output_process_id  )
      IF ( err_stat /= success ) THEN
        msg = 'Error loading VISlandCoeff data from '//trim(default_vislandcoeff_file)
        CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
        RETURN
      END IF
      ! ...VIS water
      err_stat = crtm_viswatercoeff_load( &
                   default_viswatercoeff_file, &
                   quiet             = quiet            , &
                   process_id        = process_id       , &
                   output_process_id = output_process_id  )
      IF ( err_stat /= success ) THEN
        msg = 'Error loading VISwaterCoeff data from '//trim(default_viswatercoeff_file)
        CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
        RETURN
      END IF
      ! ...VIS snow
      err_stat = crtm_vissnowcoeff_load( &
                   default_vissnowcoeff_file, &
                   quiet             = quiet            , &
                   process_id        = process_id       , &
                   output_process_id = output_process_id  )
      IF ( err_stat /= success ) THEN
        msg = 'Error loading VISsnowCoeff data from '//trim(default_vissnowcoeff_file)
        CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
        RETURN
      END IF
      ! ...VIS ice
      err_stat = crtm_visicecoeff_load( &
                   default_visicecoeff_file, &
                   quiet             = quiet            , &
                   process_id        = process_id       , &
                   output_process_id = output_process_id  )
      IF ( err_stat /= success ) THEN
        msg = 'Error loading VISiceCoeff data from '//trim(default_visicecoeff_file)
        CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
        RETURN
      END IF
    END IF visible_sensor
    
    ! ...Microwave
    microwave_sensor: IF ( any(spccoeff_ismicrowavesensor(sc)) ) THEN
      ! ...MW water
      err_stat = crtm_mwwatercoeff_load( &
                   default_mwwatercoeff_file, &
                   quiet             = quiet            , &
                   process_id        = process_id       , &
                   output_process_id = output_process_id  )
      IF ( err_stat /= success ) THEN
        msg = 'Error loading MWwaterCoeff data from '//trim(default_mwwatercoeff_file)
        CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
        RETURN
      END IF
    END IF microwave_sensor


    ! Load the ChannelInfo structure
    DO n = 1, n_sensors
      ! ...Allocate the ChannelInfo structure
      CALL crtm_channelinfo_create( channelinfo(n), sc(n)%n_Channels )
      IF ( .NOT. crtm_channelinfo_associated(channelinfo(n)) ) THEN
        err_stat = failure
        msg = 'ChannelInfo allocation failed for '//trim(sensor_id(n))//' sensor'
        CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
        RETURN
      END IF
      ! ...Fill the structure for the current sensor
      channelinfo(n)%Sensor_Index     = n
      channelinfo(n)%Channel_Index    = (/(l, l=1,sc(n)%n_Channels)/)
      channelinfo(n)%Sensor_ID        = sc(n)%Sensor_Id
      channelinfo(n)%Sensor_Type      = sc(n)%Sensor_Type
      channelinfo(n)%WMO_Satellite_ID = sc(n)%WMO_Satellite_ID
      channelinfo(n)%WMO_Sensor_ID    = sc(n)%WMO_Sensor_ID
      channelinfo(n)%Sensor_Channel   = sc(n)%Sensor_Channel
    END DO

  END FUNCTION crtm_init


!------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       CRTM_Destroy
!
! PURPOSE:
!       Function to deallocate all the shared data arrays allocated and
!       populated during the CRTM initialization.
!
! CALLING SEQUENCE:
!       Error_Status = CRTM_Destroy( ChannelInfo            , &
!                                    Process_ID = Process_ID  )
!
! OUTPUTS:
!       ChannelInfo:  Reinitialized ChannelInfo structure.
!                     UNITS:      N/A
!                     TYPE:       CRTM_ChannelInfo_type
!                     DIMENSION:  Rank-1
!                     ATTRIBUTES: INTENT(IN OUT)
!
! OPTIONAL INPUTS:
!       Process_ID:   Set this argument to the MPI process ID that this
!                     function call is running under. This value is used
!                     solely for controlling message output. If MPI is not
!                     being used, ignore this argument.
!                     UNITS:      N/A
!                     TYPE:       INTEGER
!                     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 CRTM deallocations were successful
!                        == FAILURE an unrecoverable error occurred.
!                     UNITS:      N/A
!                     TYPE:       INTEGER
!                     DIMENSION:  Scalar
!
! SIDE EFFECTS:
!       All CRTM shared data arrays and structures are deallocated.
!
! COMMENTS:
!       Note the INTENT on the output ChannelInfo argument is IN OUT rather than
!       just OUT. This is necessary because the argument may be defined upon
!       input. To prevent memory leaks, the IN OUT INTENT is a must.
!
!:sdoc-:
!------------------------------------------------------------------------------

  FUNCTION crtm_destroy( &
    ChannelInfo, &  ! Output
    Process_ID ) &  ! Optional input
  result( err_stat )
    ! Arguments
    TYPE(crtm_channelinfo_type), INTENT(IN OUT) :: channelinfo(:)
    INTEGER,           OPTIONAL, INTENT(IN)     :: process_id
    ! Function result
    INTEGER :: err_stat
    ! Local parameters
    CHARACTER(*), PARAMETER :: routine_name = 'CRTM_Destroy'
    ! Local variables
    CHARACTER(ML) :: msg, pid_msg
    INTEGER :: destroy_status

    ! Set up
    err_stat = success
    ! ...Create a process ID message tag for error messages
    IF ( PRESENT(process_id) ) THEN
      WRITE( pid_msg,'("; Process ID: ",i0)' ) process_id
    ELSE
      pid_msg = ''
    END IF


    ! Destroy all the ChannelInfo structures
    CALL crtm_channelinfo_destroy( channelinfo )
    IF ( any(crtm_channelinfo_associated(channelinfo)) ) THEN
      err_stat = failure
      msg = 'Error deallocating ChannelInfo structure(s)'
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF


    ! Destroy the shared data structure
    destroy_status = crtm_visicecoeff_destroy( process_id = process_id )
    IF ( destroy_status /= success ) THEN
      err_stat = destroy_status
      msg = 'Error deallocating shared VISiceCoeff data structure'//trim(pid_msg)
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF

    destroy_status = crtm_vissnowcoeff_destroy( process_id = process_id )
    IF ( destroy_status /= success ) THEN
      err_stat = destroy_status
      msg = 'Error deallocating shared VISsnowCoeff data structure'//trim(pid_msg)
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF

    destroy_status = crtm_viswatercoeff_destroy( process_id = process_id )
    IF ( destroy_status /= success ) THEN
      err_stat = destroy_status
      msg = 'Error deallocating shared VISwaterCoeff data structure'//trim(pid_msg)
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF

    destroy_status = crtm_vislandcoeff_destroy( process_id = process_id )
    IF ( destroy_status /= success ) THEN
      err_stat = destroy_status
      msg = 'Error deallocating shared VISlandCoeff data structure'//trim(pid_msg)
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF

    destroy_status = crtm_iricecoeff_destroy( process_id = process_id )
    IF ( destroy_status /= success ) THEN
      err_stat = destroy_status
      msg = 'Error deallocating shared IRiceCoeff data structure'//trim(pid_msg)
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF

    destroy_status = crtm_irsnowcoeff_destroy( process_id = process_id )
    IF ( destroy_status /= success ) THEN
      err_stat = destroy_status
      msg = 'Error deallocating shared IRsnowCoeff data structure'//trim(pid_msg)
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF

    destroy_status = crtm_irwatercoeff_destroy( process_id = process_id )
    IF ( destroy_status /= success ) THEN
      err_stat = destroy_status
      msg = 'Error deallocating shared IRwaterCoeff data structure'//trim(pid_msg)
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF

    destroy_status = crtm_irlandcoeff_destroy( process_id = process_id )
    IF ( destroy_status /= success ) THEN
      err_stat = destroy_status
      msg = 'Error deallocating shared IRlandCoeff data structure'//trim(pid_msg)
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF

    destroy_status = crtm_aerosolcoeff_destroy( process_id = process_id )
    IF ( destroy_status /= success ) THEN
      err_stat = destroy_status
      msg = 'Error deallocating shared AerosolCoeff data structure'//trim(pid_msg)
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF

    destroy_status = crtm_cloudcoeff_destroy( process_id = process_id )
    IF ( destroy_status /= success ) THEN
      err_stat = destroy_status
      msg = 'Error deallocating shared CloudCoeff data structure'//trim(pid_msg)
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF

    destroy_status = crtm_destroy_taucoeff( process_id = process_id )
    IF ( destroy_status /= success ) THEN
      err_stat = destroy_status
      msg = 'Error deallocating shared TauCoeff data structure(s)'//trim(pid_msg)
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF

    destroy_status = crtm_spccoeff_destroy( process_id = process_id )
    IF ( destroy_status /= success ) THEN
      err_stat = destroy_status
      msg = 'Error deallocating shared SpcCoeff data structure(s)'//trim(pid_msg)
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF

    destroy_status = crtm_mwwatercoeff_destroy( process_id = process_id )
    IF ( destroy_status /= success ) THEN
      err_stat = destroy_status
      msg = 'Error deallocating shared MWwaterCoeff data structure'//trim(pid_msg)
      CALL display_message( routine_name,trim(msg)//trim(pid_msg),err_stat )
    END IF

  END FUNCTION crtm_destroy


!--------------------------------------------------------------------------------
!:sdoc+:
!
! NAME:
!       CRTM_IsInitialized
!
! PURPOSE:
!       Logical function to test if the CRTM has been correctly initialized.
!
! CALLING SEQUENCE:
!       status = CRTM_IsInitialized( ChannelInfo )
!
! INPUTS:
!       ChannelInfo:  ChannelInfo structure array.
!                     UNITS:      N/A
!                     TYPE:       CRTM_ChannelInfo_type
!                     DIMENSION:  Rank-1
!                     ATTRIBUTES: INTENT(IN)
!
! FUNCTION RESULT:
!       Status:       The return value is a logical result indicating if the
!                     CRTM has been correctly initialised.
!                     If == .TRUE., all the ChannelInfo entries are valid.
!                        == .FALSE., any of the ChannelInfo entries are invalid.
!                     UNITS:      N/A
!                     TYPE:       LOGICAL
!                     DIMENSION:  Scalar
!
!:sdoc-:
!--------------------------------------------------------------------------------

  FUNCTION crtm_isinitialized( ChannelInfo ) RESULT( Status )
    TYPE(crtm_channelinfo_type), INTENT(IN) :: channelinfo(:)
    LOGICAL :: status
    status = all(crtm_channelinfo_associated(channelinfo))
  END FUNCTION crtm_isinitialized


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

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

END MODULE crtm_lifecycle