mod_ocean_phy_driver Module Reference

module OCEAN / Physics More...

Functions/Subroutines
subroutine, public	ocean_phy_driver_setup
	Setup. More...
subroutine, public	ocean_phy_driver_resume
	Resume. More...
subroutine, public	ocean_phy_driver (update_flag)
	Driver. More...

Detailed Description

module OCEAN / Physics

Description

ocean physics module

Author

Team SCALE

NAMELIST

• No namelist group

History Output

name	description	unit	variable
OCEAN_ALB_LW_t	tendency of OCEAN_ALB_LW	1	OCEAN_SFC_albedo_t
OCEAN_ALB_SW_t	tendency of OCEAN_ALB_SW	1	OCEAN_SFC_albedo_t
OCEAN_SFC_TEMP_t	tendency of OCEAN_SFC_TEMP	K	OCEAN_SFC_TEMP_t
OCEAN_SFC_Z0E_t	tendency of OCEAN_SFC_Z0E	m	OCEAN_SFC_Z0E_t
OCEAN_SFC_Z0H_t	tendency of OCEAN_SFC_Z0H	m	OCEAN_SFC_Z0H_t
OCEAN_SFC_Z0M_t	tendency of OCEAN_SFC_Z0M	m	OCEAN_SFC_Z0M_t
OCEAN_TEMP_t	tendency of OCEAN_TEMP	K	OCEAN_TEMP_t

Function/Subroutine Documentation

ocean_phy_driver_setup()

subroutine, public mod_ocean_phy_driver::ocean_phy_driver_setup

(

)

Setup.

Definition at line 53 of file mod_ocean_phy_driver.f90.

References scale_stdio::io_fid_log, scale_stdio::io_l, scale_ocean_phy::ocean_phy_setup(), scale_ocean_sfc::ocean_sfc_setup(), mod_ocean_admin::ocean_sw, and mod_ocean_admin::ocean_type.

Referenced by mod_ocean_driver::ocean_driver_setup().

    use scale_ocean_phy, only: &
       ocean_phy_setup
    use scale_ocean_sfc, only: &
       ocean_sfc_setup
    use mod_ocean_admin, only: &
       ocean_type, &
       ocean_sw
    implicit none
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[DRIVER] / Categ[OCEAN PHY] / Origin[SCALE-RM]'

    if ( ocean_sw ) then

       ! setup library component
       call ocean_phy_setup( ocean_type )
       call ocean_sfc_setup( ocean_type )

    else
       if( io_l ) write(io_fid_log,*) '*** this component is never called.'
    endif

    return

Here is the call graph for this function:

Here is the caller graph for this function:

ocean_phy_driver_resume()

subroutine, public mod_ocean_phy_driver::ocean_phy_driver_resume

(

)

Resume.

Definition at line 82 of file mod_ocean_phy_driver.f90.

References ocean_phy_driver(), mod_ocean_admin::ocean_sw, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), and mod_admin_restart::restart_run.

Referenced by mod_ocean_driver::ocean_driver_resume().

    use mod_admin_restart, only: &
       restart_run
    use mod_ocean_admin, only: &
       ocean_sw
    implicit none

    if ( ocean_sw ) then

       if ( .NOT. restart_run ) then ! tentative
          ! run once (only for the diagnostic value)
          call prof_rapstart('OCN_Physics', 1)
          call ocean_phy_driver( update_flag = .true. )
          call prof_rapend  ('OCN_Physics', 1)
       end if

    end if

    return

Here is the call graph for this function:

Here is the caller graph for this function:

ocean_phy_driver()

subroutine, public mod_ocean_phy_driver::ocean_phy_driver

(

logical, intent(in)

update_flag

)

Driver.

Definition at line 105 of file mod_ocean_phy_driver.f90.

References mod_ocean_vars::atmos_cossza, mod_ocean_vars::atmos_dens, mod_ocean_vars::atmos_pbl, mod_ocean_vars::atmos_pres, mod_ocean_vars::atmos_qv, mod_ocean_vars::atmos_sfc_pres, mod_ocean_vars::atmos_sflx_lw, mod_ocean_vars::atmos_sflx_prec, mod_ocean_vars::atmos_sflx_sw, mod_ocean_vars::atmos_temp, mod_ocean_vars::atmos_u, mod_ocean_vars::atmos_v, mod_ocean_vars::atmos_w, scale_grid_index::ie, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, scale_ocean_phy::ocean_phy, mod_ocean_vars::ocean_q2, scale_ocean_sfc::ocean_sfc, mod_ocean_vars::ocean_sfc_albedo, mod_ocean_vars::ocean_sfc_albedo_t, scale_ocean_sfc::ocean_sfc_simplealbedo(), mod_ocean_vars::ocean_sfc_temp, mod_ocean_vars::ocean_sfc_temp_t, mod_ocean_vars::ocean_sfc_z0e, mod_ocean_vars::ocean_sfc_z0e_t, mod_ocean_vars::ocean_sfc_z0h, mod_ocean_vars::ocean_sfc_z0h_t, mod_ocean_vars::ocean_sfc_z0m, mod_ocean_vars::ocean_sfc_z0m_t, mod_ocean_vars::ocean_sflx_evap, mod_ocean_vars::ocean_sflx_lh, mod_ocean_vars::ocean_sflx_mu, mod_ocean_vars::ocean_sflx_mv, mod_ocean_vars::ocean_sflx_mw, mod_ocean_vars::ocean_sflx_sh, mod_ocean_vars::ocean_sflx_wh, mod_ocean_vars::ocean_t2, mod_ocean_vars::ocean_temp, mod_ocean_vars::ocean_temp_t, mod_ocean_vars::ocean_u10, mod_ocean_vars::ocean_v10, scale_grid_real::real_z1, scale_roughness::roughness, scale_rm_statistics::statistics_checktotal, and scale_time::time_dtsec_ocean.

Referenced by mod_ocean_driver::ocean_driver(), and ocean_phy_driver_resume().

    use scale_atmos_hydrometeor, only: &
       hydrometeor_lhv => atmos_hydrometeor_lhv
    use scale_time, only: &
       dt => time_dtsec_ocean
    use scale_rm_statistics, only: &
       statistics_checktotal, &
       stat_total
    use scale_history, only: &
       hist_in
    use scale_grid_real, only: &
       real_z1
    use scale_roughness, only: &
       roughness
    use scale_ocean_phy, only: &
       ocean_phy
    use scale_ocean_sfc, only: &
       ocean_sfc,              &
       ocean_sfc_simplealbedo
    use mod_ocean_vars, only: &
       ocean_temp,         &
       ocean_sfc_temp,     &
       ocean_sfc_albedo,   &
       ocean_sfc_z0m,      &
       ocean_sfc_z0h,      &
       ocean_sfc_z0e,      &
       ocean_temp_t,       &
       ocean_sfc_temp_t,   &
       ocean_sfc_albedo_t, &
       ocean_sfc_z0m_t,    &
       ocean_sfc_z0h_t,    &
       ocean_sfc_z0e_t,    &
       ocean_sflx_mw,      &
       ocean_sflx_mu,      &
       ocean_sflx_mv,      &
       ocean_sflx_sh,      &
       ocean_sflx_lh,      &
       ocean_sflx_wh,      &
       ocean_sflx_evap,    &
       ocean_u10,          &
       ocean_v10,          &
       ocean_t2,           &
       ocean_q2,           &
       atmos_temp,         &
       atmos_pres,         &
       atmos_w,            &
       atmos_u,            &
       atmos_v,            &
       atmos_dens,         &
       atmos_qv,           &
       atmos_pbl,          &
       atmos_sfc_pres,     &
       atmos_sflx_lw,      &
       atmos_sflx_sw,      &
       atmos_cossza,       &
       atmos_sflx_prec
    implicit none

    logical, intent(in) :: update_flag

    real(RP) :: LHV(IA,JA) ! latent heat of vaporization [J/kg]
    real(RP) :: total      ! dummy

    integer  :: i, j
    !---------------------------------------------------------------------------

    if ( update_flag ) then

       call roughness( ocean_sfc_z0m_t(:,:), & ! [OUT]
                       ocean_sfc_z0h_t(:,:), & ! [OUT]
                       ocean_sfc_z0e_t(:,:), & ! [OUT]
                       ocean_sfc_z0m(:,:), & ! [IN]
                       ocean_sfc_z0h(:,:), & ! [IN]
                       ocean_sfc_z0e(:,:), & ! [IN]
                       atmos_u(:,:), & ! [IN]
                       atmos_v(:,:), & ! [IN]
                       real_z1(:,:), & ! [IN]
                       dt                    ) ! [IN]

       call ocean_sfc_simplealbedo( ocean_sfc_albedo_t(:,:,:), & ! [OUT]
                                    ocean_sfc_albedo(:,:,:), & ! [IN]
                                    atmos_cossza(:,:),   & ! [IN]
                                    dt                         ) ! [IN]

       call ocean_sfc( ocean_sfc_temp_t(:,:),      & ! [OUT]
                       ocean_sflx_mw(:,:),      & ! [OUT]
                       ocean_sflx_mu(:,:),      & ! [OUT]
                       ocean_sflx_mv(:,:),      & ! [OUT]
                       ocean_sflx_sh(:,:),      & ! [OUT]
                       ocean_sflx_lh(:,:),      & ! [OUT]
                       ocean_sflx_wh(:,:),      & ! [OUT]
                       ocean_u10(:,:),      & ! [OUT]
                       ocean_v10(:,:),      & ! [OUT]
                       ocean_t2(:,:),      & ! [OUT]
                       ocean_q2(:,:),      & ! [OUT]
                       atmos_temp(:,:),      & ! [IN]
                       atmos_pres(:,:),      & ! [IN]
                       atmos_w(:,:),      & ! [IN]
                       atmos_u(:,:),      & ! [IN]
                       atmos_v(:,:),      & ! [IN]
                       atmos_dens(:,:),      & ! [IN]
                       atmos_qv(:,:),      & ! [IN]
                       real_z1(:,:),      & ! [IN]
                       atmos_pbl(:,:),      & ! [IN]
                       atmos_sfc_pres(:,:),      & ! [IN]
                       atmos_sflx_lw(:,:),      & ! [IN]
                       atmos_sflx_sw(:,:),      & ! [IN]
                       ocean_temp(:,:),      & ! [IN]
                       ocean_sfc_temp(:,:),      & ! [IN]
                       ocean_sfc_albedo(:,:,i_lw), & ! [IN]
                       ocean_sfc_albedo(:,:,i_sw), & ! [IN]
                       ocean_sfc_z0m(:,:),      & ! [IN]
                       ocean_sfc_z0h(:,:),      & ! [IN]
                       ocean_sfc_z0e(:,:),      & ! [IN]
                       dt                          ) ! [IN]


       call hydrometeor_lhv( lhv(:,:), atmos_temp(:,:) )

!OCL XFILL
       do j = js, je
       do i = is, ie
          ocean_sflx_evap(i,j) = ocean_sflx_lh(i,j) / lhv(i,j)
       end do
       end do

       call ocean_phy( ocean_temp_t(:,:), & ! [OUT]
                       ocean_temp(:,:), & ! [IN]
                       ocean_sflx_wh(:,:), & ! [IN]
                       atmos_sflx_prec(:,:), & ! [IN]
                       ocean_sflx_evap(:,:), & ! [IN]
                       dt                    ) ! [IN]

       call hist_in( ocean_temp_t(:,:),      'OCEAN_TEMP_t',     'tendency of OCEAN_TEMP',     'K' )
       call hist_in( ocean_sfc_temp_t(:,:),      'OCEAN_SFC_TEMP_t', 'tendency of OCEAN_SFC_TEMP', 'K' )
       call hist_in( ocean_sfc_albedo_t(:,:,i_lw), 'OCEAN_ALB_LW_t',   'tendency of OCEAN_ALB_LW',   '1' )
       call hist_in( ocean_sfc_albedo_t(:,:,i_sw), 'OCEAN_ALB_SW_t',   'tendency of OCEAN_ALB_SW',   '1' )
       call hist_in( ocean_sfc_z0m_t(:,:),      'OCEAN_SFC_Z0M_t',  'tendency of OCEAN_SFC_Z0M',  'm' )
       call hist_in( ocean_sfc_z0h_t(:,:),      'OCEAN_SFC_Z0H_t',  'tendency of OCEAN_SFC_Z0H',  'm' )
       call hist_in( ocean_sfc_z0e_t(:,:),      'OCEAN_SFC_Z0E_t',  'tendency of OCEAN_SFC_Z0E',  'm' )

    end if

    if ( statistics_checktotal ) then
       call stat_total( total, ocean_temp_t(:,:),      'OCEAN_TEMP_t'     )
       call stat_total( total, ocean_sfc_temp_t(:,:),      'OCEAN_SFC_TEMP_t' )
       call stat_total( total, ocean_sfc_albedo_t(:,:,i_lw), 'OCEAN_ALB_LW_t'   )
       call stat_total( total, ocean_sfc_albedo_t(:,:,i_sw), 'OCEAN_ALB_SW_t'   )
       call stat_total( total, ocean_sfc_z0m_t(:,:),      'OCEAN_SFC_Z0M_t'  )
       call stat_total( total, ocean_sfc_z0h_t(:,:),      'OCEAN_SFC_Z0H_t'  )
       call stat_total( total, ocean_sfc_z0e_t(:,:),      'OCEAN_SFC_Z0E_t'  )
    end if

    return

Here is the call graph for this function:

Here is the caller graph for this function: