CRTM_Parameters.f90

Go to the documentation of this file.

!
! CRTM_Parameters
!
! Module of parameter definitions for the CRTM.
!
!
! CREATION HISTORY:
!       Written by:     Paul van Delst, 31-Jul-2000
!                       paul.vandelst@noaa.gov
!

MODULE crtm_parameters


  ! ---------------------
  ! Module use statements
  ! ---------------------
  USE type_kinds ,   ONLY: fp

  ! ------------------
  ! Default visibility
  ! ------------------
  ! Everything PRIVATE by default
  PRIVATE
  ! The MAX_N_CHANNELS methods
  PUBLIC :: crtm_set_max_nchannels    ! Subroutine
  PUBLIC :: crtm_reset_max_nchannels  ! Subroutine
  PUBLIC :: crtm_get_max_nchannels    ! Function
  PUBLIC :: crtm_isset_max_nchannels  ! Function


  !#----------------------------------------------------------------------------#
  !#                    -- ALGORITHM INDEPENDENT PARAMETERS --                  #
  !#----------------------------------------------------------------------------#

  ! ----------------------------------------------------------
  ! Number of channels (for ALL satellites - really the number
  ! of satellites x number of channels USED per satellite)
  !
  ! This is also the number of lines in the satellite 
  ! information (satinfo) file used in the NCEP GDAS.
  !
  ! The number of channels that can be used is determined,
  ! and SET, during the model initialisation.
  !
  ! In this module it is a protected variable in that it can
  ! only be set, reset, or retrieved via the *_N_CHANNELS
  ! methods.
  ! ----------------------------------------------------------
  INTEGER, PRIVATE, PARAMETER :: reset_value = -1
  INTEGER, PRIVATE, SAVE      :: max_n_channels = reset_value


  ! -----------------------------------------
  ! The maximum number of sensors that can be
  ! specified in a single initialisation call
  ! -----------------------------------------
  INTEGER, PUBLIC, PARAMETER :: max_n_sensors = 100

  ! -----------------------------------------------------
  ! The maximum number of atmospheric profiles and layers
  ! accepted. These values are arbitrary. Nothing magical
  ! -----------------------------------------------------
  INTEGER, PUBLIC, PARAMETER :: max_n_layers   = 200


  ! -----------------
  ! Literal constants
  ! -----------------
  REAL(fp), PUBLIC, PARAMETER :: zero          =  0.0_fp
  REAL(fp), PUBLIC, PARAMETER :: one           =  1.0_fp
  REAL(fp), PUBLIC, PARAMETER :: two           =  2.0_fp
  REAL(fp), PUBLIC, PARAMETER :: three         =  3.0_fp
  REAL(fp), PUBLIC, PARAMETER :: four          =  4.0_fp
  REAL(fp), PUBLIC, PARAMETER :: five          =  5.0_fp
  REAL(fp), PUBLIC, PARAMETER :: ten           = 10.0_fp
  REAL(fp), PUBLIC, PARAMETER :: point_25      =  0.25_fp
  REAL(fp), PUBLIC, PARAMETER :: point_5       =  0.5_fp
  REAL(fp), PUBLIC, PARAMETER :: point_75      =  0.75_fp
  REAL(fp), PUBLIC, PARAMETER :: onepointfive  =  1.5_fp
                                                                                                        

  ! --------------------
  ! PI-related constants
  ! --------------------
  REAL(fp), PUBLIC, PARAMETER :: pi = 3.141592653589793238462643383279_fp
  REAL(fp), PUBLIC, PARAMETER :: twopi = two * pi
  REAL(fp), PUBLIC, PARAMETER :: degrees_to_radians = pi / 180.0_fp
  REAL(fp), PUBLIC, PARAMETER :: radians_to_degrees = 180.0_fp / pi


  ! -------------------------
  ! Direction flags
  !   DOWN == From TOA to SFC
  !   UP   == From SFC to TOA
  ! -------------------------
  INTEGER, PUBLIC, PARAMETER :: down = 0
  INTEGER, PUBLIC, PARAMETER :: up   = 1


  ! ------------------------
  ! Invalid sensor ID values
  ! ------------------------
  INTEGER, PUBLIC, PARAMETER :: invalid_ncep_sensor_id   = -1
  INTEGER, PUBLIC, PARAMETER :: invalid_wmo_satellite_id = 1023
  INTEGER, PUBLIC, PARAMETER :: invalid_wmo_sensor_id    = 2047
  

  ! -------------------------
  ! Yes/No, Set/Unset flags
  ! -------------------------
  INTEGER, PUBLIC, PARAMETER :: no      = 0, yes = 1
  INTEGER, PUBLIC, PARAMETER :: not_set = 0, set = 1

  ! Default string length for SensorIDs
  INTEGER, PUBLIC, PARAMETER :: strlen = 20


  !#----------------------------------------------------------------------------#
  !#                       -- AtmAbsorption PARAMETERS --                       #
  !#----------------------------------------------------------------------------#

  ! Algorithm IDs
  INTEGER, PUBLIC, PARAMETER :: odas_algorithm     = 1
  INTEGER, PUBLIC, PARAMETER :: odps_algorithm     = 2
  INTEGER, PUBLIC, PARAMETER :: odssu_algorithm    = 3

  ! Maximum number of absorbers in the gas absorption models
  ! ---------------------------------------------------------------
  INTEGER, PUBLIC, PARAMETER :: max_n_absorbers = 7

  ! ----------------------------------------------
  ! The minimum absorber amount allowed based upon
  ! the smallest representable numbers.
  ! This value is equivalent to TINY(ONE)**0.25
  ! ----------------------------------------------
  REAL(fp), PUBLIC, PARAMETER :: minimum_absorber_amount = 1.0e-076_fp   ! TEN**(-RANGE(ONE)/4)

  ! ---------------------------------------
  ! Numerical limits for the gas absorption
  ! coefficient reconstruction
  ! ---------------------------------------
  ! Numerical limits based on precision
!  REAL(fp), PUBLIC, PARAMETER :: LIMIT_EXP = 36.0436_fp   ! ABS( LOG( TOLERANCE ) )
!  REAL(fp), PUBLIC, PARAMETER :: LIMIT_LOG = 4.5e+15_fp   ! EXP( LIMIT_EXP )
  ! Numerical limits based on experiment.
  REAL(fp), PUBLIC, PARAMETER :: limit_exp = 20.0_fp
  REAL(fp), PUBLIC, PARAMETER :: limit_log = 4.8e+08_fp   ! EXP( LIMIT_EXP )


  ! ---------------------------------------
  ! Top-Of-Atmosphere (TOA) pressure in hPa
  ! ---------------------------------------
  REAL(fp), PUBLIC, PARAMETER :: toa_pressure = 0.005_fp


  ! -------------------------------------------------------
  ! Reciprocal gravity (scaled by 100 for use with pressure
  ! in hPa) used in computing integrated absorber amounts
  ! -------------------------------------------------------
  REAL(fp), PUBLIC, PARAMETER :: reciprocal_gravity = one / 980.665_fp


  ! ----------------------------------------------------------
  ! Cosmic background temperature. Taken from
  ! Mather,J.C. et. al., 1999, "Calibrator Design for the COBE
  !    Far-Infrared Absolute Spectrophotometer (FIRAS)"
  !    Astrophysical Journal, vol 512, pp 511-520
  ! ----------------------------------------------------------
  REAL(fp), PUBLIC, PARAMETER :: tspace = 2.7253_fp



  !#----------------------------------------------------------------------------#
  !#                         -- Geometry PARAMETERS --                          #
  !#----------------------------------------------------------------------------#

  ! -----------------
  ! Distance defaults
  ! -----------------
  REAL(fp), PUBLIC, PARAMETER :: earth_radius     = 6370.0_fp  ! Mean earth radius 
  REAL(fp), PUBLIC, PARAMETER :: satellite_height = 800.0_fp


  ! ----------------
  ! Altitude extrema
  ! ----------------
  REAL(fp), PUBLIC, PARAMETER :: min_surface_altitude = -400.0_fp ! Dead sea
  REAL(fp), PUBLIC, PARAMETER :: max_surface_altitude = 8900.0_fp ! Chomolungma


  ! -----------------------------------
  ! Limits on sensor angles.
  ! - Azimuth angles that are multiples
  !   of 2pi are not accepted.
  ! -----------------------------------
  REAL(fp), PUBLIC, PARAMETER :: max_sensor_scan_angle    = 80.0_fp
  REAL(fp), PUBLIC, PARAMETER :: max_sensor_zenith_angle  = 80.0_fp
  REAL(fp), PUBLIC, PARAMETER :: max_sensor_azimuth_angle = 360.0_fp
  REAL(fp), PUBLIC, PARAMETER :: max_optran_secant  = 2.36620_fp   ! = 1/cosd(65.0)
  ! The maximum angle used to train ODAS transmittance 
  REAL(fp), PUBLIC, PARAMETER :: max_trans_zenith_angle = 63.6122_fp !corresponding to amass 2.25

  ! -------------------------------------------------
  ! Limits on source angles.
  ! - The maximum source zenith angle should
  !   be determined by the maximum angle secant
  !   used in generating the gas absorption model
  !   coefficients, i.e. a secant of 2.25 => 63.6deg.
  !   Users have requested the Value be 85deg which
  !   has a secant of ~11.47.
  ! - Azimuth angles that are multiples
  !   of 2pi are not accepted.
  ! -------------------------------------------------
  REAL(fp), PUBLIC, PARAMETER :: max_source_zenith_angle  = 85.0_fp
  REAL(fp), PUBLIC, PARAMETER :: max_secant_source_zenith = 11.473711738554476_fp
  REAL(fp), PUBLIC, PARAMETER :: max_source_azimuth_angle = 360.0_fp


  ! ----------------------------------------
  ! Default diffusivity angle and secant
  ! ACOS( 3/5 ) in degrees is (~53.13)
  ! Used to approximate the downwelling flux
  ! ----------------------------------------
  REAL(fp), PUBLIC, PARAMETER :: diffusivity_angle  = 53.130102354156_fp
  REAL(fp), PUBLIC, PARAMETER :: diffusivity_radian = 0.927295218002_fp
  REAL(fp), PUBLIC, PARAMETER :: secant_diffusivity = five / three


  ! -----------------------------------------------------------
  ! Maximum flux angle definitions. Determined by the maximum
  ! angle secant used in generating the gas absorption model
  ! coefficients, i.e. a secant of 2.25 => 63.6deg. If the user
  ! inputs a value larger than this for the Flux_Zenith_Angle,
  ! the diffusivity angles are used instead
  ! -----------------------------------------------------------
  REAL(fp), PUBLIC, PARAMETER :: max_flux_zenith_angle  = 63.612200038757_fp
  REAL(fp), PUBLIC, PARAMETER :: max_secant_flux_zenith = 2.25_fp


  !#----------------------------------------------------------------------------#
  !#            -- CloudScatter, RTSolution, AtmOptics PARAMETERS --            #
  !#----------------------------------------------------------------------------#
  INTEGER, PUBLIC, PARAMETER :: rt_ada = 56
  INTEGER, PUBLIC, PARAMETER :: rt_soi = 168
  
  INTEGER, PUBLIC, PARAMETER :: max_n_clouds   = 4 ! Max. number of clouds/profile. Needed for CSV...
  INTEGER, PUBLIC, PARAMETER :: max_n_aerosols = 20 ! Max. number of aerosols/profile. Needed for ASV
  
  REAL(fp), PUBLIC, PARAMETER :: aerosol_content_threshold = 0.000000001_fp
  REAL(fp), PUBLIC, PARAMETER :: water_content_threshold = 0.000001_fp
  REAL(fp), PUBLIC, PARAMETER :: optical_depth_threshold = 0.000001_fp

  REAL(fp), PUBLIC, PARAMETER :: bs_threshold  = 1.0e-10_fp  ! Was SCATTERING_ALBEDO_THRESHOLD
  REAL(fp), PUBLIC, PARAMETER :: scattering_albedo_threshold = bs_threshold  ! Eventually replace this with BS_THRESHOLD


  INTEGER, PUBLIC, PARAMETER :: max_n_legendre_terms = 16
  INTEGER, PUBLIC, PARAMETER :: max_n_phase_elements = 1
  INTEGER, PUBLIC, PARAMETER :: max_n_streams = 16 
  INTEGER, PUBLIC, PARAMETER :: max_n_angles = 16
  INTEGER, PUBLIC, PARAMETER :: max_n_stokes = 4
  INTEGER, PUBLIC, PARAMETER :: max_n_azimuth_fourier = 16   ! maximum number of Fourier components for azimuth angles
    
!### SOI uses HG phase function (modified by Tahara, Feb 2008)
  LOGICAL, PUBLIC, PARAMETER :: hgphase = .false.
! LOGICAL, PUBLIC, PARAMETER :: HGPHASE = .TRUE.
!### SOI uses HG phase function (modified by Tahara, Feb 2008)
                                                                                                        

  !#----------------------------------------------------------------------------#
  !#            -- Aircraft/profile pressure difference threshold --            #
  !#----------------------------------------------------------------------------#
  REAL(fp), PUBLIC, PARAMETER :: aircraft_pressure_threshold = 0.1_fp
  
  
CONTAINS


!--------------------------------------------------------------------------------
!
! NAME:
!       CRTM_Set_Max_nChannels
! 
! PURPOSE:
!       Subroutine to set the protected variable MAX_N_CHANNELS value in the
!       CRTM_Parameters module. This should *only* be done during the CRTM
!       initialisation.
!
! CALLING SEQUENCE:
!       CALL CRTM_Set_Max_nChannels( Value )
!
! INPUT ARGUMENTS:
!       Value:        The number of channels.
!                     UNITS:      N/A
!                     TYPE:       INTEGER
!                     DIMENSION:  Scalar
!                     ATTRIBUTES: INTENT(IN)
!
! SIDE EFFECTS:
!       This subroutines changes the value of the MAX_N_CHANNELS pseudo-parameter
!       in the CRTM_Parameters module.
!
! CREATION HISTORY:
!       Written by:     Paul van Delst, CIMSS/SSEC 16-Aug-2001
!                       paul.vandelst@ssec.wisc.edu
!
!--------------------------------------------------------------------------------

  SUBROUTINE crtm_set_max_nchannels( Value )
    INTEGER, INTENT(IN) :: value
    max_n_channels = Value
  END SUBROUTINE crtm_set_max_nchannels


!--------------------------------------------------------------------------------
!
! NAME:
!       CRTM_Reset_Max_nChannels
! 
! PURPOSE:
!       Subroutine to reset the protected variable MAX_N_CHANNELS value in the
!       CRTM_Parameters module to an invalid value. This should *only* be done
!       during the CRTM destruction.
!
! CALLING SEQUENCE:
!       CALL CRTM_Reset_Max_nChannels()
!
! SIDE EFFECTS:
!       This subroutines changes the value of the MAX_N_CHANNELS pseudo-parameter
!       in the CRTM_Parameters module.
!
! CREATION HISTORY:
!       Written by:     Paul van Delst, CIMSS/SSEC 16-Aug-2001
!                       paul.vandelst@ssec.wisc.edu
!
!--------------------------------------------------------------------------------

  SUBROUTINE crtm_reset_max_nchannels()
    max_n_channels = reset_value
  END SUBROUTINE crtm_reset_max_nchannels


!--------------------------------------------------------------------------------
!
! NAME:
!       CRTM_Get_Max_nChannels
! 
! PURPOSE:
!       Function to GET the protected variable MAX_N_CHANNELS value stored
!       in the CRTM_Parameters module.
!
! CALLING SEQUENCE:
!       Value = CRTM_Get_Max_nChannels()
!
! FUNCTION RESULT:
!       Value:  The number of channels.
!               UNITS:      N/A
!               TYPE:       INTEGER
!               DIMENSION:  Scalar
!
! CREATION HISTORY:
!       Written by:     Paul van Delst, CIMSS/SSEC 16-Aug-2001
!                       paul.vandelst@ssec.wisc.edu
!
!--------------------------------------------------------------------------------

  FUNCTION crtm_get_max_nchannels() RESULT(Value)
    INTEGER :: value
    Value = max_n_channels
  END FUNCTION crtm_get_max_nchannels


!--------------------------------------------------------------------------------
!
! NAME:
!       CRTM_IsSet_Max_nChannels
! 
! PURPOSE:
!       Function to determie if the MAX_N_CHANNELS protected variable in the
!       CRTM_Parameters module has been set to a valid value
!
! CALLING SEQUENCE:
!       Is_Set = CRTM_IsSet_Max_nChannels()
!
! FUNCTION RESULT:
!       Is_Set:  Logical flag for determining whether or not the
!                number of channels has been set.
!                If == .TRUE.  the MAX_N_CHANNELS protected variable is
!                              set to a valid value.
!                   == .FALSE. the MAX_N_CHANNELS protected variable
!                              value is invalid
!                UNITS:      N/A
!                TYPE:       LOGICAL
!                DIMENSION:  Scalar
!
! CREATION HISTORY:
!       Written by:     Paul van Delst, CIMSS/SSEC 16-Aug-2001
!                       paul.vandelst@ssec.wisc.edu
!
!--------------------------------------------------------------------------------

  FUNCTION crtm_isset_max_nchannels() RESULT(Is_Set)
    LOGICAL :: is_set
    is_set = (max_n_channels /= reset_value)
  END FUNCTION crtm_isset_max_nchannels

END MODULE crtm_parameters