SCALE-RM
Functions/Subroutines
scale_atmos_dyn_fvm_flux_ud1 Module Reference

module scale_atmos_dyn_fvm_flux_ud1 More...

Functions/Subroutines

subroutine, public atmos_dyn_fvm_flux_valuew_z_ud1 (valW, mflx, val, GSQRT, CDZ)
 value at XYW More...
 
subroutine, public atmos_dyn_fvm_fluxz_xyz_ud1 (flux, mflx, val, GSQRT, CDZ, IIS, IIE, JJS, JJE)
 calculation z-flux at XYZ More...
 
subroutine, public atmos_dyn_fvm_fluxx_xyz_ud1 (flux, mflx, val, GSQRT, CDZ, IIS, IIE, JJS, JJE)
 calculation X-flux at XYZ More...
 
subroutine, public atmos_dyn_fvm_fluxy_xyz_ud1 (flux, mflx, val, GSQRT, CDZ, IIS, IIE, JJS, JJE)
 calculation Y-flux at XYZ More...
 
subroutine, public atmos_dyn_fvm_fluxz_xyw_ud1 (flux, mom, val, DENS, GSQRT, J33G, CDZ, FDZ, dtrk, IIS, IIE, JJS, JJE)
 calculation z-flux at XYW More...
 
subroutine, public atmos_dyn_fvm_fluxj13_xyw_ud1 (flux, mom, val, DENS, GSQRT, J13G, MAPF, CDZ, IIS, IIE, JJS, JJE)
 calculation J13-flux at XYW More...
 
subroutine, public atmos_dyn_fvm_fluxj23_xyw_ud1 (flux, mom, val, DENS, GSQRT, J23G, MAPF, CDZ, IIS, IIE, JJS, JJE)
 calculation J23-flux at XYW More...
 
subroutine, public atmos_dyn_fvm_fluxx_xyw_ud1 (flux, mom, val, DENS, GSQRT, MAPF, CDZ, IIS, IIE, JJS, JJE)
 calculation X-flux at XYW More...
 
subroutine, public atmos_dyn_fvm_fluxy_xyw_ud1 (flux, mom, val, DENS, GSQRT, MAPF, CDZ, IIS, IIE, JJS, JJE)
 calculation Y-flux at XYW More...
 
subroutine, public atmos_dyn_fvm_fluxz_uyz_ud1 (flux, mom, val, DENS, GSQRT, J33G, CDZ, IIS, IIE, JJS, JJE)
 calculation z-flux at UY More...
 
subroutine, public atmos_dyn_fvm_fluxj13_uyz_ud1 (flux, mom, val, DENS, GSQRT, J13G, MAPF, CDZ, IIS, IIE, JJS, JJE)
 calculation J13-flux at UYZ More...
 
subroutine, public atmos_dyn_fvm_fluxj23_uyz_ud1 (flux, mom, val, DENS, GSQRT, J23G, MAPF, CDZ, IIS, IIE, JJS, JJE)
 calculation J23-flux at UYZ More...
 
subroutine, public atmos_dyn_fvm_fluxx_uyz_ud1 (flux, mom, val, DENS, GSQRT, MAPF, CDZ, IIS, IIE, JJS, JJE)
 calculation X-flux at UY More...
 
subroutine, public atmos_dyn_fvm_fluxy_uyz_ud1 (flux, mom, val, DENS, GSQRT, MAPF, CDZ, IIS, IIE, JJS, JJE)
 calculation Y-flux at UY More...
 
subroutine, public atmos_dyn_fvm_fluxz_xvz_ud1 (flux, mom, val, DENS, GSQRT, J33G, CDZ, IIS, IIE, JJS, JJE)
 calculation z-flux at XV More...
 
subroutine, public atmos_dyn_fvm_fluxj13_xvz_ud1 (flux, mom, val, DENS, GSQRT, J13G, MAPF, CDZ, IIS, IIE, JJS, JJE)
 calculation J13-flux at XVZ More...
 
subroutine, public atmos_dyn_fvm_fluxj23_xvz_ud1 (flux, mom, val, DENS, GSQRT, J23G, MAPF, CDZ, IIS, IIE, JJS, JJE)
 calculation J23-flux at XVZ More...
 
subroutine, public atmos_dyn_fvm_fluxx_xvz_ud1 (flux, mom, val, DENS, GSQRT, MAPF, CDZ, IIS, IIE, JJS, JJE)
 calculation X-flux at XV More...
 
subroutine, public atmos_dyn_fvm_fluxy_xvz_ud1 (flux, mom, val, DENS, GSQRT, MAPF, CDZ, IIS, IIE, JJS, JJE)
 calculation Y-flux at XV More...
 

Detailed Description

module scale_atmos_dyn_fvm_flux_ud1

Description
FVM flux scheme with the ud1 order
Author
Team SCALE
History
  • 2016-04-18 (S.Nishizawa) [new]

Function/Subroutine Documentation

◆ atmos_dyn_fvm_flux_valuew_z_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_flux_valuew_z_ud1 ( real(rp), dimension (ka), intent(out)  valW,
real(rp), dimension (ka), intent(in)  mflx,
real(rp), dimension (ka), intent(in)  val,
real(rp), dimension(ka), intent(in)  GSQRT,
real(rp), dimension (ka), intent(in)  CDZ 
)

value at XYW

Definition at line 101 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_debug::check(), scale_grid_index::ke, and scale_grid_index::ks.

101  implicit none
102 
103  real(RP), intent(out) :: valw (ka)
104  real(RP), intent(in) :: mflx (ka)
105  real(RP), intent(in) :: val (ka)
106  real(RP), intent(in) :: gsqrt(ka)
107  real(RP), intent(in) :: cdz (ka)
108 
109  integer :: k
110  !---------------------------------------------------------------------------
111 
112  do k = ks, ke-1
113 #ifdef DEBUG
114  call check( __line__, mflx(k) )
115 
116  call check( __line__, val(k) )
117  call check( __line__, val(k+1) )
118 
119 #endif
120  valw(k) = f1 * ( val(k+1)+val(k) ) - sign(f1,mflx(k)) * ( val(k+1)-val(k) )
121  enddo
122 #ifdef DEBUG
123  k = iundef
124 #endif
125 
126 #ifdef DEBUG
127 
128 #endif
129 
130  return
Here is the call graph for this function:

◆ atmos_dyn_fvm_fluxz_xyz_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxz_xyz_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mflx,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation z-flux at XYZ

Definition at line 141 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_debug::check(), scale_grid_index::ke, and scale_grid_index::ks.

Referenced by scale_atmos_dyn_tinteg_large_euler::atmos_dyn_tinteg_large_euler(), scale_atmos_dyn_tinteg_large_rk3::atmos_dyn_tinteg_large_rk3(), scale_atmos_dyn_tstep_large_fvm_heve::atmos_dyn_tstep_large_fvm_heve(), scale_atmos_dyn_tstep_short_fvm_heve::atmos_dyn_tstep_short_fvm_heve(), and scale_atmos_dyn_tstep_tracer_fvm_heve::atmos_dyn_tstep_tracer_fvm_heve().

141  implicit none
142 
143  real(RP), intent(out) :: flux (ka,ia,ja)
144  real(RP), intent(in) :: mflx (ka,ia,ja)
145  real(RP), intent(in) :: val (ka,ia,ja)
146  real(RP), intent(in) :: gsqrt (ka,ia,ja)
147 
148  real(RP), intent(in) :: cdz (ka)
149  integer, intent(in) :: iis, iie, jjs, jje
150 
151  real(RP) :: vel
152  integer :: k, i, j
153  !---------------------------------------------------------------------------
154 
155  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
156  do j = jjs, jje
157  do i = iis, iie
158  do k = ks, ke-1
159 #ifdef DEBUG
160  call check( __line__, mflx(k,i,j) )
161 
162  call check( __line__, val(k,i,j) )
163  call check( __line__, val(k+1,i,j) )
164 
165 #endif
166  vel = mflx(k,i,j)
167  flux(k,i,j) = vel &
168  * ( f1 * ( val(k+1,i,j)+val(k,i,j) ) - sign(f1,vel) * ( val(k+1,i,j)-val(k,i,j) ) )
169  enddo
170  enddo
171  enddo
172 #ifdef DEBUG
173  k = iundef; i = iundef; j = iundef
174 #endif
175 
176  !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
177  do j = jjs, jje
178  do i = iis, iie
179 #ifdef DEBUG
180 
181 #endif
182  flux(ks-1,i,j) = 0.0_rp
183 
184  flux(ke ,i,j) = 0.0_rp
185  enddo
186  enddo
187 #ifdef DEBUG
188  k = iundef; i = iundef; j = iundef
189 #endif
190 
191  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxx_xyz_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxx_xyz_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mflx,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension(ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation X-flux at XYZ

Definition at line 202 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_debug::check(), scale_grid_index::ke, and scale_grid_index::ks.

Referenced by scale_atmos_dyn_tinteg_large_euler::atmos_dyn_tinteg_large_euler(), scale_atmos_dyn_tinteg_large_rk3::atmos_dyn_tinteg_large_rk3(), scale_atmos_dyn_tstep_large_fvm_heve::atmos_dyn_tstep_large_fvm_heve(), scale_atmos_dyn_tstep_short_fvm_heve::atmos_dyn_tstep_short_fvm_heve(), and scale_atmos_dyn_tstep_tracer_fvm_heve::atmos_dyn_tstep_tracer_fvm_heve().

202  implicit none
203 
204  real(RP), intent(out) :: flux (ka,ia,ja)
205  real(RP), intent(in) :: mflx (ka,ia,ja)
206  real(RP), intent(in) :: val (ka,ia,ja)
207  real(RP), intent(in) :: gsqrt (ka,ia,ja)
208 
209  real(RP), intent(in) :: cdz(ka)
210  integer, intent(in) :: iis, iie, jjs, jje
211 
212  real(RP) :: vel
213  integer :: k, i, j
214  !---------------------------------------------------------------------------
215 
216  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
217  do j = jjs, jje
218  do i = iis-1, iie
219  do k = ks, ke
220 #ifdef DEBUG
221  call check( __line__, mflx(k,i,j) )
222 
223  call check( __line__, val(k,i,j) )
224  call check( __line__, val(k,i+1,j) )
225 
226 #endif
227  vel = mflx(k,i,j)
228  flux(k,i,j) = vel &
229  * ( f1 * ( val(k,i+1,j)+val(k,i,j) ) - sign(f1,vel) * ( val(k,i+1,j)-val(k,i,j) ) )
230  enddo
231  enddo
232  enddo
233 #ifdef DEBUG
234  k = iundef; i = iundef; j = iundef
235 #endif
236 
237  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxy_xyz_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxy_xyz_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mflx,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension(ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation Y-flux at XYZ

Definition at line 248 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_debug::check(), scale_grid_index::ke, and scale_grid_index::ks.

Referenced by scale_atmos_dyn_tinteg_large_euler::atmos_dyn_tinteg_large_euler(), scale_atmos_dyn_tinteg_large_rk3::atmos_dyn_tinteg_large_rk3(), scale_atmos_dyn_tstep_large_fvm_heve::atmos_dyn_tstep_large_fvm_heve(), scale_atmos_dyn_tstep_short_fvm_heve::atmos_dyn_tstep_short_fvm_heve(), and scale_atmos_dyn_tstep_tracer_fvm_heve::atmos_dyn_tstep_tracer_fvm_heve().

248  implicit none
249 
250  real(RP), intent(out) :: flux (ka,ia,ja)
251  real(RP), intent(in) :: mflx (ka,ia,ja)
252  real(RP), intent(in) :: val (ka,ia,ja)
253  real(RP), intent(in) :: gsqrt (ka,ia,ja)
254 
255  real(RP), intent(in) :: cdz(ka)
256  integer, intent(in) :: iis, iie, jjs, jje
257 
258  real(RP) :: vel
259  integer :: k, i, j
260  !---------------------------------------------------------------------------
261 
262  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
263  do j = jjs-1, jje
264  do i = iis, iie
265  do k = ks, ke
266 #ifdef DEBUG
267  call check( __line__, mflx(k,i,j) )
268 
269  call check( __line__, val(k,i,j) )
270  call check( __line__, val(k,i,j+1) )
271 
272 #endif
273  vel = mflx(k,i,j)
274  flux(k,i,j) = vel &
275  * ( f1 * ( val(k,i,j+1)+val(k,i,j) ) - sign(f1,vel) * ( val(k,i,j+1)-val(k,i,j) ) )
276  enddo
277  enddo
278  enddo
279 #ifdef DEBUG
280  k = iundef; i = iundef; j = iundef
281 #endif
282 
283  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxz_xyw_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxz_xyw_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), intent(in)  J33G,
real(rp), dimension (ka), intent(in)  CDZ,
real(rp), dimension (ka-1), intent(in)  FDZ,
real(rp), intent(in)  dtrk,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation z-flux at XYW

Definition at line 297 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_debug::check(), scale_grid_index::ke, and scale_grid_index::ks.

Referenced by scale_atmos_dyn_tstep_short_fvm_heve::atmos_dyn_tstep_short_fvm_heve().

297  implicit none
298 
299  real(RP), intent(out) :: flux (ka,ia,ja)
300  real(RP), intent(in) :: mom (ka,ia,ja)
301  real(RP), intent(in) :: val (ka,ia,ja)
302  real(RP), intent(in) :: dens (ka,ia,ja)
303  real(RP), intent(in) :: gsqrt (ka,ia,ja)
304  real(RP), intent(in) :: j33g
305 
306  real(RP), intent(in) :: cdz (ka)
307  real(RP), intent(in) :: fdz (ka-1)
308  real(RP), intent(in) :: dtrk
309  integer, intent(in) :: iis, iie, jjs, jje
310 
311  real(RP) :: vel
312  real(RP) :: sw
313  integer :: k, i, j
314  !---------------------------------------------------------------------------
315 
316  ! note than z-index is added by -1
317 
318  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
319  do j = jjs, jje
320  do i = iis, iie
321  do k = ks+1, ke-1
322 #ifdef DEBUG
323  call check( __line__, mom(k-1,i,j) )
324  call check( __line__, mom(k ,i,j) )
325 
326  call check( __line__, val(k-1,i,j) )
327  call check( __line__, val(k,i,j) )
328 
329 #endif
330  vel = ( 0.5_rp * ( mom(k-1,i,j) &
331  + mom(k,i,j) ) ) &
332  / dens(k,i,j)
333  flux(k-1,i,j) = j33g * vel &
334  * ( f1 * ( val(k,i,j)+val(k-1,i,j) ) - sign(f1,vel) * ( val(k,i,j)-val(k-1,i,j) ) )
335  enddo
336  enddo
337  enddo
338 #ifdef DEBUG
339  k = iundef; i = iundef; j = iundef
340 #endif
341 
342  !$omp parallel do private(i,j,sw) OMP_SCHEDULE_ collapse(2)
343  do j = jjs, jje
344  do i = iis, iie
345 #ifdef DEBUG
346 
347 #endif
348  flux(ks-1,i,j) = 0.0_rp ! k = KS
349 
350  ! if w>0; min(f,w*dz/dt)
351  ! else ; max(f,w*dz/dt) = -min(-f,-w*dz/dt)
352  sw = sign( 1.0_rp, mom(ks,i,j) )
353  flux(ks ,i,j) = sw * min( sw*flux(ks,i,j), sw*val(ks,i,j)*gsqrt(ks,i,j)*fdz(ks)/dtrk )
354 
355 
356  flux(ke-1,i,j) = 0.0_rp ! k = KE
357  flux(ke ,i,j) = 0.0_rp ! k = KE+1
358  enddo
359  enddo
360 
361  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxj13_xyw_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxj13_xyw_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension (ka,ia,ja), intent(in)  J13G,
real(rp), dimension ( ia,ja,2), intent(in)  MAPF,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J13-flux at XYW

Definition at line 373 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_grid_index::ke, and scale_grid_index::ks.

373  implicit none
374 
375  real(RP), intent(out) :: flux (ka,ia,ja)
376  real(RP), intent(in) :: mom (ka,ia,ja)
377  real(RP), intent(in) :: val (ka,ia,ja)
378  real(RP), intent(in) :: dens (ka,ia,ja)
379  real(RP), intent(in) :: gsqrt (ka,ia,ja)
380  real(RP), intent(in) :: j13g (ka,ia,ja)
381  real(RP), intent(in) :: mapf ( ia,ja,2)
382  real(RP), intent(in) :: cdz (ka)
383  integer, intent(in) :: iis, iie, jjs, jje
384 
385  real(RP) :: vel
386  integer :: k, i, j
387  !---------------------------------------------------------------------------
388 
389  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
390  do j = jjs, jje
391  do i = iis, iie
392  do k = ks+1, ke-1
393  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i-1,j) ) ) &
394  / dens(k,i,j)
395  flux(k-1,i,j) = j13g(k,i,j) / mapf(i,j,+2) * vel &
396  * ( f1 * ( val(k,i,j)+val(k-1,i,j) ) - sign(f1,vel) * ( val(k,i,j)-val(k-1,i,j) ) )
397  enddo
398  enddo
399  enddo
400 
401  !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
402  do j = jjs, jje
403  do i = iis, iie
404  flux(ks-1,i,j) = 0.0_rp
405 
406  flux(ke-1,i,j) = 0.0_rp
407  enddo
408  enddo
409 
410  return

◆ atmos_dyn_fvm_fluxj23_xyw_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxj23_xyw_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension (ka,ia,ja), intent(in)  J23G,
real(rp), dimension ( ia,ja,2), intent(in)  MAPF,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J23-flux at XYW

Definition at line 421 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_grid_index::ke, and scale_grid_index::ks.

421  implicit none
422 
423  real(RP), intent(out) :: flux (ka,ia,ja)
424  real(RP), intent(in) :: mom (ka,ia,ja)
425  real(RP), intent(in) :: val (ka,ia,ja)
426  real(RP), intent(in) :: dens (ka,ia,ja)
427  real(RP), intent(in) :: gsqrt (ka,ia,ja)
428  real(RP), intent(in) :: j23g (ka,ia,ja)
429  real(RP), intent(in) :: mapf ( ia,ja,2)
430  real(RP), intent(in) :: cdz (ka)
431  integer, intent(in) :: iis, iie, jjs, jje
432 
433  real(RP) :: vel
434  integer :: k, i, j
435  !---------------------------------------------------------------------------
436 
437  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
438  do j = jjs, jje
439  do i = iis, iie
440  do k = ks+1, ke-1
441  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j-1) ) ) &
442  / dens(k,i,j)
443  flux(k-1,i,j) = j23g(k,i,j) / mapf(i,j,+1) * vel &
444  * ( f1 * ( val(k,i,j)+val(k-1,i,j) ) - sign(f1,vel) * ( val(k,i,j)-val(k-1,i,j) ) )
445  enddo
446  enddo
447  enddo
448 
449  !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
450  do j = jjs, jje
451  do i = iis, iie
452  flux(ks-1,i,j) = 0.0_rp
453 
454  flux(ke-1,i,j) = 0.0_rp
455  enddo
456  enddo
457 
458  return

◆ atmos_dyn_fvm_fluxx_xyw_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxx_xyw_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension ( ia,ja,2), intent(in)  MAPF,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation X-flux at XYW

Definition at line 471 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_debug::check(), scale_grid_index::ke, and scale_grid_index::ks.

Referenced by scale_atmos_dyn_tstep_short_fvm_heve::atmos_dyn_tstep_short_fvm_heve().

471  implicit none
472 
473  real(RP), intent(out) :: flux (ka,ia,ja)
474  real(RP), intent(in) :: mom (ka,ia,ja)
475  real(RP), intent(in) :: val (ka,ia,ja)
476  real(RP), intent(in) :: dens (ka,ia,ja)
477  real(RP), intent(in) :: gsqrt (ka,ia,ja)
478  real(RP), intent(in) :: mapf ( ia,ja,2)
479 
480  real(RP), intent(in) :: cdz (ka)
481  integer, intent(in) :: iis, iie, jjs, jje
482 
483  real(RP) :: vel
484  integer :: k, i, j
485  !---------------------------------------------------------------------------
486 
487  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
488  do j = jjs, jje
489  do i = iis-1, iie
490  do k = ks, ke-1
491 #ifdef DEBUG
492  call check( __line__, mom(k ,i,j) )
493  call check( __line__, mom(k+1,i,j) )
494 
495  call check( __line__, val(k,i,j) )
496  call check( __line__, val(k,i+1,j) )
497 
498 #endif
499  vel = ( f2h(k,1,i_uyz) &
500  * mom(k+1,i,j) &
501  + f2h(k,2,i_uyz) &
502  * mom(k,i,j) ) &
503  / ( f2h(k,1,i_xyz) &
504  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
505  + f2h(k,2,i_xyz) &
506  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
507  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+2) * vel &
508  * ( f1 * ( val(k,i+1,j)+val(k,i,j) ) - sign(f1,vel) * ( val(k,i+1,j)-val(k,i,j) ) )
509  enddo
510  enddo
511  enddo
512 #ifdef DEBUG
513  k = iundef; i = iundef; j = iundef
514 #endif
515 
516  !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
517  do j = jjs, jje
518  do i = iis-1, iie
519  flux(ke,i,j) = 0.0_rp
520  enddo
521  enddo
522 #ifdef DEBUG
523  k = iundef; i = iundef; j = iundef
524 #endif
525 
526  return
scale_gridtrans::i_uyz
integer, public i_uyz
Definition: scale_gridtrans.F90:53
scale_gridtrans::i_xyz
integer, public i_xyz
Definition: scale_gridtrans.F90:49
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxy_xyw_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxy_xyw_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension ( ia,ja,2), intent(in)  MAPF,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation Y-flux at XYW

Definition at line 538 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_debug::check(), scale_grid_index::ke, and scale_grid_index::ks.

Referenced by scale_atmos_dyn_tstep_short_fvm_heve::atmos_dyn_tstep_short_fvm_heve().

538  implicit none
539 
540  real(RP), intent(out) :: flux (ka,ia,ja)
541  real(RP), intent(in) :: mom (ka,ia,ja)
542  real(RP), intent(in) :: val (ka,ia,ja)
543  real(RP), intent(in) :: dens (ka,ia,ja)
544  real(RP), intent(in) :: gsqrt (ka,ia,ja)
545  real(RP), intent(in) :: mapf ( ia,ja,2)
546 
547  real(RP), intent(in) :: cdz (ka)
548  integer, intent(in) :: iis, iie, jjs, jje
549 
550  real(RP) :: vel
551  integer :: k, i, j
552  !---------------------------------------------------------------------------
553 
554  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
555  do j = jjs-1, jje
556  do i = iis, iie
557  do k = ks, ke-1
558 #ifdef DEBUG
559  call check( __line__, mom(k ,i,j) )
560  call check( __line__, mom(k+1,i,j) )
561 
562  call check( __line__, val(k,i,j) )
563  call check( __line__, val(k,i,j+1) )
564 
565 #endif
566  vel = ( f2h(k,1,i_xvz) &
567  * mom(k+1,i,j) &
568  + f2h(k,2,i_xvz) &
569  * mom(k,i,j) ) &
570  / ( f2h(k,1,i_xyz) &
571  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
572  + f2h(k,2,i_xyz) &
573  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
574  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+1) * vel &
575  * ( f1 * ( val(k,i,j+1)+val(k,i,j) ) - sign(f1,vel) * ( val(k,i,j+1)-val(k,i,j) ) )
576  enddo
577  enddo
578  enddo
579 #ifdef DEBUG
580  k = iundef; i = iundef; j = iundef
581 #endif
582 
583  !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
584  do j = jjs-1, jje
585  do i = iis, iie
586  flux(ke,i,j) = 0.0_rp
587  enddo
588  enddo
589 #ifdef DEBUG
590  k = iundef; i = iundef; j = iundef
591 #endif
592 
593  return
scale_gridtrans::i_xvz
integer, public i_xvz
Definition: scale_gridtrans.F90:54
scale_gridtrans::i_xyz
integer, public i_xyz
Definition: scale_gridtrans.F90:49
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxz_uyz_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxz_uyz_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), intent(in)  J33G,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation z-flux at UY

Definition at line 606 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_debug::check(), scale_grid_index::ke, and scale_grid_index::ks.

Referenced by scale_atmos_dyn_tstep_short_fvm_heve::atmos_dyn_tstep_short_fvm_heve().

606  implicit none
607 
608  real(RP), intent(out) :: flux (ka,ia,ja)
609  real(RP), intent(in) :: mom (ka,ia,ja)
610  real(RP), intent(in) :: val (ka,ia,ja)
611  real(RP), intent(in) :: dens (ka,ia,ja)
612  real(RP), intent(in) :: gsqrt (ka,ia,ja)
613  real(RP), intent(in) :: j33g
614 
615  real(RP), intent(in) :: cdz (ka)
616  integer, intent(in) :: iis, iie, jjs, jje
617 
618  real(RP) :: vel
619  integer :: k, i, j
620  !---------------------------------------------------------------------------
621 
622  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
623  do j = jjs, jje
624  do i = iis, iie
625  do k = ks, ke-1
626 #ifdef DEBUG
627  call check( __line__, mom(k,i,j) )
628  call check( __line__, mom(k,i+1,j) )
629 
630  call check( __line__, val(k,i,j) )
631  call check( __line__, val(k+1,i,j) )
632 
633 #endif
634  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i+1,j) ) ) &
635  / ( f2h(k,1,i_xyz) &
636  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
637  + f2h(k,2,i_xyz) &
638  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
639  flux(k,i,j) = j33g * vel &
640  * ( f1 * ( val(k+1,i,j)+val(k,i,j) ) - sign(f1,vel) * ( val(k+1,i,j)-val(k,i,j) ) )
641  enddo
642  enddo
643  enddo
644 #ifdef DEBUG
645  k = iundef; i = iundef; j = iundef
646 #endif
647 
648  !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
649  do j = jjs, jje
650  do i = iis, iie
651 #ifdef DEBUG
652 
653 #endif
654  flux(ks-1,i,j) = 0.0_rp
655 
656  flux(ke,i,j) = 0.0_rp
657  enddo
658  enddo
659 #ifdef DEBUG
660  k = iundef; i = iundef; j = iundef
661 #endif
662 
663  return
scale_gridtrans::i_xyz
integer, public i_xyz
Definition: scale_gridtrans.F90:49
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxj13_uyz_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxj13_uyz_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension (ka,ia,ja), intent(in)  J13G,
real(rp), dimension ( ia,ja,2), intent(in)  MAPF,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J13-flux at UYZ

Definition at line 674 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_grid_index::ke, and scale_grid_index::ks.

674  implicit none
675 
676  real(RP), intent(out) :: flux (ka,ia,ja)
677  real(RP), intent(in) :: mom (ka,ia,ja)
678  real(RP), intent(in) :: val (ka,ia,ja)
679  real(RP), intent(in) :: dens (ka,ia,ja)
680  real(RP), intent(in) :: gsqrt (ka,ia,ja)
681  real(RP), intent(in) :: j13g (ka,ia,ja)
682  real(RP), intent(in) :: mapf ( ia,ja,2)
683  real(RP), intent(in) :: cdz (ka)
684  integer, intent(in) :: iis, iie, jjs, jje
685 
686  real(RP) :: vel
687  integer :: k, i, j
688  !---------------------------------------------------------------------------
689 
690  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
691  do j = jjs, jje
692  do i = iis, iie
693  do k = ks, ke-1
694  vel = ( f2h(k,1,i_uyz) &
695  * mom(k+1,i,j) &
696  + f2h(k,2,i_uyz) &
697  * mom(k,i,j) ) &
698  / ( f2h(k,1,i_xyz) &
699  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
700  + f2h(k,2,i_xyz) &
701  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
702  flux(k,i,j) = j13g(k,i,j) / mapf(i,j,+2) * vel &
703  * ( f1 * ( val(k+1,i,j)+val(k,i,j) ) - sign(f1,vel) * ( val(k+1,i,j)-val(k,i,j) ) )
704  enddo
705  enddo
706  enddo
707 
708  !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
709  do j = jjs, jje
710  do i = iis, iie
711  flux(ks-1,i,j) = 0.0_rp
712 
713  flux(ke ,i,j) = 0.0_rp
714  enddo
715  enddo
716 
717  return
scale_gridtrans::i_uyz
integer, public i_uyz
Definition: scale_gridtrans.F90:53
scale_gridtrans::i_xyz
integer, public i_xyz
Definition: scale_gridtrans.F90:49

◆ atmos_dyn_fvm_fluxj23_uyz_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxj23_uyz_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension (ka,ia,ja), intent(in)  J23G,
real(rp), dimension ( ia,ja,2), intent(in)  MAPF,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J23-flux at UYZ

Definition at line 728 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_grid_index::ke, and scale_grid_index::ks.

728  implicit none
729 
730  real(RP), intent(out) :: flux (ka,ia,ja)
731  real(RP), intent(in) :: mom (ka,ia,ja)
732  real(RP), intent(in) :: val (ka,ia,ja)
733  real(RP), intent(in) :: dens (ka,ia,ja)
734  real(RP), intent(in) :: gsqrt (ka,ia,ja)
735  real(RP), intent(in) :: j23g (ka,ia,ja)
736  real(RP), intent(in) :: mapf ( ia,ja,2)
737  real(RP), intent(in) :: cdz (ka)
738  integer, intent(in) :: iis, iie, jjs, jje
739 
740  real(RP) :: vel
741  integer :: k, i, j
742  !---------------------------------------------------------------------------
743 
744  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
745  do j = jjs, jje
746  do i = iis, iie
747  do k = ks, ke-1
748  vel = ( f2h(k,1,i_xvz) &
749  * 0.25_rp * ( mom(k+1,i,j)+mom(k+1,i+1,j)+mom(k+1,i,j-1)+mom(k+1,i+1,j-1) ) &
750  + f2h(k,2,i_xvz) &
751  * 0.25_rp * ( mom(k,i,j)+mom(k,i+1,j)+mom(k,i,j-1)+mom(k,i+1,j-1) ) ) &
752  / ( f2h(k,1,i_xyz) &
753  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
754  + f2h(k,2,i_xyz) &
755  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
756  flux(k,i,j) = j23g(k,i,j) / mapf(i,j,+1) * vel &
757  * ( f1 * ( val(k+1,i,j)+val(k,i,j) ) - sign(f1,vel) * ( val(k+1,i,j)-val(k,i,j) ) )
758  enddo
759  enddo
760  enddo
761 
762  !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
763  do j = jjs, jje
764  do i = iis, iie
765  flux(ks-1,i,j) = 0.0_rp
766 
767  flux(ke ,i,j) = 0.0_rp
768  enddo
769  enddo
770 
771  return
scale_gridtrans::i_xvz
integer, public i_xvz
Definition: scale_gridtrans.F90:54
scale_gridtrans::i_xyz
integer, public i_xyz
Definition: scale_gridtrans.F90:49

◆ atmos_dyn_fvm_fluxx_uyz_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxx_uyz_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension ( ia,ja,2), intent(in)  MAPF,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation X-flux at UY

Definition at line 783 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_debug::check(), scale_grid_index::ke, and scale_grid_index::ks.

Referenced by scale_atmos_dyn_tstep_short_fvm_heve::atmos_dyn_tstep_short_fvm_heve().

783  implicit none
784 
785  real(RP), intent(out) :: flux (ka,ia,ja)
786  real(RP), intent(in) :: mom (ka,ia,ja)
787  real(RP), intent(in) :: val (ka,ia,ja)
788  real(RP), intent(in) :: dens (ka,ia,ja)
789  real(RP), intent(in) :: gsqrt (ka,ia,ja)
790  real(RP), intent(in) :: mapf ( ia,ja,2)
791 
792  real(RP), intent(in) :: cdz (ka)
793  integer, intent(in) :: iis, iie, jjs, jje
794 
795  real(RP) :: vel
796  integer :: k, i, j
797  !---------------------------------------------------------------------------
798 
799  ! note that x-index is added by -1
800 
801  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
802  do j = jjs, jje
803  do i = iis, iie+1
804  do k = ks, ke
805 #ifdef DEBUG
806  call check( __line__, mom(k,i ,j) )
807  call check( __line__, mom(k,i-1,j) )
808 
809  call check( __line__, val(k,i-1,j) )
810  call check( __line__, val(k,i,j) )
811 
812 #endif
813  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i-1,j) ) ) &
814  / ( dens(k,i,j) )
815  flux(k,i-1,j) = gsqrt(k,i,j) / mapf(i,j,+2) * vel &
816  * ( f1 * ( val(k,i,j)+val(k,i-1,j) ) - sign(f1,vel) * ( val(k,i,j)-val(k,i-1,j) ) )
817  enddo
818  enddo
819  enddo
820 #ifdef DEBUG
821  k = iundef; i = iundef; j = iundef
822 #endif
823 
824  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxy_uyz_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxy_uyz_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension ( ia,ja,2), intent(in)  MAPF,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation Y-flux at UY

Definition at line 836 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_debug::check(), scale_grid_index::ke, and scale_grid_index::ks.

Referenced by scale_atmos_dyn_tstep_short_fvm_heve::atmos_dyn_tstep_short_fvm_heve().

836  implicit none
837 
838  real(RP), intent(out) :: flux (ka,ia,ja)
839  real(RP), intent(in) :: mom (ka,ia,ja)
840  real(RP), intent(in) :: val (ka,ia,ja)
841  real(RP), intent(in) :: dens (ka,ia,ja)
842  real(RP), intent(in) :: gsqrt (ka,ia,ja)
843  real(RP), intent(in) :: mapf ( ia,ja,2)
844 
845  real(RP), intent(in) :: cdz (ka)
846  integer, intent(in) :: iis, iie, jjs, jje
847 
848  real(RP) :: vel
849  integer :: k, i, j
850  !---------------------------------------------------------------------------
851 
852  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
853  do j = jjs-1, jje
854  do i = iis, iie
855  do k = ks, ke
856 #ifdef DEBUG
857  call check( __line__, mom(k,i ,j) )
858  call check( __line__, mom(k,i-1,j) )
859 
860  call check( __line__, val(k,i,j) )
861  call check( __line__, val(k,i,j+1) )
862 
863 #endif
864  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i+1,j) ) ) &
865  / ( 0.25_rp * ( dens(k,i,j)+dens(k,i+1,j)+dens(k,i,j+1)+dens(k,i+1,j+1) ) )
866  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+1) * vel &
867  * ( f1 * ( val(k,i,j+1)+val(k,i,j) ) - sign(f1,vel) * ( val(k,i,j+1)-val(k,i,j) ) )
868  enddo
869  enddo
870  enddo
871 #ifdef DEBUG
872  k = iundef; i = iundef; j = iundef
873 #endif
874 
875  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxz_xvz_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxz_xvz_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), intent(in)  J33G,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation z-flux at XV

Definition at line 889 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_debug::check(), scale_grid_index::ke, and scale_grid_index::ks.

Referenced by scale_atmos_dyn_tstep_short_fvm_heve::atmos_dyn_tstep_short_fvm_heve().

889  implicit none
890 
891  real(RP), intent(out) :: flux (ka,ia,ja)
892  real(RP), intent(in) :: mom (ka,ia,ja)
893  real(RP), intent(in) :: val (ka,ia,ja)
894  real(RP), intent(in) :: dens (ka,ia,ja)
895  real(RP), intent(in) :: gsqrt (ka,ia,ja)
896  real(RP), intent(in) :: j33g
897 
898  real(RP), intent(in) :: cdz (ka)
899  integer, intent(in) :: iis, iie, jjs, jje
900 
901  real(RP) :: vel
902  integer :: k, i, j
903  !---------------------------------------------------------------------------
904 
905  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
906  do j = jjs, jje
907  do i = iis, iie
908  do k = ks, ke-1
909 #ifdef DEBUG
910  call check( __line__, mom(k,i,j) )
911  call check( __line__, mom(k,i,j+1) )
912 
913  call check( __line__, val(k,i,j) )
914  call check( __line__, val(k+1,i,j) )
915 
916 #endif
917  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j+1) ) ) &
918  / ( f2h(k,1,i_xyz) &
919  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
920  + f2h(k,2,i_xyz) &
921  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
922  flux(k,i,j) = j33g * vel &
923  * ( f1 * ( val(k+1,i,j)+val(k,i,j) ) - sign(f1,vel) * ( val(k+1,i,j)-val(k,i,j) ) )
924  enddo
925  enddo
926  enddo
927 #ifdef DEBUG
928  k = iundef; i = iundef; j = iundef
929 #endif
930 
931  !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
932  do j = jjs, jje
933  do i = iis, iie
934 #ifdef DEBUG
935 
936 #endif
937  flux(ks-1,i,j) = 0.0_rp
938 
939  flux(ke,i,j) = 0.0_rp
940  enddo
941  enddo
942 #ifdef DEBUG
943  k = iundef; i = iundef; j = iundef
944 #endif
945 
946  return
scale_gridtrans::i_xyz
integer, public i_xyz
Definition: scale_gridtrans.F90:49
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxj13_xvz_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxj13_xvz_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension (ka,ia,ja), intent(in)  J13G,
real(rp), dimension ( ia,ja,2), intent(in)  MAPF,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J13-flux at XVZ

Definition at line 957 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_grid_index::ke, and scale_grid_index::ks.

957  implicit none
958 
959  real(RP), intent(out) :: flux (ka,ia,ja)
960  real(RP), intent(in) :: mom (ka,ia,ja)
961  real(RP), intent(in) :: val (ka,ia,ja)
962  real(RP), intent(in) :: dens (ka,ia,ja)
963  real(RP), intent(in) :: gsqrt (ka,ia,ja)
964  real(RP), intent(in) :: j13g (ka,ia,ja)
965  real(RP), intent(in) :: mapf ( ia,ja,2)
966  real(RP), intent(in) :: cdz (ka)
967  integer, intent(in) :: iis, iie, jjs, jje
968 
969  real(RP) :: vel
970  integer :: k, i, j
971  !---------------------------------------------------------------------------
972 
973  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
974  do j = jjs, jje
975  do i = iis, iie
976  do k = ks, ke-1
977  vel = ( f2h(k,1,i_uyz) &
978  * 0.25_rp * ( mom(k+1,i,j)+mom(k+1,i-1,j)+mom(k+1,i,j+1)+mom(k+1,i-1,j+1) ) &
979  + f2h(k,2,i_uyz) &
980  * 0.25_rp * ( mom(k,i,j)+mom(k,i-1,j)+mom(k,i,j+1)+mom(k,i-1,j+1) ) ) &
981  / ( f2h(k,1,i_xyz) &
982  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
983  + f2h(k,2,i_xyz) &
984  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
985  flux(k,i,j) = j13g(k,i,j) / mapf(i,j,+2) * vel &
986  * ( f1 * ( val(k+1,i,j)+val(k,i,j) ) - sign(f1,vel) * ( val(k+1,i,j)-val(k,i,j) ) )
987  enddo
988  enddo
989  enddo
990 
991  !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
992  do j = jjs, jje
993  do i = iis, iie
994  flux(ks-1,i,j) = 0.0_rp
995 
996  flux(ke ,i,j) = 0.0_rp
997  enddo
998  enddo
999 
1000  return
scale_gridtrans::i_uyz
integer, public i_uyz
Definition: scale_gridtrans.F90:53
scale_gridtrans::i_xyz
integer, public i_xyz
Definition: scale_gridtrans.F90:49

◆ atmos_dyn_fvm_fluxj23_xvz_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxj23_xvz_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension (ka,ia,ja), intent(in)  J23G,
real(rp), dimension ( ia,ja,2), intent(in)  MAPF,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J23-flux at XVZ

Definition at line 1011 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_grid_index::ke, and scale_grid_index::ks.

1011  implicit none
1012 
1013  real(RP), intent(out) :: flux (ka,ia,ja)
1014  real(RP), intent(in) :: mom (ka,ia,ja)
1015  real(RP), intent(in) :: val (ka,ia,ja)
1016  real(RP), intent(in) :: dens (ka,ia,ja)
1017  real(RP), intent(in) :: gsqrt (ka,ia,ja)
1018  real(RP), intent(in) :: j23g (ka,ia,ja)
1019  real(RP), intent(in) :: mapf ( ia,ja,2)
1020  real(RP), intent(in) :: cdz (ka)
1021  integer, intent(in) :: iis, iie, jjs, jje
1022 
1023  real(RP) :: vel
1024  integer :: k, i, j
1025  !---------------------------------------------------------------------------
1026 
1027  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1028  do j = jjs, jje
1029  do i = iis, iie
1030  do k = ks, ke-1
1031  vel = ( f2h(k,1,i_xvz) &
1032  * mom(k+1,i,j) &
1033  + f2h(k,2,i_xvz) &
1034  * mom(k,i,j) ) &
1035  / ( f2h(k,1,i_xyz) &
1036  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
1037  + f2h(k,2,i_xyz) &
1038  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
1039  flux(k,i,j) = j23g(k,i,j) / mapf(i,j,+1) * vel &
1040  * ( f1 * ( val(k+1,i,j)+val(k,i,j) ) - sign(f1,vel) * ( val(k+1,i,j)-val(k,i,j) ) )
1041  enddo
1042  enddo
1043  enddo
1044 
1045  !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
1046  do j = jjs, jje
1047  do i = iis, iie
1048  flux(ks-1,i,j) = 0.0_rp
1049 
1050  flux(ke ,i,j) = 0.0_rp
1051  enddo
1052  enddo
1053 
1054  return
scale_gridtrans::i_xvz
integer, public i_xvz
Definition: scale_gridtrans.F90:54
scale_gridtrans::i_xyz
integer, public i_xyz
Definition: scale_gridtrans.F90:49

◆ atmos_dyn_fvm_fluxx_xvz_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxx_xvz_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension ( ia,ja,2), intent(in)  MAPF,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation X-flux at XV

Definition at line 1066 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_debug::check(), scale_grid_index::ke, and scale_grid_index::ks.

Referenced by scale_atmos_dyn_tstep_short_fvm_heve::atmos_dyn_tstep_short_fvm_heve().

1066  implicit none
1067 
1068  real(RP), intent(out) :: flux (ka,ia,ja)
1069  real(RP), intent(in) :: mom (ka,ia,ja)
1070  real(RP), intent(in) :: val (ka,ia,ja)
1071  real(RP), intent(in) :: dens (ka,ia,ja)
1072  real(RP), intent(in) :: gsqrt (ka,ia,ja)
1073  real(RP), intent(in) :: mapf ( ia,ja,2)
1074 
1075  real(RP), intent(in) :: cdz (ka)
1076  integer, intent(in) :: iis, iie, jjs, jje
1077 
1078  real(RP) :: vel
1079  integer :: k, i, j
1080  !---------------------------------------------------------------------------
1081 
1082  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1083  do j = jjs, jje
1084  do i = iis-1, iie
1085  do k = ks, ke
1086 #ifdef DEBUG
1087  call check( __line__, mom(k,i ,j) )
1088  call check( __line__, mom(k,i,j-1) )
1089 
1090  call check( __line__, val(k,i,j) )
1091  call check( __line__, val(k,i+1,j) )
1092 
1093 #endif
1094  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j+1) ) ) &
1095  / ( 0.25_rp * ( dens(k,i,j)+dens(k,i+1,j)+dens(k,i,j+1)+dens(k,i+1,j+1) ) )
1096  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+2) * vel &
1097  * ( f1 * ( val(k,i+1,j)+val(k,i,j) ) - sign(f1,vel) * ( val(k,i+1,j)-val(k,i,j) ) )
1098  enddo
1099  enddo
1100  enddo
1101 #ifdef DEBUG
1102  k = iundef; i = iundef; j = iundef
1103 #endif
1104 
1105  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxy_xvz_ud1()

subroutine, public scale_atmos_dyn_fvm_flux_ud1::atmos_dyn_fvm_fluxy_xvz_ud1 ( real(rp), dimension (ka,ia,ja), intent(out)  flux,
real(rp), dimension (ka,ia,ja), intent(in)  mom,
real(rp), dimension (ka,ia,ja), intent(in)  val,
real(rp), dimension (ka,ia,ja), intent(in)  DENS,
real(rp), dimension (ka,ia,ja), intent(in)  GSQRT,
real(rp), dimension ( ia,ja,2), intent(in)  MAPF,
real(rp), dimension (ka), intent(in)  CDZ,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation Y-flux at XV

Definition at line 1117 of file scale_atmos_dyn_fvm_flux_ud1.F90.

References scale_debug::check(), scale_grid_index::ke, and scale_grid_index::ks.

Referenced by scale_atmos_dyn_tstep_short_fvm_heve::atmos_dyn_tstep_short_fvm_heve().

1117  implicit none
1118 
1119  real(RP), intent(out) :: flux (ka,ia,ja)
1120  real(RP), intent(in) :: mom (ka,ia,ja)
1121  real(RP), intent(in) :: val (ka,ia,ja)
1122  real(RP), intent(in) :: dens (ka,ia,ja)
1123  real(RP), intent(in) :: gsqrt (ka,ia,ja)
1124  real(RP), intent(in) :: mapf ( ia,ja,2)
1125 
1126  real(RP), intent(in) :: cdz (ka)
1127  integer, intent(in) :: iis, iie, jjs, jje
1128 
1129  real(RP) :: vel
1130  integer :: k, i, j
1131  !---------------------------------------------------------------------------
1132 
1133  ! note that y-index is added by -1
1134 
1135  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1136  do j = jjs, jje+1
1137  do i = iis, iie
1138  do k = ks, ke
1139 #ifdef DEBUG
1140  call check( __line__, mom(k,i ,j) )
1141  call check( __line__, mom(k,i,j-1) )
1142 
1143  call check( __line__, val(k,i,j-1) )
1144  call check( __line__, val(k,i,j) )
1145 
1146 #endif
1147  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j-1) ) ) &
1148  / ( dens(k,i,j) )
1149  flux(k,i,j-1) = gsqrt(k,i,j) / mapf(i,j,+1) * vel &
1150  * ( f1 * ( val(k,i,j)+val(k,i,j-1) ) - sign(f1,vel) * ( val(k,i,j)-val(k,i,j-1) ) )
1151  enddo
1152  enddo
1153  enddo
1154 #ifdef DEBUG
1155  k = iundef; i = iundef; j = iundef
1156 #endif
1157 
1158  return
Here is the call graph for this function:
Here is the caller graph for this function:
Generated by   doxygen 1.8.13