scale_ocean_phy_roughness_miller92 Module Reference

module ocean / physics / surface roughness length / miller92 More...

Functions/Subroutines
subroutine, public	ocean_phy_roughness_miller92_setup
subroutine, public	ocean_phy_roughness_miller92 (OIA, OIS, OIE, OJA, OJS, OJE, Uabs, Z0M, Z0H, Z0E)

Detailed Description

module ocean / physics / surface roughness length / miller92

Description

surface roughness length Miller92 scheme

Author

Team SCALE

NAMELIST

• PARAM_OCEAN_PHY_ROUGHNESS_MILLER92

  name	type	default value	comment
  OCEAN_PHY_ROUGHNESS_MILLER92_CM0	real(RP)	1.0E-3_RP	bulk coef. for U*
  OCEAN_PHY_ROUGHNESS_MILLER92_Z0MI	real(RP)	0.0E-0_RP	base roughness length for momentum
  OCEAN_PHY_ROUGHNESS_MILLER92_Z0MR	real(RP)	1.8E-2_RP	rough factor for momentum
  OCEAN_PHY_ROUGHNESS_MILLER92_Z0MS	real(RP)	1.1E-1_RP	smooth factor for momentum
  OCEAN_PHY_ROUGHNESS_MILLER92_Z0HI	real(RP)	1.4E-5_RP	base roughness length for heat
  OCEAN_PHY_ROUGHNESS_MILLER92_Z0HR	real(RP)	0.0E-0_RP	rough factor for heat
  OCEAN_PHY_ROUGHNESS_MILLER92_Z0HS	real(RP)	4.0E-1_RP	smooth factor for heat
  OCEAN_PHY_ROUGHNESS_MILLER92_Z0EI	real(RP)	1.3E-4_RP	base roughness length for moisture
  OCEAN_PHY_ROUGHNESS_MILLER92_Z0ER	real(RP)	0.0E-0_RP	rough factor for moisture
  OCEAN_PHY_ROUGHNESS_MILLER92_Z0ES	real(RP)	6.2E-1_RP	smooth factor for moisture

  
  

History Output

No history output

Function/Subroutine Documentation

ocean_phy_roughness_miller92_setup()

subroutine, public scale_ocean_phy_roughness_miller92::ocean_phy_roughness_miller92_setup

(

)

Definition at line 57 of file scale_ocean_phy_roughness_miller92.F90.

References scale_io::io_fid_conf, scale_ocean_phy_roughness::ocean_phy_roughness_setup(), and scale_prc::prc_abort().

Referenced by mod_ocean_driver::ocean_driver_setup().

    use scale_prc, only: &
       prc_abort
    use scale_ocean_phy_roughness, only: &
       ocean_phy_roughness_setup
    implicit none

    namelist / param_ocean_phy_roughness_miller92 / &
       ocean_phy_roughness_miller92_cm0,  &
       ocean_phy_roughness_miller92_z0mi, &
       ocean_phy_roughness_miller92_z0mr, &
       ocean_phy_roughness_miller92_z0ms, &
       ocean_phy_roughness_miller92_z0hi, &
       ocean_phy_roughness_miller92_z0hr, &
       ocean_phy_roughness_miller92_z0hs, &
       ocean_phy_roughness_miller92_z0ei, &
       ocean_phy_roughness_miller92_z0er, &
       ocean_phy_roughness_miller92_z0es

    integer :: ierr
    !---------------------------------------------------------------------------

    ! common setup
    call ocean_phy_roughness_setup

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_ocean_phy_roughness_miller92,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("OCEAN_PHY_ROUGHNESS_miller92_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("OCEAN_PHY_ROUGHNESS_miller92_setup",*) 'Not appropriate names in namelist PARAM_OCEAN_PHY_ROUGHNESS_MILLER92. Check!'
       call prc_abort
    endif
    log_nml(param_ocean_phy_roughness_miller92)

    return

Here is the call graph for this function:

Here is the caller graph for this function:

ocean_phy_roughness_miller92()

subroutine, public scale_ocean_phy_roughness_miller92::ocean_phy_roughness_miller92	(	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)	Uabs,
		real(rp), dimension (oia,oja), intent(out)	Z0M,
		real(rp), dimension (oia,oja), intent(out)	Z0H,
		real(rp), dimension (oia,oja), intent(out)	Z0E
	)

Definition at line 101 of file scale_ocean_phy_roughness_miller92.F90.

References scale_const::const_grav, scale_ocean_phy_roughness::ocean_phy_roughness_ustar_min, scale_ocean_phy_roughness::ocean_phy_roughness_visck, scale_ocean_phy_roughness::ocean_phy_roughness_z0e_min, scale_ocean_phy_roughness::ocean_phy_roughness_z0h_min, and scale_ocean_phy_roughness::ocean_phy_roughness_z0m_min.

Referenced by mod_ocean_driver::ocean_driver_calc_tendency().

    use scale_const, only: &
       grav => const_grav
    use scale_ocean_phy_roughness, only: &
       ocean_phy_roughness_visck, &
       ocean_phy_roughness_ustar_min, &
       ocean_phy_roughness_z0m_min, &
       ocean_phy_roughness_z0h_min, &
       ocean_phy_roughness_z0e_min
    implicit none

    integer,  intent(in)  :: oia, ois, oie
    integer,  intent(in)  :: oja, ojs, oje
    real(RP), intent(in)  :: uabs(oia,oja) ! velocity at the lowest atomspheric layer [m/s]
    real(RP), intent(out) :: z0m (oia,oja) ! roughness length for momentum [m]
    real(RP), intent(out) :: z0h (oia,oja) ! roughness length for heat [m]
    real(RP), intent(out) :: z0e (oia,oja) ! roughness length for vapor [m]

    real(RP) :: ustar

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

    !$omp parallel do private(Ustar)
    do j = ojs, oje
    do i = ois, oie
       ustar = max( sqrt( ocean_phy_roughness_miller92_cm0 ) * uabs(i,j), ocean_phy_roughness_ustar_min )

       z0m(i,j) = max( ocean_phy_roughness_miller92_z0mi &
                     + ocean_phy_roughness_miller92_z0mr / grav * ustar * ustar &
                     + ocean_phy_roughness_miller92_z0ms * ocean_phy_roughness_visck / ustar, &
                       ocean_phy_roughness_z0m_min )
       z0h(i,j) = max( ocean_phy_roughness_miller92_z0hi &
                     + ocean_phy_roughness_miller92_z0hr / grav * ustar * ustar &
                     + ocean_phy_roughness_miller92_z0hs * ocean_phy_roughness_visck / ustar, &
                       ocean_phy_roughness_z0h_min )
       z0e(i,j) = max( ocean_phy_roughness_miller92_z0ei &
                     + ocean_phy_roughness_miller92_z0er / grav * ustar * ustar &
                     + ocean_phy_roughness_miller92_z0es * ocean_phy_roughness_visck / ustar, &
                       ocean_phy_roughness_z0e_min )
    enddo
    enddo

    return

Here is the caller graph for this function: