SCALE-RM
Functions/Subroutines
scale_atmos_dyn_fvm_flux_cd2 Module Reference

module scale_atmos_dyn_fvm_flux_cd2 More...

Functions/Subroutines

subroutine, public atmos_dyn_fvm_flux_valuew_z_cd2 (valW, mflx, val, GSQRT, CDZ)
 value at XYW More...
 
subroutine, public atmos_dyn_fvm_fluxz_xyz_cd2 (flux, mflx, val, GSQRT, num_diff, CDZ, IIS, IIE, JJS, JJE)
 calculation z-flux at XYZ More...
 
subroutine, public atmos_dyn_fvm_fluxx_xyz_cd2 (flux, mflx, val, GSQRT, num_diff, CDZ, IIS, IIE, JJS, JJE)
 calculation X-flux at XYZ More...
 
subroutine, public atmos_dyn_fvm_fluxy_xyz_cd2 (flux, mflx, val, GSQRT, num_diff, CDZ, IIS, IIE, JJS, JJE)
 calculation Y-flux at XYZ More...
 
subroutine, public atmos_dyn_fvm_fluxz_xyw_cd2 (flux, mom, val, DENS, GSQRT, J33G, num_diff, CDZ, FDZ, dtrk, IIS, IIE, JJS, JJE)
 calculation z-flux at XYW More...
 
subroutine, public atmos_dyn_fvm_fluxj13_xyw_cd2 (flux, mom, val, DENS, GSQRT, J13G, MAPF, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation J13-flux at XYW More...
 
subroutine, public atmos_dyn_fvm_fluxj23_xyw_cd2 (flux, mom, val, DENS, GSQRT, J23G, MAPF, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation J23-flux at XYW More...
 
subroutine, public atmos_dyn_fvm_fluxx_xyw_cd2 (flux, mom, val, DENS, GSQRT, MAPF, num_diff, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation X-flux at XYW More...
 
subroutine, public atmos_dyn_fvm_fluxy_xyw_cd2 (flux, mom, val, DENS, GSQRT, MAPF, num_diff, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation Y-flux at XYW More...
 
subroutine, public atmos_dyn_fvm_fluxz_uyz_cd2 (flux, mom, val, DENS, GSQRT, J33G, num_diff, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation z-flux at UY More...
 
subroutine, public atmos_dyn_fvm_fluxj13_uyz_cd2 (flux, mom, val, DENS, GSQRT, J13G, MAPF, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation J13-flux at UYZ More...
 
subroutine, public atmos_dyn_fvm_fluxj23_uyz_cd2 (flux, mom, val, DENS, GSQRT, J23G, MAPF, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation J23-flux at UYZ More...
 
subroutine, public atmos_dyn_fvm_fluxx_uyz_cd2 (flux, mom, val, DENS, GSQRT, MAPF, num_diff, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation X-flux at UY More...
 
subroutine, public atmos_dyn_fvm_fluxy_uyz_cd2 (flux, mom, val, DENS, GSQRT, MAPF, num_diff, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation Y-flux at UY More...
 
subroutine, public atmos_dyn_fvm_fluxz_xvz_cd2 (flux, mom, val, DENS, GSQRT, J33G, num_diff, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation z-flux at XV More...
 
subroutine, public atmos_dyn_fvm_fluxj13_xvz_cd2 (flux, mom, val, DENS, GSQRT, J13G, MAPF, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation J13-flux at XVZ More...
 
subroutine, public atmos_dyn_fvm_fluxj23_xvz_cd2 (flux, mom, val, DENS, GSQRT, J23G, MAPF, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation J23-flux at XVZ More...
 
subroutine, public atmos_dyn_fvm_fluxx_xvz_cd2 (flux, mom, val, DENS, GSQRT, MAPF, num_diff, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation X-flux at XV More...
 
subroutine, public atmos_dyn_fvm_fluxy_xvz_cd2 (flux, mom, val, DENS, GSQRT, MAPF, num_diff, CDZ, TwoD, IIS, IIE, JJS, JJE)
 calculation Y-flux at XV More...
 

Detailed Description

module scale_atmos_dyn_fvm_flux_cd2

Description
FVM flux scheme with the cd2 order
Author
Team SCALE

Function/Subroutine Documentation

◆ atmos_dyn_fvm_flux_valuew_z_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_flux_valuew_z_cd2 ( 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 99 of file scale_atmos_dyn_fvm_flux_cd2.F90.

99  !$acc routine vector
100  implicit none
101 
102  real(RP), intent(out) :: valW (KA)
103  real(RP), intent(in) :: mflx (KA)
104  real(RP), intent(in) :: val (KA)
105  real(RP), intent(in) :: GSQRT(KA)
106  real(RP), intent(in) :: CDZ (KA)
107 
108  integer :: k
109  !---------------------------------------------------------------------------
110 
111  do k = ks, ke-1
112 #ifdef DEBUG
113  call check( __line__, mflx(k) )
114 
115  call check( __line__, val(k) )
116  call check( __line__, val(k+1) )
117 
118 #endif
119  valw(k) = f2 * ( val(k+1)+val(k) )
120  enddo
121 #ifdef DEBUG
122  k = iundef
123 #endif
124 
125 #ifdef DEBUG
126 
127 #endif
128 
129 
130  return

References scale_debug::check(), scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxz_xyz_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxz_xyz_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,ia,ja), intent(in)  num_diff,
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_cd2.F90.

141  use scale_const, only: &
142  eps => const_eps
143  implicit none
144 
145  real(RP), intent(inout) :: flux (KA,IA,JA)
146  real(RP), intent(in) :: mflx (KA,IA,JA)
147  real(RP), intent(in) :: val (KA,IA,JA)
148  real(RP), intent(in) :: GSQRT (KA,IA,JA)
149  real(RP), intent(in) :: num_diff(KA,IA,JA)
150  real(RP), intent(in) :: CDZ (KA)
151  integer, intent(in) :: IIS, IIE, JJS, JJE
152 
153  real(RP) :: vel
154  integer :: k, i, j
155  !---------------------------------------------------------------------------
156 
157  !$omp parallel default(none) private(i,j,k, vel) &
158  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mflx,val,flux,GSQRT,num_diff,EPS)
159 
160  !$acc data copy(flux) copyin(mflx, val, GSQRT, num_diff, CDZ)
161 
162  !$omp do OMP_SCHEDULE_ collapse(2)
163  !$acc kernels
164  do j = jjs, jje
165  do i = iis, iie
166  do k = ks, ke-1
167 #ifdef DEBUG
168  call check( __line__, mflx(k,i,j) )
169 
170  call check( __line__, val(k,i,j) )
171  call check( __line__, val(k+1,i,j) )
172 
173 #endif
174  vel = mflx(k,i,j)
175  flux(k,i,j) = vel &
176  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) ) &
177  + gsqrt(k,i,j) * num_diff(k,i,j)
178  enddo
179  enddo
180  enddo
181  !$acc end kernels
182  !$omp end do nowait
183 #ifdef DEBUG
184  k = iundef; i = iundef; j = iundef
185 #endif
186 
187  !$omp do OMP_SCHEDULE_ collapse(2)
188  !$acc kernels
189  do j = jjs, jje
190  do i = iis, iie
191 #ifdef DEBUG
192 
193 #endif
194  flux(ks-1,i,j) = 0.0_rp
195 
196  flux(ke ,i,j) = 0.0_rp
197  enddo
198  enddo
199  !$acc end kernels
200  !$omp end do nowait
201 
202  !$acc end data
203 
204  !$omp end parallel
205 #ifdef DEBUG
206  k = iundef; i = iundef; j = iundef
207 #endif
208 
209  return

References scale_debug::check(), scale_const::const_eps, scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxx_xyz_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxx_xyz_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,ia,ja), intent(in)  num_diff,
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 220 of file scale_atmos_dyn_fvm_flux_cd2.F90.

220  implicit none
221 
222  real(RP), intent(inout) :: flux (KA,IA,JA)
223  real(RP), intent(in) :: mflx (KA,IA,JA)
224  real(RP), intent(in) :: val (KA,IA,JA)
225  real(RP), intent(in) :: GSQRT (KA,IA,JA)
226  real(RP), intent(in) :: num_diff(KA,IA,JA)
227  real(RP), intent(in) :: CDZ(KA)
228  integer, intent(in) :: IIS, IIE, JJS, JJE
229 
230  real(RP) :: vel
231  integer :: k, i, j
232  !---------------------------------------------------------------------------
233 
234  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
235  !$omp private(vel) &
236  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mflx,val,flux,GSQRT,num_diff)
237  !$acc kernels
238  do j = jjs, jje
239  do i = iis-1, iie
240  do k = ks, ke
241 #ifdef DEBUG
242  call check( __line__, mflx(k,i,j) )
243 
244  call check( __line__, val(k,i,j) )
245  call check( __line__, val(k,i+1,j) )
246 
247 #endif
248  vel = mflx(k,i,j)
249  flux(k,i,j) = vel &
250  * ( f2 * ( val(k,i+1,j)+val(k,i,j) ) ) &
251  + gsqrt(k,i,j) * num_diff(k,i,j)
252  enddo
253  enddo
254  enddo
255  !$acc end kernels
256 #ifdef DEBUG
257  k = iundef; i = iundef; j = iundef
258 #endif
259 
260  return

References scale_debug::check(), scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxy_xyz_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxy_xyz_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,ia,ja), intent(in)  num_diff,
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 271 of file scale_atmos_dyn_fvm_flux_cd2.F90.

271  implicit none
272 
273  real(RP), intent(inout) :: flux (KA,IA,JA)
274  real(RP), intent(in) :: mflx (KA,IA,JA)
275  real(RP), intent(in) :: val (KA,IA,JA)
276  real(RP), intent(in) :: GSQRT (KA,IA,JA)
277  real(RP), intent(in) :: num_diff(KA,IA,JA)
278  real(RP), intent(in) :: CDZ(KA)
279  integer, intent(in) :: IIS, IIE, JJS, JJE
280 
281  real(RP) :: vel
282  integer :: k, i, j
283  !---------------------------------------------------------------------------
284 
285  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
286  !$omp private(vel) &
287  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mflx,val,flux,GSQRT,num_diff)
288  !$acc kernels
289  do j = jjs-1, jje
290  do i = iis, iie
291  do k = ks, ke
292 #ifdef DEBUG
293  call check( __line__, mflx(k,i,j) )
294 
295  call check( __line__, val(k,i,j) )
296  call check( __line__, val(k,i,j+1) )
297 
298 #endif
299  vel = mflx(k,i,j)
300  flux(k,i,j) = vel &
301  * ( f2 * ( val(k,i,j+1)+val(k,i,j) ) ) &
302  + gsqrt(k,i,j) * num_diff(k,i,j)
303  enddo
304  enddo
305  enddo
306  !$acc end kernels
307 #ifdef DEBUG
308  k = iundef; i = iundef; j = iundef
309 #endif
310 
311  return

References scale_debug::check(), scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxz_xyw_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxz_xyw_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,ia,ja), intent(in)  num_diff,
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 325 of file scale_atmos_dyn_fvm_flux_cd2.F90.

325  implicit none
326 
327  real(RP), intent(inout) :: flux (KA,IA,JA)
328  real(RP), intent(in) :: mom (KA,IA,JA)
329  real(RP), intent(in) :: val (KA,IA,JA)
330  real(RP), intent(in) :: DENS (KA,IA,JA)
331  real(RP), intent(in) :: GSQRT (KA,IA,JA)
332  real(RP), intent(in) :: J33G
333  real(RP), intent(in) :: num_diff(KA,IA,JA)
334  real(RP), intent(in) :: CDZ (KA)
335  real(RP), intent(in) :: FDZ (KA-1)
336  real(RP), intent(in) :: dtrk
337  integer, intent(in) :: IIS, IIE, JJS, JJE
338 
339  real(RP) :: vel
340  integer :: k, i, j
341  !---------------------------------------------------------------------------
342 
343  ! note than z-index is added by -1
344 
345  !$omp parallel default(none) private(i,j,k,vel) &
346  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,flux,J33G,GSQRT,num_diff,DENS,FDZ,dtrk)
347 
348  !$acc data copy(flux) copyin(mom, val, DENS, GSQRT, num_diff, CDZ, FDZ)
349 
350  !$omp do OMP_SCHEDULE_ collapse(2)
351  !$acc kernels
352  do j = jjs, jje
353  do i = iis, iie
354  do k = ks+1, ke-1
355 #ifdef DEBUG
356  call check( __line__, mom(k-1,i,j) )
357  call check( __line__, mom(k ,i,j) )
358 
359  call check( __line__, val(k-1,i,j) )
360  call check( __line__, val(k,i,j) )
361 
362 #endif
363  vel = ( 0.5_rp * ( mom(k-1,i,j) &
364  + mom(k,i,j) ) ) &
365  / dens(k,i,j)
366  flux(k-1,i,j) = j33g * vel &
367  * ( f2 * ( val(k,i,j)+val(k-1,i,j) ) ) &
368  + gsqrt(k,i,j) * num_diff(k,i,j)
369  enddo
370  enddo
371  enddo
372  !$acc end kernels
373  !$omp end do nowait
374 #ifdef DEBUG
375  k = iundef; i = iundef; j = iundef
376 #endif
377 
378  !$omp do OMP_SCHEDULE_ collapse(2)
379  !$acc kernels
380  do j = jjs, jje
381  do i = iis, iie
382 #ifdef DEBUG
383 
384 
385 #endif
386  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS.
387  ! The flux at KS can be non-zero.
388  ! To reduce calculations, all the fluxes are set to zero.
389  flux(ks-1,i,j) = 0.0_rp ! k = KS
390 
391 
392 
393  flux(ke-1,i,j) = 0.0_rp ! k = KE
394  flux(ke ,i,j) = 0.0_rp ! k = KE+1
395  enddo
396  enddo
397  !$acc end kernels
398  !$omp end do nowait
399 
400  !$acc end data
401 
402  !$omp end parallel
403 
404  return

References scale_debug::check(), scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxj13_xyw_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxj13_xyw_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J13-flux at XYW

Definition at line 416 of file scale_atmos_dyn_fvm_flux_cd2.F90.

416  implicit none
417 
418  real(RP), intent(inout) :: flux (KA,IA,JA)
419  real(RP), intent(in) :: mom (KA,IA,JA)
420  real(RP), intent(in) :: val (KA,IA,JA)
421  real(RP), intent(in) :: DENS (KA,IA,JA)
422  real(RP), intent(in) :: GSQRT (KA,IA,JA)
423  real(RP), intent(in) :: J13G (KA,IA,JA)
424  real(RP), intent(in) :: MAPF ( IA,JA,2)
425  real(RP), intent(in) :: CDZ (KA)
426  logical, intent(in) :: TwoD
427  integer, intent(in) :: IIS, IIE, JJS, JJE
428 
429  real(RP) :: vel
430  integer :: k, i, j
431  !---------------------------------------------------------------------------
432 
433  !$omp parallel default(none) private(i,j,k,vel) &
434  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J13G,MAPF)
435 
436  !$acc data copy(flux) copyin(mom, val, DENS, GSQRT, J13G, MAPF, CDZ)
437 
438  !$omp do OMP_SCHEDULE_ collapse(2)
439  !$acc kernels
440  do j = jjs, jje
441  do i = iis, iie
442  do k = ks+2, ke-1
443  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i-1,j) ) ) &
444  / dens(k,i,j)
445  vel = vel * j13g(k,i,j)
446  flux(k-1,i,j) = vel / mapf(i,j,+2) &
447  * ( f2 * ( val(k,i,j)+val(k-1,i,j) ) )
448  enddo
449  enddo
450  enddo
451  !$acc end kernels
452  !$omp end do nowait
453 
454  !$omp do OMP_SCHEDULE_ collapse(2)
455  !$acc kernels
456  do j = jjs, jje
457  do i = iis, iie
458  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS.
459  ! The flux at KS can be non-zero.
460  ! To reduce calculations, all the fluxes are set to zero.
461  flux(ks-1,i,j) = 0.0_rp ! k = KS
462 
463  ! physically incorrect but for numerical stability
464  vel = ( ( 0.5_rp * ( mom(ks+1,i,j)+mom(ks+1,i-1,j) ) ) / dens(ks+1,i,j) &
465  + ( 0.5_rp * ( mom(ks,i,j)+mom(ks,i-1,j) ) ) / dens(ks ,i,j) ) * 0.5_rp
466 ! vel = ( 0.5_RP * ( mom(KS+1,i,j)+mom(KS+1,i-1,j) ) ) &
467 ! / DENS(KS+1,i,j)
468  vel = vel * j13g(ks+1,i,j)
469  flux(ks,i,j) = vel / mapf(i,j,+2) &
470  * ( f2 * ( val(ks+1,i,j)+val(ks,i,j) ) ) ! k = KS+1
471 
472 
473  flux(ke-1,i,j) = 0.0_rp
474  enddo
475  enddo
476  !$acc end kernels
477  !$omp end do nowait
478 
479  !$acc end data
480 
481  !$omp end parallel
482 
483  return

References scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxj23_xyw_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxj23_xyw_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J23-flux at XYW

Definition at line 494 of file scale_atmos_dyn_fvm_flux_cd2.F90.

494  implicit none
495 
496  real(RP), intent(inout) :: flux (KA,IA,JA)
497  real(RP), intent(in) :: mom (KA,IA,JA)
498  real(RP), intent(in) :: val (KA,IA,JA)
499  real(RP), intent(in) :: DENS (KA,IA,JA)
500  real(RP), intent(in) :: GSQRT (KA,IA,JA)
501  real(RP), intent(in) :: J23G (KA,IA,JA)
502  real(RP), intent(in) :: MAPF ( IA,JA,2)
503  real(RP), intent(in) :: CDZ (KA)
504  logical, intent(in) :: TwoD
505  integer, intent(in) :: IIS, IIE, JJS, JJE
506 
507  real(RP) :: vel
508  integer :: k, i, j
509  !---------------------------------------------------------------------------
510 
511  !$omp parallel default(none) private(i,j,k,vel) &
512  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF)
513 
514  !$acc data copy(flux) copyin(mom, val, DENS, GSQRT, J23G, MAPF, CDZ)
515 
516  !$omp do OMP_SCHEDULE_ collapse(2)
517  !$acc kernels
518  do j = jjs, jje
519  do i = iis, iie
520  do k = ks+2, ke-1
521  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j-1) ) ) &
522  / dens(k,i,j)
523  vel = vel * j23g(k,i,j)
524  flux(k-1,i,j) = vel / mapf(i,j,+1) &
525  * ( f2 * ( val(k,i,j)+val(k-1,i,j) ) )
526  enddo
527  enddo
528  enddo
529  !$acc end kernels
530  !$omp end do nowait
531 
532  !$omp do OMP_SCHEDULE_ collapse(2)
533  !$acc kernels
534  do j = jjs, jje
535  do i = iis, iie
536  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS.
537  ! The flux at KS can be non-zero.
538  ! To reduce calculations, all the fluxes are set to zero.
539  flux(ks-1,i,j) = 0.0_rp ! k = KS
540 
541  ! physically incorrect but for numerical stability
542  vel = ( ( 0.5_rp * ( mom(ks+1,i,j)+mom(ks+1,i,j-1) ) ) / dens(ks+1,i,j) &
543  + ( 0.5_rp * ( mom(ks,i,j)+mom(ks,i,j-1) ) ) / dens(ks ,i,j) ) * 0.5_rp
544 ! vel = ( 0.5_RP * ( mom(KS+1,i,j)+mom(KS+1,i,j-1) ) ) &
545 ! / DENS(KS+1,i,j)
546  vel = vel * j23g(ks+1,i,j)
547  flux(ks,i,j) = vel / mapf(i,j,+1) &
548  * ( f2 * ( val(ks+1,i,j)+val(ks,i,j) ) ) ! k = KS+1
549 
550 
551  flux(ke-1,i,j) = 0.0_rp
552  enddo
553  enddo
554  !$acc end kernels
555  !$omp end do nowait
556 
557  !$acc end data
558 
559  !$omp end parallel
560 
561  return

References scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxx_xyw_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxx_xyw_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,ia,ja), intent(in)  num_diff,
real(rp), dimension (ka), intent(in)  CDZ,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation X-flux at XYW

Definition at line 574 of file scale_atmos_dyn_fvm_flux_cd2.F90.

574  implicit none
575 
576  real(RP), intent(inout) :: flux (KA,IA,JA)
577  real(RP), intent(in) :: mom (KA,IA,JA)
578  real(RP), intent(in) :: val (KA,IA,JA)
579  real(RP), intent(in) :: DENS (KA,IA,JA)
580  real(RP), intent(in) :: GSQRT (KA,IA,JA)
581  real(RP), intent(in) :: MAPF ( IA,JA,2)
582  real(RP), intent(in) :: num_diff(KA,IA,JA)
583  real(RP), intent(in) :: CDZ (KA)
584  logical, intent(in) :: TwoD
585  integer, intent(in) :: IIS, IIE, JJS, JJE
586 
587  real(RP) :: vel
588  integer :: k, i, j
589  !---------------------------------------------------------------------------
590 
591  !$omp parallel default(none) private(i,j,k,vel) &
592  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff) &
593  !$omp shared(CDZ)
594 
595  !$acc data copy(flux) copyin(mom, val, DENS, GSQRT, MAPF, num_diff, CDZ)
596 
597  !$omp do OMP_SCHEDULE_ collapse(2)
598  !$acc kernels
599  do j = jjs, jje
600  do i = iis-1, iie
601  do k = ks, ke-1
602 #ifdef DEBUG
603  call check( __line__, mom(k ,i,j) )
604  call check( __line__, mom(k+1,i,j) )
605 
606  call check( __line__, val(k,i,j) )
607  call check( __line__, val(k,i+1,j) )
608 
609 #endif
610  vel = ( f2h(k,1,i_uyz) &
611  * mom(k+1,i,j) &
612  + f2h(k,2,i_uyz) &
613  * mom(k,i,j) ) &
614  / ( f2h(k,1,i_uyz) &
615  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
616  + f2h(k,2,i_uyz) &
617  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
618  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+2) * vel &
619  * ( f2 * ( val(k,i+1,j)+val(k,i,j) ) ) &
620  + gsqrt(k,i,j) * num_diff(k,i,j)
621  enddo
622  enddo
623  enddo
624  !$acc end kernels
625  !$omp end do nowait
626 #ifdef DEBUG
627  k = iundef; i = iundef; j = iundef
628 #endif
629 
630  !$omp do OMP_SCHEDULE_ collapse(2)
631  !$acc kernels
632  do j = jjs, jje
633  do i = iis-1, iie
634  flux(ke,i,j) = 0.0_rp
635  enddo
636  enddo
637  !$acc end kernels
638  !$omp end do nowait
639 
640  !$acc end data
641 
642  !$omp end parallel
643 #ifdef DEBUG
644  k = iundef; i = iundef; j = iundef
645 #endif
646 
647  return

References scale_debug::check(), scale_atmos_grid_cartesc_index::i_uyz, scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxy_xyw_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxy_xyw_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,ia,ja), intent(in)  num_diff,
real(rp), dimension (ka), intent(in)  CDZ,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation Y-flux at XYW

Definition at line 659 of file scale_atmos_dyn_fvm_flux_cd2.F90.

659  implicit none
660 
661  real(RP), intent(inout) :: flux (KA,IA,JA)
662  real(RP), intent(in) :: mom (KA,IA,JA)
663  real(RP), intent(in) :: val (KA,IA,JA)
664  real(RP), intent(in) :: DENS (KA,IA,JA)
665  real(RP), intent(in) :: GSQRT (KA,IA,JA)
666  real(RP), intent(in) :: MAPF ( IA,JA,2)
667  real(RP), intent(in) :: num_diff(KA,IA,JA)
668  real(RP), intent(in) :: CDZ (KA)
669  logical, intent(in) :: TwoD
670  integer, intent(in) :: IIS, IIE, JJS, JJE
671 
672  real(RP) :: vel
673  integer :: k, i, j
674  !---------------------------------------------------------------------------
675 
676  !$omp parallel default(none) private(i,j,k,vel) &
677  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff) &
678  !$omp shared(CDZ)
679 
680  !$acc data copy(flux) copyin(mom, val, DENS, GSQRT, MAPF, num_diff, CDZ)
681 
682  !$omp do OMP_SCHEDULE_ collapse(2)
683  !$acc kernels
684  do j = jjs-1, jje
685  do i = iis, iie
686  do k = ks, ke-1
687 #ifdef DEBUG
688  call check( __line__, mom(k ,i,j) )
689  call check( __line__, mom(k+1,i,j) )
690 
691  call check( __line__, val(k,i,j) )
692  call check( __line__, val(k,i,j+1) )
693 
694 #endif
695  vel = ( f2h(k,1,i_xvz) &
696  * mom(k+1,i,j) &
697  + f2h(k,2,i_xvz) &
698  * mom(k,i,j) ) &
699  / ( f2h(k,1,i_xvz) &
700  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
701  + f2h(k,2,i_xvz) &
702  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
703  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+1) * vel &
704  * ( f2 * ( val(k,i,j+1)+val(k,i,j) ) ) &
705  + gsqrt(k,i,j) * num_diff(k,i,j)
706  enddo
707  enddo
708  enddo
709  !$acc end kernels
710  !$omp end do nowait
711 #ifdef DEBUG
712  k = iundef; i = iundef; j = iundef
713 #endif
714 
715  !$omp do OMP_SCHEDULE_ collapse(2)
716  !$acc kernels
717  do j = jjs-1, jje
718  do i = iis, iie
719  flux(ke,i,j) = 0.0_rp
720  enddo
721  enddo
722  !$acc end kernels
723  !$omp end do nowait
724 
725  !$acc end data
726 
727  !$omp end parallel
728 #ifdef DEBUG
729  k = iundef; i = iundef; j = iundef
730 #endif
731 
732  return

References scale_debug::check(), scale_atmos_grid_cartesc_index::i_xvz, scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxz_uyz_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxz_uyz_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,ia,ja), intent(in)  num_diff,
real(rp), dimension (ka), intent(in)  CDZ,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation z-flux at UY

Definition at line 745 of file scale_atmos_dyn_fvm_flux_cd2.F90.

745  implicit none
746 
747  real(RP), intent(inout) :: flux (KA,IA,JA)
748  real(RP), intent(in) :: mom (KA,IA,JA)
749  real(RP), intent(in) :: val (KA,IA,JA)
750  real(RP), intent(in) :: DENS (KA,IA,JA)
751  real(RP), intent(in) :: GSQRT (KA,IA,JA)
752  real(RP), intent(in) :: J33G
753  real(RP), intent(in) :: num_diff(KA,IA,JA)
754  real(RP), intent(in) :: CDZ (KA)
755  logical, intent(in) :: TwoD
756  integer, intent(in) :: IIS, IIE, JJS, JJE
757 
758  real(RP) :: vel
759  integer :: k, i, j
760  !---------------------------------------------------------------------------
761 
762  !$omp parallel default(none) private(i,j,k,vel) &
763  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J33G,GSQRT,num_diff) &
764  !$omp shared(CDZ,TwoD)
765 
766  !$acc data copy(flux) copyin(mom, val, DENS, GSQRT, num_diff, CDZ)
767 
768 
769  if ( twod ) then
770 
771  !$omp do OMP_SCHEDULE_
772  !$acc kernels
773  do j = jjs, jje
774  do k = ks, ke-1
775  i = iis
776 #ifdef DEBUG
777  call check( __line__, mom(k,i,j) )
778 
779  call check( __line__, val(k,i,j) )
780  call check( __line__, val(k+1,i,j) )
781 
782 #endif
783  vel = ( mom(k,i,j) ) &
784  / ( f2h(k,1,i_xyz) &
785  * dens(k+1,i,j) &
786  + f2h(k,2,i_xyz) &
787  * dens(k,i,j) )
788  flux(k,i,j) = j33g * vel &
789  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) ) &
790  + gsqrt(k,i,j) * num_diff(k,i,j)
791  enddo
792  enddo
793  !$acc end kernels
794  !$omp end do nowait
795 #ifdef DEBUG
796  k = iundef; i = iundef; j = iundef
797 #endif
798 
799  !$omp do OMP_SCHEDULE_
800  !$acc kernels
801  do j = jjs, jje
802  i = iis
803 #ifdef DEBUG
804 
805 #endif
806  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
807  ! The flux at KS-1 can be non-zero.
808  ! To reduce calculations, all the fluxes are set to zero.
809  flux(ks-1,i,j) = 0.0_rp
810 
811  flux(ke,i,j) = 0.0_rp
812  enddo
813  !$acc end kernels
814  !$omp end do nowait
815 
816  else
817 
818 
819  !$omp do OMP_SCHEDULE_ collapse(2)
820  !$acc kernels
821  do j = jjs, jje
822  do i = iis, iie
823  do k = ks, ke-1
824 #ifdef DEBUG
825  call check( __line__, mom(k,i,j) )
826  call check( __line__, mom(k,i+1,j) )
827 
828  call check( __line__, val(k,i,j) )
829  call check( __line__, val(k+1,i,j) )
830 
831 #endif
832  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i+1,j) ) ) &
833  / ( f2h(k,1,i_uyz) &
834  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
835  + f2h(k,2,i_uyz) &
836  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
837  flux(k,i,j) = j33g * vel &
838  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) ) &
839  + gsqrt(k,i,j) * num_diff(k,i,j)
840  enddo
841  enddo
842  enddo
843  !$acc end kernels
844  !$omp end do nowait
845 #ifdef DEBUG
846  k = iundef; i = iundef; j = iundef
847 #endif
848 
849  !$omp do OMP_SCHEDULE_ collapse(2)
850  !$acc kernels
851  do j = jjs, jje
852  do i = iis, iie
853 #ifdef DEBUG
854 
855 #endif
856  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
857  ! The flux at KS-1 can be non-zero.
858  ! To reduce calculations, all the fluxes are set to zero.
859  flux(ks-1,i,j) = 0.0_rp
860 
861  flux(ke,i,j) = 0.0_rp
862  enddo
863  enddo
864  !$acc end kernels
865  !$omp end do nowait
866 
867  end if
868 
869 
870  !$acc end data
871 
872  !$omp end parallel
873 #ifdef DEBUG
874  k = iundef; i = iundef; j = iundef
875 #endif
876 
877  return

References scale_debug::check(), scale_atmos_grid_cartesc_index::i_uyz, scale_atmos_grid_cartesc_index::i_xyz, scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxj13_uyz_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxj13_uyz_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J13-flux at UYZ

Definition at line 888 of file scale_atmos_dyn_fvm_flux_cd2.F90.

888  implicit none
889 
890  real(RP), intent(inout) :: flux (KA,IA,JA)
891  real(RP), intent(in) :: mom (KA,IA,JA)
892  real(RP), intent(in) :: val (KA,IA,JA)
893  real(RP), intent(in) :: DENS (KA,IA,JA)
894  real(RP), intent(in) :: GSQRT (KA,IA,JA)
895  real(RP), intent(in) :: J13G (KA,IA,JA)
896  real(RP), intent(in) :: MAPF ( IA,JA,2)
897  real(RP), intent(in) :: CDZ (KA)
898  logical, intent(in) :: TwoD
899  integer, intent(in) :: IIS, IIE, JJS, JJE
900 
901  real(RP) :: vel
902  integer :: k, i, j
903  !---------------------------------------------------------------------------
904 
905  !$omp parallel default(none) private(i,j,k,vel) &
906  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J13G,MAPF) &
907  !$omp shared(GSQRT,CDZ,TwoD)
908 
909  !$acc data copy(flux) copyin(mom, val, DENS, GSQRT, J13G, MAPF, CDZ)
910 
911 
912 
913  !$omp do OMP_SCHEDULE_ collapse(2)
914  !$acc kernels
915  do j = jjs, jje
916  do i = iis, iie
917  do k = ks, ke-1
918  vel = ( f2h(k,1,i_uyz) &
919  * mom(k+1,i,j) &
920  + f2h(k,2,i_uyz) &
921  * mom(k,i,j) ) &
922  / ( f2h(k,1,i_uyz) &
923  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
924  + f2h(k,2,i_uyz) &
925  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
926  vel = vel * j13g(k,i,j)
927  flux(k,i,j) = vel / mapf(i,j,+2) &
928  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) )
929  enddo
930  enddo
931  enddo
932  !$acc end kernels
933  !$omp end do nowait
934 
935  !$omp do OMP_SCHEDULE_ collapse(2)
936  !$acc kernels
937  do j = jjs, jje
938  do i = iis, iie
939  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
940  ! The flux at KS-1 can be non-zero.
941  ! To reduce calculations, all the fluxes are set to zero.
942  flux(ks-1,i,j) = 0.0_rp
943 
944  flux(ke ,i,j) = 0.0_rp
945  enddo
946  enddo
947  !$acc end kernels
948  !$omp end do nowait
949 
950 
951 
952  !$acc end data
953 
954  !$omp end parallel
955  return

References scale_atmos_grid_cartesc_index::i_uyz, scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxj23_uyz_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxj23_uyz_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J23-flux at UYZ

Definition at line 966 of file scale_atmos_dyn_fvm_flux_cd2.F90.

966  implicit none
967 
968  real(RP), intent(inout) :: flux (KA,IA,JA)
969  real(RP), intent(in) :: mom (KA,IA,JA)
970  real(RP), intent(in) :: val (KA,IA,JA)
971  real(RP), intent(in) :: DENS (KA,IA,JA)
972  real(RP), intent(in) :: GSQRT (KA,IA,JA)
973  real(RP), intent(in) :: J23G (KA,IA,JA)
974  real(RP), intent(in) :: MAPF ( IA,JA,2)
975  real(RP), intent(in) :: CDZ (KA)
976  logical, intent(in) :: TwoD
977  integer, intent(in) :: IIS, IIE, JJS, JJE
978 
979  real(RP) :: vel
980  integer :: k, i, j
981  !---------------------------------------------------------------------------
982 
983  !$omp parallel default(none) private(i,j,k,vel) &
984  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF) &
985  !$omp shared(GSQRT,CDZ,TwoD)
986 
987  !$acc data copy(flux) copyin(mom, val, DENS, GSQRT, J23G, MAPF, CDZ)
988 
989 
990  if ( twod ) then
991 
992  !$omp do OMP_SCHEDULE_
993  !$acc kernels
994  do j = jjs, jje
995  do k = ks, ke-1
996  i = iis
997  vel = ( f2h(k,1,i_xyz) &
998  * 0.5_rp * ( mom(k+1,i,j)+mom(k+1,i,j-1) ) &
999  + f2h(k,2,i_xyz) &
1000  * 0.5_rp * ( mom(k,i,j)+mom(k,i,j-1) ) ) &
1001  / ( f2h(k,1,i_xyz) &
1002  * dens(k+1,i,j) &
1003  + f2h(k,2,i_xyz) &
1004  * dens(k,i,j) )
1005  vel = vel * j23g(k,i,j)
1006  flux(k,i,j) = vel * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) )
1007  enddo
1008  enddo
1009  !$acc end kernels
1010  !$omp end do nowait
1011 
1012  !$omp do OMP_SCHEDULE_
1013  !$acc kernels
1014  do j = jjs, jje
1015  i = iis
1016  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
1017  ! The flux at KS-1 can be non-zero.
1018  ! To reduce calculations, all the fluxes are set to zero.
1019  flux(ks-1,i,j) = 0.0_rp
1020 
1021  flux(ke ,i,j) = 0.0_rp
1022  enddo
1023  !$acc end kernels
1024  !$omp end do nowait
1025 
1026  else
1027 
1028 
1029  !$omp do OMP_SCHEDULE_ collapse(2)
1030  !$acc kernels
1031  do j = jjs, jje
1032  do i = iis, iie
1033  do k = ks, ke-1
1034  vel = ( f2h(k,1,i_uyz) &
1035  * 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) ) &
1036  + f2h(k,2,i_uyz) &
1037  * 0.25_rp * ( mom(k,i,j)+mom(k,i+1,j)+mom(k,i,j-1)+mom(k,i+1,j-1) ) ) &
1038  / ( f2h(k,1,i_uyz) &
1039  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
1040  + f2h(k,2,i_uyz) &
1041  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
1042  vel = vel * j23g(k,i,j)
1043  flux(k,i,j) = vel / mapf(i,j,+1) &
1044  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) )
1045  enddo
1046  enddo
1047  enddo
1048  !$acc end kernels
1049  !$omp end do nowait
1050 
1051  !$omp do OMP_SCHEDULE_ collapse(2)
1052  !$acc kernels
1053  do j = jjs, jje
1054  do i = iis, iie
1055  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
1056  ! The flux at KS-1 can be non-zero.
1057  ! To reduce calculations, all the fluxes are set to zero.
1058  flux(ks-1,i,j) = 0.0_rp
1059 
1060  flux(ke ,i,j) = 0.0_rp
1061  enddo
1062  enddo
1063  !$acc end kernels
1064  !$omp end do nowait
1065 
1066 
1067  end if
1068 
1069 
1070  !$acc end data
1071 
1072  !$omp end parallel
1073  return

References scale_atmos_grid_cartesc_index::i_uyz, scale_atmos_grid_cartesc_index::i_xyz, scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxx_uyz_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxx_uyz_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,ia,ja), intent(in)  num_diff,
real(rp), dimension (ka), intent(in)  CDZ,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation X-flux at UY

Definition at line 1085 of file scale_atmos_dyn_fvm_flux_cd2.F90.

1085  implicit none
1086 
1087  real(RP), intent(inout) :: flux (KA,IA,JA)
1088  real(RP), intent(in) :: mom (KA,IA,JA)
1089  real(RP), intent(in) :: val (KA,IA,JA)
1090  real(RP), intent(in) :: DENS (KA,IA,JA)
1091  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1092  real(RP), intent(in) :: MAPF ( IA,JA,2)
1093  real(RP), intent(in) :: num_diff(KA,IA,JA)
1094  real(RP), intent(in) :: CDZ (KA)
1095  logical, intent(in) :: TwoD
1096  integer, intent(in) :: IIS, IIE, JJS, JJE
1097 
1098  real(RP) :: vel
1099  integer :: k, i, j
1100  !---------------------------------------------------------------------------
1101 
1102  ! note that x-index is added by -1
1103 
1104 
1105 
1106  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1107  !$omp private(vel) &
1108  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff)
1109  !$acc kernels
1110  do j = jjs, jje
1111  do i = iis, iie+1
1112  do k = ks, ke
1113 #ifdef DEBUG
1114  call check( __line__, mom(k,i ,j) )
1115  call check( __line__, mom(k,i-1,j) )
1116 
1117  call check( __line__, val(k,i-1,j) )
1118  call check( __line__, val(k,i,j) )
1119 
1120 #endif
1121  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i-1,j) ) ) &
1122  / ( dens(k,i,j) )
1123  flux(k,i-1,j) = gsqrt(k,i,j) / mapf(i,j,+2) * vel &
1124  * ( f2 * ( val(k,i,j)+val(k,i-1,j) ) ) &
1125  + gsqrt(k,i,j) * num_diff(k,i,j)
1126  enddo
1127  enddo
1128  enddo
1129  !$acc end kernels
1130 #ifdef DEBUG
1131  k = iundef; i = iundef; j = iundef
1132 #endif
1133 
1134 
1135 
1136  return

References scale_debug::check(), scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxy_uyz_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxy_uyz_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,ia,ja), intent(in)  num_diff,
real(rp), dimension (ka), intent(in)  CDZ,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation Y-flux at UY

Definition at line 1148 of file scale_atmos_dyn_fvm_flux_cd2.F90.

1148  implicit none
1149 
1150  real(RP), intent(inout) :: flux (KA,IA,JA)
1151  real(RP), intent(in) :: mom (KA,IA,JA)
1152  real(RP), intent(in) :: val (KA,IA,JA)
1153  real(RP), intent(in) :: DENS (KA,IA,JA)
1154  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1155  real(RP), intent(in) :: MAPF ( IA,JA,2)
1156  real(RP), intent(in) :: num_diff(KA,IA,JA)
1157  real(RP), intent(in) :: CDZ (KA)
1158  logical, intent(in) :: TwoD
1159  integer, intent(in) :: IIS, IIE, JJS, JJE
1160 
1161  real(RP) :: vel
1162  integer :: k, i, j
1163  !---------------------------------------------------------------------------
1164 
1165 
1166 
1167  if ( twod ) then
1168 
1169  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ &
1170  !$omp private(vel) &
1171  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff,TwoD)
1172  !$acc kernels
1173  do j = jjs-1, jje
1174  do k = ks, ke
1175  i = iis
1176 #ifdef DEBUG
1177  call check( __line__, mom(k,i ,j) )
1178 
1179  call check( __line__, val(k,i,j) )
1180  call check( __line__, val(k,i,j+1) )
1181 
1182 #endif
1183  vel = ( mom(k,i,j) ) &
1184  / ( 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
1185  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+1) * vel &
1186  * ( f2 * ( val(k,i,j+1)+val(k,i,j) ) ) &
1187  + gsqrt(k,i,j) * num_diff(k,i,j)
1188  enddo
1189  enddo
1190  !$acc end kernels
1191 #ifdef DEBUG
1192  k = iundef; i = iundef; j = iundef
1193 #endif
1194 
1195  else
1196 
1197 
1198  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1199  !$omp private(vel) &
1200  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff)
1201  !$acc kernels
1202  do j = jjs-1, jje
1203  do i = iis, iie
1204  do k = ks, ke
1205 #ifdef DEBUG
1206  call check( __line__, mom(k,i ,j) )
1207  call check( __line__, mom(k,i-1,j) )
1208 
1209  call check( __line__, val(k,i,j) )
1210  call check( __line__, val(k,i,j+1) )
1211 
1212 #endif
1213  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i+1,j) ) ) &
1214  / ( 0.25_rp * ( dens(k,i,j)+dens(k,i+1,j)+dens(k,i,j+1)+dens(k,i+1,j+1) ) )
1215  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+1) * vel &
1216  * ( f2 * ( val(k,i,j+1)+val(k,i,j) ) ) &
1217  + gsqrt(k,i,j) * num_diff(k,i,j)
1218  enddo
1219  enddo
1220  enddo
1221  !$acc end kernels
1222 #ifdef DEBUG
1223  k = iundef; i = iundef; j = iundef
1224 #endif
1225 
1226 
1227  end if
1228 
1229 
1230  return

References scale_debug::check(), scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxz_xvz_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxz_xvz_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,ia,ja), intent(in)  num_diff,
real(rp), dimension (ka), intent(in)  CDZ,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation z-flux at XV

Definition at line 1244 of file scale_atmos_dyn_fvm_flux_cd2.F90.

1244  implicit none
1245 
1246  real(RP), intent(inout) :: flux (KA,IA,JA)
1247  real(RP), intent(in) :: mom (KA,IA,JA)
1248  real(RP), intent(in) :: val (KA,IA,JA)
1249  real(RP), intent(in) :: DENS (KA,IA,JA)
1250  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1251  real(RP), intent(in) :: J33G
1252  real(RP), intent(in) :: num_diff(KA,IA,JA)
1253  real(RP), intent(in) :: CDZ (KA)
1254  logical, intent(in) :: TwoD
1255  integer, intent(in) :: IIS, IIE, JJS, JJE
1256 
1257  real(RP) :: vel
1258  integer :: k, i, j
1259  !---------------------------------------------------------------------------
1260 
1261  !$omp parallel default(none) private(i,j,k,vel) &
1262  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J33G,GSQRT,num_diff) &
1263  !$omp shared(CDZ,TwoD)
1264 
1265  !$acc data copy(flux) copyin(mom, val, DENS, GSQRT, num_diff, CDZ)
1266 
1267 
1268  !$omp do OMP_SCHEDULE_ collapse(2)
1269  !$acc kernels
1270  do j = jjs, jje
1271  do i = iis, iie
1272  do k = ks, ke-1
1273 #ifdef DEBUG
1274  call check( __line__, mom(k,i,j) )
1275  call check( __line__, mom(k,i,j+1) )
1276 
1277  call check( __line__, val(k,i,j) )
1278  call check( __line__, val(k+1,i,j) )
1279 
1280 #endif
1281  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j+1) ) ) &
1282  / ( f2h(k,1,i_xvz) &
1283  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
1284  + f2h(k,2,i_xvz) &
1285  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
1286  flux(k,i,j) = j33g * vel &
1287  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) ) &
1288  + gsqrt(k,i,j) * num_diff(k,i,j)
1289  enddo
1290  enddo
1291  enddo
1292  !$acc end kernels
1293  !$omp end do nowait
1294 #ifdef DEBUG
1295  k = iundef; i = iundef; j = iundef
1296 #endif
1297 
1298  !$omp do OMP_SCHEDULE_ collapse(2)
1299  !$acc kernels
1300  do j = jjs, jje
1301  do i = iis, iie
1302 #ifdef DEBUG
1303 
1304 #endif
1305  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
1306  ! The flux at KS-1 can be non-zero.
1307  ! To reduce calculations, all the fluxes are set to zero.
1308  flux(ks-1,i,j) = 0.0_rp
1309 
1310  flux(ke,i,j) = 0.0_rp
1311  enddo
1312  enddo
1313  !$acc end kernels
1314  !$omp end do nowait
1315 
1316 
1317  !$acc end data
1318 
1319  !$omp end parallel
1320 #ifdef DEBUG
1321  k = iundef; i = iundef; j = iundef
1322 #endif
1323 
1324  return

References scale_debug::check(), scale_atmos_grid_cartesc_index::i_xvz, scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxj13_xvz_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxj13_xvz_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J13-flux at XVZ

Definition at line 1335 of file scale_atmos_dyn_fvm_flux_cd2.F90.

1335  implicit none
1336 
1337  real(RP), intent(inout) :: flux (KA,IA,JA)
1338  real(RP), intent(in) :: mom (KA,IA,JA)
1339  real(RP), intent(in) :: val (KA,IA,JA)
1340  real(RP), intent(in) :: DENS (KA,IA,JA)
1341  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1342  real(RP), intent(in) :: J13G (KA,IA,JA)
1343  real(RP), intent(in) :: MAPF ( IA,JA,2)
1344  real(RP), intent(in) :: CDZ (KA)
1345  logical, intent(in) :: TwoD
1346  integer, intent(in) :: IIS, IIE, JJS, JJE
1347 
1348  real(RP) :: vel
1349  integer :: k, i, j
1350  !---------------------------------------------------------------------------
1351 
1352  !$omp parallel default(none) private(i,j,k,vel) &
1353  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J13G,MAPF) &
1354  !$omp shared(GSQRT,CDZ,TwoD)
1355 
1356  !$acc data copy(flux) copyin(mom, val, DENS, GSQRT, J13G, MAPF, CDZ)
1357 
1358 
1359 
1360  !$omp do OMP_SCHEDULE_ collapse(2)
1361  !$acc kernels
1362  do j = jjs, jje
1363  do i = iis, iie
1364  do k = ks, ke-1
1365  vel = ( f2h(k,1,i_xvz) &
1366  * 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) ) &
1367  + f2h(k,2,i_xvz) &
1368  * 0.25_rp * ( mom(k,i,j)+mom(k,i-1,j)+mom(k,i,j+1)+mom(k,i-1,j+1) ) ) &
1369  / ( f2h(k,1,i_xvz) &
1370  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
1371  + f2h(k,2,i_xvz) &
1372  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
1373  vel = vel * j13g(k,i,j)
1374  flux(k,i,j) = vel / mapf(i,j,+2) &
1375  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) )
1376  enddo
1377  enddo
1378  enddo
1379  !$acc end kernels
1380  !$omp end do nowait
1381 
1382  !$omp do OMP_SCHEDULE_ collapse(2)
1383  !$acc kernels
1384  do j = jjs, jje
1385  do i = iis, iie
1386  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
1387  ! The flux at KS-1 can be non-zero.
1388  ! To reduce calculations, all the fluxes are set to zero.
1389  flux(ks-1,i,j) = 0.0_rp
1390 
1391  flux(ke ,i,j) = 0.0_rp
1392  enddo
1393  enddo
1394  !$acc end kernels
1395  !$omp end do nowait
1396 
1397 
1398 
1399  !$acc end data
1400 
1401  !$omp end parallel
1402  return

References scale_atmos_grid_cartesc_index::i_xvz, scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxj23_xvz_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxj23_xvz_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J23-flux at XVZ

Definition at line 1413 of file scale_atmos_dyn_fvm_flux_cd2.F90.

1413  implicit none
1414 
1415  real(RP), intent(inout) :: flux (KA,IA,JA)
1416  real(RP), intent(in) :: mom (KA,IA,JA)
1417  real(RP), intent(in) :: val (KA,IA,JA)
1418  real(RP), intent(in) :: DENS (KA,IA,JA)
1419  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1420  real(RP), intent(in) :: J23G (KA,IA,JA)
1421  real(RP), intent(in) :: MAPF ( IA,JA,2)
1422  real(RP), intent(in) :: CDZ (KA)
1423  logical, intent(in) :: TwoD
1424  integer, intent(in) :: IIS, IIE, JJS, JJE
1425 
1426  real(RP) :: vel
1427  integer :: k, i, j
1428  !---------------------------------------------------------------------------
1429 
1430  !$omp parallel default(none) private(i,j,k,vel) &
1431  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF) &
1432  !$omp shared(GSQRT,CDZ,TwoD)
1433 
1434  !$acc data copy(flux) copyin(mom, val, DENS, GSQRT, J23G, MAPF, CDZ)
1435 
1436 
1437 
1438  !$omp do OMP_SCHEDULE_ collapse(2)
1439  !$acc kernels
1440  do j = jjs, jje
1441  do i = iis, iie
1442  do k = ks, ke-1
1443  vel = ( f2h(k,1,i_xvz) &
1444  * mom(k+1,i,j) &
1445  + f2h(k,2,i_xvz) &
1446  * mom(k,i,j) ) &
1447  / ( f2h(k,1,i_xvz) &
1448  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
1449  + f2h(k,2,i_xvz) &
1450  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
1451  vel = vel * j23g(k,i,j)
1452  flux(k,i,j) = vel / mapf(i,j,+1) &
1453  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) )
1454  enddo
1455  enddo
1456  enddo
1457  !$acc end kernels
1458  !$omp end do nowait
1459 
1460  !$omp do OMP_SCHEDULE_ collapse(2)
1461  !$acc kernels
1462  do j = jjs, jje
1463  do i = iis, iie
1464  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
1465  ! The flux at KS-1 can be non-zero.
1466  ! To reduce calculations, all the fluxes are set to zero.
1467  flux(ks-1,i,j) = 0.0_rp
1468 
1469  flux(ke ,i,j) = 0.0_rp
1470  enddo
1471  enddo
1472  !$acc end kernels
1473  !$omp end do nowait
1474 
1475 
1476 
1477  !$acc end data
1478 
1479  !$omp end parallel
1480  return

References scale_atmos_grid_cartesc_index::i_xvz, scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxx_xvz_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxx_xvz_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,ia,ja), intent(in)  num_diff,
real(rp), dimension (ka), intent(in)  CDZ,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation X-flux at XV

Definition at line 1492 of file scale_atmos_dyn_fvm_flux_cd2.F90.

1492  implicit none
1493 
1494  real(RP), intent(inout) :: flux (KA,IA,JA)
1495  real(RP), intent(in) :: mom (KA,IA,JA)
1496  real(RP), intent(in) :: val (KA,IA,JA)
1497  real(RP), intent(in) :: DENS (KA,IA,JA)
1498  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1499  real(RP), intent(in) :: MAPF ( IA,JA,2)
1500  real(RP), intent(in) :: num_diff(KA,IA,JA)
1501  real(RP), intent(in) :: CDZ (KA)
1502  logical, intent(in) :: TwoD
1503  integer, intent(in) :: IIS, IIE, JJS, JJE
1504 
1505  real(RP) :: vel
1506  integer :: k, i, j
1507  !---------------------------------------------------------------------------
1508 
1509 
1510 
1511  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1512  !$omp private(vel) &
1513  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff)
1514  !$acc kernels
1515  do j = jjs, jje
1516  do i = iis-1, iie
1517  do k = ks, ke
1518 #ifdef DEBUG
1519  call check( __line__, mom(k,i ,j) )
1520  call check( __line__, mom(k,i,j-1) )
1521 
1522  call check( __line__, val(k,i,j) )
1523  call check( __line__, val(k,i+1,j) )
1524 
1525 #endif
1526  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j+1) ) ) &
1527  / ( 0.25_rp * ( dens(k,i,j)+dens(k,i+1,j)+dens(k,i,j+1)+dens(k,i+1,j+1) ) )
1528  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+2) * vel &
1529  * ( f2 * ( val(k,i+1,j)+val(k,i,j) ) ) &
1530  + gsqrt(k,i,j) * num_diff(k,i,j)
1531  enddo
1532  enddo
1533  enddo
1534  !$acc end kernels
1535 #ifdef DEBUG
1536  k = iundef; i = iundef; j = iundef
1537 #endif
1538 
1539 
1540 
1541  return

References scale_debug::check(), scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxy_xvz_cd2()

subroutine, public scale_atmos_dyn_fvm_flux_cd2::atmos_dyn_fvm_fluxy_xvz_cd2 ( real(rp), dimension (ka,ia,ja), intent(inout)  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,ia,ja), intent(in)  num_diff,
real(rp), dimension (ka), intent(in)  CDZ,
logical, intent(in)  TwoD,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation Y-flux at XV

Definition at line 1553 of file scale_atmos_dyn_fvm_flux_cd2.F90.

1553  implicit none
1554 
1555  real(RP), intent(inout) :: flux (KA,IA,JA)
1556  real(RP), intent(in) :: mom (KA,IA,JA)
1557  real(RP), intent(in) :: val (KA,IA,JA)
1558  real(RP), intent(in) :: DENS (KA,IA,JA)
1559  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1560  real(RP), intent(in) :: MAPF ( IA,JA,2)
1561  real(RP), intent(in) :: num_diff(KA,IA,JA)
1562  real(RP), intent(in) :: CDZ (KA)
1563  logical, intent(in) :: TwoD
1564  integer, intent(in) :: IIS, IIE, JJS, JJE
1565 
1566  real(RP) :: vel
1567  integer :: k, i, j
1568  !---------------------------------------------------------------------------
1569 
1570  ! note that y-index is added by -1
1571 
1572 
1573 
1574  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1575  !$omp private(vel) &
1576  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff)
1577  !$acc kernels
1578  do j = jjs, jje+1
1579  do i = iis, iie
1580  do k = ks, ke
1581 #ifdef DEBUG
1582  call check( __line__, mom(k,i ,j) )
1583  call check( __line__, mom(k,i,j-1) )
1584 
1585  call check( __line__, val(k,i,j-1) )
1586  call check( __line__, val(k,i,j) )
1587 
1588 #endif
1589  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j-1) ) ) &
1590  / ( dens(k,i,j) )
1591  flux(k,i,j-1) = gsqrt(k,i,j) / mapf(i,j,+1) * vel &
1592  * ( f2 * ( val(k,i,j)+val(k,i,j-1) ) ) &
1593  + gsqrt(k,i,j) * num_diff(k,i,j)
1594  enddo
1595  enddo
1596  enddo
1597  !$acc end kernels
1598 #ifdef DEBUG
1599  k = iundef; i = iundef; j = iundef
1600 #endif
1601 
1602 
1603 
1604  return

References scale_debug::check(), scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

Here is the call graph for this function:
Here is the caller graph for this function:
scale_const::const_eps
real(rp), public const_eps
small number
Definition: scale_const.F90:35
scale_atmos_grid_cartesc_index::i_xyz
integer, public i_xyz
Definition: scale_atmos_grid_cartesC_index.F90:91
scale_const
module CONSTANT
Definition: scale_const.F90:11
Generated by   doxygen 1.8.17