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
    !---------------------------------------------------------------------------

    !$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
       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

    return
  end subroutine ocean_phy_albedo_nakajima00

end module scale_ocean_phy_albedo_nakajima00