Foam_Utility_Module.f90

Go to the documentation of this file.

!
! Helper module containing the foam-related utility routines for the
! CRTM implementation of FASTEM4 and FASTEM5
!
!
! CREATION HISTORY:
!       Written by:     Original FASTEM1-5 authors
!
!       Refactored by:  Paul van Delst, November 2011
!                       paul.vandelst@noaa.gov
!

MODULE foam_utility_module

  ! -----------------
  ! Environment setup
  ! -----------------
  ! Module use
  USE type_kinds     , ONLY: fp
  USE fitcoeff_define, ONLY: fitcoeff_1d_type
  ! Disable implicit typing
  IMPLICIT NONE


  ! ------------
  ! Visibilities
  ! ------------
  PRIVATE
  PUBLIC :: foam_coverage
  PUBLIC :: foam_coverage_tl
  PUBLIC :: foam_coverage_ad
  PUBLIC :: foam_reflectivity


  ! -----------------
  ! Module parameters
  ! -----------------
  CHARACTER(*), PARAMETER :: module_version_id = &
  '$Id: Foam_Utility_Module.f90 60152 2015-08-13 19:19:13Z paul.vandelst@noaa.gov $'

  ! Literal constants
  REAL(fp), PARAMETER :: zero = 0.0_fp
  REAL(fp), PARAMETER :: one  = 1.0_fp
  REAL(fp), PARAMETER :: two  = 2.0_fp
  

CONTAINS


  ! ===================================================================
  ! Foam coverage.
  !
  !   Monahan, E.C., and O'Muircheartaigh, I.G., (1986)
  !     Whitecaps and the passive remote sensing of the ocean surface,
  !     International Journal of Remote Sensing, 7, pp627-642.
  !
  ! The neutral stability condition is used here (i.e. the difference
  ! between the skin and air temperature is assumed to be zero) so
  ! that the form of the foam coverage equation is the same as in
  ! Tang (1974) and Liu et al. (1998)..
  !
  !   Liu, Q. et al. (1998) Monte Carlo simulations of the
  !     microwave emissivity of the sea surface.
  !     JGR, 103(C11), pp24983-24989
  !
  !   Tang, C. (1974) The effect of droplets in the air-sea
  !     transition zone on the sea brightness temperature.
  !     J. Phys. Oceanography, 4, pp579-593.
  !
  ! ===================================================================
  ! Forward model
  SUBROUTINE foam_coverage(FCCoeff, wind_speed, coverage)
    TYPE(fitcoeff_1d_type), INTENT(IN)  :: fccoeff
    REAL(fp)              , INTENT(IN)  :: wind_speed
    REAL(fp)              , INTENT(OUT) :: coverage
    IF ( wind_speed < zero ) THEN
      coverage = zero
      RETURN
    END IF
    coverage = fccoeff%C(1) * (wind_speed**fccoeff%C(2))
  END SUBROUTINE foam_coverage
  
  ! Tangent-linear model
  SUBROUTINE foam_coverage_tl(FCCoeff, wind_speed, wind_speed_TL, coverage_TL)
    TYPE(fitcoeff_1d_type), INTENT(IN)  :: fccoeff
    REAL(fp)              , INTENT(IN)  :: wind_speed
    REAL(fp)              , INTENT(IN)  :: wind_speed_tl
    REAL(fp)              , INTENT(OUT) :: coverage_tl
    IF ( wind_speed < zero ) THEN
      coverage_tl = zero
      RETURN
    END IF
    coverage_tl = fccoeff%C(1)*fccoeff%C(2) * (wind_speed**(fccoeff%C(2)-one)) * wind_speed_tl
  END SUBROUTINE foam_coverage_tl

  ! Adjoint model
  SUBROUTINE foam_coverage_ad(FCCoeff, wind_speed, coverage_AD, wind_speed_AD)
    TYPE(fitcoeff_1d_type), INTENT(IN)     :: fccoeff
    REAL(fp)              , INTENT(IN)     :: wind_speed     ! Input
    REAL(fp)              , INTENT(IN OUT) :: coverage_ad    ! Input
    REAL(fp)              , INTENT(IN OUT) :: wind_speed_ad  ! Output
    IF ( wind_speed < zero ) THEN
      coverage_ad = zero
      RETURN
    END IF
    wind_speed_ad = wind_speed_ad + &
                    fccoeff%C(1)*fccoeff%C(2) * (wind_speed**(fccoeff%C(2)-one)) * coverage_ad
    coverage_ad = zero
  END SUBROUTINE foam_coverage_ad


  ! =============================================================
  ! Foam reflectivity
  !
  ! See section d in
  !
  !   Kazumori, M. et al. (2008) Impact Study of AMSR-E Radiances
  !     in the NCEP Global Data Assimilation System,
  !     Monthly Weather Review, 136, pp541-559
  !
  ! Function dependence is on zenith angle only so no TL
  ! or AD routine.
  ! =============================================================
  SUBROUTINE foam_reflectivity( &
    FRCoeff     , &
    Zenith_Angle, &
    Frequency   , &
    Rv          , &
    Rh            )
    ! Arguments
    TYPE(fitcoeff_1d_type), INTENT(IN)  :: frcoeff
    REAL(fp)              , INTENT(IN)  :: zenith_angle
    REAL(fp)              , INTENT(IN)  :: frequency
    REAL(fp)              , INTENT(OUT) :: rv, rh
    ! Local variables
    REAL(fp) :: factor
    
    ! The vertical component is a fixed value
    rv = one - frcoeff%C(1)  ! Fixed nadir emissivity
    
    ! The horizontal component uses a regression equation
    ! to compute a factor modifying the nadir emissivity
    factor = one + zenith_angle*(frcoeff%C(2) + &
                     zenith_angle*(frcoeff%C(3) + &
                       zenith_angle*frcoeff%C(4)  )  )
    rh = one - factor*frcoeff%C(1)
    
    ! Frequency correction
    factor = frcoeff%C(5) * exp(frcoeff%C(6)*frequency)
    rv = rv * factor
    rh = rh * factor

  END SUBROUTINE foam_reflectivity  

END MODULE foam_utility_module