scale_ocean_phy_albedo_nakajima00 Module Reference

module ocean / physics / surface albedo / nakajima00 More...

Functions/Subroutines
subroutine, public	ocean_phy_albedo_nakajima00_setup
subroutine, public	ocean_phy_albedo_nakajima00 (OIA, OIS, OIE, OJA, OJS, OJE, cosSZA, SFC_albedo)

Detailed Description

module ocean / physics / surface albedo / nakajima00

Description

Ocean surface albedo nakajima00 model

Reference

Nakajima, T., M. Tsukamoto, Y. Tsushima, A. Numaguti, and T. Kimura, 2000: Modeling of the radiative process in an atmospheric general circulation model, Applied Optics, 39, 4869-4878, doi:10.1364/AO.39.004869

Author

Team SCALE

Function/Subroutine Documentation

ocean_phy_albedo_nakajima00_setup()

subroutine, public scale_ocean_phy_albedo_nakajima00::ocean_phy_albedo_nakajima00_setup

(

)

Definition at line 49 of file scale_ocean_phy_albedo_nakajima00.F90.

Referenced by mod_ocean_driver::ocean_driver_setup().

    log_newline
    log_info("OCEAN_PHY_ALBEDO_nakajima00_setup",*) 'Setup'

    return

Here is the caller graph for this function:

ocean_phy_albedo_nakajima00()

subroutine, public scale_ocean_phy_albedo_nakajima00::ocean_phy_albedo_nakajima00	(	integer, intent(in)	OIA,
		integer, intent(in)	OIS,
		integer, intent(in)	OIE,
		integer, intent(in)	OJA,
		integer, intent(in)	OJS,
		integer, intent(in)	OJE,
		real(rp), dimension (oia,oja), intent(in)	cosSZA,
		real(rp), dimension(oia,oja,n_rad_dir), intent(out)	SFC_albedo
	)

Definition at line 63 of file scale_ocean_phy_albedo_nakajima00.F90.

References scale_cpl_sfc_index::i_r_diffuse, and scale_cpl_sfc_index::i_r_direct.

Referenced by mod_ocean_driver::ocean_driver_calc_tendency().

    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

Here is the caller graph for this function: