SCALE-RM
Functions/Subroutines
scale_atmos_phy_sf_bulk Module Reference

module ATMOSPHERE / Physics Surface fluxes More...

Functions/Subroutines

subroutine, public atmos_phy_sf_bulk_setup (ATMOS_PHY_SF_TYPE)
 Setup. More...
 
subroutine, public atmos_phy_sf_bulk (ATM_TEMP, ATM_PRES, ATM_W, ATM_U, ATM_V, ATM_DENS, ATM_QTRC, ATM_Z1, dt, SFC_DENS, SFC_PRES, SFLX_LW_dn, SFLX_SW_dn, SFC_TEMP, SFC_albedo, SFC_Z0M, SFC_Z0H, SFC_Z0E, SFLX_MW, SFLX_MU, SFLX_MV, SFLX_SH, SFLX_LH, SFLX_QTRC, U10, V10, T2, Q2)
 Calculate surface flux. More...
 

Detailed Description

module ATMOSPHERE / Physics Surface fluxes

Description
Flux from/to bottom wall of atmosphere (surface) Bulk Method
Author
Team SCALE
History
  • 2011-12-03 (Y.Miyamoto) [new]
  • 2011-12-11 (H.Yashiro) [mod] integrate to SCALE-LES ver.3
  • 2012-03-23 (H.Yashiro) [mod] Explicit index parameter inclusion
  • 2012-04-10 (Y.Miyamoto) [mod] introduce coefficients for interpolation
  • 2012-09-11 (S.Nishizawa) [mod] bugfix based on the scale document
NAMELIST
  • PARAM_ATMOS_PHY_SF_BULK
    nametypedefault valuecomment
    ATMOS_PHY_SF_BETA real(RP) 1.0_RP evaporation efficiency [0-1]
    ATMOS_PHY_SF_U_MINM real(RP) 0.0_RP minimum limit of absolute velocity for momentum [m/s]
    ATMOS_PHY_SF_U_MINH real(RP) 0.0_RP minimum limit of absolute velocity for heat [m/s]
    ATMOS_PHY_SF_U_MINE real(RP) 0.0_RP minimum limit of absolute velocity for vapor [m/s]

History Output
No history output

Function/Subroutine Documentation

◆ atmos_phy_sf_bulk_setup()

subroutine, public scale_atmos_phy_sf_bulk::atmos_phy_sf_bulk_setup ( character(len=*), intent(in)  ATMOS_PHY_SF_TYPE)

Setup.

Definition at line 65 of file scale_atmos_phy_sf_bulk.F90.

References scale_atmos_phy_sf_bulkcoef::atmos_phy_sf_bulkcoef_setup(), scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_l, scale_stdio::io_lnml, and scale_process::prc_mpistop().

Referenced by scale_atmos_phy_sf::atmos_phy_sf_setup().

65  use scale_process, only: &
66  prc_mpistop
67  use scale_atmos_phy_sf_bulkcoef, only: &
68  sf_bulkcoef_setup => atmos_phy_sf_bulkcoef_setup
69  implicit none
70 
71  character(len=*), intent(in) :: atmos_phy_sf_type
72 
73  namelist / param_atmos_phy_sf_bulk / &
74  atmos_phy_sf_beta, &
75  atmos_phy_sf_u_minm, &
76  atmos_phy_sf_u_minh, &
77  atmos_phy_sf_u_mine
78 
79  integer :: ierr
80  !---------------------------------------------------------------------------
81 
82  if( io_l ) write(io_fid_log,*)
83  if( io_l ) write(io_fid_log,*) '++++++ Module[SURFACE FLUX] / Categ[ATMOS PHYSICS] / Origin[SCALElib]'
84  if( io_l ) write(io_fid_log,*) '+++ Bulk scheme'
85 
86  if ( atmos_phy_sf_type /= 'BULK' ) then
87  write(*,*) 'xxx ATMOS_PHY_SF_TYPE is not BULK. Check!'
88  call prc_mpistop
89  endif
90 
91  !--- read namelist
92  rewind(io_fid_conf)
93  read(io_fid_conf,nml=param_atmos_phy_sf_bulk,iostat=ierr)
94  if( ierr < 0 ) then !--- missing
95  if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
96  elseif( ierr > 0 ) then !--- fatal error
97  write(*,*) 'xxx Not appropriate names in namelist PARAM_ATMOS_PHY_SF_BULK. Check!'
98  call prc_mpistop
99  endif
100  if( io_lnml ) write(io_fid_log,nml=param_atmos_phy_sf_bulk)
101 
102  call sf_bulkcoef_setup
103 
104  return
scale_process::prc_mpistop
subroutine, public prc_mpistop
Abort MPI.
Definition: scale_process.F90:234
scale_atmos_phy_sf_bulkcoef::atmos_phy_sf_bulkcoef_setup
subroutine, public atmos_phy_sf_bulkcoef_setup
Definition: scale_atmos_phy_sf_bulkcoef.F90:98
scale_process
module PROCESS
Definition: scale_process.F90:14
scale_atmos_phy_sf_bulkcoef
module ATMOSPHERE / Physics Surface bulk coefficient
Definition: scale_atmos_phy_sf_bulkcoef.F90:15
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_phy_sf_bulk()

subroutine, public scale_atmos_phy_sf_bulk::atmos_phy_sf_bulk ( real(rp), dimension (ia,ja), intent(in)  ATM_TEMP,
real(rp), dimension (ia,ja), intent(in)  ATM_PRES,
real(rp), dimension (ia,ja), intent(in)  ATM_W,
real(rp), dimension (ia,ja), intent(in)  ATM_U,
real(rp), dimension (ia,ja), intent(in)  ATM_V,
real(rp), dimension (ia,ja), intent(in)  ATM_DENS,
real(rp), dimension (ia,ja,qa), intent(in)  ATM_QTRC,
real(rp), dimension (ia,ja), intent(in)  ATM_Z1,
real(dp), intent(in)  dt,
real(rp), dimension (ia,ja), intent(in)  SFC_DENS,
real(rp), dimension (ia,ja), intent(in)  SFC_PRES,
real(rp), dimension(ia,ja), intent(in)  SFLX_LW_dn,
real(rp), dimension(ia,ja), intent(in)  SFLX_SW_dn,
real(rp), dimension (ia,ja), intent(in)  SFC_TEMP,
real(rp), dimension(ia,ja,2), intent(in)  SFC_albedo,
real(rp), dimension (ia,ja), intent(inout)  SFC_Z0M,
real(rp), dimension (ia,ja), intent(inout)  SFC_Z0H,
real(rp), dimension (ia,ja), intent(inout)  SFC_Z0E,
real(rp), dimension (ia,ja), intent(out)  SFLX_MW,
real(rp), dimension (ia,ja), intent(out)  SFLX_MU,
real(rp), dimension (ia,ja), intent(out)  SFLX_MV,
real(rp), dimension (ia,ja), intent(out)  SFLX_SH,
real(rp), dimension (ia,ja), intent(out)  SFLX_LH,
real(rp), dimension (ia,ja,qa), intent(out)  SFLX_QTRC,
real(rp), dimension (ia,ja), intent(out)  U10,
real(rp), dimension (ia,ja), intent(out)  V10,
real(rp), dimension (ia,ja), intent(out)  T2,
real(rp), dimension (ia,ja), intent(out)  Q2 
)

Calculate surface flux.

Definition at line 121 of file scale_atmos_phy_sf_bulk.F90.

References scale_atmos_phy_sf_bulkcoef::atmos_phy_sf_bulkcoef, scale_bulkflux::bulkflux, scale_const::const_cpdry, scale_const::const_rdry, scale_tracer::i_qv, scale_grid_index::ie, scale_stdio::io_fid_log, scale_stdio::io_l, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, dc_log::log(), and scale_roughness::roughness.

Referenced by scale_atmos_phy_sf::atmos_phy_sf_setup().

121  use scale_grid_index
122  use scale_tracer
123  use scale_const, only: &
124  cpdry => const_cpdry, &
125  rdry => const_rdry
126  use scale_atmos_phy_sf_bulkcoef, only: &
127  sf_bulkcoef => atmos_phy_sf_bulkcoef
128  use scale_atmos_saturation, only: &
129  saturation_pres2qsat_all => atmos_saturation_pres2qsat_all
130  use scale_atmos_thermodyn, only: &
131  atmos_thermodyn_templhv
132  use scale_roughness, only: &
133  roughness
134  use scale_bulkflux, only: &
135  bulkflux
136  implicit none
137 
138  real(RP), intent(in) :: atm_temp (ia,ja) ! temperature at the lowermost layer (cell center) [K]
139  real(RP), intent(in) :: atm_pres (ia,ja) ! pressure at the lowermost layer (cell center) [Pa]
140  real(RP), intent(in) :: atm_w (ia,ja) ! velocity w at the lowermost layer (cell center) [m/s]
141  real(RP), intent(in) :: atm_u (ia,ja) ! velocity u at the lowermost layer (cell center) [m/s]
142  real(RP), intent(in) :: atm_v (ia,ja) ! velocity v at the lowermost layer (cell center) [m/s]
143  real(RP), intent(in) :: atm_dens (ia,ja) ! density at the lowermost layer (cell center) [kg/m3]
144  real(RP), intent(in) :: atm_qtrc (ia,ja,qa) ! tracer at the lowermost layer (cell center) [kg/kg]
145  real(RP), intent(in) :: atm_z1 (ia,ja) ! height of the lowermost grid from surface (cell center) [m]
146  real(DP), intent(in) :: dt ! delta time
147  real(RP), intent(in) :: sfc_dens (ia,ja) ! density at the surface atmosphere [kg/m3]
148  real(RP), intent(in) :: sfc_pres (ia,ja) ! pressure at the surface atmosphere [Pa]
149  real(RP), intent(in) :: sflx_lw_dn(ia,ja) ! downward longwave radiation flux at the surface [J/m2/s]
150  real(RP), intent(in) :: sflx_sw_dn(ia,ja) ! downward shortwave radiation flux at the surface [J/m2/s]
151  real(RP), intent(in) :: sfc_temp (ia,ja) ! temperature at the surface skin [K]
152  real(RP), intent(in) :: sfc_albedo(ia,ja,2) ! surface albedo (LW/SW) [0-1]
153  real(RP), intent(inout) :: sfc_z0m (ia,ja) ! surface roughness length (momentum) [m]
154  real(RP), intent(inout) :: sfc_z0h (ia,ja) ! surface roughness length (heat) [m]
155  real(RP), intent(inout) :: sfc_z0e (ia,ja) ! surface roughness length (vapor) [m]
156  real(RP), intent(out) :: sflx_mw (ia,ja) ! surface flux for z-momentum (area center) [m/s*kg/m2/s]
157  real(RP), intent(out) :: sflx_mu (ia,ja) ! surface flux for x-momentum (area center) [m/s*kg/m2/s]
158  real(RP), intent(out) :: sflx_mv (ia,ja) ! surface flux for y-momentum (area center) [m/s*kg/m2/s]
159  real(RP), intent(out) :: sflx_sh (ia,ja) ! surface flux for sensible heat (area center) [J/m2/s]
160  real(RP), intent(out) :: sflx_lh (ia,ja) ! surface flux for latent heat (area center) [J/m2/s]
161  real(RP), intent(out) :: sflx_qtrc (ia,ja,qa) ! surface flux for tracer mass (area center) [kg/m2/s]
162  real(RP), intent(out) :: u10 (ia,ja) ! velocity u at 10m height
163  real(RP), intent(out) :: v10 (ia,ja) ! velocity v at 10m height
164  real(RP), intent(out) :: t2 (ia,ja) ! temperature t at 2m height
165  real(RP), intent(out) :: q2 (ia,ja) ! water vapor q at 2m height
166 
167  real(RP) :: sfc_z0m_t(ia,ja)
168  real(RP) :: sfc_z0h_t(ia,ja)
169  real(RP) :: sfc_z0e_t(ia,ja)
170 
171  real(RP) :: sfc_qsat(ia,ja) ! saturatad water vapor mixing ratio [kg/kg]
172 
173  real(RP) :: pbl(ia,ja)
174  real(RP) :: lhv (ia,ja)
175  real(RP) :: ustar
176  real(RP) :: tstar
177  real(RP) :: qstar
178  real(RP) :: uabs
179 
180  integer :: i, j
181  !---------------------------------------------------------------------------
182 
183  if( io_l ) write(io_fid_log,*) '*** Physics step: Surface flux(bulk)'
184 
185  call roughness( sfc_z0m_t(:,:), & ! [OUT]
186  sfc_z0h_t(:,:), & ! [OUT]
187  sfc_z0e_t(:,:), & ! [OUT]
188  sfc_z0m(:,:), & ! [IN]
189  sfc_z0h(:,:), & ! [IN]
190  sfc_z0e(:,:), & ! [IN]
191  atm_u(:,:), & ! [IN]
192  atm_v(:,:), & ! [IN]
193  atm_z1(:,:), & ! [IN]
194  dt ) ! [IN]
195 
196  sfc_z0m(:,:) = sfc_z0m(:,:) + sfc_z0m_t(:,:) * dt
197  sfc_z0h(:,:) = sfc_z0h(:,:) + sfc_z0h_t(:,:) * dt
198  sfc_z0e(:,:) = sfc_z0e(:,:) + sfc_z0e_t(:,:) * dt
199 
200  call saturation_pres2qsat_all( sfc_qsat(:,:), & ! [OUT]
201  sfc_temp(:,:), & ! [IN]
202  sfc_pres(:,:) ) ! [IN]
203 
204  call atmos_thermodyn_templhv( lhv, atm_temp )
205 
206  sflx_qtrc(:,:,:) = 0.0_rp
207  pbl(:,:) = 100.0_rp ! tentative
208  do j = js, je
209  do i = is, ie
210 
211  call bulkflux( &
212  ustar, & ! [OUT]
213  tstar, & ! [OUT]
214  qstar, & ! [OUT]
215  uabs, & ! [OUT]
216  atm_temp(i,j), & ! [IN]
217  sfc_temp(i,j), & ! [IN]
218  atm_pres(i,j), & ! [IN]
219  sfc_pres(i,j), & ! [IN]
220  atm_qtrc(i,j,i_qv), & ! [IN]
221  sfc_qsat(i,j), & ! [IN]
222  atm_u(i,j), & ! [IN]
223  atm_v(i,j), & ! [IN]
224  atm_z1(i,j), & ! [IN]
225  pbl(i,j), & ! [IN]
226  sfc_z0m(i,j), & ! [IN]
227  sfc_z0h(i,j), & ! [IN]
228  sfc_z0e(i,j) ) ! [IN]
229 
230  !-----< momentum >-----
231  sflx_mw(i,j) = -atm_dens(i,j) * ustar**2 / uabs * atm_w(i,j)
232  sflx_mu(i,j) = -atm_dens(i,j) * ustar**2 / uabs * atm_u(i,j)
233  sflx_mv(i,j) = -atm_dens(i,j) * ustar**2 / uabs * atm_v(i,j)
234 
235  !-----< heat flux >-----
236  sflx_sh(i,j) = -cpdry * atm_dens(i,j) * ustar * tstar
237  sflx_lh(i,j) = -lhv(i,j) * atm_dens(i,j) * ustar * qstar * atmos_phy_sf_beta
238 
239  !-----< mass flux >-----
240  sflx_qtrc(i,j,i_qv) = sflx_lh(i,j) / lhv(i,j)
241  enddo
242  enddo
243 
244  !-----< U10, T2, q2 >-----
245 
246  do j = js, je
247  do i = is, ie
248  u10(i,j) = atm_u(i,j) * log( 10.0_rp / sfc_z0m(i,j) ) / log( atm_z1(i,j) / sfc_z0m(i,j) )
249  v10(i,j) = atm_v(i,j) * log( 10.0_rp / sfc_z0m(i,j) ) / log( atm_z1(i,j) / sfc_z0m(i,j) )
250  t2(i,j) = sfc_temp(i,j) + ( atm_temp(i,j) - sfc_temp(i,j) ) &
251  * ( log( 2.0_rp / sfc_z0m(i,j) ) * log( 2.0_rp / sfc_z0h(i,j) ) ) &
252  / ( log( atm_z1(i,j) / sfc_z0m(i,j) ) * log( atm_z1(i,j) / sfc_z0h(i,j) ) )
253  q2(i,j) = sfc_qsat(i,j) + ( atm_qtrc(i,j,i_qv) - sfc_qsat(i,j) ) &
254  * ( log( 2.0_rp / sfc_z0m(i,j) ) * log( 2.0_rp / sfc_z0e(i,j) ) ) &
255  / ( log( atm_z1(i,j) / sfc_z0m(i,j) ) * log( atm_z1(i,j) / sfc_z0e(i,j) ) )
256  enddo
257  enddo
258 
259  return
scale_grid_index::is
integer, public is
start point of inner domain: x, local
Definition: scale_grid_index.F90:60
scale_grid_index::je
integer, public je
end point of inner domain: y, local
Definition: scale_grid_index.F90:63
scale_const::const_cpdry
real(rp), public const_cpdry
specific heat (dry air,constant pressure) [J/kg/K]
Definition: scale_const.F90:58
scale_atmos_saturation
module ATMOSPHERE / Saturation adjustment
Definition: scale_atmos_sub_saturation.F90:16
scale_tracer::qa
integer, public qa
Definition: scale_tracer.F90:44
scale_const::const_rdry
real(rp), public const_rdry
specific gas constant (dry air) [J/kg/K]
Definition: scale_const.F90:57
scale_grid_index
module grid index
Definition: scale_grid_index.F90:14
scale_tracer
module TRACER
Definition: scale_tracer.F90:18
scale_grid_index::ia
integer, public ia
of x whole cells (local, with HALO)
Definition: scale_grid_index.F90:55
scale_tracer::i_qv
integer, public i_qv
Definition: scale_tracer.F90:47
scale_bulkflux::bulkflux
procedure(bc), pointer, public bulkflux
Definition: scale_bulkflux.F90:71
scale_grid_index::js
integer, public js
start point of inner domain: y, local
Definition: scale_grid_index.F90:62
scale_roughness::roughness
procedure(rl), pointer, public roughness
Definition: scale_roughness.F90:61
scale_atmos_phy_sf_bulkcoef::atmos_phy_sf_bulkcoef
procedure(bc), pointer, public atmos_phy_sf_bulkcoef
Definition: scale_atmos_phy_sf_bulkcoef.F90:62
dc_log::log
subroutine, public log(type, message)
Definition: dc_log.f90:133
scale_const
module CONSTANT
Definition: scale_const.F90:14
scale_grid_index::ie
integer, public ie
end point of inner domain: x, local
Definition: scale_grid_index.F90:61
scale_bulkflux
module Surface bulk flux
Definition: scale_bulkflux.F90:12
scale_atmos_thermodyn
module ATMOSPHERE / Thermodynamics
Definition: scale_atmos_sub_thermodyn.F90:19
scale_roughness
module Surface roughness length
Definition: scale_roughness.F90:13
scale_atmos_phy_sf_bulkcoef
module ATMOSPHERE / Physics Surface bulk coefficient
Definition: scale_atmos_phy_sf_bulkcoef.F90:15
scale_grid_index::ja
integer, public ja
of y whole cells (local, with HALO)
Definition: scale_grid_index.F90:56
Here is the call graph for this function:
Here is the caller graph for this function:
Generated by   doxygen 1.8.13