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