scale_bulkflux Module Reference

module Surface bulk flux More...

Functions/Subroutines
subroutine, public	bulkflux_setup (dx)

Variables
procedure(bc), pointer, public	bulkflux => NULL()

Detailed Description

module Surface bulk flux

Description

calculation of surface bulk flux

Author

Team SCALE

History

NAMELIST

• PARAM_BULKFLUX

  name	type	default value	comment
  BULKFLUX_TYPE	character(len=H_SHORT)	'B91W01'
  BULKFLUX_ITR_SA_MAX	integer	5	maximum iteration number for successive approximation
  BULKFLUX_ITR_NR_MAX	integer	10	maximum iteration number for Newton-Raphson method
  BULKFLUX_ERR_MIN	real(RP)	1.0E-3_RP	minimum value of error
  BULKFLUX_WSCF	real(RP)		empirical scaling factor of Wstar (Beljaars 1994)
  BULKFLUX_UABS_MIN	real(RP)	1.0E-2_RP	minimum of Uabs [m/s]
  BULKFLUX_WSTAR_MIN	real(RP)	1.0E-4_RP	minimum of W* [m/s]

  
  

History Output

No history output

Function/Subroutine Documentation

bulkflux_setup()

subroutine, public scale_bulkflux::bulkflux_setup

(

real(rp), intent(in)

dx

)

Definition at line 118 of file scale_bulkflux.F90.

References bulkflux, scale_const::const_cpdry, scale_const::const_eps, scale_const::const_epstvap, scale_const::const_grav, scale_const::const_karman, scale_const::const_pre00, scale_const::const_rdry, scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_fid_nml, scale_stdio::io_l, scale_stdio::io_nml, and scale_process::prc_mpistop().

Referenced by mod_rm_driver::scalerm().

    use scale_process, only: &
       prc_mpistop
    implicit none
    real(RP), intent(in) :: dx

    namelist / param_bulkflux / &
       bulkflux_type,       &
       bulkflux_itr_sa_max, &
       bulkflux_itr_nr_max, &
       bulkflux_err_min,    &
       bulkflux_wscf,       &
       bulkflux_uabs_min,   &
       bulkflux_wstar_min

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

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[BULKFLUX] / Categ[COUPLER] / Origin[SCALElib]'

    ! WSCF = 1.2 for dx > 1 km (Beljaars 1994)
    ! lim_{dx->0} WSCF = 0  for LES (Fig. 6 Kitamura and Ito 2016 BLM)
    bulkflux_wscf = 1.2_rp * min(dx*1.e-3_rp, 1.0_rp)


    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_bulkflux,iostat=ierr)

    if( ierr < 0 ) then !--- missing
       if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       write(*,*) 'xxx Not appropriate names in namelist PARAM_BULKFLUX. Check!'
       call prc_mpistop
    endif
    if( io_nml ) write(io_fid_nml,nml=param_bulkflux)

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '*** Scheme for surface bulk flux : ', trim(bulkflux_type)
    select case(bulkflux_type)
    case('U95')
       if( io_l ) write(io_fid_log,*) '*** => Uno et al.(1995)'
       bulkflux => bulkflux_u95
    case('B91W01')
       if( io_l ) write(io_fid_log,*) '*** => Beljaars (1991) and Wilson (2001)'
       bulkflux => bulkflux_b91w01
    case default
       write(*,*) 'xxx Unsupported BULKFLUX_type. STOP'
       call prc_mpistop
    end select

    return

Here is the call graph for this function:

Here is the caller graph for this function:

Variable Documentation

bulkflux

procedure(bc), pointer, public scale_bulkflux::bulkflux => NULL()

Definition at line 77 of file scale_bulkflux.F90.

Referenced by scale_atmos_phy_sf_bulk::atmos_phy_sf_bulk(), bulkflux_setup(), scale_land_sfc_const::land_sfc_const(), scale_land_sfc_thick_slab::land_sfc_thick_slab(), scale_land_sfc_thin_slab::land_sfc_thin_slab(), scale_ocean_sfc_const::ocean_sfc_const(), scale_ocean_sfc_slab::ocean_sfc_slab(), and scale_urban_phy_slc::urban_phy_slc().

  procedure(bc), pointer :: BULKFLUX => null()