SCALE-RM
Functions/Subroutines
scale_atmos_dyn_fvm_flux_ud3koren1993 Module Reference

module scale_atmos_dyn_fvm_flux_ud3Koren1993 More...

Functions/Subroutines

subroutine, public atmos_dyn_fvm_flux_valuew_z_ud3koren1993 (valW, mflx, val, GSQRT, CDZ)
 value at XYW More...
 
subroutine, public atmos_dyn_fvm_fluxz_xyz_ud3koren1993 (flux, mflx, val, GSQRT, num_diff, CDZ, IIS, IIE, JJS, JJE)
 calculation z-flux at XYZ More...
 
subroutine, public atmos_dyn_fvm_fluxx_xyz_ud3koren1993 (flux, mflx, val, GSQRT, num_diff, CDZ, IIS, IIE, JJS, JJE)
 calculation X-flux at XYZ More...
 
subroutine, public atmos_dyn_fvm_fluxy_xyz_ud3koren1993 (flux, mflx, val, GSQRT, num_diff, CDZ, IIS, IIE, JJS, JJE)
 calculation Y-flux at XYZ More...
 
subroutine, public atmos_dyn_fvm_fluxz_xyw_ud3koren1993 (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_ud3koren1993 (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_ud3koren1993 (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_ud3koren1993 (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_ud3koren1993 (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_ud3koren1993 (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_ud3koren1993 (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_ud3koren1993 (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_ud3koren1993 (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_ud3koren1993 (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_ud3koren1993 (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_ud3koren1993 (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_ud3koren1993 (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_ud3koren1993 (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_ud3koren1993 (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_ud3Koren1993

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

Function/Subroutine Documentation

◆ atmos_dyn_fvm_flux_valuew_z_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_flux_valuew_z_ud3koren1993 ( 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 105 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.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().

105  implicit none
106 
107  real(RP), intent(out) :: valW (KA)
108  real(RP), intent(in) :: mflx (KA)
109  real(RP), intent(in) :: val (KA)
110  real(RP), intent(in) :: GSQRT(KA)
111  real(RP), intent(in) :: CDZ (KA)
112 
113  integer :: k
114  !---------------------------------------------------------------------------
115 
116  do k = ks+1, ke-2
117 #ifdef DEBUG
118  call check( __line__, mflx(k) )
119 
120  call check( __line__, val(k) )
121  call check( __line__, val(k+1) )
122 
123  call check( __line__, val(k-1) )
124  call check( __line__, val(k+2) )
125 
126 #endif
127  valw(k) = ( val(k) &
128  + 0.5_rp * phi(val(k+1),val(k),val(k-1)) * ( val(k)-val(k-1) ) ) &
129  * ( 0.5_rp + sign(0.5_rp,mflx(k)) ) &
130  + ( val(k+1) &
131  + 0.5_rp * phi(val(k),val(k+1),val(k+2)) * ( val(k+1)-val(k+2) ) ) &
132  * ( 0.5_rp - sign(0.5_rp,mflx(k)) )
133  enddo
134 #ifdef DEBUG
135  k = iundef
136 #endif
137 
138 #ifdef DEBUG
139 
140  call check( __line__, mflx(ks) )
141  call check( __line__, val(ks ) )
142  call check( __line__, val(ks+1) )
143  call check( __line__, mflx(ke-1) )
144  call check( __line__, val(ke ) )
145  call check( __line__, val(ke-1) )
146 
147 #endif
148 
149 
150 
151  valw(ks) = f2 * ( val(ks+1)+val(ks) )
152  valw(ke-1) = f2 * ( val(ke)+val(ke-1) )
153 
154 
155  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxz_xyz_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxz_xyz_ud3koren1993 ( 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 166 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.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().

166  implicit none
167 
168  real(RP), intent(inout) :: flux (KA,IA,JA)
169  real(RP), intent(in) :: mflx (KA,IA,JA)
170  real(RP), intent(in) :: val (KA,IA,JA)
171  real(RP), intent(in) :: GSQRT (KA,IA,JA)
172  real(RP), intent(in) :: num_diff(KA,IA,JA)
173  real(RP), intent(in) :: CDZ (KA)
174  integer, intent(in) :: IIS, IIE, JJS, JJE
175 
176  real(RP) :: vel
177  integer :: k, i, j
178  !---------------------------------------------------------------------------
179 
180  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
181  !$omp private(vel) &
182  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mflx,val,flux,GSQRT,num_diff)
183  do j = jjs, jje
184  do i = iis, iie
185  do k = ks+1, ke-2
186 #ifdef DEBUG
187  call check( __line__, mflx(k,i,j) )
188 
189  call check( __line__, val(k,i,j) )
190  call check( __line__, val(k+1,i,j) )
191 
192  call check( __line__, val(k-1,i,j) )
193  call check( __line__, val(k+2,i,j) )
194 
195 #endif
196  vel = mflx(k,i,j)
197  flux(k,i,j) = vel &
198  * ( ( val(k,i,j) &
199  + 0.5_rp * phi(val(k+1,i,j),val(k,i,j),val(k-1,i,j)) * ( val(k,i,j)-val(k-1,i,j) ) ) &
200  * ( 0.5_rp + sign(0.5_rp,vel) ) &
201  + ( val(k+1,i,j) &
202  + 0.5_rp * phi(val(k,i,j),val(k+1,i,j),val(k+2,i,j)) * ( val(k+1,i,j)-val(k+2,i,j) ) ) &
203  * ( 0.5_rp - sign(0.5_rp,vel) ) ) &
204  + gsqrt(k,i,j) * num_diff(k,i,j)
205  enddo
206  enddo
207  enddo
208 #ifdef DEBUG
209  k = iundef; i = iundef; j = iundef
210 #endif
211 
212  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
213  !$omp private(vel) &
214  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mflx,val,flux,GSQRT,num_diff)
215  do j = jjs, jje
216  do i = iis, iie
217 #ifdef DEBUG
218 
219  call check( __line__, mflx(ks,i,j) )
220  call check( __line__, val(ks ,i,j) )
221  call check( __line__, val(ks+1,i,j) )
222  call check( __line__, mflx(ke-1,i,j) )
223  call check( __line__, val(ke ,i,j) )
224  call check( __line__, val(ke-1,i,j) )
225 
226 #endif
227  flux(ks-1,i,j) = 0.0_rp
228 
229  vel = mflx(ks,i,j)
230  flux(ks,i,j) = vel &
231  * ( f2 * ( val(ks+1,i,j)+val(ks,i,j) ) ) &
232  + gsqrt(ks,i,j) * num_diff(ks,i,j)
233  vel = mflx(ke-1,i,j)
234  flux(ke-1,i,j) = vel &
235  * ( f2 * ( val(ke,i,j)+val(ke-1,i,j) ) ) &
236  + gsqrt(ke-1,i,j) * num_diff(ke-1,i,j)
237 
238  flux(ke ,i,j) = 0.0_rp
239  enddo
240  enddo
241 #ifdef DEBUG
242  k = iundef; i = iundef; j = iundef
243 #endif
244 
245  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxx_xyz_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxx_xyz_ud3koren1993 ( 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 256 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.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().

256  implicit none
257 
258  real(RP), intent(inout) :: flux (KA,IA,JA)
259  real(RP), intent(in) :: mflx (KA,IA,JA)
260  real(RP), intent(in) :: val (KA,IA,JA)
261  real(RP), intent(in) :: GSQRT (KA,IA,JA)
262  real(RP), intent(in) :: num_diff(KA,IA,JA)
263  real(RP), intent(in) :: CDZ(KA)
264  integer, intent(in) :: IIS, IIE, JJS, JJE
265 
266  real(RP) :: vel
267  integer :: k, i, j
268  !---------------------------------------------------------------------------
269 
270  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
271  !$omp private(vel) &
272  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mflx,val,flux,GSQRT,num_diff)
273  do j = jjs, jje
274  do i = iis-1, iie
275  do k = ks, ke
276 #ifdef DEBUG
277  call check( __line__, mflx(k,i,j) )
278 
279  call check( __line__, val(k,i,j) )
280  call check( __line__, val(k,i+1,j) )
281 
282  call check( __line__, val(k,i-1,j) )
283  call check( __line__, val(k,i+2,j) )
284 
285 #endif
286  vel = mflx(k,i,j)
287  flux(k,i,j) = vel &
288  * ( ( val(k,i,j) &
289  + 0.5_rp * phi(val(k,i+1,j),val(k,i,j),val(k,i-1,j)) * ( val(k,i,j)-val(k,i-1,j) ) ) &
290  * ( 0.5_rp + sign(0.5_rp,vel) ) &
291  + ( val(k,i+1,j) &
292  + 0.5_rp * phi(val(k,i,j),val(k,i+1,j),val(k,i+2,j)) * ( val(k,i+1,j)-val(k,i+2,j) ) ) &
293  * ( 0.5_rp - sign(0.5_rp,vel) ) ) &
294  + gsqrt(k,i,j) * num_diff(k,i,j)
295  enddo
296  enddo
297  enddo
298 #ifdef DEBUG
299  k = iundef; i = iundef; j = iundef
300 #endif
301 
302  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxy_xyz_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxy_xyz_ud3koren1993 ( 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 313 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.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().

313  implicit none
314 
315  real(RP), intent(inout) :: flux (KA,IA,JA)
316  real(RP), intent(in) :: mflx (KA,IA,JA)
317  real(RP), intent(in) :: val (KA,IA,JA)
318  real(RP), intent(in) :: GSQRT (KA,IA,JA)
319  real(RP), intent(in) :: num_diff(KA,IA,JA)
320  real(RP), intent(in) :: CDZ(KA)
321  integer, intent(in) :: IIS, IIE, JJS, JJE
322 
323  real(RP) :: vel
324  integer :: k, i, j
325  !---------------------------------------------------------------------------
326 
327  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
328  !$omp private(vel) &
329  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mflx,val,flux,GSQRT,num_diff)
330  do j = jjs-1, jje
331  do i = iis, iie
332  do k = ks, ke
333 #ifdef DEBUG
334  call check( __line__, mflx(k,i,j) )
335 
336  call check( __line__, val(k,i,j) )
337  call check( __line__, val(k,i,j+1) )
338 
339  call check( __line__, val(k,i,j-1) )
340  call check( __line__, val(k,i,j+2) )
341 
342 #endif
343  vel = mflx(k,i,j)
344  flux(k,i,j) = vel &
345  * ( ( val(k,i,j) &
346  + 0.5_rp * phi(val(k,i,j+1),val(k,i,j),val(k,i,j-1)) * ( val(k,i,j)-val(k,i,j-1) ) ) &
347  * ( 0.5_rp + sign(0.5_rp,vel) ) &
348  + ( val(k,i,j+1) &
349  + 0.5_rp * phi(val(k,i,j),val(k,i,j+1),val(k,i,j+2)) * ( val(k,i,j+1)-val(k,i,j+2) ) ) &
350  * ( 0.5_rp - sign(0.5_rp,vel) ) ) &
351  + gsqrt(k,i,j) * num_diff(k,i,j)
352  enddo
353  enddo
354  enddo
355 #ifdef DEBUG
356  k = iundef; i = iundef; j = iundef
357 #endif
358 
359  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxz_xyw_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxz_xyw_ud3koren1993 ( 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 373 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.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().

373  implicit none
374 
375  real(RP), intent(inout) :: flux (KA,IA,JA)
376  real(RP), intent(in) :: mom (KA,IA,JA)
377  real(RP), intent(in) :: val (KA,IA,JA)
378  real(RP), intent(in) :: DENS (KA,IA,JA)
379  real(RP), intent(in) :: GSQRT (KA,IA,JA)
380  real(RP), intent(in) :: J33G
381  real(RP), intent(in) :: num_diff(KA,IA,JA)
382  real(RP), intent(in) :: CDZ (KA)
383  real(RP), intent(in) :: FDZ (KA-1)
384  real(RP), intent(in) :: dtrk
385  integer, intent(in) :: IIS, IIE, JJS, JJE
386 
387  real(RP) :: vel
388  integer :: k, i, j
389  !---------------------------------------------------------------------------
390 
391  ! note than z-index is added by -1
392 
393  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
394  !$omp private(vel) &
395  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,flux,J33G,GSQRT,num_diff,DENS)
396  do j = jjs, jje
397  do i = iis, iie
398  do k = ks+2, ke-1
399 #ifdef DEBUG
400  call check( __line__, mom(k-1,i,j) )
401  call check( __line__, mom(k ,i,j) )
402 
403  call check( __line__, val(k-1,i,j) )
404  call check( __line__, val(k,i,j) )
405 
406  call check( __line__, val(k-2,i,j) )
407  call check( __line__, val(k+1,i,j) )
408 
409 #endif
410  vel = ( 0.5_rp * ( mom(k-1,i,j) &
411  + mom(k,i,j) ) ) &
412  / dens(k,i,j)
413  flux(k-1,i,j) = j33g * vel &
414  * ( ( val(k-1,i,j) &
415  + 0.5_rp * phi(val(k,i,j),val(k-1,i,j),val(k-2,i,j)) * ( val(k-1,i,j)-val(k-2,i,j) ) ) &
416  * ( 0.5_rp + sign(0.5_rp,vel) ) &
417  + ( val(k,i,j) &
418  + 0.5_rp * phi(val(k-1,i,j),val(k,i,j),val(k+1,i,j)) * ( val(k,i,j)-val(k+1,i,j) ) ) &
419  * ( 0.5_rp - sign(0.5_rp,vel) ) ) &
420  + gsqrt(k,i,j) * num_diff(k,i,j)
421  enddo
422  enddo
423  enddo
424 #ifdef DEBUG
425  k = iundef; i = iundef; j = iundef
426 #endif
427 
428  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
429  !$omp private(vel) &
430  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J33G,GSQRT,num_diff,FDZ,dtrk)
431  do j = jjs, jje
432  do i = iis, iie
433 #ifdef DEBUG
434 
435  call check( __line__, val(ks ,i,j) )
436  call check( __line__, val(ks+1,i,j) )
437 
438 
439 #endif
440  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS.
441  ! The flux at KS can be non-zero.
442  ! To reduce calculations, all the fluxes are set to zero.
443  flux(ks-1,i,j) = 0.0_rp ! k = KS
444 
445  vel = ( 0.5_rp * ( mom(ks,i,j) &
446  + mom(ks+1,i,j) ) ) &
447  / dens(ks+1,i,j)
448  flux(ks,i,j) = j33g * vel &
449  * ( f2 * ( val(ks+1,i,j)+val(ks,i,j) ) ) &
450  + gsqrt(ks+1,i,j) * num_diff(ks+1,i,j) ! k = KS+1
451 
452 
453 
454  flux(ke-1,i,j) = 0.0_rp ! k = KE
455  flux(ke ,i,j) = 0.0_rp ! k = KE+1
456  enddo
457  enddo
458 
459  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxj13_xyw_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxj13_xyw_ud3koren1993 ( 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 471 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

471  implicit none
472 
473  real(RP), intent(inout) :: flux (KA,IA,JA)
474  real(RP), intent(in) :: mom (KA,IA,JA)
475  real(RP), intent(in) :: val (KA,IA,JA)
476  real(RP), intent(in) :: DENS (KA,IA,JA)
477  real(RP), intent(in) :: GSQRT (KA,IA,JA)
478  real(RP), intent(in) :: J13G (KA,IA,JA)
479  real(RP), intent(in) :: MAPF ( IA,JA,2)
480  real(RP), intent(in) :: CDZ (KA)
481  integer, intent(in) :: IIS, IIE, JJS, JJE
482 
483  real(RP) :: vel
484  integer :: k, i, j
485  !---------------------------------------------------------------------------
486 
487  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
488  !$omp private(vel) &
489  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J13G,MAPF)
490  do j = jjs, jje
491  do i = iis, iie
492  do k = ks+2, ke-1
493  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i-1,j) ) ) &
494  / dens(k,i,j)
495  flux(k-1,i,j) = j13g(k,i,j) / mapf(i,j,+2) * vel &
496  * ( ( val(k-1,i,j) &
497  + 0.5_rp * phi(val(k,i,j),val(k-1,i,j),val(k-2,i,j)) * ( val(k-1,i,j)-val(k-2,i,j) ) ) &
498  * ( 0.5_rp + sign(0.5_rp,vel) ) &
499  + ( val(k,i,j) &
500  + 0.5_rp * phi(val(k-1,i,j),val(k,i,j),val(k+1,i,j)) * ( val(k,i,j)-val(k+1,i,j) ) ) &
501  * ( 0.5_rp - sign(0.5_rp,vel) ) )
502  enddo
503  enddo
504  enddo
505 
506  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
507  !$omp private(vel) &
508  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J13G,MAPF)
509  do j = jjs, jje
510  do i = iis, iie
511  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS.
512  ! The flux at KS can be non-zero.
513  ! To reduce calculations, all the fluxes are set to zero.
514  flux(ks-1,i,j) = 0.0_rp ! k = KS
515 
516  ! physically incorrect but for numerical stability
517  vel = ( ( 0.5_rp * ( mom(ks+1,i,j)+mom(ks+1,i-1,j) ) ) / dens(ks+1,i,j) &
518  + ( 0.5_rp * ( mom(ks,i,j)+mom(ks,i-1,j) ) ) / dens(ks ,i,j) ) * 0.5_rp
519 ! vel = ( 0.5_RP * ( mom(KS+1,i,j)+mom(KS+1,i-1,j) ) ) &
520 ! / DENS(KS+1,i,j)
521  flux(ks,i,j) = j13g(ks+1,i,j) / mapf(i,j,+2) * vel &
522  * ( f2 * ( val(ks+1,i,j)+val(ks,i,j) ) ) ! k = KS+1
523 
524 
525  flux(ke-1,i,j) = 0.0_rp
526  enddo
527  enddo
528 
529  return
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxj23_xyw_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxj23_xyw_ud3koren1993 ( 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 540 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

540  implicit none
541 
542  real(RP), intent(inout) :: flux (KA,IA,JA)
543  real(RP), intent(in) :: mom (KA,IA,JA)
544  real(RP), intent(in) :: val (KA,IA,JA)
545  real(RP), intent(in) :: DENS (KA,IA,JA)
546  real(RP), intent(in) :: GSQRT (KA,IA,JA)
547  real(RP), intent(in) :: J23G (KA,IA,JA)
548  real(RP), intent(in) :: MAPF ( IA,JA,2)
549  real(RP), intent(in) :: CDZ (KA)
550  integer, intent(in) :: IIS, IIE, JJS, JJE
551 
552  real(RP) :: vel
553  integer :: k, i, j
554  !---------------------------------------------------------------------------
555 
556  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
557  !$omp private(vel) &
558  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF)
559  do j = jjs, jje
560  do i = iis, iie
561  do k = ks+2, ke-1
562  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j-1) ) ) &
563  / dens(k,i,j)
564  flux(k-1,i,j) = j23g(k,i,j) / mapf(i,j,+1) * vel &
565  * ( ( val(k-1,i,j) &
566  + 0.5_rp * phi(val(k,i,j),val(k-1,i,j),val(k-2,i,j)) * ( val(k-1,i,j)-val(k-2,i,j) ) ) &
567  * ( 0.5_rp + sign(0.5_rp,vel) ) &
568  + ( val(k,i,j) &
569  + 0.5_rp * phi(val(k-1,i,j),val(k,i,j),val(k+1,i,j)) * ( val(k,i,j)-val(k+1,i,j) ) ) &
570  * ( 0.5_rp - sign(0.5_rp,vel) ) )
571  enddo
572  enddo
573  enddo
574 
575  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
576  !$omp private(vel) &
577  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF)
578  do j = jjs, jje
579  do i = iis, iie
580  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS.
581  ! The flux at KS can be non-zero.
582  ! To reduce calculations, all the fluxes are set to zero.
583  flux(ks-1,i,j) = 0.0_rp ! k = KS
584 
585  ! physically incorrect but for numerical stability
586  vel = ( ( 0.5_rp * ( mom(ks+1,i,j)+mom(ks+1,i,j-1) ) ) / dens(ks+1,i,j) &
587  + ( 0.5_rp * ( mom(ks,i,j)+mom(ks,i,j-1) ) ) / dens(ks ,i,j) ) * 0.5_rp
588 ! vel = ( 0.5_RP * ( mom(KS+1,i,j)+mom(KS+1,i,j-1) ) ) &
589 ! / DENS(KS+1,i,j)
590  flux(ks,i,j) = j23g(ks+1,i,j) / mapf(i,j,+1) * vel &
591  * ( f2 * ( val(ks+1,i,j)+val(ks,i,j) ) ) ! k = KS+1
592 
593 
594  flux(ke-1,i,j) = 0.0_rp
595  enddo
596  enddo
597 
598  return
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxx_xyw_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxx_xyw_ud3koren1993 ( 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 611 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.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().

611  implicit none
612 
613  real(RP), intent(inout) :: flux (KA,IA,JA)
614  real(RP), intent(in) :: mom (KA,IA,JA)
615  real(RP), intent(in) :: val (KA,IA,JA)
616  real(RP), intent(in) :: DENS (KA,IA,JA)
617  real(RP), intent(in) :: GSQRT (KA,IA,JA)
618  real(RP), intent(in) :: MAPF ( IA,JA,2)
619  real(RP), intent(in) :: num_diff(KA,IA,JA)
620  real(RP), intent(in) :: CDZ (KA)
621  integer, intent(in) :: IIS, IIE, JJS, JJE
622 
623  real(RP) :: vel
624  integer :: k, i, j
625  !---------------------------------------------------------------------------
626 
627  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
628  !$omp private(vel) &
629  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff) &
630  !$omp shared(CDZ)
631  do j = jjs, jje
632  do i = iis-1, iie
633  do k = ks, ke-1
634 #ifdef DEBUG
635  call check( __line__, mom(k ,i,j) )
636  call check( __line__, mom(k+1,i,j) )
637 
638  call check( __line__, val(k,i,j) )
639  call check( __line__, val(k,i+1,j) )
640 
641  call check( __line__, val(k,i-1,j) )
642  call check( __line__, val(k,i+2,j) )
643 
644 #endif
645  vel = ( f2h(k,1,i_uyz) &
646  * mom(k+1,i,j) &
647  + f2h(k,2,i_uyz) &
648  * mom(k,i,j) ) &
649  / ( f2h(k,1,i_uyz) &
650  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
651  + f2h(k,2,i_uyz) &
652  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
653  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+2) * vel &
654  * ( ( val(k,i,j) &
655  + 0.5_rp * phi(val(k,i+1,j),val(k,i,j),val(k,i-1,j)) * ( val(k,i,j)-val(k,i-1,j) ) ) &
656  * ( 0.5_rp + sign(0.5_rp,vel) ) &
657  + ( val(k,i+1,j) &
658  + 0.5_rp * phi(val(k,i,j),val(k,i+1,j),val(k,i+2,j)) * ( val(k,i+1,j)-val(k,i+2,j) ) ) &
659  * ( 0.5_rp - sign(0.5_rp,vel) ) ) &
660  + gsqrt(k,i,j) * num_diff(k,i,j)
661  enddo
662  enddo
663  enddo
664 #ifdef DEBUG
665  k = iundef; i = iundef; j = iundef
666 #endif
667 
668  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
669  !$omp private(vel) &
670  !$omp shared(JJS,JJE,IIS,IIE,KE,flux)
671  do j = jjs, jje
672  do i = iis-1, iie
673  flux(ke,i,j) = 0.0_rp
674  enddo
675  enddo
676 #ifdef DEBUG
677  k = iundef; i = iundef; j = iundef
678 #endif
679 
680  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_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxy_xyw_ud3koren1993 ( 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 692 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.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().

692  implicit none
693 
694  real(RP), intent(inout) :: flux (KA,IA,JA)
695  real(RP), intent(in) :: mom (KA,IA,JA)
696  real(RP), intent(in) :: val (KA,IA,JA)
697  real(RP), intent(in) :: DENS (KA,IA,JA)
698  real(RP), intent(in) :: GSQRT (KA,IA,JA)
699  real(RP), intent(in) :: MAPF ( IA,JA,2)
700  real(RP), intent(in) :: num_diff(KA,IA,JA)
701  real(RP), intent(in) :: CDZ (KA)
702  integer, intent(in) :: IIS, IIE, JJS, JJE
703 
704  real(RP) :: vel
705  integer :: k, i, j
706  !---------------------------------------------------------------------------
707 
708  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
709  !$omp private(vel) &
710  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff) &
711  !$omp shared(CDZ)
712  do j = jjs-1, jje
713  do i = iis, iie
714  do k = ks, ke-1
715 #ifdef DEBUG
716  call check( __line__, mom(k ,i,j) )
717  call check( __line__, mom(k+1,i,j) )
718 
719  call check( __line__, val(k,i,j) )
720  call check( __line__, val(k,i,j+1) )
721 
722  call check( __line__, val(k,i,j-1) )
723  call check( __line__, val(k,i,j+2) )
724 
725 #endif
726  vel = ( f2h(k,1,i_xvz) &
727  * mom(k+1,i,j) &
728  + f2h(k,2,i_xvz) &
729  * mom(k,i,j) ) &
730  / ( f2h(k,1,i_xvz) &
731  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
732  + f2h(k,2,i_xvz) &
733  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
734  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+1) * vel &
735  * ( ( val(k,i,j) &
736  + 0.5_rp * phi(val(k,i,j+1),val(k,i,j),val(k,i,j-1)) * ( val(k,i,j)-val(k,i,j-1) ) ) &
737  * ( 0.5_rp + sign(0.5_rp,vel) ) &
738  + ( val(k,i,j+1) &
739  + 0.5_rp * phi(val(k,i,j),val(k,i,j+1),val(k,i,j+2)) * ( val(k,i,j+1)-val(k,i,j+2) ) ) &
740  * ( 0.5_rp - sign(0.5_rp,vel) ) ) &
741  + gsqrt(k,i,j) * num_diff(k,i,j)
742  enddo
743  enddo
744  enddo
745 #ifdef DEBUG
746  k = iundef; i = iundef; j = iundef
747 #endif
748 
749  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
750  !$omp private(vel) &
751  !$omp shared(JJS,JJE,IIS,IIE,KE,flux)
752  do j = jjs-1, jje
753  do i = iis, iie
754  flux(ke,i,j) = 0.0_rp
755  enddo
756  enddo
757 #ifdef DEBUG
758  k = iundef; i = iundef; j = iundef
759 #endif
760 
761  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_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxz_uyz_ud3koren1993 ( 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 774 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.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().

774  implicit none
775 
776  real(RP), intent(inout) :: flux (KA,IA,JA)
777  real(RP), intent(in) :: mom (KA,IA,JA)
778  real(RP), intent(in) :: val (KA,IA,JA)
779  real(RP), intent(in) :: DENS (KA,IA,JA)
780  real(RP), intent(in) :: GSQRT (KA,IA,JA)
781  real(RP), intent(in) :: J33G
782  real(RP), intent(in) :: num_diff(KA,IA,JA)
783  real(RP), intent(in) :: CDZ (KA)
784  integer, intent(in) :: IIS, IIE, JJS, JJE
785 
786  real(RP) :: vel
787  integer :: k, i, j
788  !---------------------------------------------------------------------------
789 
790  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
791  !$omp private(vel) &
792  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J33G,GSQRT,num_diff) &
793  !$omp shared(CDZ)
794  do j = jjs, jje
795  do i = iis, iie
796  do k = ks+1, ke-2
797 #ifdef DEBUG
798  call check( __line__, mom(k,i,j) )
799  call check( __line__, mom(k,i+1,j) )
800 
801  call check( __line__, val(k,i,j) )
802  call check( __line__, val(k+1,i,j) )
803 
804  call check( __line__, val(k-1,i,j) )
805  call check( __line__, val(k+2,i,j) )
806 
807 #endif
808  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i+1,j) ) ) &
809  / ( f2h(k,1,i_uyz) &
810  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
811  + f2h(k,2,i_uyz) &
812  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
813  flux(k,i,j) = j33g * vel &
814  * ( ( val(k,i,j) &
815  + 0.5_rp * phi(val(k+1,i,j),val(k,i,j),val(k-1,i,j)) * ( val(k,i,j)-val(k-1,i,j) ) ) &
816  * ( 0.5_rp + sign(0.5_rp,vel) ) &
817  + ( val(k+1,i,j) &
818  + 0.5_rp * phi(val(k,i,j),val(k+1,i,j),val(k+2,i,j)) * ( val(k+1,i,j)-val(k+2,i,j) ) ) &
819  * ( 0.5_rp - sign(0.5_rp,vel) ) ) &
820  + gsqrt(k,i,j) * num_diff(k,i,j)
821  enddo
822  enddo
823  enddo
824 #ifdef DEBUG
825  k = iundef; i = iundef; j = iundef
826 #endif
827 
828  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
829  !$omp private(vel) &
830  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,flux,J33G,GSQRT,num_diff,DENS,CDZ)
831  do j = jjs, jje
832  do i = iis, iie
833 #ifdef DEBUG
834 
835  call check( __line__, mom(ks,i ,j) )
836  call check( __line__, mom(ks,i+1,j) )
837  call check( __line__, val(ks+1,i,j) )
838  call check( __line__, val(ks,i,j) )
839 
840 #endif
841  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
842  ! The flux at KS-1 can be non-zero.
843  ! To reduce calculations, all the fluxes are set to zero.
844  flux(ks-1,i,j) = 0.0_rp
845 
846  vel = ( 0.5_rp * ( mom(ks,i,j)+mom(ks,i+1,j) ) ) &
847  / ( f2h(ks,1,i_uyz) &
848  * 0.5_rp * ( dens(ks+1,i,j)+dens(ks+1,i+1,j) ) &
849  + f2h(ks,2,i_uyz) &
850  * 0.5_rp * ( dens(ks,i,j)+dens(ks,i+1,j) ) )
851  flux(ks,i,j) = j33g * vel &
852  * ( f2 * ( val(ks+1,i,j)+val(ks,i,j) ) ) &
853  + gsqrt(ks,i,j) * num_diff(ks,i,j)
854  vel = ( 0.5_rp * ( mom(ke-1,i,j)+mom(ke-1,i+1,j) ) ) &
855  / ( f2h(ke-1,1,i_uyz) &
856  * 0.5_rp * ( dens(ke,i,j)+dens(ke,i+1,j) ) &
857  + f2h(ke-1,2,i_uyz) &
858  * 0.5_rp * ( dens(ke-1,i,j)+dens(ke-1,i+1,j) ) )
859  flux(ke-1,i,j) = j33g * vel &
860  * ( f2 * ( val(ke,i,j)+val(ke-1,i,j) ) ) &
861  + gsqrt(ke-1,i,j) * num_diff(ke-1,i,j)
862 
863  flux(ke,i,j) = 0.0_rp
864  enddo
865  enddo
866 #ifdef DEBUG
867  k = iundef; i = iundef; j = iundef
868 #endif
869 
870  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_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxj13_uyz_ud3koren1993 ( 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 881 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

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

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxj23_uyz_ud3koren1993 ( 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 971 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

971  implicit none
972 
973  real(RP), intent(inout) :: flux (KA,IA,JA)
974  real(RP), intent(in) :: mom (KA,IA,JA)
975  real(RP), intent(in) :: val (KA,IA,JA)
976  real(RP), intent(in) :: DENS (KA,IA,JA)
977  real(RP), intent(in) :: GSQRT (KA,IA,JA)
978  real(RP), intent(in) :: J23G (KA,IA,JA)
979  real(RP), intent(in) :: MAPF ( IA,JA,2)
980  real(RP), intent(in) :: CDZ (KA)
981  integer, intent(in) :: IIS, IIE, JJS, JJE
982 
983  real(RP) :: vel
984  integer :: k, i, j
985  !---------------------------------------------------------------------------
986 
987  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
988  !$omp private(vel) &
989  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF) &
990  !$omp shared(GSQRT,CDZ)
991  do j = jjs, jje
992  do i = iis, iie
993  do k = ks+1, ke-2
994  vel = ( f2h(k,1,i_uyz) &
995  * 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) ) &
996  + f2h(k,2,i_uyz) &
997  * 0.25_rp * ( mom(k,i,j)+mom(k,i+1,j)+mom(k,i,j-1)+mom(k,i+1,j-1) ) ) &
998  / ( f2h(k,1,i_uyz) &
999  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i+1,j) ) &
1000  + f2h(k,2,i_uyz) &
1001  * 0.5_rp * ( dens(k,i,j)+dens(k,i+1,j) ) )
1002  flux(k,i,j) = j23g(k,i,j) / mapf(i,j,+1) * vel &
1003  * ( ( val(k,i,j) &
1004  + 0.5_rp * phi(val(k+1,i,j),val(k,i,j),val(k-1,i,j)) * ( val(k,i,j)-val(k-1,i,j) ) ) &
1005  * ( 0.5_rp + sign(0.5_rp,vel) ) &
1006  + ( val(k+1,i,j) &
1007  + 0.5_rp * phi(val(k,i,j),val(k+1,i,j),val(k+2,i,j)) * ( val(k+1,i,j)-val(k+2,i,j) ) ) &
1008  * ( 0.5_rp - sign(0.5_rp,vel) ) )
1009  enddo
1010  enddo
1011  enddo
1012 
1013  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
1014  !$omp private(vel) &
1015  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF) &
1016  !$omp shared(GSQRT,CDZ)
1017  do j = jjs, jje
1018  do i = iis, iie
1019  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
1020  ! The flux at KS-1 can be non-zero.
1021  ! To reduce calculations, all the fluxes are set to zero.
1022  flux(ks-1,i,j) = 0.0_rp
1023 
1024  vel = ( f2h(ks,1,i_uyz) &
1025  * 0.25_rp * ( mom(ks+1,i,j)+mom(ks+1,i+1,j)+mom(ks+1,i,j-1)+mom(ks+1,i+1,j-1) ) &
1026  + f2h(ks,2,i_uyz) &
1027  * 0.25_rp * ( mom(ks,i,j)+mom(ks,i+1,j)+mom(ks,i,j-1)+mom(ks,i+1,j-1) ) ) &
1028  / ( f2h(ks,1,i_uyz) &
1029  * 0.5_rp * ( dens(ks+1,i,j)+dens(ks+1,i+1,j) ) &
1030  + f2h(ks,2,i_uyz) &
1031  * 0.5_rp * ( dens(ks,i,j)+dens(ks,i+1,j) ) )
1032  flux(ks,i,j) = j23g(ks,i,j) / mapf(i,j,+1) * vel &
1033  * ( f2 * ( val(ks+1,i,j)+val(ks,i,j) ) )
1034 
1035  vel = ( f2h(ke-1,1,i_uyz) &
1036  * 0.25_rp * ( mom(ke,i,j)+mom(ke,i+1,j)+mom(ke,i,j-1)+mom(ke,i+1,j-1) ) &
1037  + f2h(ke-1,2,i_uyz) &
1038  * 0.25_rp * ( mom(ke-1,i,j)+mom(ke-1,i+1,j)+mom(ke-1,i,j-1)+mom(ke-1,i+1,j-1) ) ) &
1039  / ( f2h(ke-1,1,i_uyz) &
1040  * 0.5_rp * ( dens(ke,i,j)+dens(ke,i+1,j) ) &
1041  + f2h(ke-1,2,i_uyz) &
1042  * 0.5_rp * ( dens(ke-1,i,j)+dens(ke-1,i+1,j) ) )
1043  flux(ke-1,i,j) = j23g(ke-1,i,j) / mapf(i,j,+1) * vel &
1044  * ( f2 * ( val(ke,i,j)+val(ke-1,i,j) ) )
1045 
1046  flux(ke ,i,j) = 0.0_rp
1047  enddo
1048  enddo
1049 
1050  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_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxx_uyz_ud3koren1993 ( 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 1062 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.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().

1062  implicit none
1063 
1064  real(RP), intent(inout) :: flux (KA,IA,JA)
1065  real(RP), intent(in) :: mom (KA,IA,JA)
1066  real(RP), intent(in) :: val (KA,IA,JA)
1067  real(RP), intent(in) :: DENS (KA,IA,JA)
1068  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1069  real(RP), intent(in) :: MAPF ( IA,JA,2)
1070  real(RP), intent(in) :: num_diff(KA,IA,JA)
1071  real(RP), intent(in) :: CDZ (KA)
1072  integer, intent(in) :: IIS, IIE, JJS, JJE
1073 
1074  real(RP) :: vel
1075  integer :: k, i, j
1076  !---------------------------------------------------------------------------
1077 
1078  ! note that x-index is added by -1
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,GSQRT,MAPF,num_diff)
1083  do j = jjs, jje
1084  do i = iis, iie+1
1085  do k = ks, ke
1086 #ifdef DEBUG
1087  call check( __line__, mom(k,i ,j) )
1088  call check( __line__, mom(k,i-1,j) )
1089 
1090  call check( __line__, val(k,i-1,j) )
1091  call check( __line__, val(k,i,j) )
1092 
1093  call check( __line__, val(k,i-2,j) )
1094  call check( __line__, val(k,i+1,j) )
1095 
1096 #endif
1097  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i-1,j) ) ) &
1098  / ( dens(k,i,j) )
1099  flux(k,i-1,j) = gsqrt(k,i,j) / mapf(i,j,+2) * vel &
1100  * ( ( val(k,i-1,j) &
1101  + 0.5_rp * phi(val(k,i,j),val(k,i-1,j),val(k,i-2,j)) * ( val(k,i-1,j)-val(k,i-2,j) ) ) &
1102  * ( 0.5_rp + sign(0.5_rp,vel) ) &
1103  + ( val(k,i,j) &
1104  + 0.5_rp * phi(val(k,i-1,j),val(k,i,j),val(k,i+1,j)) * ( val(k,i,j)-val(k,i+1,j) ) ) &
1105  * ( 0.5_rp - sign(0.5_rp,vel) ) ) &
1106  + gsqrt(k,i,j) * num_diff(k,i,j)
1107  enddo
1108  enddo
1109  enddo
1110 #ifdef DEBUG
1111  k = iundef; i = iundef; j = iundef
1112 #endif
1113 
1114  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxy_uyz_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxy_uyz_ud3koren1993 ( 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 1126 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.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().

1126  implicit none
1127 
1128  real(RP), intent(inout) :: flux (KA,IA,JA)
1129  real(RP), intent(in) :: mom (KA,IA,JA)
1130  real(RP), intent(in) :: val (KA,IA,JA)
1131  real(RP), intent(in) :: DENS (KA,IA,JA)
1132  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1133  real(RP), intent(in) :: MAPF ( IA,JA,2)
1134  real(RP), intent(in) :: num_diff(KA,IA,JA)
1135  real(RP), intent(in) :: CDZ (KA)
1136  integer, intent(in) :: IIS, IIE, JJS, JJE
1137 
1138  real(RP) :: vel
1139  integer :: k, i, j
1140  !---------------------------------------------------------------------------
1141 
1142  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1143  !$omp private(vel) &
1144  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff)
1145  do j = jjs-1, jje
1146  do i = iis, iie
1147  do k = ks, ke
1148 #ifdef DEBUG
1149  call check( __line__, mom(k,i ,j) )
1150  call check( __line__, mom(k,i-1,j) )
1151 
1152  call check( __line__, val(k,i,j) )
1153  call check( __line__, val(k,i,j+1) )
1154 
1155  call check( __line__, val(k,i,j-1) )
1156  call check( __line__, val(k,i,j+2) )
1157 
1158 #endif
1159  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i+1,j) ) ) &
1160  / ( 0.25_rp * ( dens(k,i,j)+dens(k,i+1,j)+dens(k,i,j+1)+dens(k,i+1,j+1) ) )
1161  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+1) * vel &
1162  * ( ( val(k,i,j) &
1163  + 0.5_rp * phi(val(k,i,j+1),val(k,i,j),val(k,i,j-1)) * ( val(k,i,j)-val(k,i,j-1) ) ) &
1164  * ( 0.5_rp + sign(0.5_rp,vel) ) &
1165  + ( val(k,i,j+1) &
1166  + 0.5_rp * phi(val(k,i,j),val(k,i,j+1),val(k,i,j+2)) * ( val(k,i,j+1)-val(k,i,j+2) ) ) &
1167  * ( 0.5_rp - sign(0.5_rp,vel) ) ) &
1168  + gsqrt(k,i,j) * num_diff(k,i,j)
1169  enddo
1170  enddo
1171  enddo
1172 #ifdef DEBUG
1173  k = iundef; i = iundef; j = iundef
1174 #endif
1175 
1176  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxz_xvz_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxz_xvz_ud3koren1993 ( 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 1190 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.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().

1190  implicit none
1191 
1192  real(RP), intent(inout) :: flux (KA,IA,JA)
1193  real(RP), intent(in) :: mom (KA,IA,JA)
1194  real(RP), intent(in) :: val (KA,IA,JA)
1195  real(RP), intent(in) :: DENS (KA,IA,JA)
1196  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1197  real(RP), intent(in) :: J33G
1198  real(RP), intent(in) :: num_diff(KA,IA,JA)
1199  real(RP), intent(in) :: CDZ (KA)
1200  integer, intent(in) :: IIS, IIE, JJS, JJE
1201 
1202  real(RP) :: vel
1203  integer :: k, i, j
1204  !---------------------------------------------------------------------------
1205 
1206  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1207  !$omp private(vel) &
1208  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J33G,GSQRT,num_diff) &
1209  !$omp shared(CDZ)
1210  do j = jjs, jje
1211  do i = iis, iie
1212  do k = ks+1, ke-2
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+1,i,j) )
1219 
1220  call check( __line__, val(k-1,i,j) )
1221  call check( __line__, val(k+2,i,j) )
1222 
1223 #endif
1224  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j+1) ) ) &
1225  / ( f2h(k,1,i_xvz) &
1226  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
1227  + f2h(k,2,i_xvz) &
1228  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
1229  flux(k,i,j) = j33g * vel &
1230  * ( ( val(k,i,j) &
1231  + 0.5_rp * phi(val(k+1,i,j),val(k,i,j),val(k-1,i,j)) * ( val(k,i,j)-val(k-1,i,j) ) ) &
1232  * ( 0.5_rp + sign(0.5_rp,vel) ) &
1233  + ( val(k+1,i,j) &
1234  + 0.5_rp * phi(val(k,i,j),val(k+1,i,j),val(k+2,i,j)) * ( val(k+1,i,j)-val(k+2,i,j) ) ) &
1235  * ( 0.5_rp - sign(0.5_rp,vel) ) ) &
1236  + gsqrt(k,i,j) * num_diff(k,i,j)
1237  enddo
1238  enddo
1239  enddo
1240 #ifdef DEBUG
1241  k = iundef; i = iundef; j = iundef
1242 #endif
1243 
1244  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
1245  !$omp private(vel) &
1246  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,flux,J33G,GSQRT,num_diff,DENS,CDZ)
1247  do j = jjs, jje
1248  do i = iis, iie
1249 #ifdef DEBUG
1250 
1251  call check( __line__, mom(ks,i ,j) )
1252  call check( __line__, mom(ks,i,j+1) )
1253  call check( __line__, val(ks+1,i,j) )
1254  call check( __line__, val(ks,i,j) )
1255 
1256 #endif
1257  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
1258  ! The flux at KS-1 can be non-zero.
1259  ! To reduce calculations, all the fluxes are set to zero.
1260  flux(ks-1,i,j) = 0.0_rp
1261 
1262  vel = ( 0.5_rp * ( mom(ks,i,j)+mom(ks,i,j+1) ) ) &
1263  / ( f2h(ks,1,i_xvz) &
1264  * 0.5_rp * ( dens(ks+1,i,j)+dens(ks+1,i,j+1) ) &
1265  + f2h(ks,2,i_xvz) &
1266  * 0.5_rp * ( dens(ks,i,j)+dens(ks,i,j+1) ) )
1267  flux(ks,i,j) = j33g * vel &
1268  * ( f2 * ( val(ks+1,i,j)+val(ks,i,j) ) ) &
1269  + gsqrt(ks,i,j) * num_diff(ks,i,j)
1270  vel = ( 0.5_rp * ( mom(ke-1,i,j)+mom(ke-1,i,j+1) ) ) &
1271  / ( f2h(ke-1,1,i_xvz) &
1272  * 0.5_rp * ( dens(ke,i,j)+dens(ke,i,j+1) ) &
1273  + f2h(ke-1,2,i_xvz) &
1274  * 0.5_rp * ( dens(ke-1,i,j)+dens(ke-1,i,j+1) ) )
1275  flux(ke-1,i,j) = j33g * vel &
1276  * ( f2 * ( val(ke,i,j)+val(ke-1,i,j) ) ) &
1277  + gsqrt(ke-1,i,j) * num_diff(ke-1,i,j)
1278 
1279  flux(ke,i,j) = 0.0_rp
1280  enddo
1281  enddo
1282 #ifdef DEBUG
1283  k = iundef; i = iundef; j = iundef
1284 #endif
1285 
1286  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_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxj13_xvz_ud3koren1993 ( 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 1297 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

1297  implicit none
1298 
1299  real(RP), intent(inout) :: flux (KA,IA,JA)
1300  real(RP), intent(in) :: mom (KA,IA,JA)
1301  real(RP), intent(in) :: val (KA,IA,JA)
1302  real(RP), intent(in) :: DENS (KA,IA,JA)
1303  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1304  real(RP), intent(in) :: J13G (KA,IA,JA)
1305  real(RP), intent(in) :: MAPF ( IA,JA,2)
1306  real(RP), intent(in) :: CDZ (KA)
1307  integer, intent(in) :: IIS, IIE, JJS, JJE
1308 
1309  real(RP) :: vel
1310  integer :: k, i, j
1311  !---------------------------------------------------------------------------
1312 
1313  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1314  !$omp private(vel) &
1315  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J13G,MAPF) &
1316  !$omp shared(GSQRT,CDZ)
1317  do j = jjs, jje
1318  do i = iis, iie
1319  do k = ks+1, ke-2
1320  vel = ( f2h(k,1,i_xvz) &
1321  * 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) ) &
1322  + f2h(k,2,i_xvz) &
1323  * 0.25_rp * ( mom(k,i,j)+mom(k,i-1,j)+mom(k,i,j+1)+mom(k,i-1,j+1) ) ) &
1324  / ( f2h(k,1,i_xvz) &
1325  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
1326  + f2h(k,2,i_xvz) &
1327  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
1328  flux(k,i,j) = j13g(k,i,j) / mapf(i,j,+2) * vel &
1329  * ( ( val(k,i,j) &
1330  + 0.5_rp * phi(val(k+1,i,j),val(k,i,j),val(k-1,i,j)) * ( val(k,i,j)-val(k-1,i,j) ) ) &
1331  * ( 0.5_rp + sign(0.5_rp,vel) ) &
1332  + ( val(k+1,i,j) &
1333  + 0.5_rp * phi(val(k,i,j),val(k+1,i,j),val(k+2,i,j)) * ( val(k+1,i,j)-val(k+2,i,j) ) ) &
1334  * ( 0.5_rp - sign(0.5_rp,vel) ) )
1335  enddo
1336  enddo
1337  enddo
1338 
1339  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
1340  !$omp private(vel) &
1341  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J13G,MAPF) &
1342  !$omp shared(GSQRT,CDZ)
1343  do j = jjs, jje
1344  do i = iis, iie
1345  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
1346  ! The flux at KS-1 can be non-zero.
1347  ! To reduce calculations, all the fluxes are set to zero.
1348  flux(ks-1,i,j) = 0.0_rp
1349 
1350  vel = ( f2h(ks,1,i_xvz) &
1351  * 0.25_rp * ( mom(ks+1,i,j)+mom(ks+1,i-1,j)+mom(ks+1,i,j+1)+mom(ks+1,i-1,j+1) ) &
1352  + f2h(ks,2,i_xvz) &
1353  * 0.25_rp * ( mom(ks,i,j)+mom(ks,i-1,j)+mom(ks,i,j+1)+mom(ks,i-1,j+1) ) ) &
1354  / ( f2h(ks,1,i_xvz) &
1355  * 0.5_rp * ( dens(ks+1,i,j)+dens(ks+1,i,j+1) ) &
1356  + f2h(ks,2,i_xvz) &
1357  * 0.5_rp * ( dens(ks,i,j)+dens(ks,i,j+1) ) )
1358  flux(ks,i,j) = j13g(ks,i,j) / mapf(i,j,+2) * vel &
1359  * ( f2 * ( val(ks+1,i,j)+val(ks,i,j) ) )
1360 
1361  vel = ( f2h(ke-1,1,i_xvz) &
1362  * 0.25_rp * ( mom(ke,i,j)+mom(ke,i-1,j)+mom(ke,i,j+1)+mom(ke,i-1,j+1) ) &
1363  + f2h(ke-1,2,i_xvz) &
1364  * 0.25_rp * ( mom(ke-1,i,j)+mom(ke-1,i-1,j)+mom(ke-1,i,j+1)+mom(ke-1,i-1,j+1) ) ) &
1365  / ( f2h(ke-1,1,i_xvz) &
1366  * 0.5_rp * ( dens(ke,i,j)+dens(ke,i,j+1) ) &
1367  + f2h(ke-1,2,i_xvz) &
1368  * 0.5_rp * ( dens(ke-1,i,j)+dens(ke-1,i,j+1) ) )
1369  flux(ke-1,i,j) = j13g(ke-1,i,j) / mapf(i,j,+2) * vel &
1370  * ( f2 * ( val(ke,i,j)+val(ke-1,i,j) ) )
1371 
1372  flux(ke ,i,j) = 0.0_rp
1373  enddo
1374  enddo
1375 
1376  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_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxj23_xvz_ud3koren1993 ( 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 1387 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

1387  implicit none
1388 
1389  real(RP), intent(inout) :: flux (KA,IA,JA)
1390  real(RP), intent(in) :: mom (KA,IA,JA)
1391  real(RP), intent(in) :: val (KA,IA,JA)
1392  real(RP), intent(in) :: DENS (KA,IA,JA)
1393  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1394  real(RP), intent(in) :: J23G (KA,IA,JA)
1395  real(RP), intent(in) :: MAPF ( IA,JA,2)
1396  real(RP), intent(in) :: CDZ (KA)
1397  integer, intent(in) :: IIS, IIE, JJS, JJE
1398 
1399  real(RP) :: vel
1400  integer :: k, i, j
1401  !---------------------------------------------------------------------------
1402 
1403  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1404  !$omp private(vel) &
1405  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF) &
1406  !$omp shared(GSQRT,CDZ)
1407  do j = jjs, jje
1408  do i = iis, iie
1409  do k = ks+1, ke-2
1410  vel = ( f2h(k,1,i_xvz) &
1411  * mom(k+1,i,j) &
1412  + f2h(k,2,i_xvz) &
1413  * mom(k,i,j) ) &
1414  / ( f2h(k,1,i_xvz) &
1415  * 0.5_rp * ( dens(k+1,i,j)+dens(k+1,i,j+1) ) &
1416  + f2h(k,2,i_xvz) &
1417  * 0.5_rp * ( dens(k,i,j)+dens(k,i,j+1) ) )
1418  flux(k,i,j) = j23g(k,i,j) / mapf(i,j,+1) * vel &
1419  * ( ( val(k,i,j) &
1420  + 0.5_rp * phi(val(k+1,i,j),val(k,i,j),val(k-1,i,j)) * ( val(k,i,j)-val(k-1,i,j) ) ) &
1421  * ( 0.5_rp + sign(0.5_rp,vel) ) &
1422  + ( val(k+1,i,j) &
1423  + 0.5_rp * phi(val(k,i,j),val(k+1,i,j),val(k+2,i,j)) * ( val(k+1,i,j)-val(k+2,i,j) ) ) &
1424  * ( 0.5_rp - sign(0.5_rp,vel) ) )
1425  enddo
1426  enddo
1427  enddo
1428 
1429  !$omp parallel do default(none) private(i,j) OMP_SCHEDULE_ collapse(2) &
1430  !$omp private(vel) &
1431  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,J23G,MAPF) &
1432  !$omp shared(GSQRT,CDZ)
1433  do j = jjs, jje
1434  do i = iis, iie
1435  ! The boundary condition is qflx_hi + qflxJ13 + qfluxJ23 = 0 at KS-1.
1436  ! The flux at KS-1 can be non-zero.
1437  ! To reduce calculations, all the fluxes are set to zero.
1438  flux(ks-1,i,j) = 0.0_rp
1439 
1440  vel = ( f2h(ks,1,i_xvz) &
1441  * mom(ks+1,i,j) &
1442  + f2h(ks,2,i_xvz) &
1443  * mom(ks,i,j) ) &
1444  / ( f2h(ks,1,i_xvz) &
1445  * 0.5_rp * ( dens(ks+1,i,j)+dens(ks+1,i,j+1) ) &
1446  + f2h(ks,2,i_xvz) &
1447  * 0.5_rp * ( dens(ks,i,j)+dens(ks,i,j+1) ) )
1448  flux(ks,i,j) = j23g(ks,i,j) / mapf(i,j,+1) * vel &
1449  * ( f2 * ( val(ks+1,i,j)+val(ks,i,j) ) )
1450 
1451  vel = ( f2h(ke-1,1,i_xvz) &
1452  * mom(ke,i,j) &
1453  + f2h(ke-1,2,i_xvz) &
1454  * mom(ke-1,i,j) ) &
1455  / ( f2h(ke-1,1,i_xvz) &
1456  * 0.5_rp * ( dens(ke,i,j)+dens(ke,i,j+1) ) &
1457  + f2h(ke-1,2,i_xvz) &
1458  * 0.5_rp * ( dens(ke-1,i,j)+dens(ke-1,i,j+1) ) )
1459  flux(ke-1,i,j) = j23g(ke-1,i,j) / mapf(i,j,+1) * vel &
1460  * ( f2 * ( val(ke,i,j)+val(ke-1,i,j) ) )
1461 
1462  flux(ke ,i,j) = 0.0_rp
1463  enddo
1464  enddo
1465 
1466  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_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxx_xvz_ud3koren1993 ( 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 1478 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.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().

1478  implicit none
1479 
1480  real(RP), intent(inout) :: flux (KA,IA,JA)
1481  real(RP), intent(in) :: mom (KA,IA,JA)
1482  real(RP), intent(in) :: val (KA,IA,JA)
1483  real(RP), intent(in) :: DENS (KA,IA,JA)
1484  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1485  real(RP), intent(in) :: MAPF ( IA,JA,2)
1486  real(RP), intent(in) :: num_diff(KA,IA,JA)
1487  real(RP), intent(in) :: CDZ (KA)
1488  integer, intent(in) :: IIS, IIE, JJS, JJE
1489 
1490  real(RP) :: vel
1491  integer :: k, i, j
1492  !---------------------------------------------------------------------------
1493 
1494  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1495  !$omp private(vel) &
1496  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff)
1497  do j = jjs, jje
1498  do i = iis-1, iie
1499  do k = ks, ke
1500 #ifdef DEBUG
1501  call check( __line__, mom(k,i ,j) )
1502  call check( __line__, mom(k,i,j-1) )
1503 
1504  call check( __line__, val(k,i,j) )
1505  call check( __line__, val(k,i+1,j) )
1506 
1507  call check( __line__, val(k,i-1,j) )
1508  call check( __line__, val(k,i+2,j) )
1509 
1510 #endif
1511  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j+1) ) ) &
1512  / ( 0.25_rp * ( dens(k,i,j)+dens(k,i+1,j)+dens(k,i,j+1)+dens(k,i+1,j+1) ) )
1513  flux(k,i,j) = gsqrt(k,i,j) / mapf(i,j,+2) * vel &
1514  * ( ( val(k,i,j) &
1515  + 0.5_rp * phi(val(k,i+1,j),val(k,i,j),val(k,i-1,j)) * ( val(k,i,j)-val(k,i-1,j) ) ) &
1516  * ( 0.5_rp + sign(0.5_rp,vel) ) &
1517  + ( val(k,i+1,j) &
1518  + 0.5_rp * phi(val(k,i,j),val(k,i+1,j),val(k,i+2,j)) * ( val(k,i+1,j)-val(k,i+2,j) ) ) &
1519  * ( 0.5_rp - sign(0.5_rp,vel) ) ) &
1520  + gsqrt(k,i,j) * num_diff(k,i,j)
1521  enddo
1522  enddo
1523  enddo
1524 #ifdef DEBUG
1525  k = iundef; i = iundef; j = iundef
1526 #endif
1527 
1528  return
Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_dyn_fvm_fluxy_xvz_ud3koren1993()

subroutine, public scale_atmos_dyn_fvm_flux_ud3koren1993::atmos_dyn_fvm_fluxy_xvz_ud3koren1993 ( 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 1540 of file scale_atmos_dyn_fvm_flux_ud3Koren1993.F90.

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

Referenced by scale_atmos_dyn_fvm_flux::atmos_dyn_fvm_flux_setup().

1540  implicit none
1541 
1542  real(RP), intent(inout) :: flux (KA,IA,JA)
1543  real(RP), intent(in) :: mom (KA,IA,JA)
1544  real(RP), intent(in) :: val (KA,IA,JA)
1545  real(RP), intent(in) :: DENS (KA,IA,JA)
1546  real(RP), intent(in) :: GSQRT (KA,IA,JA)
1547  real(RP), intent(in) :: MAPF ( IA,JA,2)
1548  real(RP), intent(in) :: num_diff(KA,IA,JA)
1549  real(RP), intent(in) :: CDZ (KA)
1550  integer, intent(in) :: IIS, IIE, JJS, JJE
1551 
1552  real(RP) :: vel
1553  integer :: k, i, j
1554  !---------------------------------------------------------------------------
1555 
1556  ! note that y-index is added by -1
1557 
1558  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1559  !$omp private(vel) &
1560  !$omp shared(JJS,JJE,IIS,IIE,KS,KE,mom,val,DENS,flux,GSQRT,MAPF,num_diff)
1561  do j = jjs, jje+1
1562  do i = iis, iie
1563  do k = ks, ke
1564 #ifdef DEBUG
1565  call check( __line__, mom(k,i ,j) )
1566  call check( __line__, mom(k,i,j-1) )
1567 
1568  call check( __line__, val(k,i,j-1) )
1569  call check( __line__, val(k,i,j) )
1570 
1571  call check( __line__, val(k,i,j-2) )
1572  call check( __line__, val(k,i,j+1) )
1573 
1574 #endif
1575  vel = ( 0.5_rp * ( mom(k,i,j)+mom(k,i,j-1) ) ) &
1576  / ( dens(k,i,j) )
1577  flux(k,i,j-1) = gsqrt(k,i,j) / mapf(i,j,+1) * vel &
1578  * ( ( val(k,i,j-1) &
1579  + 0.5_rp * phi(val(k,i,j),val(k,i,j-1),val(k,i,j-2)) * ( val(k,i,j-1)-val(k,i,j-2) ) ) &
1580  * ( 0.5_rp + sign(0.5_rp,vel) ) &
1581  + ( val(k,i,j) &
1582  + 0.5_rp * phi(val(k,i,j-1),val(k,i,j),val(k,i,j+1)) * ( val(k,i,j)-val(k,i,j+1) ) ) &
1583  * ( 0.5_rp - sign(0.5_rp,vel) ) ) &
1584  + gsqrt(k,i,j) * num_diff(k,i,j)
1585  enddo
1586  enddo
1587  enddo
1588 #ifdef DEBUG
1589  k = iundef; i = iundef; j = iundef
1590 #endif
1591 
1592  return
Here is the call graph for this function:
Here is the caller graph for this function:
Generated by   doxygen 1.8.14