scale_ocean_phy_roughness_moon07 Module Reference

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

Functions/Subroutines
subroutine, public	ocean_phy_roughness_moon07_setup
subroutine, public	ocean_phy_roughness_moon07 (OIA, OIS, OIE, OJA, OJS, OJE, Uabs, Z1, mask, Z0M, Z0H, Z0E)

Detailed Description

module ocean / physics / surface roughness length / moon07

Description

surface roughness length Moon07 scheme IL-JU MOON, ISAAC GINIS, TETSU HARA, AND BIJU THOMAS, 2007: A Physics-Based Parameterization of Air-Sea Momentum Flux at High Wind Speeds and Its Impact on Hurricane Intensity Predictions, Mon. Wea. Rev., 135, 2869-2878

Author

Team SCALE

NAMELIST

• PARAM_OCEAN_PHY_ROUGHNESS_MOON07

  name	type	default value	comment
  OCEAN_PHY_ROUGHNESS_MOON07_ITELIM	integer	10	maximum iteration number

  
  

History Output

No history output

Function/Subroutine Documentation

ocean_phy_roughness_moon07_setup()

subroutine, public scale_ocean_phy_roughness_moon07::ocean_phy_roughness_moon07_setup

Definition at line 51 of file scale_ocean_phy_roughness_moon07.F90.

    use scale_prc, only: &
       prc_abort
    use scale_ocean_phy_roughness, only: &
       ocean_phy_roughness_setup
    implicit none
 
    namelist / param_ocean_phy_roughness_moon07 / &
       ocean_phy_roughness_moon07_itelim
 
    integer :: ierr
    !---------------------------------------------------------------------------
 
    ! common setup
    call ocean_phy_roughness_setup
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_ocean_phy_roughness_moon07,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("OCEAN_PHY_ROUGHNESS_moon07_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("OCEAN_PHY_ROUGHNESS_moon07_setup",*) 'Not appropriate names in namelist PARAM_OCEAN_PHY_ROUGHNESS_MOON07. Check!'
       call prc_abort
    endif
    log_nml(param_ocean_phy_roughness_moon07)
 
    return

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().

Here is the call graph for this function:

Here is the caller graph for this function:

ocean_phy_roughness_moon07()

subroutine, public scale_ocean_phy_roughness_moon07::ocean_phy_roughness_moon07	(	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(in)	Z1,
		logical, dimension(oia,oja), intent(in)	mask,
		real(rp), dimension(oia,oja), intent(inout)	Z0M,
		real(rp), dimension(oia,oja), intent(out)	Z0H,
		real(rp), dimension(oia,oja), intent(out)	Z0E
	)

Definition at line 86 of file scale_ocean_phy_roughness_moon07.F90.

    use scale_const, only: &
       undef  => const_undef,  &
       grav   => const_grav,   &
       karman => const_karman
    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(in) :: Z1  (OIA,OJA) ! cell center height at the lowest atmospheric layer [m]
    logical,  intent(in) :: mask(OIA,OJA)
 
    real(RP), intent(inout) :: 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]
 
    ! works
    real(RP) :: Ustar
    real(RP) :: U10M
 
    integer  :: ite
    integer  :: i, j
    !---------------------------------------------------------------------------
    !$acc data copyin(Uabs,Z1,mask) copyout(Z0H,Z0E) copy(Z0M)
 
    !$omp parallel do default(none) OMP_SCHEDULE_ &
    !$omp shared(OJS,OJE,OIS,OIE, &
    !$omp        GRAV,UNDEF, &
    !$omp        OCEAN_PHY_ROUGHNESS_Z0M_min,OCEAN_PHY_ROUGHNESS_Z0H_min,OCEAN_PHY_ROUGHNESS_Z0E_min, &
    !$omp        OCEAN_PHY_ROUGHNESS_Ustar_min,OCEAN_PHY_ROUGHNESS_visck,OCEAN_PHY_ROUGHNESS_moon07_itelim, &
    !$omp        Z0M,Z0H,Z0E,Uabs,Z1,mask) &
    !$omp private(i,j,ite,U10M,Ustar)
    !$acc kernels
    do j = ojs, oje
    !$acc loop private(ite,U10M,Ustar)
    do i = ois, oie
       if ( mask(i,j) ) then
 
          z0m(i,j) = max( z0m(i,j), ocean_phy_roughness_z0m_min )
 
          !$acc loop seq
          do ite = 1, ocean_phy_roughness_moon07_itelim
             ustar = max( karman * uabs(i,j) / log( z1(i,j)/z0m(i,j) ), ocean_phy_roughness_ustar_min )
             u10m = ustar / karman * log( 10.0_rp/z0m(i,j) )
 
             if ( u10m <= 12.5_rp ) then
                z0m(i,j) = 0.0185_rp * ustar**2 / grav
             else
                z0m(i,j) = ( 0.085_rp * ( -  0.56_rp  * ustar**2   &
                                          + 20.255_rp * ustar      &
                                          +  2.458_rp            ) &
                           - 0.58_rp ) *  1.0e-3_rp
             end if
             z0m(i,j) = max( z0m(i,j), ocean_phy_roughness_z0m_min )
          enddo
 
          !  Fairall et al. TOGA V3.0
          !  Fairall et al. (2003) JCLI, vol. 16, 571-591. Eq. (28)
          z0h(i,j) = min( 5.5e-5_rp / ( z0m(i,j) * ustar / ocean_phy_roughness_visck )**0.6_rp, &
                          1.1e-4_rp )
          z0e(i,j) = z0h(i,j)
          z0h(i,j) = max( z0h(i,j), ocean_phy_roughness_z0h_min )
          z0e(i,j) = max( z0e(i,j), ocean_phy_roughness_z0e_min )
 
       else
          z0h(i,j) = undef
          z0e(i,j) = undef
       end if
    enddo
    enddo
    !$acc end kernels
 
    !$acc end data
    return

References scale_const::const_grav, scale_const::const_karman, scale_const::const_undef, 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().

Here is the caller graph for this function: