radiation_cloud_overlap.F90

 
 
module module_radiation_cloud_overlap
  use machine, only : kind_phys
  implicit none
  public :: cmp_dcorr_lgth
  
  interface cmp_dcorr_lgth
     module procedure cmp_dcorr_lgth_hogan
     module procedure cmp_dcorr_lgth_oreopoulos
  end interface
  
contains
 
  ! ######################################################################################
  ! Hogan et al. (2010)
  ! "Effect of improving representation of horizontal and vertical cloud structure on the 
  ! Earth's global radiation budget. Part I: Review and parametrization"
  ! https://rmets.onlinelibrary.wiley.com/doi/full/10.1002/qj.647
  ! ######################################################################################
  subroutine cmp_dcorr_lgth_hogan(nCol, lat, con_pi, dcorr_lgth)
    ! Inputs
    integer, intent(in) :: &
         nCol       ! Number of horizontal grid-points
    real(kind_phys), intent(in) :: &
         con_pi     ! Physical constant: Pi
    real(kind_phys), dimension(:), intent(in) :: &
         lat        ! Latitude  
    ! Outputs
    real(kind_phys), dimension(:),intent(out) :: &
         dcorr_lgth ! Decorrelation length  
    
    ! Local variables
    integer :: iCol
    
    ! Parameters
    real(kind_phys),parameter :: min_dcorr = 0.6 ! (see section 2.3)
    
    do icol =1,ncol
       dcorr_lgth(icol) = max(min_dcorr, 2.78-4.6*abs(lat(icol)/con_pi))        ! (eq. 13)
    enddo
    
  end subroutine cmp_dcorr_lgth_hogan
  ! ######################################################################################
  ! Oreopoulos et al. (2012)
  ! "Radiative impacts of cloud heterogeneity and overlap in an
  ! atmospheric General Circulation Model"
  ! 10.5194/acp-12-9097-2012
  ! ######################################################################################
  subroutine cmp_dcorr_lgth_oreopoulos(nCol, lat, juldat, yearlength, dcorr_lgth)
    ! Inputs
    integer, intent(in) :: &
         nCol,        & ! Number of horizontal grid-points
         yearlength     ! Number of days in year
 
    real(kind_phys), intent(in) :: &
         juldat         ! Julian date
    real(kind_phys), dimension(:), intent(in) :: &
         lat            ! Latitude  
    
    ! Outputs
    real(kind_phys), dimension(:),intent(out) :: &
         dcorr_lgth    ! Decorrelation length (km)
    
    ! Parameters for the Gaussian fits per Eqs. (10) and (11) (See Table 1)
    real(kind_phys), parameter ::  & !  
         am1 = 1.4315_kind_phys,   & ! 
         am2 = 2.1219_kind_phys,   & !
         am4 = -25.584_kind_phys,  & !
         amr = 7.0_kind_phys         !
    
    ! Local variables
    integer :: iCol     
    real(kind_phys) :: am3
    
    do icol = 1, ncol
       if (juldat .gt. 181._kind_phys) then
          am3 = -4._kind_phys * amr * (juldat - 272._kind_phys) / yearlength  ! (eq. 11a)
       else
          am3 =  4._kind_phys * amr * (juldat - 91._kind_phys)  / yearlength  ! (eq. 11b)
       endif
       dcorr_lgth(icol) = am1 + am2 * exp( -(lat(icol) - am3)**2 / am4**2)    ! (eq. 10)
    enddo
    
  end subroutine cmp_dcorr_lgth_oreopoulos
  
  subroutine get_alpha_exper(nCol, nLay, iovr, iovr_exprand, dzlay,    &
                             dcorr_lgth, cld_frac, alpha)
    
    ! Inputs
    integer, intent(in) :: &
         nCol,     & ! Number of horizontal grid points
         nLay        ! Number of vertical grid points
    integer, intent(in) :: &
         iovr,     &
         iovr_exprand
    real(kind_phys), dimension(:), intent(in) :: &
         dcorr_lgth  ! Decorrelation length (km)
    real(kind_phys), dimension(:,:), intent(in) :: &
         dzlay       !
    real(kind_phys), dimension(:,:), intent(in) ::  &
         cld_frac
    
    ! Outputs
    real(kind_phys), dimension(:,:) :: &
         alpha       ! Cloud overlap parameter
    
    ! Local variables
    integer :: iCol,iLay
    
    do icol = 1, ncol
       alpha(icol,1) = 0.0d0
       do ilay = 2, nlay
          alpha(icol,ilay) = exp( -(dzlay(icol,ilay)) / dcorr_lgth(icol))
       enddo
    enddo
   
    ! Revise alpha for exponential-random cloud overlap
    ! Decorrelate layers when a clear layer follows a cloudy layer to enforce
    ! random correlation between non-adjacent blocks of cloudy layers
    if (iovr == iovr_exprand) then
      do ilay = 2, nlay
        do icol = 1, ncol
          if (cld_frac(icol,ilay) == 0.0 .and. cld_frac(icol,ilay-1) > 0.0) then
            alpha(icol,ilay) = 0.0
          endif
        enddo
      enddo
    endif
 
    return
    
  end subroutine get_alpha_exper
end module module_radiation_cloud_overlap