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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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
History
  • 2016-04-18 (S.Nishizawa) [new]

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 101 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

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) = f2 * ( val(k+1)+val(k) )
121  enddo
122 #ifdef DEBUG
123  k = iundef
124 #endif
125 
126 #ifdef DEBUG
127 
128 #endif
129 
130 
131 
132  valw(ks) = f2 * ( val(ks+1)+val(ks) )
133  valw(ke-1) = f2 * ( val(ke)+val(ke-1) )
134 
135 
136  return
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 147 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

147  implicit none
148 
149  real(RP), intent(inout) :: flux (KA,IA,JA)
150  real(RP), intent(in) :: mflx (KA,IA,JA)
151  real(RP), intent(in) :: val (KA,IA,JA)
152  real(RP), intent(in) :: GSQRT (KA,IA,JA)
153  real(RP), intent(in) :: num_diff(KA,IA,JA)
154  real(RP), intent(in) :: CDZ (KA)
155  integer, intent(in) :: IIS, IIE, JJS, JJE
156 
157  real(RP) :: vel
158  integer :: k, i, j
159  !---------------------------------------------------------------------------
160 
161  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
162  !$omp private(vel) &
163  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mflx,val,flux,GSQRT,num_diff)
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 #ifdef DEBUG
182  k = iundef; i = iundef; j = iundef
183 #endif
184 
185  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
186  !$omp private(vel) &
187  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mflx,val,flux,GSQRT,num_diff)
188  do j = jjs, jje
189  do i = iis, iie
190 #ifdef DEBUG
191 
192 #endif
193  flux(ks-1,i,j) = 0.0_rp
194 
195  flux(ke ,i,j) = 0.0_rp
196  enddo
197  enddo
198 #ifdef DEBUG
199  k = iundef; i = iundef; j = iundef
200 #endif
201 
202  return
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 213 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

213  implicit none
214 
215  real(RP), intent(inout) :: flux (KA,IA,JA)
216  real(RP), intent(in) :: mflx (KA,IA,JA)
217  real(RP), intent(in) :: val (KA,IA,JA)
218  real(RP), intent(in) :: GSQRT (KA,IA,JA)
219  real(RP), intent(in) :: num_diff(KA,IA,JA)
220  real(RP), intent(in) :: CDZ(KA)
221  integer, intent(in) :: IIS, IIE, JJS, JJE
222 
223  real(RP) :: vel
224  integer :: k, i, j
225  !---------------------------------------------------------------------------
226 
227  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
228  !$omp private(vel) &
229  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mflx,val,flux,GSQRT,num_diff)
230  do j = jjs, jje
231  do i = iis-1, iie
232  do k = ks, ke
233 #ifdef DEBUG
234  call check( __line__, mflx(k,i,j) )
235 
236  call check( __line__, val(k,i,j) )
237  call check( __line__, val(k,i+1,j) )
238 
239 #endif
240  vel = mflx(k,i,j)
241  flux(k,i,j) = vel &
242  * ( f2 * ( val(k,i+1,j)+val(k,i,j) ) ) &
243  + gsqrt(k,i,j) * num_diff(k,i,j)
244  enddo
245  enddo
246  enddo
247 #ifdef DEBUG
248  k = iundef; i = iundef; j = iundef
249 #endif
250 
251  return
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 262 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

262  implicit none
263 
264  real(RP), intent(inout) :: flux (KA,IA,JA)
265  real(RP), intent(in) :: mflx (KA,IA,JA)
266  real(RP), intent(in) :: val (KA,IA,JA)
267  real(RP), intent(in) :: GSQRT (KA,IA,JA)
268  real(RP), intent(in) :: num_diff(KA,IA,JA)
269  real(RP), intent(in) :: CDZ(KA)
270  integer, intent(in) :: IIS, IIE, JJS, JJE
271 
272  real(RP) :: vel
273  integer :: k, i, j
274  !---------------------------------------------------------------------------
275 
276  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
277  !$omp private(vel) &
278  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mflx,val,flux,GSQRT,num_diff)
279  do j = jjs-1, jje
280  do i = iis, iie
281  do k = ks, ke
282 #ifdef DEBUG
283  call check( __line__, mflx(k,i,j) )
284 
285  call check( __line__, val(k,i,j) )
286  call check( __line__, val(k,i,j+1) )
287 
288 #endif
289  vel = mflx(k,i,j)
290  flux(k,i,j) = vel &
291  * ( f2 * ( val(k,i,j+1)+val(k,i,j) ) ) &
292  + gsqrt(k,i,j) * num_diff(k,i,j)
293  enddo
294  enddo
295  enddo
296 #ifdef DEBUG
297  k = iundef; i = iundef; j = iundef
298 #endif
299 
300  return
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 314 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

314  implicit none
315 
316  real(RP), intent(inout) :: flux (KA,IA,JA)
317  real(RP), intent(in) :: mom (KA,IA,JA)
318  real(RP), intent(in) :: val (KA,IA,JA)
319  real(RP), intent(in) :: DENS (KA,IA,JA)
320  real(RP), intent(in) :: GSQRT (KA,IA,JA)
321  real(RP), intent(in) :: J33G
322  real(RP), intent(in) :: num_diff(KA,IA,JA)
323  real(RP), intent(in) :: CDZ (KA)
324  real(RP), intent(in) :: FDZ (KA-1)
325  real(RP), intent(in) :: dtrk
326  integer, intent(in) :: IIS, IIE, JJS, JJE
327 
328  real(RP) :: vel
329  integer :: k, i, j
330  !---------------------------------------------------------------------------
331 
332  ! note than z-index is added by -1
333 
334  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
335  !$omp private(vel) &
336  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,flux,J33G,GSQRT,num_diff,DENS)
337  do j = jjs, jje
338  do i = iis, iie
339  do k = ks+1, ke-1
340 #ifdef DEBUG
341  call check( __line__, mom(k-1,i,j) )
342  call check( __line__, mom(k ,i,j) )
343 
344  call check( __line__, val(k-1,i,j) )
345  call check( __line__, val(k,i,j) )
346 
347 #endif
348  vel = ( 0.5_rp * ( mom(k-1,i,j) &
349  + mom(k,i,j) ) ) &
350  / dens(k,i,j)
351  flux(k-1,i,j) = j33g * vel &
352  * ( f2 * ( val(k,i,j)+val(k-1,i,j) ) ) &
353  + gsqrt(k,i,j) * num_diff(k,i,j)
354  enddo
355  enddo
356  enddo
357 #ifdef DEBUG
358  k = iundef; i = iundef; j = iundef
359 #endif
360 
361  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
362  !$omp private(vel) &
363  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J33G,GSQRT,num_diff,FDZ,dtrk)
364  do j = jjs, jje
365  do i = iis, iie
366 #ifdef DEBUG
367 
368 
369 #endif
370  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS.
371  ! The flux at KS can be non-zero.
372  ! To reduce calculations, all the fluxes are set to zero.
373  flux(ks-1,i,j) = 0.0_rp ! k = KS
374 
375 
376 
377  flux(ke-1,i,j) = 0.0_rp ! k = KE
378  flux(ke ,i,j) = 0.0_rp ! k = KE+1
379  enddo
380  enddo
381 
382  return
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J13-flux at XYW

Definition at line 394 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

394  implicit none
395 
396  real(RP), intent(inout) :: flux (KA,IA,JA)
397  real(RP), intent(in) :: mom (KA,IA,JA)
398  real(RP), intent(in) :: val (KA,IA,JA)
399  real(RP), intent(in) :: DENS (KA,IA,JA)
400  real(RP), intent(in) :: GSQRT (KA,IA,JA)
401  real(RP), intent(in) :: J13G (KA,IA,JA)
402  real(RP), intent(in) :: MAPF ( IA,JA,2)
403  real(RP), intent(in) :: CDZ (KA)
404  integer, intent(in) :: IIS, IIE, JJS, JJE
405 
406  real(RP) :: vel
407  integer :: k, i, j
408  !---------------------------------------------------------------------------
409 
410  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
411  !$omp private(vel) &
412  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J13G,MAPF)
413  do j = jjs, jje
414  do i = iis, iie
415  do k = ks+2, ke-1
416  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i-1,j) ) ) &
417  / dens(k,i,j)
418  flux(k-1,i,j) = j13g(k,i,j) / mapf(i,j,+2) * vel &
419  * ( f2 * ( val(k,i,j)+val(k-1,i,j) ) )
420  enddo
421  enddo
422  enddo
423 
424  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
425  !$omp private(vel) &
426  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J13G,MAPF)
427  do j = jjs, jje
428  do i = iis, iie
429  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS.
430  ! The flux at KS can be non-zero.
431  ! To reduce calculations, all the fluxes are set to zero.
432  flux(ks-1,i,j) = 0.0_rp ! k = KS
433 
434  ! physically incorrect but for numerical stability
435  vel = ( ( 0.5_rp * ( mom(ks+1,i,j)+mom(ks+1,i-1,j) ) ) / dens(ks+1,i,j) &
436  + ( 0.5_rp * ( mom(ks,i,j)+mom(ks,i-1,j) ) ) / dens(ks ,i,j) ) * 0.5_rp
437 ! vel = ( 0.5_RP * ( mom(KS+1,i,j)+mom(KS+1,i-1,j) ) ) &
438 ! / DENS(KS+1,i,j)
439  flux(ks,i,j) = j13g(ks+1,i,j) / mapf(i,j,+2) * vel &
440  * ( f2 * ( val(ks+1,i,j)+val(ks,i,j) ) ) ! k = KS+1
441 
442 
443  flux(ke-1,i,j) = 0.0_rp
444  enddo
445  enddo
446 
447  return
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J23-flux at XYW

Definition at line 458 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

458  implicit none
459 
460  real(RP), intent(inout) :: flux (KA,IA,JA)
461  real(RP), intent(in) :: mom (KA,IA,JA)
462  real(RP), intent(in) :: val (KA,IA,JA)
463  real(RP), intent(in) :: DENS (KA,IA,JA)
464  real(RP), intent(in) :: GSQRT (KA,IA,JA)
465  real(RP), intent(in) :: J23G (KA,IA,JA)
466  real(RP), intent(in) :: MAPF ( IA,JA,2)
467  real(RP), intent(in) :: CDZ (KA)
468  integer, intent(in) :: IIS, IIE, JJS, JJE
469 
470  real(RP) :: vel
471  integer :: k, i, j
472  !---------------------------------------------------------------------------
473 
474  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
475  !$omp private(vel) &
476  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF)
477  do j = jjs, jje
478  do i = iis, iie
479  do k = ks+2, ke-1
480  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j-1) ) ) &
481  / dens(k,i,j)
482  flux(k-1,i,j) = j23g(k,i,j) / mapf(i,j,+1) * vel &
483  * ( f2 * ( val(k,i,j)+val(k-1,i,j) ) )
484  enddo
485  enddo
486  enddo
487 
488  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
489  !$omp private(vel) &
490  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF)
491  do j = jjs, jje
492  do i = iis, iie
493  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS.
494  ! The flux at KS can be non-zero.
495  ! To reduce calculations, all the fluxes are set to zero.
496  flux(ks-1,i,j) = 0.0_rp ! k = KS
497 
498  ! physically incorrect but for numerical stability
499  vel = ( ( 0.5_rp * ( mom(ks+1,i,j)+mom(ks+1,i,j-1) ) ) / dens(ks+1,i,j) &
500  + ( 0.5_rp * ( mom(ks,i,j)+mom(ks,i,j-1) ) ) / dens(ks ,i,j) ) * 0.5_rp
501 ! vel = ( 0.5_RP * ( mom(KS+1,i,j)+mom(KS+1,i,j-1) ) ) &
502 ! / DENS(KS+1,i,j)
503  flux(ks,i,j) = j23g(ks+1,i,j) / mapf(i,j,+1) * vel &
504  * ( f2 * ( val(ks+1,i,j)+val(ks,i,j) ) ) ! k = KS+1
505 
506 
507  flux(ke-1,i,j) = 0.0_rp
508  enddo
509  enddo
510 
511  return
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation X-flux at XYW

Definition at line 524 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

524  implicit none
525 
526  real(RP), intent(inout) :: flux (KA,IA,JA)
527  real(RP), intent(in) :: mom (KA,IA,JA)
528  real(RP), intent(in) :: val (KA,IA,JA)
529  real(RP), intent(in) :: DENS (KA,IA,JA)
530  real(RP), intent(in) :: GSQRT (KA,IA,JA)
531  real(RP), intent(in) :: MAPF ( IA,JA,2)
532  real(RP), intent(in) :: num_diff(KA,IA,JA)
533  real(RP), intent(in) :: CDZ (KA)
534  integer, intent(in) :: IIS, IIE, JJS, JJE
535 
536  real(RP) :: vel
537  integer :: k, i, j
538  !---------------------------------------------------------------------------
539 
540  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
541  !$omp private(vel) &
542  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff) &
543  !$omp shared(CDZ)
544  do j = jjs, jje
545  do i = iis-1, iie
546  do k = ks, ke-1
547 #ifdef DEBUG
548  call check( __line__, mom(k ,i,j) )
549  call check( __line__, mom(k+1,i,j) )
550 
551  call check( __line__, val(k,i,j) )
552  call check( __line__, val(k,i+1,j) )
553 
554 #endif
555  vel = ( f2h(k,1,i_uyz) &
556  * mom(k+1,i,j) &
557  + f2h(k,2,i_uyz) &
558  * mom(k,i,j) ) &
559  / ( f2h(k,1,i_uyz) &
560  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
561  + f2h(k,2,i_uyz) &
562  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
563  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+2) * vel &
564  * ( f2 * ( val(k,i+1,j)+val(k,i,j) ) ) &
565  + gsqrt(k,i,j) * num_diff(k,i,j)
566  enddo
567  enddo
568  enddo
569 #ifdef DEBUG
570  k = iundef; i = iundef; j = iundef
571 #endif
572 
573  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
574  !$omp private(vel) &
575  !$omp shared(JJS,JJE,IIS,IIE,KE,flux)
576  do j = jjs, jje
577  do i = iis-1, iie
578  flux(ke,i,j) = 0.0_rp
579  enddo
580  enddo
581 #ifdef DEBUG
582  k = iundef; i = iundef; j = iundef
583 #endif
584 
585  return
scale_gridtrans::i_uyz
integer, public i_uyz
Definition: scale_gridtrans.F90:53
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation Y-flux at XYW

Definition at line 597 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

597  implicit none
598 
599  real(RP), intent(inout) :: flux (KA,IA,JA)
600  real(RP), intent(in) :: mom (KA,IA,JA)
601  real(RP), intent(in) :: val (KA,IA,JA)
602  real(RP), intent(in) :: DENS (KA,IA,JA)
603  real(RP), intent(in) :: GSQRT (KA,IA,JA)
604  real(RP), intent(in) :: MAPF ( IA,JA,2)
605  real(RP), intent(in) :: num_diff(KA,IA,JA)
606  real(RP), intent(in) :: CDZ (KA)
607  integer, intent(in) :: IIS, IIE, JJS, JJE
608 
609  real(RP) :: vel
610  integer :: k, i, j
611  !---------------------------------------------------------------------------
612 
613  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
614  !$omp private(vel) &
615  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff) &
616  !$omp shared(CDZ)
617  do j = jjs-1, jje
618  do i = iis, iie
619  do k = ks, ke-1
620 #ifdef DEBUG
621  call check( __line__, mom(k ,i,j) )
622  call check( __line__, mom(k+1,i,j) )
623 
624  call check( __line__, val(k,i,j) )
625  call check( __line__, val(k,i,j+1) )
626 
627 #endif
628  vel = ( f2h(k,1,i_xvz) &
629  * mom(k+1,i,j) &
630  + f2h(k,2,i_xvz) &
631  * mom(k,i,j) ) &
632  / ( f2h(k,1,i_xvz) &
633  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
634  + f2h(k,2,i_xvz) &
635  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
636  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+1) * vel &
637  * ( f2 * ( val(k,i,j+1)+val(k,i,j) ) ) &
638  + gsqrt(k,i,j) * num_diff(k,i,j)
639  enddo
640  enddo
641  enddo
642 #ifdef DEBUG
643  k = iundef; i = iundef; j = iundef
644 #endif
645 
646  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
647  !$omp private(vel) &
648  !$omp shared(JJS,JJE,IIS,IIE,KE,flux)
649  do j = jjs-1, jje
650  do i = iis, iie
651  flux(ke,i,j) = 0.0_rp
652  enddo
653  enddo
654 #ifdef DEBUG
655  k = iundef; i = iundef; j = iundef
656 #endif
657 
658  return
scale_gridtrans::i_xvz
integer, public i_xvz
Definition: scale_gridtrans.F90:54
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation z-flux at UY

Definition at line 671 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

671  implicit none
672 
673  real(RP), intent(inout) :: flux (KA,IA,JA)
674  real(RP), intent(in) :: mom (KA,IA,JA)
675  real(RP), intent(in) :: val (KA,IA,JA)
676  real(RP), intent(in) :: DENS (KA,IA,JA)
677  real(RP), intent(in) :: GSQRT (KA,IA,JA)
678  real(RP), intent(in) :: J33G
679  real(RP), intent(in) :: num_diff(KA,IA,JA)
680  real(RP), intent(in) :: CDZ (KA)
681  integer, intent(in) :: IIS, IIE, JJS, JJE
682 
683  real(RP) :: vel
684  integer :: k, i, j
685  !---------------------------------------------------------------------------
686 
687  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
688  !$omp private(vel) &
689  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J33G,GSQRT,num_diff) &
690  !$omp shared(CDZ)
691  do j = jjs, jje
692  do i = iis, iie
693  do k = ks, ke-1
694 #ifdef DEBUG
695  call check( __line__, mom(k,i,j) )
696  call check( __line__, mom(k,i+1,j) )
697 
698  call check( __line__, val(k,i,j) )
699  call check( __line__, val(k+1,i,j) )
700 
701 #endif
702  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i+1,j) ) ) &
703  / ( f2h(k,1,i_uyz) &
704  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
705  + f2h(k,2,i_uyz) &
706  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
707  flux(k,i,j) = j33g * vel &
708  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) ) &
709  + gsqrt(k,i,j) * num_diff(k,i,j)
710  enddo
711  enddo
712  enddo
713 #ifdef DEBUG
714  k = iundef; i = iundef; j = iundef
715 #endif
716 
717  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
718  !$omp private(vel) &
719  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,flux,J33G,GSQRT,num_diff,DENS,CDZ)
720  do j = jjs, jje
721  do i = iis, iie
722 #ifdef DEBUG
723 
724 #endif
725  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
726  ! The flux at KS-1 can be non-zero.
727  ! To reduce calculations, all the fluxes are set to zero.
728  flux(ks-1,i,j) = 0.0_rp
729 
730  flux(ke,i,j) = 0.0_rp
731  enddo
732  enddo
733 #ifdef DEBUG
734  k = iundef; i = iundef; j = iundef
735 #endif
736 
737  return
scale_gridtrans::i_uyz
integer, public i_uyz
Definition: scale_gridtrans.F90:53
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J13-flux at UYZ

Definition at line 748 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

748  implicit none
749 
750  real(RP), intent(inout) :: flux (KA,IA,JA)
751  real(RP), intent(in) :: mom (KA,IA,JA)
752  real(RP), intent(in) :: val (KA,IA,JA)
753  real(RP), intent(in) :: DENS (KA,IA,JA)
754  real(RP), intent(in) :: GSQRT (KA,IA,JA)
755  real(RP), intent(in) :: J13G (KA,IA,JA)
756  real(RP), intent(in) :: MAPF ( IA,JA,2)
757  real(RP), intent(in) :: CDZ (KA)
758  integer, intent(in) :: IIS, IIE, JJS, JJE
759 
760  real(RP) :: vel
761  integer :: k, i, j
762  !---------------------------------------------------------------------------
763 
764  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
765  !$omp private(vel) &
766  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J13G,MAPF) &
767  !$omp shared(GSQRT,CDZ)
768  do j = jjs, jje
769  do i = iis, iie
770  do k = ks, ke-1
771  vel = ( f2h(k,1,i_uyz) &
772  * mom(k+1,i,j) &
773  + f2h(k,2,i_uyz) &
774  * mom(k,i,j) ) &
775  / ( f2h(k,1,i_uyz) &
776  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
777  + f2h(k,2,i_uyz) &
778  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
779  flux(k,i,j) = j13g(k,i,j) / mapf(i,j,+2) * vel &
780  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) )
781  enddo
782  enddo
783  enddo
784 
785  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
786  !$omp private(vel) &
787  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J13G,MAPF) &
788  !$omp shared(GSQRT,CDZ)
789  do j = jjs, jje
790  do i = iis, iie
791  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
792  ! The flux at KS-1 can be non-zero.
793  ! To reduce calculations, all the fluxes are set to zero.
794  flux(ks-1,i,j) = 0.0_rp
795 
796  flux(ke ,i,j) = 0.0_rp
797  enddo
798  enddo
799 
800  return
scale_gridtrans::i_uyz
integer, public i_uyz
Definition: scale_gridtrans.F90:53
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J23-flux at UYZ

Definition at line 811 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

811  implicit none
812 
813  real(RP), intent(inout) :: flux (KA,IA,JA)
814  real(RP), intent(in) :: mom (KA,IA,JA)
815  real(RP), intent(in) :: val (KA,IA,JA)
816  real(RP), intent(in) :: DENS (KA,IA,JA)
817  real(RP), intent(in) :: GSQRT (KA,IA,JA)
818  real(RP), intent(in) :: J23G (KA,IA,JA)
819  real(RP), intent(in) :: MAPF ( IA,JA,2)
820  real(RP), intent(in) :: CDZ (KA)
821  integer, intent(in) :: IIS, IIE, JJS, JJE
822 
823  real(RP) :: vel
824  integer :: k, i, j
825  !---------------------------------------------------------------------------
826 
827  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
828  !$omp private(vel) &
829  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF) &
830  !$omp shared(GSQRT,CDZ)
831  do j = jjs, jje
832  do i = iis, iie
833  do k = ks, ke-1
834  vel = ( f2h(k,1,i_uyz) &
835  * 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) ) &
836  + f2h(k,2,i_uyz) &
837  * 0.25_rp * ( mom(k,i,j)+mom(k,i+1,j)+mom(k,i,j-1)+mom(k,i+1,j-1) ) ) &
838  / ( f2h(k,1,i_uyz) &
839  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
840  + f2h(k,2,i_uyz) &
841  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
842  flux(k,i,j) = j23g(k,i,j) / mapf(i,j,+1) * vel &
843  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) )
844  enddo
845  enddo
846  enddo
847 
848  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
849  !$omp private(vel) &
850  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF) &
851  !$omp shared(GSQRT,CDZ)
852  do j = jjs, jje
853  do i = iis, iie
854  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
855  ! The flux at KS-1 can be non-zero.
856  ! To reduce calculations, all the fluxes are set to zero.
857  flux(ks-1,i,j) = 0.0_rp
858 
859  flux(ke ,i,j) = 0.0_rp
860  enddo
861  enddo
862 
863  return
scale_gridtrans::i_uyz
integer, public i_uyz
Definition: scale_gridtrans.F90:53
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation X-flux at UY

Definition at line 875 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

875  implicit none
876 
877  real(RP), intent(inout) :: flux (KA,IA,JA)
878  real(RP), intent(in) :: mom (KA,IA,JA)
879  real(RP), intent(in) :: val (KA,IA,JA)
880  real(RP), intent(in) :: DENS (KA,IA,JA)
881  real(RP), intent(in) :: GSQRT (KA,IA,JA)
882  real(RP), intent(in) :: MAPF ( IA,JA,2)
883  real(RP), intent(in) :: num_diff(KA,IA,JA)
884  real(RP), intent(in) :: CDZ (KA)
885  integer, intent(in) :: IIS, IIE, JJS, JJE
886 
887  real(RP) :: vel
888  integer :: k, i, j
889  !---------------------------------------------------------------------------
890 
891  ! note that x-index is added by -1
892 
893  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
894  !$omp private(vel) &
895  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff)
896  do j = jjs, jje
897  do i = iis, iie+1
898  do k = ks, ke
899 #ifdef DEBUG
900  call check( __line__, mom(k,i ,j) )
901  call check( __line__, mom(k,i-1,j) )
902 
903  call check( __line__, val(k,i-1,j) )
904  call check( __line__, val(k,i,j) )
905 
906 #endif
907  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i-1,j) ) ) &
908  / ( dens(k,i,j) )
909  flux(k,i-1,j) = gsqrt(k,i,j) / mapf(i,j,+2) * vel &
910  * ( f2 * ( val(k,i,j)+val(k,i-1,j) ) ) &
911  + gsqrt(k,i,j) * num_diff(k,i,j)
912  enddo
913  enddo
914  enddo
915 #ifdef DEBUG
916  k = iundef; i = iundef; j = iundef
917 #endif
918 
919  return
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation Y-flux at UY

Definition at line 931 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

931  implicit none
932 
933  real(RP), intent(inout) :: flux (KA,IA,JA)
934  real(RP), intent(in) :: mom (KA,IA,JA)
935  real(RP), intent(in) :: val (KA,IA,JA)
936  real(RP), intent(in) :: DENS (KA,IA,JA)
937  real(RP), intent(in) :: GSQRT (KA,IA,JA)
938  real(RP), intent(in) :: MAPF ( IA,JA,2)
939  real(RP), intent(in) :: num_diff(KA,IA,JA)
940  real(RP), intent(in) :: CDZ (KA)
941  integer, intent(in) :: IIS, IIE, JJS, JJE
942 
943  real(RP) :: vel
944  integer :: k, i, j
945  !---------------------------------------------------------------------------
946 
947  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
948  !$omp private(vel) &
949  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff)
950  do j = jjs-1, jje
951  do i = iis, iie
952  do k = ks, ke
953 #ifdef DEBUG
954  call check( __line__, mom(k,i ,j) )
955  call check( __line__, mom(k,i-1,j) )
956 
957  call check( __line__, val(k,i,j) )
958  call check( __line__, val(k,i,j+1) )
959 
960 #endif
961  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i+1,j) ) ) &
962  / ( 0.25_rp * ( dens(k,i,j)+dens(k,i+1,j)+dens(k,i,j+1)+dens(k,i+1,j+1) ) )
963  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+1) * vel &
964  * ( f2 * ( val(k,i,j+1)+val(k,i,j) ) ) &
965  + gsqrt(k,i,j) * num_diff(k,i,j)
966  enddo
967  enddo
968  enddo
969 #ifdef DEBUG
970  k = iundef; i = iundef; j = iundef
971 #endif
972 
973  return
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation z-flux at XV

Definition at line 987 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

987  implicit none
988 
989  real(RP), intent(inout) :: flux (KA,IA,JA)
990  real(RP), intent(in) :: mom (KA,IA,JA)
991  real(RP), intent(in) :: val (KA,IA,JA)
992  real(RP), intent(in) :: DENS (KA,IA,JA)
993  real(RP), intent(in) :: GSQRT (KA,IA,JA)
994  real(RP), intent(in) :: J33G
995  real(RP), intent(in) :: num_diff(KA,IA,JA)
996  real(RP), intent(in) :: CDZ (KA)
997  integer, intent(in) :: IIS, IIE, JJS, JJE
998 
999  real(RP) :: vel
1000  integer :: k, i, j
1001  !---------------------------------------------------------------------------
1002 
1003  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1004  !$omp private(vel) &
1005  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J33G,GSQRT,num_diff) &
1006  !$omp shared(CDZ)
1007  do j = jjs, jje
1008  do i = iis, iie
1009  do k = ks, ke-1
1010 #ifdef DEBUG
1011  call check( __line__, mom(k,i,j) )
1012  call check( __line__, mom(k,i,j+1) )
1013 
1014  call check( __line__, val(k,i,j) )
1015  call check( __line__, val(k+1,i,j) )
1016 
1017 #endif
1018  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j+1) ) ) &
1019  / ( f2h(k,1,i_xvz) &
1020  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
1021  + f2h(k,2,i_xvz) &
1022  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
1023  flux(k,i,j) = j33g * vel &
1024  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) ) &
1025  + gsqrt(k,i,j) * num_diff(k,i,j)
1026  enddo
1027  enddo
1028  enddo
1029 #ifdef DEBUG
1030  k = iundef; i = iundef; j = iundef
1031 #endif
1032 
1033  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
1034  !$omp private(vel) &
1035  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,flux,J33G,GSQRT,num_diff,DENS,CDZ)
1036  do j = jjs, jje
1037  do i = iis, iie
1038 #ifdef DEBUG
1039 
1040 #endif
1041  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
1042  ! The flux at KS-1 can be non-zero.
1043  ! To reduce calculations, all the fluxes are set to zero.
1044  flux(ks-1,i,j) = 0.0_rp
1045 
1046  flux(ke,i,j) = 0.0_rp
1047  enddo
1048  enddo
1049 #ifdef DEBUG
1050  k = iundef; i = iundef; j = iundef
1051 #endif
1052 
1053  return
scale_gridtrans::i_xvz
integer, public i_xvz
Definition: scale_gridtrans.F90:54
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J13-flux at XVZ

Definition at line 1064 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

1064  implicit none
1065 
1066  real(RP), intent(inout) :: flux (KA,IA,JA)
1067  real(RP), intent(in) :: mom (KA,IA,JA)
1068  real(RP), intent(in) :: val (KA,IA,JA)
1069  real(RP), intent(in) :: DENS (KA,IA,JA)
1070  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1071  real(RP), intent(in) :: J13G (KA,IA,JA)
1072  real(RP), intent(in) :: MAPF ( IA,JA,2)
1073  real(RP), intent(in) :: CDZ (KA)
1074  integer, intent(in) :: IIS, IIE, JJS, JJE
1075 
1076  real(RP) :: vel
1077  integer :: k, i, j
1078  !---------------------------------------------------------------------------
1079 
1080  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1081  !$omp private(vel) &
1082  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J13G,MAPF) &
1083  !$omp shared(GSQRT,CDZ)
1084  do j = jjs, jje
1085  do i = iis, iie
1086  do k = ks, ke-1
1087  vel = ( f2h(k,1,i_xvz) &
1088  * 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) ) &
1089  + f2h(k,2,i_xvz) &
1090  * 0.25_rp * ( mom(k,i,j)+mom(k,i-1,j)+mom(k,i,j+1)+mom(k,i-1,j+1) ) ) &
1091  / ( f2h(k,1,i_xvz) &
1092  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
1093  + f2h(k,2,i_xvz) &
1094  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
1095  flux(k,i,j) = j13g(k,i,j) / mapf(i,j,+2) * vel &
1096  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) )
1097  enddo
1098  enddo
1099  enddo
1100 
1101  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
1102  !$omp private(vel) &
1103  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J13G,MAPF) &
1104  !$omp shared(GSQRT,CDZ)
1105  do j = jjs, jje
1106  do i = iis, iie
1107  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
1108  ! The flux at KS-1 can be non-zero.
1109  ! To reduce calculations, all the fluxes are set to zero.
1110  flux(ks-1,i,j) = 0.0_rp
1111 
1112  flux(ke ,i,j) = 0.0_rp
1113  enddo
1114  enddo
1115 
1116  return
scale_gridtrans::i_xvz
integer, public i_xvz
Definition: scale_gridtrans.F90:54
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation J23-flux at XVZ

Definition at line 1127 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

1127  implicit none
1128 
1129  real(RP), intent(inout) :: flux (KA,IA,JA)
1130  real(RP), intent(in) :: mom (KA,IA,JA)
1131  real(RP), intent(in) :: val (KA,IA,JA)
1132  real(RP), intent(in) :: DENS (KA,IA,JA)
1133  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1134  real(RP), intent(in) :: J23G (KA,IA,JA)
1135  real(RP), intent(in) :: MAPF ( IA,JA,2)
1136  real(RP), intent(in) :: CDZ (KA)
1137  integer, intent(in) :: IIS, IIE, JJS, JJE
1138 
1139  real(RP) :: vel
1140  integer :: k, i, j
1141  !---------------------------------------------------------------------------
1142 
1143  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1144  !$omp private(vel) &
1145  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF) &
1146  !$omp shared(GSQRT,CDZ)
1147  do j = jjs, jje
1148  do i = iis, iie
1149  do k = ks, ke-1
1150  vel = ( f2h(k,1,i_xvz) &
1151  * mom(k+1,i,j) &
1152  + f2h(k,2,i_xvz) &
1153  * mom(k,i,j) ) &
1154  / ( f2h(k,1,i_xvz) &
1155  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
1156  + f2h(k,2,i_xvz) &
1157  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
1158  flux(k,i,j) = j23g(k,i,j) / mapf(i,j,+1) * vel &
1159  * ( f2 * ( val(k+1,i,j)+val(k,i,j) ) )
1160  enddo
1161  enddo
1162  enddo
1163 
1164  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
1165  !$omp private(vel) &
1166  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF) &
1167  !$omp shared(GSQRT,CDZ)
1168  do j = jjs, jje
1169  do i = iis, iie
1170  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
1171  ! The flux at KS-1 can be non-zero.
1172  ! To reduce calculations, all the fluxes are set to zero.
1173  flux(ks-1,i,j) = 0.0_rp
1174 
1175  flux(ke ,i,j) = 0.0_rp
1176  enddo
1177  enddo
1178 
1179  return
scale_gridtrans::i_xvz
integer, public i_xvz
Definition: scale_gridtrans.F90:54
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation X-flux at XV

Definition at line 1191 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

1191  implicit none
1192 
1193  real(RP), intent(inout) :: flux (KA,IA,JA)
1194  real(RP), intent(in) :: mom (KA,IA,JA)
1195  real(RP), intent(in) :: val (KA,IA,JA)
1196  real(RP), intent(in) :: DENS (KA,IA,JA)
1197  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1198  real(RP), intent(in) :: MAPF ( IA,JA,2)
1199  real(RP), intent(in) :: num_diff(KA,IA,JA)
1200  real(RP), intent(in) :: CDZ (KA)
1201  integer, intent(in) :: IIS, IIE, JJS, JJE
1202 
1203  real(RP) :: vel
1204  integer :: k, i, j
1205  !---------------------------------------------------------------------------
1206 
1207  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1208  !$omp private(vel) &
1209  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff)
1210  do j = jjs, jje
1211  do i = iis-1, iie
1212  do k = ks, ke
1213 #ifdef DEBUG
1214  call check( __line__, mom(k,i ,j) )
1215  call check( __line__, mom(k,i,j-1) )
1216 
1217  call check( __line__, val(k,i,j) )
1218  call check( __line__, val(k,i+1,j) )
1219 
1220 #endif
1221  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j+1) ) ) &
1222  / ( 0.25_rp * ( dens(k,i,j)+dens(k,i+1,j)+dens(k,i,j+1)+dens(k,i+1,j+1) ) )
1223  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+2) * vel &
1224  * ( f2 * ( val(k,i+1,j)+val(k,i,j) ) ) &
1225  + gsqrt(k,i,j) * num_diff(k,i,j)
1226  enddo
1227  enddo
1228  enddo
1229 #ifdef DEBUG
1230  k = iundef; i = iundef; j = iundef
1231 #endif
1232 
1233  return
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,
integer, intent(in)  IIS,
integer, intent(in)  IIE,
integer, intent(in)  JJS,
integer, intent(in)  JJE 
)

calculation Y-flux at XV

Definition at line 1245 of file scale_atmos_dyn_fvm_flux_cd2.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

1245  implicit none
1246 
1247  real(RP), intent(inout) :: flux (KA,IA,JA)
1248  real(RP), intent(in) :: mom (KA,IA,JA)
1249  real(RP), intent(in) :: val (KA,IA,JA)
1250  real(RP), intent(in) :: DENS (KA,IA,JA)
1251  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1252  real(RP), intent(in) :: MAPF ( IA,JA,2)
1253  real(RP), intent(in) :: num_diff(KA,IA,JA)
1254  real(RP), intent(in) :: CDZ (KA)
1255  integer, intent(in) :: IIS, IIE, JJS, JJE
1256 
1257  real(RP) :: vel
1258  integer :: k, i, j
1259  !---------------------------------------------------------------------------
1260 
1261  ! note that y-index is added by -1
1262 
1263  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1264  !$omp private(vel) &
1265  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff)
1266  do j = jjs, jje+1
1267  do i = iis, iie
1268  do k = ks, ke
1269 #ifdef DEBUG
1270  call check( __line__, mom(k,i ,j) )
1271  call check( __line__, mom(k,i,j-1) )
1272 
1273  call check( __line__, val(k,i,j-1) )
1274  call check( __line__, val(k,i,j) )
1275 
1276 #endif
1277  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j-1) ) ) &
1278  / ( dens(k,i,j) )
1279  flux(k,i,j-1) = gsqrt(k,i,j) / mapf(i,j,+1) * vel &
1280  * ( f2 * ( val(k,i,j)+val(k,i,j-1) ) ) &
1281  + gsqrt(k,i,j) * num_diff(k,i,j)
1282  enddo
1283  enddo
1284  enddo
1285 #ifdef DEBUG
1286  k = iundef; i = iundef; j = iundef
1287 #endif
1288 
1289  return
Here is the call graph for this function:
Here is the caller graph for this function:
Generated by   doxygen 1.8.14