SCALE-RM
scale_grid_cartesian.F90
Go to the documentation of this file.
1 !-------------------------------------------------------------------------------
16 !-------------------------------------------------------------------------------
17 module scale_grid
18  !-----------------------------------------------------------------------------
19  !
20  !++ used modules
21  !
22  use scale_precision
23  use scale_stdio
24  use scale_prof
25  use scale_grid_index
26  !-----------------------------------------------------------------------------
27  implicit none
28  private
29  !-----------------------------------------------------------------------------
30  !
31  !++ Public procedure
32  !
33  public :: grid_setup
34  public :: grid_allocate
35  public :: grid_generate
36 
37  !-----------------------------------------------------------------------------
38  !
39  !++ Public parameters & variables
40  !
41  real(RP), public :: dz = 500.0_rp
42  real(RP), public :: dx = 500.0_rp
43  real(RP), public :: dy = 500.0_rp
44 
45  real(RP), public :: buffer_dz = 0.0_rp
46  real(RP), public :: buffer_dx = 0.0_rp
47  real(RP), public :: buffer_dy = 0.0_rp
48  real(RP), public :: bufffact = 1.0_rp
49  real(RP), public :: bufffact_x = -1.0_rp
50  real(RP), public :: bufffact_y = -1.0_rp
51  real(RP), public :: bufffact_z = -1.0_rp
52 
53  real(RP), public :: grid_domain_center_x
54  real(RP), public :: grid_domain_center_y
55 
56  real(RP), public, allocatable :: grid_cz (:)
57  real(RP), public, allocatable :: grid_cx (:)
58  real(RP), public, allocatable :: grid_cy (:)
59  real(RP), public, allocatable :: grid_cdz (:)
60  real(RP), public, allocatable :: grid_cdx (:)
61  real(RP), public, allocatable :: grid_cdy (:)
62  real(RP), public, allocatable :: grid_rcdz(:)
63  real(RP), public, allocatable :: grid_rcdx(:)
64  real(RP), public, allocatable :: grid_rcdy(:)
65 
66  real(RP), public, allocatable :: grid_fz (:)
67  real(RP), public, allocatable :: grid_fx (:)
68  real(RP), public, allocatable :: grid_fy (:)
69  real(RP), public, allocatable :: grid_fdz (:)
70  real(RP), public, allocatable :: grid_fdx (:)
71  real(RP), public, allocatable :: grid_fdy (:)
72  real(RP), public, allocatable :: grid_rfdz(:)
73  real(RP), public, allocatable :: grid_rfdx(:)
74  real(RP), public, allocatable :: grid_rfdy(:)
75 
76  real(RP), public, allocatable :: grid_cbfz(:)
77  real(RP), public, allocatable :: grid_cbfx(:)
78  real(RP), public, allocatable :: grid_cbfy(:)
79  real(RP), public, allocatable :: grid_fbfz(:)
80  real(RP), public, allocatable :: grid_fbfx(:)
81  real(RP), public, allocatable :: grid_fbfy(:)
82 
83  real(RP), public, allocatable :: grid_fxg (:)
84  real(RP), public, allocatable :: grid_fyg (:)
85  real(RP), public, allocatable :: grid_cxg (:)
86  real(RP), public, allocatable :: grid_cyg (:)
87  real(RP), public, allocatable :: grid_fdxg (:)
88  real(RP), public, allocatable :: grid_fdyg (:)
89  real(RP), public, allocatable :: grid_cdxg (:)
90  real(RP), public, allocatable :: grid_cdyg (:)
91  real(RP), public, allocatable :: grid_fbfxg(:)
92  real(RP), public, allocatable :: grid_fbfyg(:)
93  real(RP), public, allocatable :: grid_cbfxg(:)
94  real(RP), public, allocatable :: grid_cbfyg(:)
95 
96  !-----------------------------------------------------------------------------
97  !
98  !++ Private procedure
99  !
100  private :: grid_read
101 
102  !-----------------------------------------------------------------------------
103  !
104  !++ Private parameters & variables
105  !
106  character(len=H_LONG), private :: grid_in_basename = ''
107  character(len=H_LONG), private :: grid_out_basename = ''
108  real(RP), private :: grid_offset_x = 0.0_rp
109  real(RP), private :: grid_offset_y = 0.0_rp
110 
111  integer, private :: buffer_nz = -1
112  integer, private :: buffer_nx = -1
113  integer, private :: buffer_ny = -1
114 
115  integer, private, parameter :: kmax_user_lim = 300
116  real(RP), private :: fz(kmax_user_lim)
117 
118  logical, private :: debug = .false.
119 
120  !-----------------------------------------------------------------------------
121 contains
122  !-----------------------------------------------------------------------------
124  subroutine grid_setup
125  use scale_process, only: &
126  prc_mpistop
127  implicit none
128 
129  namelist / param_grid / &
130  grid_in_basename, &
131  grid_out_basename, &
132  grid_offset_x, &
133  grid_offset_y, &
134  dx, &
135  dy, &
136  dz, &
137  buffer_dz, &
138  buffer_dx, &
139  buffer_dy, &
140  buffer_nz, &
141  buffer_nx, &
142  buffer_ny, &
143  bufffact, &
144  bufffact_x, &
145  bufffact_y, &
146  bufffact_z, &
147  fz, &
148  debug
149 
150  integer :: ierr
151  !---------------------------------------------------------------------------
152 
153  if( io_l ) write(io_fid_log,*)
154  if( io_l ) write(io_fid_log,*) '++++++ Module[GRID] / Categ[ATMOS-RM GRID] / Origin[SCALElib]'
155 
156  call grid_allocate
157 
158  !--- read namelist
159  rewind(io_fid_conf)
160  read(io_fid_conf,nml=param_grid,iostat=ierr)
161  if( ierr < 0 ) then !--- missing
162  if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
163  elseif( ierr > 0 ) then !--- fatal error
164  write(*,*) 'xxx Not appropriate names in namelist PARAM_GRID. Check!'
165  call prc_mpistop
166  endif
167  if( io_nml ) write(io_fid_nml,nml=param_grid)
168 
169  if ( bufffact_x < 0.0_rp ) bufffact_x = bufffact
170  if ( bufffact_y < 0.0_rp ) bufffact_y = bufffact
171  if ( bufffact_z < 0.0_rp ) bufffact_z = bufffact
172 
173  if( io_l ) write(io_fid_log,*)
174  if( io_l ) write(io_fid_log,*) '*** Atmosphere grid information ***'
175  if( io_l ) write(io_fid_log,'(1x,A,3(1x,F9.3))') '*** delta Z, X, Y [m] :', dz, dx, dy
176 
177  if ( grid_in_basename /= '' ) then
178  call grid_read
179  else
180  if( io_l ) write(io_fid_log,*) '*** Not found input grid file. Grid position is calculated.'
181 
182  call grid_generate
183  endif
184 
185  if( io_l ) write(io_fid_log,*)
186  if( io_l ) write(io_fid_log,*) '*** Domain size [km] (local) :'
187  if( io_l ) write(io_fid_log,'(1x,6(A,F9.3))') '*** X:', &
188  grid_fx(0) *1.e-3_rp, ' -HALO- ', grid_fx(is-1)*1.e-3_rp, ' | ', &
189  grid_cx(is)*1.e-3_rp, ' - ', grid_cx(ie) *1.e-3_rp, ' | ', &
190  grid_fx(ie)*1.e-3_rp, ' -HALO- ', grid_fx(ia) *1.e-3_rp
191  if( io_l ) write(io_fid_log,'(1x,6(A,F9.3))') '*** Y:', &
192  grid_fy(0) *1.e-3_rp, ' -HALO- ', grid_fy(js-1)*1.e-3_rp, ' | ', &
193  grid_cy(js)*1.e-3_rp, ' - ', grid_cy(je) *1.e-3_rp, ' | ', &
194  grid_fy(je)*1.e-3_rp, ' -HALO- ', grid_fy(ja) *1.e-3_rp
195 
196  return
197  end subroutine grid_setup
198 
199  !-----------------------------------------------------------------------------
201  subroutine grid_allocate
202  use scale_rm_process, only: &
203  prc_num_x, &
204  prc_num_y
205  implicit none
206  !---------------------------------------------------------------------------
207 
208  ! working
209  fz(:) = -1.0_rp
210 
211  ! local domain
212  allocate( grid_cz(ka) )
213  allocate( grid_cx(ia) )
214  allocate( grid_cy(ja) )
215  allocate( grid_cdz(ka) )
216  allocate( grid_cdx(ia) )
217  allocate( grid_cdy(ja) )
218  allocate( grid_rcdz(ka) )
219  allocate( grid_rcdx(ia) )
220  allocate( grid_rcdy(ja) )
221 
222  allocate( grid_fz(0:ka) )
223  allocate( grid_fx(0:ia) )
224  allocate( grid_fy(0:ja) )
225  allocate( grid_fdz(ka-1) )
226  allocate( grid_fdx(ia-1) )
227  allocate( grid_fdy(ja-1) )
228  allocate( grid_rfdz(ka-1) )
229  allocate( grid_rfdx(ia-1) )
230  allocate( grid_rfdy(ja-1) )
231 
232  allocate( grid_cbfz(ka) )
233  allocate( grid_cbfx(ia) )
234  allocate( grid_cbfy(ja) )
235  allocate( grid_fbfz(ka) )
236  allocate( grid_fbfx(ia) )
237  allocate( grid_fbfy(ja) )
238 
239  ! global domain
240  allocate( grid_fxg(0:iag) )
241  allocate( grid_fyg(0:jag) )
242  allocate( grid_cxg( iag) )
243  allocate( grid_cyg( jag) )
244  allocate( grid_fdxg(iag-1) )
245  allocate( grid_fdyg(jag-1) )
246  allocate( grid_cdxg( iag) )
247  allocate( grid_cdyg( jag) )
248  allocate( grid_cbfxg( iag) )
249  allocate( grid_cbfyg( jag) )
250  allocate( grid_fbfxg( iag) )
251  allocate( grid_fbfyg( jag) )
252 
253  return
254  end subroutine grid_allocate
255 
256  !-----------------------------------------------------------------------------
258  subroutine grid_read
259  use gtool_file, only: &
260  fileread
261  use scale_process, only: &
262  prc_myrank
263  implicit none
264 
265  character(len=H_LONG) :: bname
266  !---------------------------------------------------------------------------
267 
268  if( io_l ) write(io_fid_log,*)
269  if( io_l ) write(io_fid_log,*) '*** Input grid file ***'
270 
271  write(bname,'(A,A,F15.3)') trim(grid_in_basename)
272 
273  call fileread( grid_cz(:), bname, 'CZ', 1, prc_myrank )
274  call fileread( grid_cx(:), bname, 'CX', 1, prc_myrank )
275  call fileread( grid_cy(:), bname, 'CY', 1, prc_myrank )
276  call fileread( grid_fz(:), bname, 'FZ', 1, prc_myrank )
277  call fileread( grid_fx(:), bname, 'FX', 1, prc_myrank )
278  call fileread( grid_fy(:), bname, 'FY', 1, prc_myrank )
279 
280  call fileread( grid_cxg(:), bname, 'CXG', 1, prc_myrank )
281  call fileread( grid_cyg(:), bname, 'CYG', 1, prc_myrank )
282  call fileread( grid_fxg(:), bname, 'FXG', 1, prc_myrank )
283  call fileread( grid_fyg(:), bname, 'FYG', 1, prc_myrank )
284 
285  call fileread( grid_cdz(:), bname, 'CDZ', 1, prc_myrank )
286  call fileread( grid_cdx(:), bname, 'CDX', 1, prc_myrank )
287  call fileread( grid_cdy(:), bname, 'CDY', 1, prc_myrank )
288  call fileread( grid_fdz(:), bname, 'FDZ', 1, prc_myrank )
289  call fileread( grid_fdx(:), bname, 'FDX', 1, prc_myrank )
290  call fileread( grid_fdy(:), bname, 'FDY', 1, prc_myrank )
291 
292  grid_rcdz(:) = 1.0_rp / grid_cdz(:)
293  grid_rcdx(:) = 1.0_rp / grid_cdx(:)
294  grid_rcdy(:) = 1.0_rp / grid_cdy(:)
295  grid_rfdz(:) = 1.0_rp / grid_fdz(:)
296  grid_rfdx(:) = 1.0_rp / grid_fdx(:)
297  grid_rfdy(:) = 1.0_rp / grid_fdy(:)
298 
299  call fileread( grid_cbfz(:), bname, 'CBFZ', 1, prc_myrank )
300  call fileread( grid_cbfx(:), bname, 'CBFX', 1, prc_myrank )
301  call fileread( grid_cbfy(:), bname, 'CBFY', 1, prc_myrank )
302  call fileread( grid_fbfz(:), bname, 'FBFZ', 1, prc_myrank )
303  call fileread( grid_fbfx(:), bname, 'FBFX', 1, prc_myrank )
304  call fileread( grid_fbfy(:), bname, 'FBFY', 1, prc_myrank )
305 
306  return
307  end subroutine grid_read
308 
309  !-----------------------------------------------------------------------------
311  subroutine grid_generate
312  use scale_process, only: &
313  prc_mpistop, &
314  prc_myrank
315  use scale_rm_process, only: &
316  prc_2drank, &
317  prc_num_x, &
318  prc_num_y
319  implicit none
320 
321  real(RP), allocatable :: buffz(:), buffx(:), buffy(:)
322  real(RP) :: bufftotz, bufftotx, bufftoty
323 
324  integer :: kbuff, ibuff, jbuff
325  integer :: kmain, imain, jmain
326 
327  logical :: use_user_input
328 
329  integer :: k, i, j, ii, jj
330  !---------------------------------------------------------------------------
331 
332  !##### coordinate in global domain #####
333 
334  allocate( buffx(0:iag) )
335  allocate( buffy(0:jag) )
336 
337  ! X-direction
338  ! calculate buffer grid size
339 
340  if ( buffer_nx > 0 ) then
341  if ( 2*buffer_nx > imaxg ) then
342  write(*,*) 'xxx Buffer grid size (', buffer_nx, &
343  'x2) must be smaller than global domain size (X). Use smaller BUFFER_NX!'
344  call prc_mpistop
345  endif
346 
347  buffx(0) = dx
348  bufftotx = 0.0_rp
349  do i = 1, buffer_nx
350  buffx(i) = buffx(i-1) * bufffact_x
351  bufftotx = bufftotx + buffx(i)
352  enddo
353  ibuff = buffer_nx
354  imain = imaxg - 2*buffer_nx
355 
356  buffer_dx = bufftotx
357  else
358  buffx(0) = dx
359  bufftotx = 0.0_rp
360  do i = 1, iag
361  if( bufftotx >= buffer_dx ) exit
362  buffx(i) = buffx(i-1) * bufffact_x
363  bufftotx = bufftotx + buffx(i)
364  enddo
365  ibuff = i - 1
366  imain = imaxg - 2*ibuff
367 
368  if ( imain < 0 ) then
369  write(*,*) 'xxx Buffer length (', bufftotx, &
370  'x2[m]) must be smaller than global domain size (X). Use smaller BUFFER_DX!'
371  call prc_mpistop
372  endif
373  endif
374 
375  ! horizontal coordinate (global domain)
376  grid_fxg(ihalo) = grid_offset_x
377  do i = ihalo-1, 0, -1
378  grid_fxg(i) = grid_fxg(i+1) - buffx(ibuff)
379  enddo
380 
381  do i = 1, ihalo
382  grid_cxg(i) = 0.5_rp * ( grid_fxg(i)+grid_fxg(i-1) )
383  enddo
384 
385  if ( ibuff > 0 ) then
386  do i = ihalo+1, ihalo+ibuff
387  grid_fxg(i) = grid_fxg(i-1) + buffx(ibuff+ihalo+1-i)
388  grid_cxg(i) = 0.5_rp * ( grid_fxg(i)+grid_fxg(i-1) )
389  enddo
390  endif
391 
392  do i = ihalo+ibuff+1, ihalo+ibuff+imain
393  grid_fxg(i) = grid_fxg(i-1) + dx
394  grid_cxg(i) = 0.5_rp * ( grid_fxg(i)+grid_fxg(i-1) )
395  enddo
396 
397  if ( ibuff > 0 ) then
398  do i = ihalo+ibuff+imain+1, ihalo+ibuff+imain+ibuff
399  grid_fxg(i) = grid_fxg(i-1) + buffx(i-ihalo-ibuff-imain)
400  grid_cxg(i) = 0.5_rp * ( grid_fxg(i)+grid_fxg(i-1) )
401  enddo
402  endif
403 
404  do i = ihalo+ibuff+imain+ibuff+1, ihalo+ibuff+imain+ibuff+ihalo
405  grid_fxg(i) = grid_fxg(i-1) + buffx(ibuff)
406  grid_cxg(i) = 0.5_rp * ( grid_fxg(i)+grid_fxg(i-1) )
407  enddo
408 
409  do i = 1, iag
410  grid_cdxg(i) = grid_fxg(i) - grid_fxg(i-1)
411  end do
412  do i = 1, iag-1
413  grid_fdxg(i) = grid_cxg(i+1)-grid_cxg(i)
414  end do
415 
416  ! calc buffer factor (global domain)
417  grid_cbfxg(:) = 0.0_rp
418  grid_fbfxg(:) = 0.0_rp
419  do i = 1, ihalo
420  grid_cbfxg(i) = 1.0_rp
421  grid_fbfxg(i) = 1.0_rp
422  enddo
423 
424  if ( ibuff > 0 ) then
425  do i = ihalo+1, ihalo+ibuff
426  grid_cbfxg(i) = (bufftotx+grid_fxg(ihalo)-grid_cxg(i)) / bufftotx
427  grid_fbfxg(i) = (bufftotx+grid_fxg(ihalo)-grid_fxg(i)) / bufftotx
428  enddo
429 
430  do i = ihalo+ibuff+imain+1, ihalo+ibuff+imain+ibuff
431  grid_cbfxg(i) = (bufftotx-grid_fxg(iag-ihalo)+grid_cxg(i)) / bufftotx
432  grid_fbfxg(i) = (bufftotx-grid_fxg(iag-ihalo)+grid_fxg(i)) / bufftotx
433  enddo
434  endif
435 
436  do i = ihalo+ibuff+imain+ibuff+1, ihalo+ibuff+imain+ibuff+ihalo
437  grid_cbfxg(i) = 1.0_rp
438  grid_fbfxg(i) = 1.0_rp
439  enddo
440 
441  grid_cbfxg(:) = max( min( grid_cbfxg(:), 1.0_rp ), 0.0_rp )
442  grid_fbfxg(:) = max( min( grid_fbfxg(:), 1.0_rp ), 0.0_rp )
443 
444  ! Y-direction
445  ! calculate buffer grid size
446 
447  if ( buffer_ny > 0 ) then
448  if ( 2*buffer_ny > jmaxg ) then
449  write(*,*) 'xxx Buffer grid size (', buffer_ny, &
450  'x2) must be smaller than global domain size (Y). Use smaller BUFFER_NY!'
451  call prc_mpistop
452  endif
453 
454  buffy(0) = dy
455  bufftoty = 0.0_rp
456  do j = 1, buffer_ny
457  buffy(j) = buffy(j-1) * bufffact_y
458  bufftoty = bufftoty + buffy(j)
459  enddo
460  jbuff = buffer_ny
461  jmain = jmaxg - 2*buffer_ny
462 
463  buffer_dy = bufftoty
464  else
465  buffy(0) = dy
466  bufftoty = 0.0_rp
467  do j = 1, jag
468  if( bufftoty >= buffer_dy ) exit
469  buffy(j) = buffy(j-1) * bufffact_y
470  bufftoty = bufftoty + buffy(j)
471  enddo
472  jbuff = j - 1
473  jmain = jmaxg - 2*jbuff
474 
475  if ( jmain < 0 ) then
476  write(*,*) 'xxx Buffer length (', bufftoty, &
477  'x2[m]) must be smaller than global domain size (Y). Use smaller BUFFER_DY!'
478  call prc_mpistop
479  endif
480  endif
481 
482  ! horizontal coordinate (global domain)
483  grid_fyg(jhalo) = grid_offset_y
484  do j = jhalo-1, 0, -1
485  grid_fyg(j) = grid_fyg(j+1) - buffy(jbuff)
486  enddo
487 
488  do j = 1, jhalo
489  grid_cyg(j) = 0.5_rp * ( grid_fyg(j)+grid_fyg(j-1) )
490  enddo
491 
492  if ( jbuff > 0 ) then
493  do j = jhalo+1, jhalo+jbuff
494  grid_fyg(j) = grid_fyg(j-1) + buffy(jbuff+jhalo+1-j)
495  grid_cyg(j) = 0.5_rp * ( grid_fyg(j)+grid_fyg(j-1) )
496  enddo
497  endif
498 
499  do j = jhalo+jbuff+1, jhalo+jbuff+jmain
500  grid_fyg(j) = grid_fyg(j-1) + dy
501  grid_cyg(j) = 0.5_rp * ( grid_fyg(j)+grid_fyg(j-1) )
502  enddo
503 
504  if ( jbuff > 0 ) then
505  do j = jhalo+jbuff+jmain+1, jhalo+jbuff+jmain+jbuff
506  grid_fyg(j) = grid_fyg(j-1) + buffy(j-jhalo-jbuff-jmain)
507  grid_cyg(j) = 0.5_rp * ( grid_fyg(j)+grid_fyg(j-1) )
508  enddo
509  endif
510 
511  do j = jhalo+jbuff+jmain+jbuff+1, jhalo+jbuff+jmain+jbuff+jhalo
512  grid_fyg(j) = grid_fyg(j-1) + buffy(jbuff)
513  grid_cyg(j) = 0.5_rp * ( grid_fyg(j)+grid_fyg(j-1) )
514  enddo
515 
516  do j = 1, jag
517  grid_cdyg(j) = grid_fyg(j) - grid_fyg(j-1)
518  end do
519  do j = 1, jag-1
520  grid_fdyg(j) = grid_cyg(j+1)-grid_cyg(j)
521  end do
522 
523  ! calc buffer factor (global domain)
524  grid_cbfyg(:) = 0.0_rp
525  grid_fbfyg(:) = 0.0_rp
526  do j = 1, jhalo
527  grid_cbfyg(j) = 1.0_rp
528  grid_fbfyg(j) = 1.0_rp
529  enddo
530 
531  if ( jbuff > 0 ) then
532  do j = jhalo+1, jhalo+jbuff
533  grid_cbfyg(j) = (bufftoty+grid_fyg(jhalo)-grid_cyg(j)) / bufftoty
534  grid_fbfyg(j) = (bufftoty+grid_fyg(jhalo)-grid_fyg(j)) / bufftoty
535  enddo
536 
537  do j = jhalo+jbuff+jmain+1, jhalo+jbuff+jmain+jbuff
538  grid_cbfyg(j) = (bufftoty-grid_fyg(jag-jhalo)+grid_cyg(j)) / bufftoty
539  grid_fbfyg(j) = (bufftoty-grid_fyg(jag-jhalo)+grid_fyg(j)) / bufftoty
540  enddo
541  endif
542 
543  do j = jhalo+jbuff+jmain+jbuff+1, jhalo+jbuff+jmain+jbuff+jhalo
544  grid_cbfyg(j) = 1.0_rp
545  grid_fbfyg(j) = 1.0_rp
546  enddo
547  grid_cbfyg(:) = max( min( grid_cbfyg(:), 1.0_rp ), 0.0_rp )
548  grid_fbfyg(:) = max( min( grid_fbfyg(:), 1.0_rp ), 0.0_rp )
549 
550  deallocate( buffx )
551  deallocate( buffy )
552 
553  !##### coordinate in local domain #####
554 
555  allocate( buffz(0:ka) )
556 
557  use_user_input = .false.
558  if ( maxval(fz(1:kmax_user_lim)) > 0.0_rp ) then ! try to use input from namelist
559  if( io_l ) write(io_fid_log,*) '*** Z coordinate is given from NAMELIST.'
560 
561  if ( kmax < 2 ) then
562  write(*,*) 'xxx KMAX must be larger than 1. Check!', kmax
563  call prc_mpistop
564  endif
565 
566  if ( kmax > kmax_user_lim ) then
567  write(*,*) 'xxx KMAX must be smaller than ', kmax_user_lim, '. Check!', kmax
568  call prc_mpistop
569  endif
570 
571  if ( minval(fz(1:kmax)) <= 0.0_rp ) then
572  write(*,*) 'xxx FZ must be positive. Check! minval(FZ(1:KMAX))=', minval(fz(1:kmax))
573  call prc_mpistop
574  endif
575 
576  use_user_input = .true.
577  endif
578 
579  if ( use_user_input ) then ! input from namelist
580 
581  ! Z-direction
582  ! calculate buffer grid size
583 
584  if ( buffer_nz > 0 ) then
585  if ( buffer_nz > kmax ) then
586  write(*,*) 'xxx Buffer grid size (', buffer_nz, &
587  ') must be smaller than global domain size (Z). Use smaller BUFFER_NZ!'
588  call prc_mpistop
589  endif
590 
591  bufftotz = 0.0_rp
592  do k = kmax, kmax-buffer_nz+1, -1
593  bufftotz = bufftotz + ( fz(k) - fz(k-1) )
594  enddo
595  kbuff = buffer_nz
596  kmain = kmax - buffer_nz
597 
598  buffer_dz = bufftotz
599  else
600  if ( buffer_dz > fz(kmax) ) then
601  write(*,*) 'xxx Buffer length (', buffer_dz, &
602  '[m]) must be smaller than global domain size (Z). Use smaller BUFFER_DZ!'
603  call prc_mpistop
604  endif
605 
606  bufftotz = 0.0_rp
607  do k = kmax, 2, -1
608  if( bufftotz >= buffer_dz ) exit
609  bufftotz = bufftotz + ( fz(k) - fz(k-1) )
610  enddo
611  kbuff = kmax - k
612  kmain = k
613  endif
614 
615  ! vertical coordinate (local=global domain)
616  grid_fz(ks-1) = 0.0_rp
617 
618  dz = fz(1)
619  do k = ks-2, 0, -1
620  grid_fz(k) = grid_fz(k+1) - dz
621  enddo
622 
623  do k = ks, ke
624  grid_fz(k) = fz(k-ks+1)
625  enddo
626 
627  dz = fz(kmax) - fz(kmax-1)
628  do k = ke+1, ka
629  grid_fz(k) = grid_fz(k-1) + dz
630  enddo
631 
632  do k = 1, ka
633  grid_cz(k) = 0.5_rp * ( grid_fz(k)+grid_fz(k-1) )
634  enddo
635 
636  else ! calc using DZ
637 
638  ! Z-direction
639  ! calculate buffer grid size
640 
641  if ( buffer_nz > 0 ) then
642  if ( buffer_nz > kmax ) then
643  write(*,*) 'xxx Buffer grid size (', buffer_nz, &
644  ') must be smaller than global domain size (Z). Use smaller BUFFER_NZ!'
645  call prc_mpistop
646  endif
647 
648  buffz(0) = dz
649  bufftotz = 0.0_rp
650  do k = 1, buffer_nz
651  buffz(k) = buffz(k-1) * bufffact_z
652  bufftotz = bufftotz + buffz(k)
653  enddo
654  kbuff = buffer_nz
655  kmain = kmax - buffer_nz
656 
657  buffer_dz = bufftotz
658  else
659  buffz(0) = dz
660  bufftotz = 0.0_rp
661  do k = 1, ka
662  if( bufftotz >= buffer_dz ) exit
663  buffz(k) = buffz(k-1) * bufffact_z
664  bufftotz = bufftotz + buffz(k)
665  enddo
666  kbuff = k - 1
667  kmain = kmax - kbuff
668 
669  if ( kmain < 0 ) then
670  write(*,*) 'xxx Buffer length (', bufftotz, &
671  '[m]) must be smaller than global domain size (Z). Use smaller BUFFER_DZ!'
672  call prc_mpistop
673  endif
674  endif
675 
676  ! vertical coordinate (local=global domain)
677  grid_fz(ks-1) = 0.0_rp
678  do k = ks-2, 0, -1
679  grid_fz(k) = grid_fz(k+1) - dz
680  enddo
681 
682  do k = 1, ks-1
683  grid_cz(k) = 0.5_rp * ( grid_fz(k)+grid_fz(k-1) )
684  enddo
685 
686  do k = ks, ks+kmain-1
687  grid_fz(k) = grid_fz(k-1) + dz
688  grid_cz(k) = 0.5_rp * ( grid_fz(k)+grid_fz(k-1) )
689  enddo
690 
691  if ( kbuff > 0 ) then
692  do k = ks+kmain, ke
693  grid_fz(k) = grid_fz(k-1) + buffz(k-ks-kmain+1)
694  grid_cz(k) = 0.5_rp * ( grid_fz(k)+grid_fz(k-1) )
695  enddo
696  endif
697 
698  do k = ke+1, ka
699  grid_fz(k) = grid_fz(k-1) + buffz(kbuff)
700  grid_cz(k) = 0.5_rp * ( grid_fz(k)+grid_fz(k-1) )
701  enddo
702 
703  endif
704 
705  ! calc buffer factor (global domain)
706  grid_cbfz(:) = 0.0_rp
707  grid_fbfz(:) = 0.0_rp
708  if ( kbuff > 0 ) then
709  do k = ks+kmain, ke
710  grid_cbfz(k) = (bufftotz-grid_fz(ke)+grid_cz(k)) / bufftotz
711  grid_fbfz(k) = (bufftotz-grid_fz(ke)+grid_fz(k)) / bufftotz
712  enddo
713  endif
714 
715  do k = ke+1, ka
716  grid_cbfz(k) = 1.0_rp
717  grid_fbfz(k) = 1.0_rp
718  enddo
719  grid_cbfz(:) = max( min( grid_cbfz(:), 1.0_rp ), 0.0_rp )
720  grid_fbfz(:) = max( min( grid_fbfz(:), 1.0_rp ), 0.0_rp )
721 
722  deallocate( buffz )
723 
724  ! vertical coordinate (local domain)
725  do k = 1, ka
726  grid_cdz(k) = grid_fz(k) - grid_fz(k-1)
727  grid_rcdz(k) = 1.0_rp / grid_cdz(k)
728  enddo
729 
730  do k = 1, ka-1
731  grid_fdz(k) = grid_cz(k+1)-grid_cz(k)
732  grid_rfdz(k) = 1.0_rp / grid_fdz(k)
733  enddo
734 
735  ! X-direction
736  ! horizontal coordinate (local domain)
737  do i = 0, ia
738  ii = i + prc_2drank(prc_myrank,1) * imax
739 
740  grid_fx(i) = grid_fxg(ii)
741  enddo
742 
743  do i = 1, ia
744  ii = i + prc_2drank(prc_myrank,1) * imax
745 
746  grid_cx(i) = grid_cxg(ii)
747  grid_cbfx(i) = grid_cbfxg(ii)
748  grid_fbfx(i) = grid_fbfxg(ii)
749 
750  grid_cdx(i) = grid_fx(i) - grid_fx(i-1)
751  grid_rcdx(i) = 1.0_rp / grid_cdx(i)
752  enddo
753 
754  do i = 1, ia-1
755  grid_fdx(i) = grid_cx(i+1)-grid_cx(i)
756  grid_rfdx(i) = 1.0_rp / grid_fdx(i)
757  enddo
758 
759  ! Y-direction
760  ! horizontal coordinate (local domain)
761  do j = 0, ja
762  jj = j + prc_2drank(prc_myrank,2) * jmax
763 
764  grid_fy(j) = grid_fyg(jj)
765  enddo
766 
767  do j = 1, ja
768  jj = j + prc_2drank(prc_myrank,2) * jmax
769 
770  grid_cy(j) = grid_cyg(jj)
771  grid_cbfy(j) = grid_cbfyg(jj)
772  grid_fbfy(j) = grid_fbfyg(jj)
773 
774  grid_cdy(j) = grid_fy(j) - grid_fy(j-1)
775  grid_rcdy(j) = 1.0_rp / grid_cdy(j)
776  enddo
777 
778  do j = 1, ja-1
779  grid_fdy(j) = grid_cy(j+1)-grid_cy(j)
780  grid_rfdy(j) = 1.0_rp / grid_fdy(j)
781  enddo
782 
783  grid_domain_center_x = 0.5_rp * ( grid_fxg(ihalo) + grid_fxg(iag-ihalo) )
784  grid_domain_center_y = 0.5_rp * ( grid_fyg(jhalo) + grid_fyg(jag-jhalo) )
785 
786  ! report
787  if( io_l ) write(io_fid_log,*)
788  if( io_l ) write(io_fid_log,*) '*** Main/buffer Grid (global) :'
789  if( io_l ) write(io_fid_log,'(1x,2(A,I6))') '*** Z: buffer = ', kbuff,' x 1, main = ',kmain
790  if( io_l ) write(io_fid_log,'(1x,2(A,I6))') '*** X: buffer = ', ibuff,' x 2, main = ',imain
791  if( io_l ) write(io_fid_log,'(1x,2(A,I6))') '*** Y: buffer = ', jbuff,' x 2, main = ',jmain
792  if( io_l ) write(io_fid_log,*)
793  if( io_l ) write(io_fid_log,*) '*** Domain size [km] (global) :'
794  if( io_l ) write(io_fid_log,'(1x,7(A,F9.3))') '*** Z:', &
795  grid_fz(0) *1.e-3_rp, ' -HALO- ', &
796  grid_fz(ks-1) *1.e-3_rp, ' | ', &
797  grid_cz(ks) *1.e-3_rp, ' - ', &
798  grid_cz(ke-kbuff)*1.e-3_rp, ' | ', &
799  grid_fz(ke-kbuff)*1.e-3_rp, ' -buffer- ', &
800  grid_fz(ke) *1.e-3_rp, ' -HALO- ', &
801  grid_fz(ka) *1.e-3_rp
802  if( io_l ) write(io_fid_log,'(1x,8(A,F9.3))') '*** X:', &
803  grid_fxg(0) *1.e-3_rp, ' -HALO- ', &
804  grid_fxg(ihalo) *1.e-3_rp, ' -buffer- ', &
805  grid_fxg(ihalo+ibuff) *1.e-3_rp, ' | ', &
806  grid_cxg(ihalo+ibuff+1) *1.e-3_rp, ' - ', &
807  grid_cxg(iag-ihalo-ibuff)*1.e-3_rp, ' | ', &
808  grid_fxg(iag-ihalo-ibuff)*1.e-3_rp, ' -buffer- ', &
809  grid_fxg(iag-ihalo) *1.e-3_rp, ' -HALO- ', &
810  grid_fxg(iag) *1.e-3_rp
811  if( io_l ) write(io_fid_log,'(1x,8(A,F9.3))') '*** Y:', &
812  grid_fyg(0) *1.e-3_rp, ' -HALO- ', &
813  grid_fyg(jhalo) *1.e-3_rp, ' -buffer- ', &
814  grid_fyg(jhalo+jbuff) *1.e-3_rp, ' | ', &
815  grid_cyg(jhalo+jbuff+1) *1.e-3_rp, ' - ', &
816  grid_cyg(jag-jhalo-jbuff)*1.e-3_rp, ' | ', &
817  grid_fyg(jag-jhalo-jbuff)*1.e-3_rp, ' -buffer- ', &
818  grid_fyg(jag-jhalo) *1.e-3_rp, ' -HALO- ', &
819  grid_fyg(jag) *1.e-3_rp
820  if( io_l ) write(io_fid_log,*)
821  if( io_l ) write(io_fid_log,*) '*** Center Position of Grid (global) :'
822  if( io_l ) write(io_fid_log,'(1x,A,F12.3)') '*** X: ', grid_domain_center_x
823  if( io_l ) write(io_fid_log,'(1x,A,F12.3)') '*** Y: ', grid_domain_center_y
824 
825  if( io_l ) write(io_fid_log,*)
826  if( io_l ) write(io_fid_log,'(1x,A)') &
827  '|============= Vertical Coordinate =============|'
828  if( io_l ) write(io_fid_log,'(1x,A)') &
829  '| k z zh dz buffer k |'
830  if( io_l ) write(io_fid_log,'(1x,A)') &
831  '| [m] [m] [m] factor |'
832 
833  do k = ka, ke+1, -1
834  if( io_l ) write(io_fid_log,'(1x,A,F9.2,A,F9.2,I5,A)') &
835  '| ',grid_fz(k),' ', grid_fbfz(k),k,' |'
836  if( io_l ) write(io_fid_log,'(1x,A,I5,F9.2,A,2F9.2,A)') &
837  '|',k,grid_cz(k),' ',grid_cdz(k), grid_cbfz(k),' |'
838  enddo
839 
840  k = ke
841  if( io_l ) write(io_fid_log,'(1x,A,F9.2,A,F9.2,I5,A)') &
842  '| ',grid_fz(k),' ', grid_fbfz(k),k,' | KE = TOA'
843  if( io_l ) write(io_fid_log,'(1x,A,I5,F9.2,A,2F9.2,A)') &
844  '|',k,grid_cz(k),' ',grid_cdz(k), grid_cbfz(k),' |'
845 
846  do k = ke-1, ks, -1
847  if( io_l ) write(io_fid_log,'(1x,A,F9.2,A,F9.2,I5,A)') &
848  '| ',grid_fz(k),' ', grid_fbfz(k),k,' |'
849  if( io_l ) write(io_fid_log,'(1x,A,I5,F9.2,A,2F9.2,A)') &
850  '|',k,grid_cz(k),' ',grid_cdz(k), grid_cbfz(k),' |'
851  enddo
852 
853  k = ks-1
854  if( io_l ) write(io_fid_log,'(1x,A,F9.2,A,F9.2,I5,A)') &
855  '| ',grid_fz(k),' ', grid_fbfz(k),k,' | KS-1 = surface'
856  if( io_l ) write(io_fid_log,'(1x,A,I5,F9.2,A,2F9.2,A)') &
857  '|',k,grid_cz(k),' ',grid_cdz(k), grid_cbfz(k),' |'
858 
859  do k = ks-2, 1, -1
860  if( io_l ) write(io_fid_log,'(1x,A,F9.2,A,F9.2,I5,A)') &
861  '| ',grid_fz(k),' ', grid_fbfz(k),k,' |'
862  if( io_l ) write(io_fid_log,'(1x,A,I5,F9.2,A,2F9.2,A)') &
863  '|',k,grid_cz(k),' ',grid_cdz(k), grid_cbfz(k),' |'
864  enddo
865 
866  k = 0
867  if( io_l ) write(io_fid_log,'(1x,A,F9.2,A,I5,A)') &
868  '| ',grid_fz(k),' ',k,' |'
869 
870  if( io_l ) write(io_fid_log,'(1x,A)') &
871  '|===============================================|'
872 
873  if ( debug ) then
874  if( io_l ) write(io_fid_log,*)
875  if( io_l ) write(io_fid_log,*) ' ', 0, grid_fx(0)
876  do i = 1, ia-1
877  if( io_l ) write(io_fid_log,*) i, grid_cx(i), grid_cbfx(i), grid_cdx(i)
878  if( io_l ) write(io_fid_log,*) ' ', i, grid_fx(i), grid_fbfx(i), grid_fdx(i)
879  enddo
880  i = ia
881  if( io_l ) write(io_fid_log,*) i, grid_cx(i), grid_cbfx(i), grid_cdx(i)
882  if( io_l ) write(io_fid_log,*) ' ', i, grid_fx(i), grid_fbfx(i)
883 
884  if( io_l ) write(io_fid_log,*)
885  if( io_l ) write(io_fid_log,*) ' ', 0, grid_fy(0)
886  do j = 1, ja-1
887  if( io_l ) write(io_fid_log,*) j, grid_cy(j), grid_cbfy(j), grid_cdy(j)
888  if( io_l ) write(io_fid_log,*) ' ', j, grid_fy(j), grid_fbfy(j), grid_fdy(j)
889  enddo
890  j = ja
891  if( io_l ) write(io_fid_log,*) j, grid_cy(j), grid_cbfy(j), grid_cdy(j)
892  if( io_l ) write(io_fid_log,*) ' ', j, grid_fy(j), grid_fbfy(j)
893  endif
894 
895  return
896  end subroutine grid_generate
897 
898 end module scale_grid
scale_grid_index::imax
integer, public imax
of computational cells: x, local
Definition: scale_grid_index.F90:44
scale_rm_process::prc_num_x
integer, public prc_num_x
x length of 2D processor topology
Definition: scale_rm_process.F90:43
scale_grid_index::is
integer, public is
start point of inner domain: x, local
Definition: scale_grid_index.F90:60
scale_grid::buffer_dx
real(rp), public buffer_dx
thickness of buffer region [m]: x
Definition: scale_grid_cartesian.F90:46
scale_grid::bufffact_z
real(rp), public bufffact_z
strech factor for dz of buffer region
Definition: scale_grid_cartesian.F90:51
gtool_file
module GTOOL_FILE
Definition: gtool_file.f90:17
scale_grid::grid_domain_center_x
real(rp), public grid_domain_center_x
center position of global domain [m]: x
Definition: scale_grid_cartesian.F90:53
scale_grid_index::je
integer, public je
end point of inner domain: y, local
Definition: scale_grid_index.F90:63
scale_grid::grid_rcdy
real(rp), dimension(:), allocatable, public grid_rcdy
reciprocal of center-dy
Definition: scale_grid_cartesian.F90:64
scale_grid::dy
real(rp), public dy
length in the main region [m]: y
Definition: scale_grid_cartesian.F90:43
scale_grid::grid_cyg
real(rp), dimension(:), allocatable, public grid_cyg
center coordinate [m]: y, global
Definition: scale_grid_cartesian.F90:86
scale_process::prc_mpistop
subroutine, public prc_mpistop
Abort MPI.
Definition: scale_process.F90:267
scale_grid::grid_cbfyg
real(rp), dimension(:), allocatable, public grid_cbfyg
center buffer factor (0-1): y, global
Definition: scale_grid_cartesian.F90:94
scale_grid::grid_cxg
real(rp), dimension(:), allocatable, public grid_cxg
center coordinate [m]: x, global
Definition: scale_grid_cartesian.F90:85
scale_grid::grid_fdy
real(rp), dimension(:), allocatable, public grid_fdy
y-length of grid(j+1) to grid(j) [m]
Definition: scale_grid_cartesian.F90:71
scale_grid::dx
real(rp), public dx
length in the main region [m]: x
Definition: scale_grid_cartesian.F90:42
scale_stdio::io_l
logical, public io_l
output log or not? (this process)
Definition: scale_stdio.F90:61
scale_grid::grid_cdyg
real(rp), dimension(:), allocatable, public grid_cdyg
center coordinate [m]: y, global
Definition: scale_grid_cartesian.F90:90
scale_grid::grid_rcdx
real(rp), dimension(:), allocatable, public grid_rcdx
reciprocal of center-dx
Definition: scale_grid_cartesian.F90:63
scale_grid::grid_cz
real(rp), dimension(:), allocatable, public grid_cz
center coordinate [m]: z, local=global
Definition: scale_grid_cartesian.F90:56
scale_grid::grid_fxg
real(rp), dimension(:), allocatable, public grid_fxg
face coordinate [m]: x, global
Definition: scale_grid_cartesian.F90:83
scale_grid::grid_fbfy
real(rp), dimension(:), allocatable, public grid_fbfy
face buffer factor (0-1): y
Definition: scale_grid_cartesian.F90:81
scale_stdio
module STDIO
Definition: scale_stdio.F90:12
scale_grid_index::ke
integer, public ke
end point of inner domain: z, local
Definition: scale_grid_index.F90:59
scale_grid::bufffact_x
real(rp), public bufffact_x
strech factor for dx of buffer region
Definition: scale_grid_cartesian.F90:49
scale_grid::grid_rfdy
real(rp), dimension(:), allocatable, public grid_rfdy
reciprocal of face-dy
Definition: scale_grid_cartesian.F90:74
scale_grid::grid_cdxg
real(rp), dimension(:), allocatable, public grid_cdxg
center coordinate [m]: x, global
Definition: scale_grid_cartesian.F90:89
scale_grid::grid_rcdz
real(rp), dimension(:), allocatable, public grid_rcdz
reciprocal of center-dz
Definition: scale_grid_cartesian.F90:62
scale_grid::dz
real(rp), public dz
length in the main region [m]: z
Definition: scale_grid_cartesian.F90:41
scale_grid::grid_fx
real(rp), dimension(:), allocatable, public grid_fx
face coordinate [m]: x, local
Definition: scale_grid_cartesian.F90:67
scale_rm_process::prc_num_y
integer, public prc_num_y
y length of 2D processor topology
Definition: scale_rm_process.F90:44
scale_grid_index
module grid index
Definition: scale_grid_index.F90:14
scale_stdio::io_nml
logical, public io_nml
output log or not? (for namelist, this process)
Definition: scale_stdio.F90:62
scale_grid_index::ia
integer, public ia
of whole cells: x, local, with HALO
Definition: scale_grid_index.F90:55
scale_grid_index::jag
integer, public jag
of computational grids
Definition: scale_grid_index.F90:82
scale_grid::grid_fbfx
real(rp), dimension(:), allocatable, public grid_fbfx
face buffer factor (0-1): x
Definition: scale_grid_cartesian.F90:80
scale_grid::grid_fdz
real(rp), dimension(:), allocatable, public grid_fdz
z-length of grid(k+1) to grid(k) [m]
Definition: scale_grid_cartesian.F90:69
scale_grid_index::ka
integer, public ka
of whole cells: z, local, with HALO
Definition: scale_grid_index.F90:54
scale_grid::grid_fbfxg
real(rp), dimension(:), allocatable, public grid_fbfxg
face buffer factor (0-1): x, global
Definition: scale_grid_cartesian.F90:91
scale_grid::grid_fdyg
real(rp), dimension(:), allocatable, public grid_fdyg
center coordinate [m]: y, global
Definition: scale_grid_cartesian.F90:88
scale_grid_index::kmax
integer, public kmax
of computational cells: z, local
Definition: scale_grid_index.F90:43
scale_grid::grid_fz
real(rp), dimension(:), allocatable, public grid_fz
face coordinate [m]: z, local=global
Definition: scale_grid_cartesian.F90:66
scale_grid::grid_fbfz
real(rp), dimension(:), allocatable, public grid_fbfz
face buffer factor (0-1): z
Definition: scale_grid_cartesian.F90:79
scale_grid_index::jhalo
integer, public jhalo
of halo cells: y
Definition: scale_grid_index.F90:52
scale_grid_index::js
integer, public js
start point of inner domain: y, local
Definition: scale_grid_index.F90:62
scale_grid_index::iag
integer, public iag
of computational grids
Definition: scale_grid_index.F90:81
scale_grid::grid_cbfx
real(rp), dimension(:), allocatable, public grid_cbfx
center buffer factor (0-1): x
Definition: scale_grid_cartesian.F90:77
scale_process
module PROCESS
Definition: scale_process.F90:14
scale_grid::bufffact_y
real(rp), public bufffact_y
strech factor for dy of buffer region
Definition: scale_grid_cartesian.F90:50
scale_grid::grid_fbfyg
real(rp), dimension(:), allocatable, public grid_fbfyg
face buffer factor (0-1): y, global
Definition: scale_grid_cartesian.F90:92
scale_grid::grid_generate
subroutine, public grid_generate
Generate horizontal&vertical grid.
Definition: scale_grid_cartesian.F90:312
scale_grid_index::ks
integer, public ks
start point of inner domain: z, local
Definition: scale_grid_index.F90:58
scale_grid::grid_cbfz
real(rp), dimension(:), allocatable, public grid_cbfz
center buffer factor (0-1): z
Definition: scale_grid_cartesian.F90:76
scale_process::prc_myrank
integer, public prc_myrank
process num in local communicator
Definition: scale_process.F90:79
scale_grid::grid_cx
real(rp), dimension(:), allocatable, public grid_cx
center coordinate [m]: x, local
Definition: scale_grid_cartesian.F90:57
scale_grid
module GRID (cartesian)
Definition: scale_grid_cartesian.F90:17
scale_rm_process
module RM PROCESS
Definition: scale_rm_process.F90:10
scale_grid::buffer_dy
real(rp), public buffer_dy
thickness of buffer region [m]: y
Definition: scale_grid_cartesian.F90:47
scale_grid::grid_domain_center_y
real(rp), public grid_domain_center_y
center position of global domain [m]: y
Definition: scale_grid_cartesian.F90:54
scale_prof
module profiler
Definition: scale_prof.F90:10
scale_grid_index::ie
integer, public ie
end point of inner domain: x, local
Definition: scale_grid_index.F90:61
scale_grid_index::imaxg
integer, public imaxg
of computational cells: x, global
Definition: scale_grid_index.F90:79
scale_grid::buffer_dz
real(rp), public buffer_dz
thickness of buffer region [m]: z
Definition: scale_grid_cartesian.F90:45
scale_grid::grid_fdxg
real(rp), dimension(:), allocatable, public grid_fdxg
center coordinate [m]: x, global
Definition: scale_grid_cartesian.F90:87
scale_precision
module PRECISION
Definition: scale_precision.F90:13
scale_rm_process::prc_2drank
integer, dimension(:,:), allocatable, public prc_2drank
node index in 2D topology
Definition: scale_rm_process.F90:46
scale_grid::grid_cdz
real(rp), dimension(:), allocatable, public grid_cdz
z-length of control volume [m]
Definition: scale_grid_cartesian.F90:59
scale_grid::grid_fdx
real(rp), dimension(:), allocatable, public grid_fdx
x-length of grid(i+1) to grid(i) [m]
Definition: scale_grid_cartesian.F90:70
scale_grid::grid_cdy
real(rp), dimension(:), allocatable, public grid_cdy
y-length of control volume [m]
Definition: scale_grid_cartesian.F90:61
scale_grid::grid_cbfy
real(rp), dimension(:), allocatable, public grid_cbfy
center buffer factor (0-1): y
Definition: scale_grid_cartesian.F90:78
scale_grid::grid_rfdx
real(rp), dimension(:), allocatable, public grid_rfdx
reciprocal of face-dx
Definition: scale_grid_cartesian.F90:73
scale_grid::grid_rfdz
real(rp), dimension(:), allocatable, public grid_rfdz
reciprocal of face-dz
Definition: scale_grid_cartesian.F90:72
scale_stdio::io_fid_conf
integer, public io_fid_conf
Config file ID.
Definition: scale_stdio.F90:55
scale_stdio::io_fid_log
integer, public io_fid_log
Log file ID.
Definition: scale_stdio.F90:56
scale_grid::grid_allocate
subroutine, public grid_allocate
Allocate arrays.
Definition: scale_grid_cartesian.F90:202
scale_grid_index::jmaxg
integer, public jmaxg
of computational cells: y, global
Definition: scale_grid_index.F90:80
scale_grid::bufffact
real(rp), public bufffact
default strech factor for dx/dy/dz of buffer region
Definition: scale_grid_cartesian.F90:48
scale_grid::grid_cdx
real(rp), dimension(:), allocatable, public grid_cdx
x-length of control volume [m]
Definition: scale_grid_cartesian.F90:60
scale_grid_index::jmax
integer, public jmax
of computational cells: y, local
Definition: scale_grid_index.F90:45
scale_stdio::io_fid_nml
integer, public io_fid_nml
Log file ID (only for output namelist)
Definition: scale_stdio.F90:57
scale_grid::grid_fyg
real(rp), dimension(:), allocatable, public grid_fyg
face coordinate [m]: y, global
Definition: scale_grid_cartesian.F90:84
scale_grid::grid_cy
real(rp), dimension(:), allocatable, public grid_cy
center coordinate [m]: y, local
Definition: scale_grid_cartesian.F90:58
scale_grid::grid_setup
subroutine, public grid_setup
Setup.
Definition: scale_grid_cartesian.F90:125
scale_grid::grid_cbfxg
real(rp), dimension(:), allocatable, public grid_cbfxg
center buffer factor (0-1): x, global
Definition: scale_grid_cartesian.F90:93
scale_grid_index::ihalo
integer, public ihalo
of halo cells: x
Definition: scale_grid_index.F90:51
scale_grid_index::ja
integer, public ja
of whole cells: y, local, with HALO
Definition: scale_grid_index.F90:56
scale_grid::grid_fy
real(rp), dimension(:), allocatable, public grid_fy
face coordinate [m]: y, local
Definition: scale_grid_cartesian.F90:68
Generated by   doxygen 1.8.14