scale_ocean_sfc.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
module scale_ocean_sfc
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  use scale_grid_index
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: ocean_sfc_setup

  abstract interface
     subroutine ocnsfc( &
           SST_t,  &
           ZMFLX,  &
           XMFLX,  &
           YMFLX,  &
           SHFLX,  &
           LHFLX,  &
           WHFLX,  &
           U10,    &
           V10,    &
           T2,     &
           Q2,     &
           TMPA,   &
           PRSA,   &
           WA,     &
           UA,     &
           VA,     &
           RHOA,   &
           QVA,    &
           Z1,     &
           PBL,    &
           PRSS,   &
           LWD,    &
           SWD,    &
           TW,     &
           SST,    &
           ALB_LW, &
           ALB_SW, &
           Z0M,    &
           Z0H,    &
           Z0E,    &
           dt      )
       use scale_precision
       use scale_grid_index
       implicit none

       real(RP), intent(out) :: sst_t(ia,ja) ! tendency of sea surface temperature
       real(RP), intent(out) :: zmflx(ia,ja) ! z-momentum flux at the surface [kg/m2/s]
       real(RP), intent(out) :: xmflx(ia,ja) ! x-momentum flux at the surface [kg/m2/s]
       real(RP), intent(out) :: ymflx(ia,ja) ! y-momentum flux at the surface [kg/m2/s]
       real(RP), intent(out) :: shflx(ia,ja) ! sensible heat flux at the surface [W/m2]
       real(RP), intent(out) :: lhflx(ia,ja) ! latent heat flux at the surface [W/m2]
       real(RP), intent(out) :: whflx(ia,ja) ! water heat flux at the surface [W/m2]
       real(RP), intent(out) :: u10  (ia,ja) ! velocity u at 10m [m/s]
       real(RP), intent(out) :: v10  (ia,ja) ! velocity v at 10m [m/s]
       real(RP), intent(out) :: t2   (ia,ja) ! temperature at 2m [K]
       real(RP), intent(out) :: q2   (ia,ja) ! water vapor at 2m [kg/kg]

       real(RP), intent(in) :: tmpa(ia,ja) ! temperature at the lowest atmospheric layer [K]
       real(RP), intent(in) :: prsa(ia,ja) ! pressure at the lowest atmospheric layer [Pa]
       real(RP), intent(in) :: wa  (ia,ja) ! velocity w at the lowest atmospheric layer [m/s]
       real(RP), intent(in) :: ua  (ia,ja) ! velocity u at the lowest atmospheric layer [m/s]
       real(RP), intent(in) :: va  (ia,ja) ! velocity v at the lowest atmospheric layer [m/s]
       real(RP), intent(in) :: rhoa(ia,ja) ! density at the lowest atmospheric layer [kg/m3]
       real(RP), intent(in) :: qva (ia,ja) ! ratio of water vapor mass to total mass at the lowest atmospheric layer [kg/kg]
       real(RP), intent(in) :: z1  (ia,ja) ! cell center height at the lowest atmospheric layer [m]
       real(RP), intent(in) :: pbl (ia,ja) ! the top of atmospheric mixing layer [m]
       real(RP), intent(in) :: prss(ia,ja) ! pressure at the surface [Pa]
       real(RP), intent(in) :: lwd (ia,ja) ! downward long-wave radiation flux at the surface (upward positive) [W/m2]
       real(RP), intent(in) :: swd (ia,ja) ! downward short-wave radiation flux at the surface (upward positive) [W/m2]

       real(RP), intent(in) :: tw    (ia,ja) ! water temperature [K]
       real(RP), intent(in) :: sst   (ia,ja) ! sea surface temperature [K]
       real(RP), intent(in) :: alb_lw(ia,ja) ! surface albedo for LW [0-1]
       real(RP), intent(in) :: alb_sw(ia,ja) ! surface albedo for SW [0-1]
       real(RP), intent(in) :: z0m   (ia,ja) ! roughness length for momentum [m]
       real(RP), intent(in) :: z0h   (ia,ja) ! roughness length for heat [m]
       real(RP), intent(in) :: z0e   (ia,ja) ! roughness length for vapor [m]
       real(DP), intent(in) :: dt            ! delta time
     end subroutine ocnsfc

     subroutine alb( &
           SFC_albedo_t, &
           SFC_albedo,   &
           cosSZA,       &
           dt            )
       use scale_precision
       use scale_grid_index
       implicit none

       real(RP), intent(out) :: sfc_albedo_t(ia,ja,2) ! tendency of sea surface albedo [0-1]
       real(RP), intent(in)  :: sfc_albedo  (ia,ja,2) ! sea surface                    [0-1]
       real(RP), intent(in)  :: cossza      (ia,ja)   ! cos(solar zenith angle)        [0-1]
       real(DP), intent(in)  :: dt                    ! delta time
     end subroutine alb
  end interface
  procedure(ocnsfc), pointer :: ocean_sfc              => null()
  procedure(alb),    pointer :: ocean_sfc_simplealbedo => null()
  public :: ocean_sfc
  public :: ocean_sfc_simplealbedo

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains

  subroutine ocean_sfc_setup( OCEAN_TYPE )
    use scale_process, only: &
       prc_mpistop
    use scale_ocean_sfc_slab, only: &
       ocean_sfc_slab_setup, &
       ocean_sfc_slab,       &
       ocean_sfc_slab_simplealbedo
    implicit none

    character(len=*), intent(in) :: OCEAN_TYPE
    !---------------------------------------------------------------------------

    select case( ocean_type )
    case ( 'CONST', 'FILE' )
       call ocean_sfc_slab_setup( ocean_type )
       ocean_sfc              => ocean_sfc_slab
       ocean_sfc_simplealbedo => ocean_sfc_slab_simplealbedo
    case ( 'SLAB' )
       call ocean_sfc_slab_setup( ocean_type )
       ocean_sfc              => ocean_sfc_slab
       ocean_sfc_simplealbedo => ocean_sfc_slab_simplealbedo
    end select

    return
  end subroutine ocean_sfc_setup

end module scale_ocean_sfc