scale_ocean_phy_albedo_nakajima00.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "scalelib.h"
module scale_ocean_phy_albedo_nakajima00
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  use scale_prof
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: ocean_phy_albedo_nakajima00_setup
  public :: ocean_phy_albedo_nakajima00
 
  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
 
  subroutine ocean_phy_albedo_nakajima00_setup
 
    log_newline
    log_info("OCEAN_PHY_ALBEDO_nakajima00_setup",*) 'Setup'
 
    return
  end subroutine ocean_phy_albedo_nakajima00_setup
 
  !-----------------------------------------------------------------------------
  subroutine ocean_phy_albedo_nakajima00( &
       OIA, OIS, OIE, &
       OJA, OJS, OJE, &
       cosSZA,        &
!      tau,           &
       SFC_albedo     )
    use scale_cpl_sfc_index
    implicit none
 
    integer,  intent(in)  :: oia, ois, oie
    integer,  intent(in)  :: oja, ojs, oje
    real(rp), intent(in)  :: cossza    (oia,oja)           ! cos(solar zenith angle)         (0-1)
!   real(RP), intent(in)  :: tau       (OIA,OJA)           ! cloud optical thickness
    real(rp), intent(out) :: sfc_albedo(oia,oja,n_rad_dir) ! sea surface albedo (short wave) (0-1)
 
    real(rp), parameter :: c_ocean_albedo(5,3) = reshape( &
         (/ -2.8108_rp   , -1.3651_rp,  2.9210e1_rp, -4.3907e1_rp,  1.8125e1_rp, &
             6.5626e-1_rp, -8.7253_rp, -2.7749e1_rp,  4.9486e1_rp, -1.8345e1_rp, &
            -6.5423e-1_rp,  9.9967_rp,  2.7769_rp  , -1.7620e1_rp,  7.0838_rp    /), &
         (/ 5,3 /) )
 
    real(rp) :: am1
!   real(RP) :: tr1
    real(rp) :: s
 
    integer :: i, j, n
    !---------------------------------------------------------------------------
    !$acc data copyin(cosSZA) copyout(SFC_albedo)
 
    !$acc kernels
    !$omp parallel do private(am1,s)
    do j = ojs, oje
    do i = ois, oie
       am1 = max( min( cossza(i,j), 0.961_rp ), 0.0349_rp )
!      tr1 = max( min( cosSZA(i,j) / ( 4.0_RP * tau(i,j) ), 1.0_RP ), 0.05_RP )
 
       s = 0.0_rp
!OCL UNROLL
       !$acc loop seq
       do n = 1, 5
          ! power of am1 is different form the paper, in which the number is wrong.
!         s = s + c_ocean_albedo(n,1) * tr1**(n-1) &
!               + c_ocean_albedo(n,2) * tr1**(n-1) * am1 &
!               + c_ocean_albedo(n,3) * tr1**(n-1) * am1**2
          s = s + c_ocean_albedo(n,1) &
                + c_ocean_albedo(n,2) * am1 &
                + c_ocean_albedo(n,3) * am1**2
       enddo
 
       sfc_albedo(i,j,i_r_direct ) = exp(s)
       sfc_albedo(i,j,i_r_diffuse) = sfc_albedo(i,j,i_r_direct)
    enddo
    enddo
    !$acc end kernels 
 
    !$acc end data
    return
  end subroutine ocean_phy_albedo_nakajima00
 
end module scale_ocean_phy_albedo_nakajima00