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_fbfxg(:)
88  real(RP), public, allocatable :: grid_fbfyg(:)
89  real(RP), public, allocatable :: grid_cbfxg(:)
90  real(RP), public, allocatable :: grid_cbfyg(:)
91 
92  !-----------------------------------------------------------------------------
93  !
94  !++ Private procedure
95  !
96  private :: grid_read
97 
98  !-----------------------------------------------------------------------------
99  !
100  !++ Private parameters & variables
101  !
102  character(len=H_LONG), private :: grid_in_basename = ''
103  character(len=H_LONG), private :: grid_out_basename = ''
104  real(RP), private :: grid_offset_x = 0.0_rp
105  real(RP), private :: grid_offset_y = 0.0_rp
106 
107  integer, private, parameter :: kmax_user_lim = 300
108  real(RP), private :: fz(kmax_user_lim)
109 
110  logical, private :: debug = .false.
111 
112  !-----------------------------------------------------------------------------
113 contains
114  !-----------------------------------------------------------------------------
116  subroutine grid_setup
117  use scale_process, only: &
118  prc_mpistop
119  implicit none
120 
121  namelist / param_grid / &
122  grid_in_basename, &
123  grid_out_basename, &
124  grid_offset_x, &
125  grid_offset_y, &
126  dx, &
127  dy, &
128  dz, &
129  buffer_dz, &
130  buffer_dx, &
131  buffer_dy, &
132  bufffact, &
133  bufffact_x, &
134  bufffact_y, &
135  bufffact_z, &
136  fz, &
137  debug
138 
139  integer :: ierr
140  !---------------------------------------------------------------------------
141 
142  if( io_l ) write(io_fid_log,*)
143  if( io_l ) write(io_fid_log,*) '++++++ Module[GRID] / Categ[ATMOS-RM GRID] / Origin[SCALElib]'
144 
145  call grid_allocate
146 
147  !--- read namelist
148  rewind(io_fid_conf)
149  read(io_fid_conf,nml=param_grid,iostat=ierr)
150  if( ierr < 0 ) then !--- missing
151  if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
152  elseif( ierr > 0 ) then !--- fatal error
153  write(*,*) 'xxx Not appropriate names in namelist PARAM_GRID. Check!'
154  call prc_mpistop
155  endif
156  if( io_lnml ) write(io_fid_log,nml=param_grid)
157 
158  if ( bufffact_x < 0.0_rp ) bufffact_x = bufffact
159  if ( bufffact_y < 0.0_rp ) bufffact_y = bufffact
160  if ( bufffact_z < 0.0_rp ) bufffact_z = bufffact
161 
162  if( io_l ) write(io_fid_log,*)
163  if( io_l ) write(io_fid_log,*) '*** Atmosphere grid information ***'
164  if( io_l ) write(io_fid_log,'(1x,A,3(F7.1))') '*** delta Z, X, Y [m] :', dz, dx, dy
165 
166  if ( grid_in_basename /= '' ) then
167  call grid_read
168  else
169  if( io_l ) write(io_fid_log,*) '*** Not found input grid file. Grid position is calculated.'
170 
171  call grid_generate
172  endif
173 
174  if( io_l ) write(io_fid_log,*)
175  if( io_l ) write(io_fid_log,*) '*** Domain size [km] (local) :'
176  if( io_l ) write(io_fid_log,'(1x,6(A,F8.3))') ' X:', &
177  grid_fx(0) *1.e-3_rp, ' -HALO- ', grid_fx(is-1)*1.e-3_rp, ' | ', &
178  grid_cx(is)*1.e-3_rp, ' - ', grid_cx(ie) *1.e-3_rp, ' | ', &
179  grid_fx(ie)*1.e-3_rp, ' -HALO- ', grid_fx(ia) *1.e-3_rp
180  if( io_l ) write(io_fid_log,'(1x,6(A,F8.3))') ' Y:', &
181  grid_fy(0) *1.e-3_rp, ' -HALO- ', grid_fy(js-1)*1.e-3_rp, ' | ', &
182  grid_cy(js)*1.e-3_rp, ' - ', grid_cy(je) *1.e-3_rp, ' | ', &
183  grid_fy(je)*1.e-3_rp, ' -HALO- ', grid_fy(ja) *1.e-3_rp
184 
185  return
186  end subroutine grid_setup
187 
188  !-----------------------------------------------------------------------------
190  subroutine grid_allocate
191  use scale_rm_process, only: &
192  prc_num_x, &
193  prc_num_y
194  implicit none
195 
196  integer :: IAG, JAG
197  !---------------------------------------------------------------------------
198 
199  ! working
200  fz(:) = -1.0_rp
201 
202  ! local domain
203  allocate( grid_cz(ka) )
204  allocate( grid_cx(ia) )
205  allocate( grid_cy(ja) )
206  allocate( grid_cdz(ka) )
207  allocate( grid_cdx(ia) )
208  allocate( grid_cdy(ja) )
209  allocate( grid_rcdz(ka) )
210  allocate( grid_rcdx(ia) )
211  allocate( grid_rcdy(ja) )
212 
213  allocate( grid_fz(0:ka) )
214  allocate( grid_fx(0:ia) )
215  allocate( grid_fy(0:ja) )
216  allocate( grid_fdz(ka-1) )
217  allocate( grid_fdx(ia-1) )
218  allocate( grid_fdy(ja-1) )
219  allocate( grid_rfdz(ka-1) )
220  allocate( grid_rfdx(ia-1) )
221  allocate( grid_rfdy(ja-1) )
222 
223  allocate( grid_cbfz(ka) )
224  allocate( grid_cbfx(ia) )
225  allocate( grid_cbfy(ja) )
226  allocate( grid_fbfz(ka) )
227  allocate( grid_fbfx(ia) )
228  allocate( grid_fbfy(ja) )
229 
230  ! array size (global domain)
231  iag = ihalo + imax*prc_num_x + ihalo
232  jag = jhalo + jmax*prc_num_y + jhalo
233 
234  ! global domain
235  allocate( grid_fxg(0:iag) )
236  allocate( grid_fyg(0:jag) )
237  allocate( grid_cxg( iag) )
238  allocate( grid_cyg( jag) )
239  allocate( grid_cbfxg( iag) )
240  allocate( grid_cbfyg( jag) )
241  allocate( grid_fbfxg( iag) )
242  allocate( grid_fbfyg( jag) )
243 
244  return
245  end subroutine grid_allocate
246 
247  !-----------------------------------------------------------------------------
249  subroutine grid_read
250  use gtool_file, only: &
251  fileread
252  use scale_process, only: &
253  prc_myrank
254  implicit none
255 
256  character(len=H_LONG) :: bname
257  !---------------------------------------------------------------------------
258 
259  if( io_l ) write(io_fid_log,*)
260  if( io_l ) write(io_fid_log,*) '*** Input grid file ***'
261 
262  write(bname,'(A,A,F15.3)') trim(grid_in_basename)
263 
264  call fileread( grid_cz(:), bname, 'CZ', 1, prc_myrank )
265  call fileread( grid_cx(:), bname, 'CX', 1, prc_myrank )
266  call fileread( grid_cy(:), bname, 'CY', 1, prc_myrank )
267  call fileread( grid_fz(:), bname, 'FZ', 1, prc_myrank )
268  call fileread( grid_fx(:), bname, 'FX', 1, prc_myrank )
269  call fileread( grid_fy(:), bname, 'FY', 1, prc_myrank )
270 
271  call fileread( grid_cxg(:), bname, 'CXG', 1, prc_myrank )
272  call fileread( grid_cyg(:), bname, 'CYG', 1, prc_myrank )
273  call fileread( grid_fxg(:), bname, 'FXG', 1, prc_myrank )
274  call fileread( grid_fyg(:), bname, 'FYG', 1, prc_myrank )
275 
276  call fileread( grid_cdz(:), bname, 'CDZ', 1, prc_myrank )
277  call fileread( grid_cdx(:), bname, 'CDX', 1, prc_myrank )
278  call fileread( grid_cdy(:), bname, 'CDY', 1, prc_myrank )
279  call fileread( grid_fdz(:), bname, 'FDZ', 1, prc_myrank )
280  call fileread( grid_fdx(:), bname, 'FDX', 1, prc_myrank )
281  call fileread( grid_fdy(:), bname, 'FDY', 1, prc_myrank )
282 
283  grid_rcdz(:) = 1.0_rp / grid_cdz(:)
284  grid_rcdx(:) = 1.0_rp / grid_cdx(:)
285  grid_rcdy(:) = 1.0_rp / grid_cdy(:)
286  grid_rfdz(:) = 1.0_rp / grid_fdz(:)
287  grid_rfdx(:) = 1.0_rp / grid_fdx(:)
288  grid_rfdy(:) = 1.0_rp / grid_fdy(:)
289 
290  call fileread( grid_cbfz(:), bname, 'CBFZ', 1, prc_myrank )
291  call fileread( grid_cbfx(:), bname, 'CBFX', 1, prc_myrank )
292  call fileread( grid_cbfy(:), bname, 'CBFY', 1, prc_myrank )
293  call fileread( grid_fbfz(:), bname, 'FBFZ', 1, prc_myrank )
294  call fileread( grid_fbfx(:), bname, 'FBFX', 1, prc_myrank )
295  call fileread( grid_fbfy(:), bname, 'FBFY', 1, prc_myrank )
296 
297  return
298  end subroutine grid_read
299 
300  !-----------------------------------------------------------------------------
302  subroutine grid_generate
303  use scale_process, only: &
304  prc_mpistop, &
305  prc_myrank
306  use scale_rm_process, only: &
307  prc_2drank, &
308  prc_num_x, &
309  prc_num_y
310  implicit none
311 
312  integer :: IAG ! # of x whole cells (global, with HALO)
313  integer :: JAG ! # of y whole cells (global, with HALO)
314 
315  real(RP), allocatable :: buffz(:), buffx(:), buffy(:)
316  real(RP) :: bufftotz, bufftotx, bufftoty
317 
318  integer :: kbuff, ibuff, jbuff
319  integer :: kmain, imain, jmain
320 
321  integer :: k, i, j, ii, jj
322  !---------------------------------------------------------------------------
323 
324  !##### coordinate in global domain #####
325 
326  ! array size (global domain)
327  iag = ihalo + imax*prc_num_x + ihalo
328  jag = jhalo + jmax*prc_num_y + jhalo
329 
330  allocate( buffx(0:iag) )
331  allocate( buffy(0:jag) )
332 
333  ! X-direction
334  ! calculate buffer grid size
335  buffx(0) = dx
336  bufftotx = 0.0_rp
337 
338  do i = 1, iag
339  if( bufftotx >= buffer_dx ) exit
340 
341  buffx(i) = buffx(i-1) * bufffact_x
342  bufftotx = bufftotx + buffx(i)
343  enddo
344  ibuff = i - 1
345  imain = iag - 2*ibuff - 2*ihalo
346 
347  if ( imain < 0 ) then
348  write(*,*) 'xxx Buffer size (', bufftotx*2.0_rp, ') must be smaller than global domain size (X). Use smaller BUFFER_DX!'
349  call prc_mpistop
350  endif
351 
352  ! horizontal coordinate (global domaim)
353  grid_fxg(ihalo) = grid_offset_x
354  do i = ihalo-1, 0, -1
355  grid_fxg(i) = grid_fxg(i+1) - buffx(ibuff)
356  enddo
357 
358  do i = 1, ihalo
359  grid_cxg(i) = 0.5_rp * ( grid_fxg(i)+grid_fxg(i-1) )
360  enddo
361 
362  if ( ibuff > 0 ) then
363  do i = ihalo+1, ihalo+ibuff
364  grid_fxg(i) = grid_fxg(i-1) + buffx(ibuff+ihalo+1-i)
365  grid_cxg(i) = 0.5_rp * ( grid_fxg(i)+grid_fxg(i-1) )
366  enddo
367  endif
368 
369  do i = ihalo+ibuff+1, ihalo+ibuff+imain
370  grid_fxg(i) = grid_fxg(i-1) + dx
371  grid_cxg(i) = 0.5_rp * ( grid_fxg(i)+grid_fxg(i-1) )
372  enddo
373 
374  if ( ibuff > 0 ) then
375  do i = ihalo+ibuff+imain+1, ihalo+ibuff+imain+ibuff
376  grid_fxg(i) = grid_fxg(i-1) + buffx(i-ihalo-ibuff-imain)
377  grid_cxg(i) = 0.5_rp * ( grid_fxg(i)+grid_fxg(i-1) )
378  enddo
379  endif
380 
381  do i = ihalo+ibuff+imain+ibuff+1, ihalo+ibuff+imain+ibuff+ihalo
382  grid_fxg(i) = grid_fxg(i-1) + buffx(ibuff)
383  grid_cxg(i) = 0.5_rp * ( grid_fxg(i)+grid_fxg(i-1) )
384  enddo
385 
386  ! calc buffer factor (global domaim)
387  grid_cbfxg(:) = 0.0_rp
388  grid_fbfxg(:) = 0.0_rp
389  do i = 1, ihalo
390  grid_cbfxg(i) = 1.0_rp
391  grid_fbfxg(i) = 1.0_rp
392  enddo
393 
394  if ( ibuff > 0 ) then
395  do i = ihalo+1, ihalo+ibuff
396  grid_cbfxg(i) = (bufftotx+grid_fxg(ihalo)-grid_cxg(i)) / bufftotx
397  grid_fbfxg(i) = (bufftotx+grid_fxg(ihalo)-grid_fxg(i)) / bufftotx
398  enddo
399 
400  do i = ihalo+ibuff+imain+1, ihalo+ibuff+imain+ibuff
401  grid_cbfxg(i) = (bufftotx-grid_fxg(iag-ihalo)+grid_cxg(i)) / bufftotx
402  grid_fbfxg(i) = (bufftotx-grid_fxg(iag-ihalo)+grid_fxg(i)) / bufftotx
403  enddo
404  endif
405 
406  do i = ihalo+ibuff+imain+ibuff+1, ihalo+ibuff+imain+ibuff+ihalo
407  grid_cbfxg(i) = 1.0_rp
408  grid_fbfxg(i) = 1.0_rp
409  enddo
410 
411  grid_cbfxg(:) = max( min( grid_cbfxg(:), 1.0_rp ), 0.0_rp )
412  grid_fbfxg(:) = max( min( grid_fbfxg(:), 1.0_rp ), 0.0_rp )
413 
414  ! Y-direction
415  ! calculate buffer grid size
416  buffy(0) = dy
417  bufftoty = 0.0_rp
418 
419  do j = 1, jag
420  if( bufftoty >= buffer_dy ) exit
421 
422  buffy(j) = buffy(j-1) * bufffact_y
423  bufftoty = bufftoty + buffy(j)
424  enddo
425  jbuff = j - 1
426  jmain = jag - 2*jbuff - 2*jhalo
427 
428  if ( jmain < 0 ) then
429  write(*,*) 'xxx Buffer size (', bufftoty*2.0_rp, ') must be smaller than global domain size (Y). Use smaller BUFFER_DY!'
430  call prc_mpistop
431  endif
432 
433  ! horizontal coordinate (global domaim)
434  grid_fyg(jhalo) = grid_offset_y
435  do j = jhalo-1, 0, -1
436  grid_fyg(j) = grid_fyg(j+1) - buffy(jbuff)
437  enddo
438 
439  do j = 1, jhalo
440  grid_cyg(j) = 0.5_rp * ( grid_fyg(j)+grid_fyg(j-1) )
441  enddo
442 
443  if ( jbuff > 0 ) then
444  do j = jhalo+1, jhalo+jbuff
445  grid_fyg(j) = grid_fyg(j-1) + buffy(jbuff+jhalo+1-j)
446  grid_cyg(j) = 0.5_rp * ( grid_fyg(j)+grid_fyg(j-1) )
447  enddo
448  endif
449 
450  do j = jhalo+jbuff+1, jhalo+jbuff+jmain
451  grid_fyg(j) = grid_fyg(j-1) + dy
452  grid_cyg(j) = 0.5_rp * ( grid_fyg(j)+grid_fyg(j-1) )
453  enddo
454 
455  if ( jbuff > 0 ) then
456  do j = jhalo+jbuff+jmain+1, jhalo+jbuff+jmain+jbuff
457  grid_fyg(j) = grid_fyg(j-1) + buffy(j-jhalo-jbuff-jmain)
458  grid_cyg(j) = 0.5_rp * ( grid_fyg(j)+grid_fyg(j-1) )
459  enddo
460  endif
461 
462  do j = jhalo+jbuff+jmain+jbuff+1, jhalo+jbuff+jmain+jbuff+jhalo
463  grid_fyg(j) = grid_fyg(j-1) + buffy(jbuff)
464  grid_cyg(j) = 0.5_rp * ( grid_fyg(j)+grid_fyg(j-1) )
465  enddo
466 
467  ! calc buffer factor (global domaim)
468  grid_cbfyg(:) = 0.0_rp
469  grid_fbfyg(:) = 0.0_rp
470  do j = 1, jhalo
471  grid_cbfyg(j) = 1.0_rp
472  grid_fbfyg(j) = 1.0_rp
473  enddo
474 
475  if ( jbuff > 0 ) then
476  do j = jhalo+1, jhalo+jbuff
477  grid_cbfyg(j) = (bufftoty+grid_fyg(jhalo)-grid_cyg(j)) / bufftoty
478  grid_fbfyg(j) = (bufftoty+grid_fyg(jhalo)-grid_fyg(j)) / bufftoty
479  enddo
480 
481  do j = jhalo+jbuff+jmain+1, jhalo+jbuff+jmain+jbuff
482  grid_cbfyg(j) = (bufftoty-grid_fyg(jag-jhalo)+grid_cyg(j)) / bufftoty
483  grid_fbfyg(j) = (bufftoty-grid_fyg(jag-jhalo)+grid_fyg(j)) / bufftoty
484  enddo
485  endif
486 
487  do j = jhalo+jbuff+jmain+jbuff+1, jhalo+jbuff+jmain+jbuff+jhalo
488  grid_cbfyg(j) = 1.0_rp
489  grid_fbfyg(j) = 1.0_rp
490  enddo
491  grid_cbfyg(:) = max( min( grid_cbfyg(:), 1.0_rp ), 0.0_rp )
492  grid_fbfyg(:) = max( min( grid_fbfyg(:), 1.0_rp ), 0.0_rp )
493 
494  deallocate( buffx )
495  deallocate( buffy )
496 
497  !##### coordinate in local domain #####
498 
499  allocate( buffz(0:ka) )
500 
501  if ( minval(fz(1:kmax)) > 0.0_rp ) then ! input from namelist
502  if( io_l ) write(io_fid_log,*) '*** Z coordinate is given from NAMELIST.'
503 
504  if ( kmax < 2 ) then
505  write(*,*) 'xxx If you use FZ, KMAX must be larger than 1. Check!', kmax
506  call prc_mpistop
507  endif
508 
509  ! Z-direction
510  ! calculate buffer grid size
511  bufftotz = 0.0_rp
512 
513  do k = kmax, 2, -1
514  if( bufftotz >= buffer_dz ) exit
515 
516  bufftotz = bufftotz + ( fz(k) - fz(k-1) )
517  enddo
518  kbuff = kmax - k
519  kmain = k
520 
521  if ( kmain < 0 ) then
522  write(*,*) 'xxx Buffer size (', bufftotz, ') must be smaller than domain size (z). Use smaller BUFFER_DZ!'
523  call prc_mpistop
524  endif
525 
526  ! vartical coordinate (local=global domaim)
527  grid_fz(ks-1) = 0.0_rp
528 
529  dz = fz(1)
530  do k = ks-2, 0, -1
531  grid_fz(k) = grid_fz(k+1) - dz
532  enddo
533 
534  do k = ks, ke
535  grid_fz(k) = fz(k-ks+1)
536  enddo
537 
538  dz = fz(kmax) - fz(kmax-1)
539  do k = ke+1, ka
540  grid_fz(k) = grid_fz(k-1) + dz
541  enddo
542 
543  do k = 1, ka
544  grid_cz(k) = 0.5_rp * ( grid_fz(k)+grid_fz(k-1) )
545  enddo
546 
547  else ! calc using DZ
548 
549  ! Z-direction
550  ! calculate buffer grid size
551  buffz(0) = dz
552  bufftotz = 0.0_rp
553 
554  do k = 1, ka
555  if( bufftotz >= buffer_dz ) exit
556 
557  buffz(k) = buffz(k-1) * bufffact_z
558  bufftotz = bufftotz + buffz(k)
559  enddo
560  kbuff = k - 1
561  kmain = ke - ks + 1 - kbuff
562 
563  if ( kmain < 0 ) then
564  write(*,*) 'xxx Buffer size (', bufftotz, ') must be smaller than domain size (z). Use smaller BUFFER_DZ!'
565  call prc_mpistop
566  endif
567 
568  ! vartical coordinate (local=global domaim)
569  grid_fz(ks-1) = 0.0_rp
570  do k = ks-2, 0, -1
571  grid_fz(k) = grid_fz(k+1) - dz
572  enddo
573 
574  do k = 1, ks-1
575  grid_cz(k) = 0.5_rp * ( grid_fz(k)+grid_fz(k-1) )
576  enddo
577 
578  do k = ks, ks+kmain-1
579  grid_fz(k) = grid_fz(k-1) + dz
580  grid_cz(k) = 0.5_rp * ( grid_fz(k)+grid_fz(k-1) )
581  enddo
582 
583  if ( kbuff > 0 ) then
584  do k = ks+kmain, ke
585  grid_fz(k) = grid_fz(k-1) + buffz(k-ks-kmain+1)
586  grid_cz(k) = 0.5_rp * ( grid_fz(k)+grid_fz(k-1) )
587  enddo
588  endif
589 
590  do k = ke+1, ka
591  grid_fz(k) = grid_fz(k-1) + buffz(kbuff)
592  grid_cz(k) = 0.5_rp * ( grid_fz(k)+grid_fz(k-1) )
593  enddo
594 
595  endif
596 
597  ! calc buffer factor (global domaim)
598  grid_cbfz(:) = 0.0_rp
599  grid_fbfz(:) = 0.0_rp
600  if ( kbuff > 0 ) then
601  do k = ks+kmain, ke
602  grid_cbfz(k) = (bufftotz-grid_fz(ke)+grid_cz(k)) / bufftotz
603  grid_fbfz(k) = (bufftotz-grid_fz(ke)+grid_fz(k)) / bufftotz
604  enddo
605  endif
606 
607  do k = ke+1, ka
608  grid_cbfz(k) = 1.0_rp
609  grid_fbfz(k) = 1.0_rp
610  enddo
611  grid_cbfz(:) = max( min( grid_cbfz(:), 1.0_rp ), 0.0_rp )
612  grid_fbfz(:) = max( min( grid_fbfz(:), 1.0_rp ), 0.0_rp )
613 
614  deallocate( buffz )
615 
616  ! vartical coordinate (local domaim)
617  do k = 1, ka
618  grid_cdz(k) = grid_fz(k) - grid_fz(k-1)
619  grid_rcdz(k) = 1.0_rp / grid_cdz(k)
620  enddo
621 
622  do k = 1, ka-1
623  grid_fdz(k) = grid_cz(k+1)-grid_cz(k)
624  grid_rfdz(k) = 1.0_rp / grid_fdz(k)
625  enddo
626 
627  ! X-direction
628  ! horizontal coordinate (local domaim)
629  do i = 0, ia
630  ii = i + prc_2drank(prc_myrank,1) * imax
631 
632  grid_fx(i) = grid_fxg(ii)
633  enddo
634 
635  do i = 1, ia
636  ii = i + prc_2drank(prc_myrank,1) * imax
637 
638  grid_cx(i) = grid_cxg(ii)
639  grid_cbfx(i) = grid_cbfxg(ii)
640  grid_fbfx(i) = grid_fbfxg(ii)
641 
642  grid_cdx(i) = grid_fx(i) - grid_fx(i-1)
643  grid_rcdx(i) = 1.0_rp / grid_cdx(i)
644  enddo
645 
646  do i = 1, ia-1
647  grid_fdx(i) = grid_cx(i+1)-grid_cx(i)
648  grid_rfdx(i) = 1.0_rp / grid_fdx(i)
649  enddo
650 
651  ! Y-direction
652  ! horizontal coordinate (local domaim)
653  do j = 0, ja
654  jj = j + prc_2drank(prc_myrank,2) * jmax
655 
656  grid_fy(j) = grid_fyg(jj)
657  enddo
658 
659  do j = 1, ja
660  jj = j + prc_2drank(prc_myrank,2) * jmax
661 
662  grid_cy(j) = grid_cyg(jj)
663  grid_cbfy(j) = grid_cbfyg(jj)
664  grid_fbfy(j) = grid_fbfyg(jj)
665 
666  grid_cdy(j) = grid_fy(j) - grid_fy(j-1)
667  grid_rcdy(j) = 1.0_rp / grid_cdy(j)
668  enddo
669 
670  do j = 1, ja-1
671  grid_fdy(j) = grid_cy(j+1)-grid_cy(j)
672  grid_rfdy(j) = 1.0_rp / grid_fdy(j)
673  enddo
674 
675  grid_domain_center_x = 0.5_rp * ( grid_fxg(0) + grid_fxg(iag) )
676  grid_domain_center_y = 0.5_rp * ( grid_fyg(0) + grid_fyg(jag) )
677 
678  ! report
679  if( io_l ) write(io_fid_log,*)
680  if( io_l ) write(io_fid_log,*) '*** Main/buffer Grid (global) :'
681  if( io_l ) write(io_fid_log,'(1x,2(A,I6))') ' Z: buffer = ', kbuff,' x 1, main = ',kmain
682  if( io_l ) write(io_fid_log,'(1x,2(A,I6))') ' X: buffer = ', ibuff,' x 2, main = ',imain
683  if( io_l ) write(io_fid_log,'(1x,2(A,I6))') ' Y: buffer = ', jbuff,' x 2, main = ',jmain
684  if( io_l ) write(io_fid_log,*)
685  if( io_l ) write(io_fid_log,*) '*** Domain size [km] (global) :'
686  if( io_l ) write(io_fid_log,'(1x,7(A,F8.3))') ' Z:', &
687  grid_fz(0) *1.e-3_rp, ' -HALO- ', &
688  grid_fz(ks-1) *1.e-3_rp, ' | ', &
689  grid_cz(ks) *1.e-3_rp, ' - ', &
690  grid_cz(ke-kbuff)*1.e-3_rp, ' | ', &
691  grid_fz(ke-kbuff)*1.e-3_rp, ' -buffer- ', &
692  grid_fz(ke) *1.e-3_rp, ' -HALO- ', &
693  grid_fz(ka) *1.e-3_rp
694  if( io_l ) write(io_fid_log,'(1x,8(A,F8.3))') ' X:', &
695  grid_fxg(0) *1.e-3_rp, ' -HALO- ', &
696  grid_fxg(ihalo) *1.e-3_rp, ' -buffer- ', &
697  grid_fxg(ihalo+ibuff) *1.e-3_rp, ' | ', &
698  grid_cxg(ihalo+ibuff+1) *1.e-3_rp, ' - ', &
699  grid_cxg(iag-ihalo-ibuff)*1.e-3_rp, ' | ', &
700  grid_fxg(iag-ihalo-ibuff)*1.e-3_rp, ' -buffer- ', &
701  grid_fxg(iag-ihalo) *1.e-3_rp, ' -HALO- ', &
702  grid_fxg(iag) *1.e-3_rp
703  if( io_l ) write(io_fid_log,'(1x,8(A,F8.3))') ' Y:', &
704  grid_fyg(0) *1.e-3_rp, ' -HALO- ', &
705  grid_fyg(jhalo) *1.e-3_rp, ' -buffer- ', &
706  grid_fyg(jhalo+jbuff) *1.e-3_rp, ' | ', &
707  grid_cyg(jhalo+jbuff+1) *1.e-3_rp, ' - ', &
708  grid_cyg(jag-jhalo-jbuff)*1.e-3_rp, ' | ', &
709  grid_fyg(jag-jhalo-jbuff)*1.e-3_rp, ' -buffer- ', &
710  grid_fyg(jag-jhalo) *1.e-3_rp, ' -HALO- ', &
711  grid_fyg(jag) *1.e-3_rp
712  if( io_l ) write(io_fid_log,*)
713  if( io_l ) write(io_fid_log,*) '*** Center Position of Grid (global) :'
714  if( io_l ) write(io_fid_log,'(1x,A,F12.3)') ' X: ', grid_domain_center_x
715  if( io_l ) write(io_fid_log,'(1x,A,F12.3)') ' Y: ', grid_domain_center_y
716 
717  if ( debug ) then
718  if( io_l ) write(io_fid_log,*)
719  if( io_l ) write(io_fid_log,'(1x,A)') &
720  '|============= Vertical Coordinate =============|'
721  if( io_l ) write(io_fid_log,'(1x,A)') &
722  '| k z zh dz buffer k |'
723  if( io_l ) write(io_fid_log,'(1x,A)') &
724  '| [m] [m] [m] factor |'
725 
726  do k = ka, ke+1, -1
727  if( io_l ) write(io_fid_log,'(1x,A,F9.2,A,F9.2,I5,A)') &
728  '| ',grid_fz(k),' ', grid_fbfz(k),k,' |'
729  if( io_l ) write(io_fid_log,'(1x,A,I5,F9.2,A,2F9.2,A)') &
730  '|',k,grid_cz(k),' ',grid_cdz(k), grid_cbfz(k),' |'
731  enddo
732 
733  k = ke
734  if( io_l ) write(io_fid_log,'(1x,A,F9.2,A,F9.2,I5,A)') &
735  '| ',grid_fz(k),' ', grid_fbfz(k),k,' | KE = TOA'
736  if( io_l ) write(io_fid_log,'(1x,A,I5,F9.2,A,2F9.2,A)') &
737  '|',k,grid_cz(k),' ',grid_cdz(k), grid_cbfz(k),' |'
738 
739  do k = ke-1, ks, -1
740  if( io_l ) write(io_fid_log,'(1x,A,F9.2,A,F9.2,I5,A)') &
741  '| ',grid_fz(k),' ', grid_fbfz(k),k,' |'
742  if( io_l ) write(io_fid_log,'(1x,A,I5,F9.2,A,2F9.2,A)') &
743  '|',k,grid_cz(k),' ',grid_cdz(k), grid_cbfz(k),' |'
744  enddo
745 
746  k = ks-1
747  if( io_l ) write(io_fid_log,'(1x,A,F9.2,A,F9.2,I5,A)') &
748  '| ',grid_fz(k),' ', grid_fbfz(k),k,' | KS-1 = surface'
749  if( io_l ) write(io_fid_log,'(1x,A,I5,F9.2,A,2F9.2,A)') &
750  '|',k,grid_cz(k),' ',grid_cdz(k), grid_cbfz(k),' |'
751 
752  do k = ks-2, 1, -1
753  if( io_l ) write(io_fid_log,'(1x,A,F9.2,A,F9.2,I5,A)') &
754  '| ',grid_fz(k),' ', grid_fbfz(k),k,' |'
755  if( io_l ) write(io_fid_log,'(1x,A,I5,F9.2,A,2F9.2,A)') &
756  '|',k,grid_cz(k),' ',grid_cdz(k), grid_cbfz(k),' |'
757  enddo
758 
759  k = 0
760  if( io_l ) write(io_fid_log,'(1x,A,F9.2,A,F9.2,I5,A)') &
761  '| ',grid_fz(k),' ', grid_fbfz(k),k,' |'
762 
763  if( io_l ) write(io_fid_log,'(1x,A)') &
764  '|===============================================|'
765 
766 ! if( IO_L ) write(IO_FID_LOG,*)
767 ! if( IO_L ) write(IO_FID_LOG,*) ' ', 0, GRID_FZ(0)
768 ! do k = 1, KA-1
769 ! if( IO_L ) write(IO_FID_LOG,*) k, GRID_CZ(k), GRID_CBFZ(k), GRID_CDZ(k)
770 ! if( IO_L ) write(IO_FID_LOG,*) ' ', k, GRID_FZ(k), GRID_FBFZ(k), GRID_FDZ(k)
771 ! enddo
772 ! k = KA
773 ! if( IO_L ) write(IO_FID_LOG,*) k, GRID_CZ(k), GRID_CBFZ(k), GRID_CDZ(k)
774 ! if( IO_L ) write(IO_FID_LOG,*) ' ', k, GRID_FZ(k), GRID_FBFZ(k)
775 !
776 ! if( IO_L ) write(IO_FID_LOG,*)
777 ! if( IO_L ) write(IO_FID_LOG,*) ' ', 0, GRID_FX(0)
778 ! do i = 1, IA-1
779 ! if( IO_L ) write(IO_FID_LOG,*) i, GRID_CX(i), GRID_CBFX(i), GRID_CDX(i)
780 ! if( IO_L ) write(IO_FID_LOG,*) ' ', i, GRID_FX(i), GRID_FBFX(i), GRID_FDX(i)
781 ! enddo
782 ! i = IA
783 ! if( IO_L ) write(IO_FID_LOG,*) i, GRID_CX(i), GRID_CBFX(i), GRID_CDX(i)
784 ! if( IO_L ) write(IO_FID_LOG,*) ' ', i, GRID_FX(i), GRID_FBFX(i)
785 !
786 ! if( IO_L ) write(IO_FID_LOG,*)
787 ! if( IO_L ) write(IO_FID_LOG,*) ' ', 0, GRID_FY(0)
788 ! do j = 1, JA-1
789 ! if( IO_L ) write(IO_FID_LOG,*) j, GRID_CY(j), GRID_CBFY(j), GRID_CDY(j)
790 ! if( IO_L ) write(IO_FID_LOG,*) ' ', j, GRID_FY(j), GRID_FBFY(j), GRID_FDY(j)
791 ! enddo
792 ! j = JA
793 ! if( IO_L ) write(IO_FID_LOG,*) j, GRID_CY(j), GRID_CBFY(j), GRID_CDY(j)
794 ! if( IO_L ) write(IO_FID_LOG,*) ' ', j, GRID_FY(j), GRID_FBFY(j)
795  endif
796 
797  return
798  end subroutine grid_generate
799 
800 end module scale_grid
scale_grid_index::imax
integer, public imax
of computational cells: x
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:234
scale_grid::grid_cbfyg
real(rp), dimension(:), allocatable, public grid_cbfyg
center buffer factor [0-1]: y, global
Definition: scale_grid_cartesian.F90:90
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:59
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_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_grid_index::ia
integer, public ia
of x whole cells (local, with HALO)
Definition: scale_grid_index.F90:55
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 z whole cells (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:87
scale_grid_index::kmax
integer, public kmax
of computational cells: z
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::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:88
scale_grid::grid_generate
subroutine, public grid_generate
Generate horizontal&vertical grid.
Definition: scale_grid_cartesian.F90:303
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:78
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::buffer_dz
real(rp), public buffer_dz
thickness of buffer region [m]: z
Definition: scale_grid_cartesian.F90:45
scale_stdio::io_lnml
logical, public io_lnml
output log or not? (for namelist, this process)
Definition: scale_stdio.F90:60
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:191
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
Definition: scale_grid_index.F90:45
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:117
scale_grid::grid_cbfxg
real(rp), dimension(:), allocatable, public grid_cbfxg
center buffer factor [0-1]: x, global
Definition: scale_grid_cartesian.F90:89
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 y whole cells (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.13