SCALE-RM
scale_atmos_dyn_tstep_large_fvm_heve.F90
Go to the documentation of this file.
1 !-------------------------------------------------------------------------------
10 !-------------------------------------------------------------------------------
11 #include "scalelib.h"
12 module scale_atmos_dyn_tstep_large_fvm_heve
13  !-----------------------------------------------------------------------------
14  !
15  !++ used modules
16  !
17  use scale_precision
18  use scale_io
19  use scale_prof
20  use scale_atmos_grid_cartesc_index
21  use scale_index
22  use scale_tracer
23 
24 #ifdef DEBUG
25  use scale_debug, only: &
26  check
27  use scale_const, only: &
28  undef => const_undef, &
29  iundef => const_undef2
30 #endif
31  !-----------------------------------------------------------------------------
32  implicit none
33  private
34  !-----------------------------------------------------------------------------
35  !
36  !++ Public procedure
37  !
38  public :: atmos_dyn_tstep_large_fvm_heve_setup
39  public :: atmos_dyn_tstep_large_fvm_heve
40 
41  !-----------------------------------------------------------------------------
42  !
43  !++ Public parameters & variables
44  !
45  !-----------------------------------------------------------------------------
46  !
47  !++ Private procedure
48  !
49  !-----------------------------------------------------------------------------
50  !
51  !++ Private parameters & variables
52  !
53  ! tendency
54  real(RP), private, allocatable :: DENS_t(:,:,:)
55  real(RP), private, allocatable :: MOMZ_t(:,:,:)
56  real(RP), private, allocatable :: MOMX_t(:,:,:)
57  real(RP), private, allocatable :: MOMY_t(:,:,:)
58  real(RP), private, allocatable :: RHOT_t(:,:,:)
59  real(RP), private, allocatable, target :: RHOQ_t(:,:,:,:)
60  real(RP), private, allocatable, target :: ZERO(:,:,:)
61  real(RP), private, pointer :: RHOQ_tn(:,:,:)
62 
63  real(RP), private, allocatable :: DENS_damp(:,:,:)
64 
65  ! flux
66  real(RP), private, allocatable, target :: mflx(:,:,:,:) ! rho * vel(x,y,z) * GSQRT / mapf
67 
68  ! for communication
69  integer :: I_COMM_DENS = 1
70  integer :: I_COMM_MOMZ = 2
71  integer :: I_COMM_MOMX = 3
72  integer :: I_COMM_MOMY = 4
73  integer :: I_COMM_RHOT = 5
74  integer, allocatable :: I_COMM_PROG(:)
75 
76  integer :: I_COMM_DENS_t = 1
77  integer :: I_COMM_MOMZ_t = 2
78  integer :: I_COMM_MOMX_t = 3
79  integer :: I_COMM_MOMY_t = 4
80  integer :: I_COMM_RHOT_t = 5
81 
82  integer :: I_COMM_DENS_damp = 6
83 
84  integer, allocatable :: I_COMM_RHOQ_t(:)
85  integer, allocatable :: I_COMM_QTRC(:)
86 
87  integer :: I_COMM_mflx_z = 1
88  integer :: I_COMM_mflx_x = 2
89  integer :: I_COMM_mflx_y = 3
90 
91  ! for history
92  integer :: HIST_mflx(3)
93  integer :: HIST_tflx(3)
94  integer :: HIST_phys(5)
95  integer :: HIST_damp(5)
96  integer, allocatable :: HIST_qflx(:,:)
97  integer, allocatable :: HIST_phys_QTRC(:)
98  integer, allocatable :: HIST_damp_QTRC(:)
99 
100  ! for monitor
101  real(RP), allocatable, target :: zero_x(:,:)
102  real(RP), allocatable, target :: zero_y(:,:)
103  integer :: MONIT_damp_mass
104  integer :: MONIT_damp_qtot
105  integer :: MONIT_mflx_west
106  integer :: MONIT_mflx_east
107  integer :: MONIT_mflx_south
108  integer :: MONIT_mflx_north
109  integer :: MONIT_qflx_west
110  integer :: MONIT_qflx_east
111  integer :: MONIT_qflx_south
112  integer :: MONIT_qflx_north
113 
114  character(len=H_SHORT) :: EVAL_TYPE_NUMFILTER = 'TENDENCY'
115 
116  !-----------------------------------------------------------------------------
117 contains
118 
119  !-----------------------------------------------------------------------------
121  subroutine atmos_dyn_tstep_large_fvm_heve_setup( &
122  DENS, MOMZ, MOMX, MOMY, RHOT, QTRC, PROG )
123  use scale_prc, only: &
124  prc_abort
125  use scale_prc_cartesc, only: &
126  prc_twod, &
127  prc_has_e, &
128  prc_has_w, &
129  prc_has_n, &
130  prc_has_s
131  use scale_const, only: &
132  ohm => const_ohm, &
133  undef => const_undef
134  use scale_comm_cartesc, only: &
135  comm_vars8_init
136  use scale_file_history, only: &
137  file_history_reg, &
138  file_history_put
139  use scale_monitor, only: &
140  monitor_reg, &
141  monitor_put
142  implicit none
143 
144  ! MPI_RECV_INIT requires intent(inout)
145  real(rp), intent(inout) :: dens(ka,ia,ja)
146  real(rp), intent(inout) :: momz(ka,ia,ja)
147  real(rp), intent(inout) :: momx(ka,ia,ja)
148  real(rp), intent(inout) :: momy(ka,ia,ja)
149  real(rp), intent(inout) :: rhot(ka,ia,ja)
150  real(rp), intent(inout) :: qtrc(ka,ia,ja,qa)
151  real(rp), intent(inout) :: prog(ka,ia,ja,va)
152 
153  integer :: iv, iq
154 
155  namelist /atmos_dyn_tstep_large_fvm_heve/ &
156  eval_type_numfilter
157 
158  integer :: ierr
159  !---------------------------------------------------------------------------
160 
161  !--- read namelist
162  rewind(io_fid_conf)
163  read(io_fid_conf,nml=atmos_dyn_tstep_large_fvm_heve,iostat=ierr)
164  if( ierr < 0 ) then !--- missing
165  log_info("ATMOS_DYN_Tstep_large_fvm_heve_setup",*) 'Not found namelist. Default used.'
166  elseif( ierr > 0 ) then !--- fatal error
167  log_error("ATMOS_DYN_Tstep_large_fvm_heve_setup",*) 'Not appropriate names in namelist ATMOS_DYN_TSTEP_LARGE_FVM_HEVE. Check!'
168  call prc_abort
169  endif
170  log_nml(atmos_dyn_tstep_large_fvm_heve)
171 
172  select case( eval_type_numfilter )
173  case( 'TENDENCY', 'FILTER' )
174  case default
175  log_error("ATMOS_DYN_Tstep_large_fvm_heve_setup",*) 'The specfied value of EVAL_TYPE_NUMFILTER is not appropriate. Check!'
176  call prc_abort
177  end select
178  !--
179 
180  allocate( dens_t(ka,ia,ja) )
181  allocate( momz_t(ka,ia,ja) )
182  allocate( momx_t(ka,ia,ja) )
183  allocate( momy_t(ka,ia,ja) )
184  allocate( rhot_t(ka,ia,ja) )
185  allocate( rhoq_t(ka,ia,ja,qa) )
186 
187  allocate( dens_damp(ka,ia,ja) )
188 
189  allocate( mflx(ka,ia,ja,3) )
190 
191  allocate( i_comm_prog(max(va,1)) )
192  allocate( i_comm_qtrc(qa) )
193  allocate( i_comm_rhoq_t(qa) )
194 
195  call comm_vars8_init( 'DENS', dens, i_comm_dens )
196  call comm_vars8_init( 'MOMZ', momz, i_comm_momz )
197  call comm_vars8_init( 'MOMX', momx, i_comm_momx )
198  call comm_vars8_init( 'MOMY', momy, i_comm_momy )
199  call comm_vars8_init( 'RHOT', rhot, i_comm_rhot )
200  do iv = 1, va
201  i_comm_prog(iv) = 5 + iv
202  call comm_vars8_init( 'PROG', prog(:,:,:,iv), i_comm_prog(iv) )
203  end do
204 
205  call comm_vars8_init( 'DENS_t', dens_t, i_comm_dens_t )
206  call comm_vars8_init( 'MOMZ_t', momz_t, i_comm_momz_t )
207  call comm_vars8_init( 'MOMX_t', momx_t, i_comm_momx_t )
208  call comm_vars8_init( 'MOMY_t', momy_t, i_comm_momy_t )
209  call comm_vars8_init( 'RHOT_t', rhot_t, i_comm_rhot_t )
210 
211  do iq = 1, qa
212  i_comm_rhoq_t(iq) = 5 + va + iq
213  i_comm_qtrc(iq) = 5 + va + iq
214 
215  call comm_vars8_init( 'RHOQ_t', rhoq_t(:,:,:,iq), i_comm_rhoq_t(iq) )
216  call comm_vars8_init( 'QTRC', qtrc(:,:,:,iq), i_comm_qtrc(iq) )
217  end do
218 
219  call comm_vars8_init( 'DENS_t', dens_t, i_comm_dens_damp )
220 
221  call comm_vars8_init( 'mflx_Z', mflx(:,:,:,zdir), i_comm_mflx_z )
222  call comm_vars8_init( 'mflx_X', mflx(:,:,:,xdir), i_comm_mflx_x )
223  call comm_vars8_init( 'mflx_Y', mflx(:,:,:,ydir), i_comm_mflx_y )
224 
225  allocate( zero(ka,ia,ja) )
226  zero(:,:,:) = 0.0_rp
227 
228  mflx(:,:,:,:) = undef
229  if ( prc_twod ) then
230  mflx(:,:,:,xdir) = 0.0_rp
231  end if
232 
233 
234  ! history
235  call file_history_reg( 'ZFLX_MOM', 'momentum flux of z-direction', 'kg/m2/s', & ! [IN]
236  hist_mflx(1), & ! [OUT]
237  dim_type='ZHXY' ) ! [IN]
238  call file_history_reg( 'XFLX_MOM', 'momentum flux of x-direction', 'kg/m2/s', & ! [IN]
239  hist_mflx(2), & ! [OUT]
240  dim_type='ZXHY' ) ! [IN]
241  call file_history_reg( 'YFLX_MOM', 'momentum flux of y-direction', 'kg/m2/s', & ! [IN]
242  hist_mflx(3), & ! [OUT]
243  dim_type='ZXYH' ) ! [IN]
244 
245  call file_history_reg( 'ZFLX_RHOT', 'potential temperature flux of z-direction', 'K*kg/m2/s', & ! [IN]
246  hist_tflx(1), & ! [OUT]
247  dim_type='ZHXY' ) ! [IN]
248  call file_history_reg( 'XFLX_RHOT', 'potential temperature flux of x-direction', 'K*kg/m2/s', & ! [IN]
249  hist_tflx(2), & ! [OUT]
250  dim_type='ZXHY' ) ! [IN]
251  call file_history_reg( 'YFLX_RHOT', 'potential temperature flux of y-direction', 'K*kg/m2/s', & ! [IN]
252  hist_tflx(3), & ! [OUT]
253  dim_type='ZXYH' ) ! [IN]
254 
255  call file_history_reg( 'DENS_t_phys', 'tendency of dencity due to physics', 'kg/m3/s', & ! [IN]
256  hist_phys(1) ) ! [OUT]
257  call file_history_reg( 'MOMZ_t_phys', 'tendency of momentum z due to physics', 'kg/m2/s2', & ! [IN]
258  hist_phys(2), & ! [OUT]
259  dim_type='ZHXY' ) ! [IN]
260  call file_history_reg( 'MOMX_t_phys', 'tendency of momentum x due to physics', 'kg/m2/s2', & ! [IN]
261  hist_phys(3), & ! [OUT]
262  dim_type='ZXHY' ) ! [IN]
263  call file_history_reg( 'MOMY_t_phys', 'tendency of momentum y due to physics', 'kg/m2/s2', & ! [IN]
264  hist_phys(4), & ! [OUT]
265  dim_type='ZXYH' ) ! [IN]
266  call file_history_reg( 'RHOT_t_phys', 'tendency of rho*theta temperature due to physics', 'K*kg/m3/s', & ! [IN]
267  hist_phys(5) ) ! [OUT]
268 
269  call file_history_reg( 'DENS_t_damp', 'tendency of dencity due to damping', 'kg/m3/s', & ! [IN]
270  hist_damp(1) ) ! [OUT]
271  call file_history_reg( 'MOMZ_t_damp', 'tendency of momentum z due to damping', 'kg/m2/s2', & ! [IN]
272  hist_damp(2), & ! [OUT]
273  dim_type='ZHXY' ) ! [IN]
274  call file_history_reg( 'MOMX_t_damp', 'tendency of momentum x due to damping', 'kg/m2/s2', & ! [IN]
275  hist_damp(3), & ! [OUT]
276  dim_type='ZXHY' ) ! [IN]
277  call file_history_reg( 'MOMY_t_damp', 'tendency of momentum y due to damping', 'kg/m2/s2', & ! [IN]
278  hist_damp(4), & ! [OUT]
279  dim_type='ZXYH' ) ! [IN]
280  call file_history_reg( 'RHOT_t_damp', 'tendency of rho*theta temperature due to damping', 'K kg/m3/s', & ! [IN]
281  hist_damp(5) ) ! [OUT]
282 
283  allocate( hist_qflx(3,qa) )
284  allocate( hist_phys_qtrc(qa) )
285  allocate( hist_damp_qtrc(qa) )
286  do iq = 1, qa
287  call file_history_reg( 'ZFLX_'//trim(tracer_name(iq)), trim(tracer_name(iq))//' flux of z-direction', 'kg/m2/s', & ! [IN]
288  hist_qflx(1,iq), & ! [OUT]
289  dim_type='ZHXY' ) ! [IN]
290  call file_history_reg( 'XFLX_'//trim(tracer_name(iq)), trim(tracer_name(iq))//' flux of x-direction', 'kg/m2/s', & ! [IN]
291  hist_qflx(2,iq), & ! [OUT]
292  dim_type='ZXHY' ) ! [IN]
293  call file_history_reg( 'YFLX_'//trim(tracer_name(iq)), trim(tracer_name(iq))//' flux of y-direction', 'kg/m2/s', & ! [IN]
294  hist_qflx(3,iq), & ! [OUT]
295  dim_type='ZXYH' ) ! [IN]
296  call file_history_reg( trim(tracer_name(iq))//'_t_phys', 'tendency of '//trim(tracer_name(iq))//' due to physics', 'kg/m3/s', &
297  hist_phys_qtrc(iq) )
298  call file_history_reg( trim(tracer_name(iq))//'_t_damp', 'tendency of '//trim(tracer_name(iq))//' due to damping', 'kg/m3/s', &
299  hist_damp_qtrc(iq) )
300  end do
301 
302  ! for history at t=0
303  do iv = 1, 3
304  call file_history_put( hist_mflx(iv), zero(:,:,:) )
305  call file_history_put( hist_tflx(iv), zero(:,:,:) )
306  end do
307  do iv = 1, 5
308  call file_history_put( hist_phys(iv), zero(:,:,:) )
309  call file_history_put( hist_damp(iv), zero(:,:,:) )
310  end do
311  do iq = 1, qa
312  do iv = 1, 3
313  call file_history_put( hist_qflx(iv,iq), zero(:,:,:) )
314  end do
315  call file_history_put( hist_phys_qtrc(iq), zero(:,:,:) )
316  call file_history_put( hist_damp_qtrc(iq), zero(:,:,:) )
317  end do
318 
319 
320  ! for monitor
321  allocate( zero_x(ka,ja) )
322  allocate( zero_y(ka,ia) )
323  zero_x(:,:) = 0.0_rp
324  zero_y(:,:) = 0.0_rp
325 
326  call monitor_reg( "MASSTND_DAMP", "mass tendency by the damping", "kg", & ! [IN]
327  monit_damp_mass, & ! [OUT]
328  is_tendency=.true. ) ! [IN]
329 
330  call monitor_reg( "MASSFLX_WEST", "mass flux at the western boundary", "kg", & ! [IN]
331  monit_mflx_west, & ! [OUT]
332  dim_type="ZY-W", is_tendency=.true. ) ! [IN]
333  call monitor_reg( "MASSFLX_EAST", "mass flux at the eastern boundary", "kg", & ! [IN]
334  monit_mflx_east, & ! [OUT]
335  dim_type="ZY-E", is_tendency=.true. ) ! [IN]
336  call monitor_reg( "MASSFLX_SOUTH", "mass flux at the southern boundary", "kg", & ! [IN]
337  monit_mflx_south, & ! [OUT]
338  dim_type="ZX-S", is_tendency=.true. ) ! [IN]
339  call monitor_reg( "MASSFLX_NORTH", "mass flux at the northern boundary", "kg", & ! [IN]
340  monit_mflx_north, & ! [OUT]
341  dim_type="ZX-N", is_tendency=.true. ) ! [IN]
342 
343  call monitor_reg( "QTOTTND_DAMP", "water mass tendency by the damping", "kg", & ! [IN]
344  monit_damp_qtot, & ! [OUT]
345  is_tendency=.true. ) ! [IN]
346 
347  call monitor_reg( "QTOTFLX_WEST", "water mass flux at the western boundary", "kg", & ! [IN]
348  monit_qflx_west, & ! [OUT]
349  dim_type="ZY-W", is_tendency=.true. ) ! [IN]
350  call monitor_reg( "QTOTFLX_EAST", "water mass flux at the eastern boundary", "kg", & ! [IN]
351  monit_qflx_east, & ! [OUT]
352  dim_type="ZY-E", is_tendency=.true. ) ! [IN]
353  call monitor_reg( "QTOTFLX_SOUTH", "water mass flux at the southern boundary", "kg", & ! [IN]
354  monit_qflx_south, & ! [OUT]
355  dim_type="ZX-S", is_tendency=.true. ) ! [IN]
356  call monitor_reg( "QTOTFLX_NORTH", "water mass flux at the northern boundary", "kg", & ! [IN]
357  monit_qflx_north, & ! [OUT]
358  dim_type="ZX-N", is_tendency=.true. ) ! [IN]
359 
360  ! at t=0
361  call monitor_put( monit_damp_mass, zero(:,:,:) )
362  call monitor_put( monit_mflx_west, zero_x(:,:) )
363  call monitor_put( monit_mflx_east, zero_x(:,:) )
364  call monitor_put( monit_mflx_south, zero_y(:,:) )
365  call monitor_put( monit_mflx_north, zero_y(:,:) )
366 
367  call monitor_put( monit_damp_qtot, zero(:,:,:) )
368  call monitor_put( monit_qflx_west, zero_x(:,:) )
369  call monitor_put( monit_qflx_east, zero_x(:,:) )
370  call monitor_put( monit_qflx_south, zero_y(:,:) )
371  call monitor_put( monit_qflx_north, zero_y(:,:) )
372 
373  return
374  end subroutine atmos_dyn_tstep_large_fvm_heve_setup
375 
376  !-----------------------------------------------------------------------------
378  subroutine atmos_dyn_tstep_large_fvm_heve( &
379  DENS, MOMZ, MOMX, MOMY, RHOT, QTRC, PROG, &
380  DENS_av, MOMZ_av, MOMX_av, MOMY_av, RHOT_av, QTRC_av, &
381  num_diff, num_diff_q, &
382  QTRC0, &
383  DENS_tp, MOMZ_tp, MOMX_tp, MOMY_tp, RHOT_tp, RHOQ_tp, &
384  CORIOLI, &
385  CDZ, CDX, CDY, FDZ, FDX, FDY, &
386  RCDZ, RCDX, RCDY, RFDZ, RFDX, RFDY, &
387  PHI, GSQRT, &
388  J13G, J23G, J33G, MAPF, &
389  AQ_R, AQ_CV, AQ_CP, AQ_MASS, &
390  REF_dens, REF_pott, REF_qv, REF_pres, &
391  BND_W, BND_E, BND_S, BND_N, TwoD, &
392  ND_COEF, ND_COEF_Q, ND_LAPLACIAN_NUM, &
393  ND_SFC_FACT, ND_USE_RS, &
394  BND_QA, BND_IQ, BND_SMOOTHER_FACT, &
395  DAMP_DENS, DAMP_VELZ, DAMP_VELX, &
396  DAMP_VELY, DAMP_POTT, DAMP_QTRC, &
397  DAMP_alpha_DENS, DAMP_alpha_VELZ, DAMP_alpha_VELX, &
398  DAMP_alpha_VELY, DAMP_alpha_POTT, DAMP_alpha_QTRC, &
399  MFLUX_OFFSET_X, MFLUX_OFFSET_Y, &
400  wdamp_coef, divdmp_coef, &
401  FLAG_TRACER_SPLIT_TEND, &
402  FLAG_FCT_MOMENTUM, FLAG_FCT_T, FLAG_FCT_TRACER, &
403  FLAG_FCT_ALONG_STREAM, &
404  USE_AVERAGE, &
405  I_QV, &
406  DTLS, DTSS, Llast )
407  use scale_const, only: &
408  eps => const_eps
409  use scale_comm_cartesc, only: &
410  comm_vars8, &
411  comm_wait
412  use scale_atmos_dyn_common, only: &
413  atmos_dyn_fill_halo, &
414  atmos_dyn_prep_pres_linearization, &
415  atmos_dyn_divergence
416  use scale_atmos_dyn_fvm_numfilter, only: &
417  atmos_dyn_fvm_numfilter_flux, &
418  atmos_dyn_fvm_numfilter_flux_q, &
419  atmos_dyn_fvm_apply_numfilter
420  use scale_file_history, only: &
421  file_history_query, &
422  file_history_put, &
423  file_history_set_disable
424  use scale_monitor, only: &
425  monitor_put
426  use scale_atmos_dyn_tinteg_short, only: &
427  atmos_dyn_tinteg_short
428  use scale_atmos_dyn_tinteg_tracer, only: &
429  atmos_dyn_tinteg_tracer
430  use scale_spnudge, only: &
431  spnudge_uv, &
432  spnudge_uv_divfree, &
433  spnudge_pt, &
434  spnudge_qv, &
435  spnudge_u_alpha, &
436  spnudge_v_alpha, &
437  spnudge_pt_alpha, &
438  spnudge_qv_alpha, &
439  spnudge_uv_lm, &
440  spnudge_uv_mm, &
441  spnudge_pt_lm, &
442  spnudge_pt_mm, &
443  spnudge_qv_lm, &
444  spnudge_qv_mm
445  use scale_dft, only: &
446  dft_g2g_divfree, &
447  dft_g2g
448  implicit none
449 
450  real(rp), intent(inout) :: dens(ka,ia,ja)
451  real(rp), intent(inout) :: momz(ka,ia,ja)
452  real(rp), intent(inout) :: momx(ka,ia,ja)
453  real(rp), intent(inout) :: momy(ka,ia,ja)
454  real(rp), intent(inout) :: rhot(ka,ia,ja)
455  real(rp), intent(inout) :: qtrc(ka,ia,ja,qa)
456  real(rp), intent(inout) :: prog(ka,ia,ja,va)
457 
458  real(rp), intent(inout) :: dens_av(ka,ia,ja)
459  real(rp), intent(inout) :: momz_av(ka,ia,ja)
460  real(rp), intent(inout) :: momx_av(ka,ia,ja)
461  real(rp), intent(inout) :: momy_av(ka,ia,ja)
462  real(rp), intent(inout) :: rhot_av(ka,ia,ja)
463  real(rp), intent(inout) :: qtrc_av(ka,ia,ja,qa)
464 
465  real(rp), intent(out) :: num_diff(ka,ia,ja,5,3)
466  real(rp), intent(out) :: num_diff_q(ka,ia,ja,3)
467 
468  real(rp), intent(in) :: qtrc0(ka,ia,ja,qa)
469 
470  real(rp), intent(in) :: dens_tp(ka,ia,ja)
471  real(rp), intent(in) :: momz_tp(ka,ia,ja)
472  real(rp), intent(in) :: momx_tp(ka,ia,ja)
473  real(rp), intent(in) :: momy_tp(ka,ia,ja)
474  real(rp), intent(in) :: rhot_tp(ka,ia,ja)
475  real(rp), intent(in) :: rhoq_tp(ka,ia,ja,qa)
476 
477  real(rp), intent(in) :: corioli(ia,ja)
478 
479  real(rp), intent(in) :: cdz (ka)
480  real(rp), intent(in) :: cdx (ia)
481  real(rp), intent(in) :: cdy (ja)
482  real(rp), intent(in) :: fdz (ka-1)
483  real(rp), intent(in) :: fdx (ia-1)
484  real(rp), intent(in) :: fdy (ja-1)
485  real(rp), intent(in) :: rcdz(ka)
486  real(rp), intent(in) :: rcdx(ia)
487  real(rp), intent(in) :: rcdy(ja)
488  real(rp), intent(in) :: rfdz(ka-1)
489  real(rp), intent(in) :: rfdx(ia-1)
490  real(rp), intent(in) :: rfdy(ja-1)
491 
492  real(rp), intent(in) :: phi (ka,ia,ja)
493  real(rp), intent(in) :: gsqrt(ka,ia,ja,7)
494  real(rp), intent(in) :: j13g (ka,ia,ja,7)
495  real(rp), intent(in) :: j23g (ka,ia,ja,7)
496  real(rp), intent(in) :: j33g
497  real(rp), intent(in) :: mapf (ia,ja,2,4)
498 
499  real(rp), intent(in) :: aq_r (qa)
500  real(rp), intent(in) :: aq_cv (qa)
501  real(rp), intent(in) :: aq_cp (qa)
502  real(rp), intent(in) :: aq_mass(qa)
503 
504  real(rp), intent(in) :: ref_dens(ka,ia,ja)
505  real(rp), intent(in) :: ref_pott(ka,ia,ja)
506  real(rp), intent(in) :: ref_qv (ka,ia,ja)
507  real(rp), intent(in) :: ref_pres(ka,ia,ja)
508 
509  logical, intent(in) :: bnd_w
510  logical, intent(in) :: bnd_e
511  logical, intent(in) :: bnd_s
512  logical, intent(in) :: bnd_n
513  logical, intent(in) :: twod
514 
515  real(rp), intent(in) :: nd_coef
516  real(rp), intent(in) :: nd_coef_q
517  integer, intent(in) :: nd_laplacian_num
518  real(rp), intent(in) :: nd_sfc_fact
519  logical, intent(in) :: nd_use_rs
520 
521  integer, intent(in) :: bnd_qa
522  integer, intent(in) :: bnd_iq(qa)
523  real(rp), intent(in) :: bnd_smoother_fact
524 
525  real(rp), intent(in) :: damp_dens(ka,ia,ja)
526  real(rp), intent(in) :: damp_velz(ka,ia,ja)
527  real(rp), intent(in) :: damp_velx(ka,ia,ja)
528  real(rp), intent(in) :: damp_vely(ka,ia,ja)
529  real(rp), intent(in) :: damp_pott(ka,ia,ja)
530  real(rp), intent(in) :: damp_qtrc(ka,ia,ja,bnd_qa)
531 
532  real(rp), intent(in) :: damp_alpha_dens(ka,ia,ja)
533  real(rp), intent(in) :: damp_alpha_velz(ka,ia,ja)
534  real(rp), intent(in) :: damp_alpha_velx(ka,ia,ja)
535  real(rp), intent(in) :: damp_alpha_vely(ka,ia,ja)
536  real(rp), intent(in) :: damp_alpha_pott(ka,ia,ja)
537  real(rp), intent(in) :: damp_alpha_qtrc(ka,ia,ja,bnd_qa)
538  real(rp), intent(in) :: mflux_offset_x(ka,ja,2)
539  real(rp), intent(in) :: mflux_offset_y(ka,ia,2)
540 
541  real(rp), intent(in) :: wdamp_coef(ka)
542  real(rp), intent(in) :: divdmp_coef
543 
544  logical, intent(in) :: flag_tracer_split_tend
545  logical, intent(in) :: flag_fct_momentum
546  logical, intent(in) :: flag_fct_t
547  logical, intent(in) :: flag_fct_tracer
548  logical, intent(in) :: flag_fct_along_stream
549 
550  logical, intent(in) :: use_average
551 
552  integer, intent(in) :: i_qv
553 
554  real(dp), intent(in) :: dtls
555  real(dp), intent(in) :: dtss
556  logical , intent(in) :: llast
557 
558  ! for time integartion
559  real(rp) :: dens00 (ka,ia,ja) ! saved density before small step loop
560  real(rp) :: qflx (ka,ia,ja,3)
561 
562  ! diagnostic variables
563  real(rp) :: ddiv (ka,ia,ja) ! 3 dimensional divergence
564  real(rp) :: dpres0 (ka,ia,ja) ! pressure deviation
565  real(rp) :: rt2p (ka,ia,ja) ! factor of RHOT to PRES
566  real(rp) :: ref_rhot(ka,ia,ja) ! reference of RHOT
567 
568  real(rp) :: dens_tq(ka,ia,ja)
569  real(rp) :: diff(ka,ia,ja), diff2(ka,ia,ja), diff3(ka,ia,ja)
570  real(rp) :: damp
571  real(rp) :: damp_t_dens(ka,ia,ja)
572  real(rp) :: damp_t_momz(ka,ia,ja)
573  real(rp) :: damp_t_momx(ka,ia,ja)
574  real(rp) :: damp_t_momy(ka,ia,ja)
575  real(rp) :: damp_t_rhot(ka,ia,ja)
576  real(rp) :: damp_t_qtrc(ka,ia,ja)
577 
578  real(rp) :: tflx(ka,ia,ja,3)
579 
580  ! For tracer advection
581  real(rp) :: mflx_av(ka,ia,ja,3) ! rho * vel(x,y,z) @ (u,v,w)-face average
582 
583  real(rp) :: dtl
584  real(rp) :: dts
585  integer :: nstep
586 
587  ! for history
588  logical :: do_put
589 
590  ! for monitor
591  real(rp), target :: qflx_west (ka,ja)
592  real(rp), target :: qflx_east (ka,ja)
593  real(rp), target :: qflx_south(ka,ia)
594  real(rp), target :: qflx_north(ka,ia)
595  logical :: monit_lateral_flag(3)
596 
597  integer :: i, j, k, iq, iqb, step
598  integer :: iv
599  integer :: n
600  !---------------------------------------------------------------------------
601 
602  call prof_rapstart("DYN_Large_Preparation", 2)
603 
604  dts = real(dtss, kind=rp) ! short time step
605  dtl = real(dtls, kind=rp) ! large time step
606  nstep = ceiling( ( dtl - eps ) / dts )
607  dts = dtl / nstep ! dts is divisor of dtl and smaller or equal to dtss
608 
609  monit_lateral_flag(zdir) = .false.
610  monit_lateral_flag(xdir) = monit_mflx_west > 0 .or. monit_mflx_east > 0
611  monit_lateral_flag(ydir) = monit_mflx_south > 0 .or. monit_mflx_north > 0
612 
613 #ifdef DEBUG
614  log_info("ATMOS_DYN_Tstep_large_fvm_heve",*) 'Dynamics large time step'
615  log_info_cont('(1x,A,F0.2,A,F0.2,A,I0)') &
616  '-> DT_large, DT_small, DT_large/DT_small : ', dtl, ', ', dts, ', ', nstep
617 
618  dens00(:,:,:) = undef
619  num_diff(:,:,:,:,:) = undef
620 #endif
621 
622 !OCL XFILL
623  dens00(:,:,:) = dens(:,:,:)
624 
625  if ( use_average ) then
626 !OCL XFILL
627  !$omp parallel do
628  do j = jsb, jeb
629  do i = isb, ieb
630  do k = ks, ke
631  dens_av(k,i,j) = 0.0_rp
632  momz_av(k,i,j) = 0.0_rp
633  momx_av(k,i,j) = 0.0_rp
634  momy_av(k,i,j) = 0.0_rp
635  rhot_av(k,i,j) = 0.0_rp
636  end do
637  end do
638  end do
639  endif
640 
641 !OCL XFILL
642  mflx_av(:,:,:,:) = 0.0_rp
643 
644 
645 !OCL XFILL
646  !$omp parallel do
647  do j = js, je
648  do k = ks, ke
649  qflx_west(k,j) = 0.0_rp
650  qflx_east(k,j) = 0.0_rp
651  end do
652  end do
653 !OCL XFILL
654  !$omp parallel do
655  do i = is, ie
656  do k = ks, ke
657  qflx_south(k,i) = 0.0_rp
658  qflx_north(k,i) = 0.0_rp
659  end do
660  end do
661 
662  !- prepare some variables for pressure linearization
663  call atmos_dyn_prep_pres_linearization( &
664  dpres0, rt2p, ref_rhot, & ! (out)
665  rhot, qtrc, ref_pres, aq_r, aq_cv, aq_cp, aq_mass ) ! (in)
666 
667  call prof_rapend ("DYN_Large_Preparation", 2)
668 
669  !###########################################################################
670  ! Update DENS,MONZ,MOMX,MOMY,MOMZ,RHOT
671  !###########################################################################
672 
673  call prof_rapstart("DYN_Large_Tendency", 2)
674 
675 !OCL XFILL
676  dens_tq(:,:,:) = 0.0_rp
677 
678  do iq = 1, qa
679 
680  iqb = bnd_iq(iq)
681 
682  if ( iqb > 0 ) then
683 
684  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
685 !OCL XFILL
686  do j = js-1, je+1
687  do i = max(is-1,1), min(ie+1,ia)
688  do k = ks, ke
689  diff(k,i,j) = qtrc(k,i,j,iq) - damp_qtrc(k,i,j,iqb)
690  enddo
691  enddo
692  enddo
693 
694  call file_history_query( hist_damp_qtrc(iq), do_put )
695  if ( twod ) then
696  !$omp parallel do default(none) OMP_SCHEDULE_ &
697  !$omp private(j,k,damp) &
698  !$omp shared(IS,JS,JE,KS,KE,iq,iqb) &
699  !$omp shared(RHOQ_t,RHOQ_tp,DENS_tq,DAMP_alpha_QTRC,diff,BND_SMOOTHER_FACT,DENS00,TRACER_MASS,I_QV) &
700  !$omp shared(damp_t_QTRC,do_put)
701 !OCL XFILL
702  do j = js, je
703  do k = ks, ke
704  damp = - damp_alpha_qtrc(k,is,j,iqb) &
705  * ( diff(k,is,j) & ! rayleigh damping
706  - ( diff(k,is,j-1) + diff(k,is,j+1) - diff(k,is,j)*2.0_rp ) &
707  * 0.25_rp * bnd_smoother_fact ) ! horizontal smoother
708  damp = damp * dens00(k,is,j)
709  if ( do_put ) damp_t_qtrc(k,is,j) = damp
710  rhoq_t(k,is,j,iq) = rhoq_tp(k,is,j,iq) + damp
711  dens_tq(k,is,j) = dens_tq(k,is,j) + damp * tracer_mass(iq) ! only for mass tracer
712  enddo
713  enddo
714  else
715  if( iq == i_qv .and. spnudge_qv ) then
716 
717  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
718 !OCL XFILL
719  do j = js, je
720  do i = is, ie
721  do k = ks, ke
722  diff2(k,i,j) = diff(k,i,j)
723  enddo
724  enddo
725  enddo
726 
727  call dft_g2g(ka,ks,ke,ia,is,ie,ja,js,je,spnudge_qv_lm,spnudge_qv_mm,diff2)
728 
729  else
730 
731  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
732 !OCL XFILL
733  do j = js, je
734  do i = is, ie
735  do k = ks, ke
736  diff2(k,i,j) = 0
737  enddo
738  enddo
739  enddo
740 
741  endif
742 
743  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
744  !$omp private(i,j,k,damp) &
745  !$omp shared(JS,JE,IS,IE,KS,KE,iq,iqb) &
746  !$omp shared(RHOQ_t,RHOQ_tp,DENS_tq,DAMP_alpha_QTRC,diff,diff2,BND_SMOOTHER_FACT,SPNUDGE_qv_alpha,DENS00,TRACER_MASS) &
747  !$omp shared(damp_t_QTRC,do_put)
748 !OCL XFILL
749  do j = js, je
750  do i = is, ie
751  do k = ks, ke
752  damp = - damp_alpha_qtrc(k,i,j,iqb) &
753  * ( diff(k,i,j) & ! rayleigh damping
754  - ( diff(k,i-1,j) + diff(k,i+1,j) + diff(k,i,j-1) + diff(k,i,j+1) - diff(k,i,j)*4.0_rp ) &
755  * 0.125_rp * bnd_smoother_fact ) & ! horizontal smoother
756  - spnudge_qv_alpha(k,i,j) * diff2(k,i,j)
757  damp = damp * dens00(k,i,j)
758  if ( do_put ) damp_t_qtrc(k,i,j) = damp
759  rhoq_t(k,i,j,iq) = rhoq_tp(k,i,j,iq) + damp
760  dens_tq(k,i,j) = dens_tq(k,i,j) + damp * tracer_mass(iq) ! only for mass tracer
761  enddo
762  enddo
763  enddo
764  end if
765 
766  if ( llast ) then
767  if ( do_put ) call file_history_put( hist_damp_qtrc(iq), damp_t_qtrc(:,:,:) )
768  call file_history_put( hist_phys_qtrc(iq), rhoq_tp(:,:,:,iq) )
769  end if
770 
771  call atmos_dyn_fill_halo( rhoq_t(:,:,:,iq), 0.0_rp, .false., .true. )
772  call comm_vars8( rhoq_t(:,:,:,iq), i_comm_rhoq_t(iq) )
773  call comm_wait ( rhoq_t(:,:,:,iq), i_comm_rhoq_t(iq), .false. )
774 
775  else
776 
777  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(3)
778 !OCL XFILL
779  do j = 1, ja
780  do i = 1, ia
781  do k = 1, ka
782  rhoq_t(k,i,j,iq) = rhoq_tp(k,i,j,iq)
783  enddo
784  enddo
785  enddo
786 
787  end if
788 
789  end do
790 
791  call prof_rapend ("DYN_Large_Tendency", 2)
792 
793  call prof_rapstart("DYN_Large_Boundary", 2)
794 
795  if ( bnd_w .and. (.not. twod) ) then
796  !$omp parallel do private(j,k) OMP_SCHEDULE_
797  do j = js, je
798  do k = ks, ke
799  mflx(k,is-1,j,xdir) = ( momx(k,is-1,j) + mflux_offset_x(k,j,1) ) &
800  * gsqrt(k,is-1,j,i_uyz) / mapf(is-1,j,2,i_uy)
801  enddo
802  enddo
803  end if
804  if ( bnd_e .and. (.not. twod) ) then
805  !$omp parallel do private(j,k) OMP_SCHEDULE_
806  do j = js, je
807  do k = ks, ke
808  mflx(k,ie,j,xdir) = ( momx(k,ie,j) + mflux_offset_x(k,j,2) ) &
809  * gsqrt(k,ie,j,i_uyz) / mapf(ie,j,2,i_uy)
810  enddo
811  enddo
812  end if
813  if ( bnd_s ) then
814  !$omp parallel do private(i,k) OMP_SCHEDULE_
815  do i = is, ie
816  do k = ks, ke
817  mflx(k,i,js-1,ydir) = ( momy(k,i,js-1) + mflux_offset_y(k,i,1) ) &
818  * gsqrt(k,i,js-1,i_xvz) / mapf(i,js-1,1,i_xv)
819  enddo
820  enddo
821  end if
822  if ( bnd_n ) then
823  !$omp parallel do private(i,k) OMP_SCHEDULE_
824  do i = is, ie
825  do k = ks, ke
826  mflx(k,i,je,ydir) = ( momy(k,i,je) + mflux_offset_y(k,i,2) ) &
827  * gsqrt(k,i,je,i_xvz) / mapf(i,je,1,i_xv)
828  enddo
829  enddo
830  end if
831 
832  call prof_rapend ("DYN_Large_Boundary", 2)
833 
834 !OCL XFILL
835  !$omp parallel do collapse(2)
836  do j = 1, ja
837  do i = 1, ia
838  do k = 1, ka
839  damp_t_dens(k,i,j) = 0.0_rp
840  damp_t_momz(k,i,j) = 0.0_rp
841  damp_t_momx(k,i,j) = 0.0_rp
842  damp_t_momy(k,i,j) = 0.0_rp
843  damp_t_rhot(k,i,j) = 0.0_rp
844  end do
845  end do
846  end do
847 
848  call atmos_dyn_fill_halo( dens_damp, 0.0_rp, .true., .true. )
849  call atmos_dyn_fill_halo( dens_t, 0.0_rp, .false., .true. )
850 
851  do step = 1, nstep
852 
853  !-----< prepare tendency >-----
854 
855  call prof_rapstart("DYN_Large_Tendency", 2)
856 
857  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
858 !OCL XFILL
859  do j = js-1, je+1
860  do i = max(is-1,1), min(ie+1,ia)
861  do k = ks, ke
862  diff(k,i,j) = dens(k,i,j) - damp_dens(k,i,j)
863  enddo
864  enddo
865  enddo
866 
867  call file_history_query( hist_damp(1), do_put )
868  if ( twod ) then
869  !$omp parallel do default(none) OMP_SCHEDULE_ &
870  !$omp private(j,k) &
871  !$omp shared(IS,JS,JE,KS,KE) &
872  !$omp shared(DAMP_alpha_DENS,diff,DENS_tq,DENS_t,DENS_tp,BND_SMOOTHER_FACT,EPS) &
873  !$omp shared(damp_t_DENS,DENS_damp,do_put,nstep)
874 !OCL XFILL
875  do j = js, je
876  do k = ks, ke
877  dens_damp(k,is,j) = - damp_alpha_dens(k,is,j) &
878  * ( diff(k,is,j) & ! rayleigh damping
879  - ( diff(k,is,j-1) + diff(k,is,j+1) - diff(k,is,j)*2.0_rp ) &
880  * 0.25_rp * bnd_smoother_fact ) & ! horizontal smoother
881  + dens_tq(k,is,j) * ( 0.5_rp - sign( 0.5_rp, damp_alpha_dens(k,is,j)-eps ) ) ! dencity change due to rayleigh damping for tracers
882  dens_t(k,is,j) = dens_tp(k,is,j) & ! tendency from physical step
883  + dens_damp(k,is,j)
884  if ( do_put ) damp_t_dens(k,is,j) = damp_t_dens(k,is,j) + dens_damp(k,is,j) / nstep
885  enddo
886  enddo
887  else
888  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
889  !$omp private(i,j,k) &
890  !$omp shared(JS,JE,IS,IE,KS,KE) &
891  !$omp shared(DAMP_alpha_DENS,diff,DENS_tq,DENS_t,DENS_tp,BND_SMOOTHER_FACT,EPS) &
892  !$omp shared(damp_t_DENS,DENS_damp,do_put,nstep)
893 !OCL XFILL
894  do j = js, je
895  do i = is, ie
896  do k = ks, ke
897  dens_damp(k,i,j) = - damp_alpha_dens(k,i,j) &
898  * ( diff(k,i,j) & ! rayleigh damping
899  - ( diff(k,i-1,j) + diff(k,i+1,j) + diff(k,i,j-1) + diff(k,i,j+1) - diff(k,i,j)*4.0_rp ) &
900  * 0.125_rp * bnd_smoother_fact ) & ! horizontal smoother
901  + dens_tq(k,i,j) * ( 0.5_rp - sign( 0.5_rp, damp_alpha_dens(k,i,j)-eps ) ) ! dencity change due to rayleigh damping for tracers
902  dens_t(k,i,j) = dens_tp(k,i,j) & ! tendency from physical step
903  + dens_damp(k,i,j)
904  if ( do_put ) damp_t_dens(k,i,j) = damp_t_dens(k,i,j) + dens_damp(k,i,j) / nstep
905  enddo
906  enddo
907  enddo
908  end if
909 
910  call comm_vars8( dens_damp(:,:,:), i_comm_dens_damp )
911  call comm_vars8( dens_t(:,:,:), i_comm_dens_t )
912 
913  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
914 !OCL XFILL
915  do j = js-1, je+1
916  do i = max(is-1,1), min(ie+1,ia)
917  do k = ks, ke-1
918  diff(k,i,j) = momz(k,i,j) - damp_velz(k,i,j) * ( dens(k,i,j)+dens(k+1,i,j) ) * 0.5_rp
919  enddo
920  enddo
921  enddo
922 
923  call file_history_query( hist_damp(2), do_put )
924  if ( twod ) then
925  !$omp parallel do default(none) OMP_SCHEDULE_ &
926  !$omp private(j,k,damp) &
927  !$omp shared(IS,JS,JE,KS,KE) &
928  !$omp shared(DENS_damp,DENS,MOMZ) &
929  !$omp shared(DAMP_alpha_VELZ,diff,BND_SMOOTHER_FACT,MOMZ_t,MOMZ_tp) &
930  !$omp shared(damp_t_MOMZ,do_put,nstep)
931 !OCL XFILL
932  do j = js, je
933  do k = ks, ke-1
934  damp = - damp_alpha_velz(k,is,j) &
935  * ( diff(k,is,j) & ! rayleigh damping
936  - ( diff(k,is,j-1) + diff(k,is,j+1) - diff(k,is,j)*2.0_rp ) &
937  * 0.25_rp * bnd_smoother_fact ) ! horizontal smoother
938 
939  momz_t(k,is,j) = momz_tp(k,is,j) & ! tendency from physical step
940  + damp &
941  + ( dens_damp(k,is,j) + dens_damp(k+1,is,j) ) * momz(k,is,j) / ( dens(k,is,j) + dens(k,is,j) )
942  if ( do_put ) damp_t_momz(k,is,j) = damp_t_momz(k,is,j) + damp / nstep
943  enddo
944  enddo
945  else
946  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
947  !$omp private(i,j,k,damp) &
948  !$omp shared(JS,JE,IS,IE,KS,KE) &
949  !$omp shared(DENS_damp,DENS,MOMZ) &
950  !$omp shared(DAMP_alpha_VELZ,diff,BND_SMOOTHER_FACT,MOMZ_t,MOMZ_tp) &
951  !$omp shared(damp_t_MOMZ,do_put,nstep)
952 !OCL XFILL
953  do j = js, je
954  do i = is, ie
955  do k = ks, ke-1
956  damp = - damp_alpha_velz(k,i,j) &
957  * ( diff(k,i,j) & ! rayleigh damping
958  - ( diff(k,i-1,j) + diff(k,i+1,j) + diff(k,i,j-1) + diff(k,i,j+1) - diff(k,i,j)*4.0_rp ) &
959  * 0.125_rp * bnd_smoother_fact ) ! horizontal smoother
960 
961  momz_t(k,i,j) = momz_tp(k,i,j) & ! tendency from physical step
962  + damp &
963  + ( dens_damp(k,i,j) + dens_damp(k+1,i,j) ) * momz(k,i,j) / ( dens(k,i,j) + dens(k,i,j) )
964  if ( do_put ) damp_t_momz(k,i,j) = damp_t_momz(k,i,j) + damp / nstep
965  enddo
966  enddo
967  enddo
968  end if
969 !OCL XFILL
970  do j = js, je
971  do i = is, ie
972  momz_t( 1:ks-1,i,j) = 0.0_rp
973  momz_t(ke:ka ,i,j) = 0.0_rp
974  enddo
975  enddo
976  call comm_vars8( momz_t(:,:,:), i_comm_momz_t )
977 
978  call comm_wait( dens_damp(:,:,:), i_comm_dens_damp )
979 
980  if ( twod ) then
981  !$omp parallel do private(j,k) OMP_SCHEDULE_
982 !OCL XFILL
983  do j = js-1, je+1
984  do k = ks, ke
985  diff(k,is,j) = momx(k,is,j) - damp_velx(k,is,j) * dens(k,is,j)
986  enddo
987  enddo
988  else
989  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
990 !OCL XFILL
991  do j = js-1, je+1
992  do i = is-1, ie+1
993  do k = ks, ke
994  diff(k,i,j) = momx(k,i,j) - damp_velx(k,i,j) * ( dens(k,i,j)+dens(k,i+1,j) ) * 0.5_rp
995  enddo
996  enddo
997  enddo
998  end if
999 
1000 
1001  call file_history_query( hist_damp(3), do_put )
1002  if ( twod ) then
1003 !OCL XFILL
1004  !$omp parallel do default(none) OMP_SCHEDULE_ &
1005  !$omp private(j,k,damp) &
1006  !$omp shared(IS,JS,JE,KS,KE) &
1007  !$omp shared(DENS_damp,DENS,MOMX) &
1008  !$omp shared(DAMP_alpha_VELX,diff,BND_SMOOTHER_FACT,MOMX_tp,MOMX_t) &
1009  !$omp shared(damp_t_MOMX,do_put,nstep)
1010  do j = js, je
1011  do k = ks, ke
1012  damp = - damp_alpha_velx(k,is,j) &
1013  * ( diff(k,is,j) & ! rayleigh damping
1014  - ( diff(k,is,j-1) + diff(k,is,j+1) - diff(k,is,j)*2.0_rp ) &
1015  * 0.25_rp * bnd_smoother_fact ) ! horizontal smoother
1016  momx_t(k,is,j) = momx_tp(k,is,j) & ! tendency from physical step
1017  + damp &
1018  + dens_damp(k,is,j) * momx(k,is,j) / dens(k,is,j)
1019  if ( do_put ) damp_t_momx(k,is,j) = damp_t_momx(k,is,j) + damp / nstep
1020  enddo
1021  enddo
1022  else
1023  if( spnudge_uv ) then
1024  if( spnudge_uv_divfree ) then
1025 
1026  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1027 !OCL XFILL
1028  do j = js-1, je+1
1029  do i = is-1, ie+1
1030  do k = ks, ke
1031  diff3(k,i,j) = momy(k,i,j) - damp_vely(k,i,j) * ( dens(k,i,j)+dens(k,i,j+1) ) * 0.5_rp
1032  enddo
1033  enddo
1034  enddo
1035 
1036  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1037 !OCL XFILL
1038  do j = js, je
1039  do i = is, ie
1040  do k = ks, ke
1041  diff2(k,i,j) = diff(k,i,j) / ( ( dens(k,i,j)+dens(k,i+1,j) ) * 0.5_rp )
1042  enddo
1043  enddo
1044  enddo
1045 
1046  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1047 !OCL XFILL
1048  do j = js, je
1049  do i = is, ie
1050  do k = ks, ke
1051  diff3(k,i,j) = diff3(k,i,j) / ( ( dens(k,i,j)+dens(k,i,j+1) ) * 0.5_rp )
1052  enddo
1053  enddo
1054  enddo
1055 
1056  call dft_g2g_divfree(ka,ks,ke,ia,is,ie,ja,js,je,spnudge_uv_lm,spnudge_uv_mm,diff2,diff3)
1057 
1058  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1059 !OCL XFILL
1060  do j = js, je
1061  do i = is, ie
1062  do k = ks, ke
1063  diff2(k,i,j) = diff2(k,i,j) * ( dens(k,i,j)+dens(k,i+1,j) ) * 0.5_rp
1064  enddo
1065  enddo
1066  enddo
1067 
1068  else
1069 
1070  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1071 !OCL XFILL
1072  do j = js, je
1073  do i = is, ie
1074  do k = ks, ke
1075  diff2(k,i,j) = diff(k,i,j) / ( ( dens(k,i,j)+dens(k,i+1,j) ) * 0.5_rp )
1076  enddo
1077  enddo
1078  enddo
1079 
1080  call dft_g2g(ka,ks,ke,ia,is,ie,ja,js,je,spnudge_uv_lm,spnudge_uv_mm,diff2)
1081 
1082  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1083  do j = js, je
1084  do i = is, ie
1085  do k = ks, ke
1086  diff2(k,i,j) = diff2(k,i,j) * ( dens(k,i,j)+dens(k,i+1,j) ) * 0.5_rp
1087  enddo
1088  enddo
1089  enddo
1090 
1091  endif
1092  else
1093 
1094  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1095 !OCL XFILL
1096  do j = js, je
1097  do i = is, ie
1098  do k = ks, ke
1099  diff2(k,i,j) = 0
1100  enddo
1101  enddo
1102  enddo
1103 
1104  endif
1105 
1106 !OCL XFILL
1107  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
1108  !$omp private(i,j,k,damp) &
1109  !$omp shared(JS,JE,IS,IE,KS,KE) &
1110  !$omp shared(DENS_damp,DENS,MOMX) &
1111  !$omp shared(DAMP_alpha_VELX,diff,diff2,BND_SMOOTHER_FACT,SPNUDGE_u_alpha,MOMX_tp,MOMX_t) &
1112  !$omp shared(damp_t_MOMX,do_put,nstep)
1113  do j = js, je
1114  do i = is, ie
1115  do k = ks, ke
1116  damp = - damp_alpha_velx(k,i,j) &
1117  * ( diff(k,i,j) & ! rayleigh damping
1118  - ( diff(k,i-1,j) + diff(k,i+1,j) + diff(k,i,j-1) + diff(k,i,j+1) - diff(k,i,j)*4.0_rp ) &
1119  * 0.125_rp * bnd_smoother_fact ) & ! horizontal smoother
1120  - spnudge_u_alpha(k,i,j) * diff2(k,i,j)
1121  momx_t(k,i,j) = momx_tp(k,i,j) & ! tendency from physical step
1122  + damp &
1123  + ( dens_damp(k,i,j) + dens_damp(k,i+1,j) ) * momx(k,i,j) / ( dens(k,i,j) + dens(k,i+1,j) )
1124  if ( do_put ) damp_t_momx(k,i,j) = damp_t_momx(k,i,j) + damp / nstep
1125  enddo
1126  enddo
1127  enddo
1128  end if
1129 
1130  call atmos_dyn_fill_halo( momx_t, 0.0_rp, .false., .true. )
1131  call comm_vars8( momx_t(:,:,:), i_comm_momx_t )
1132 
1133  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1134 !OCL XFILL
1135  do j = js-1, je+1
1136  do i = max(is-1,1), min(ie+1,ia)
1137  do k = ks, ke
1138  diff(k,i,j) = momy(k,i,j) - damp_vely(k,i,j) * ( dens(k,i,j)+dens(k,i,j+1) ) * 0.5_rp
1139  enddo
1140  enddo
1141  enddo
1142 
1143  call file_history_query( hist_damp(4), do_put )
1144  if ( twod ) then
1145 !OCL XFILL
1146  !$omp parallel do default(none) OMP_SCHEDULE_ &
1147  !$omp private(j,k,damp) &
1148  !$omp shared(IS,JS,JE,KS,KE) &
1149  !$omp shared(DENS_damp,DENS,MOMY) &
1150  !$omp shared(DAMP_alpha_VELY,diff,BND_SMOOTHER_FACT,MOMY_tp,MOMY_t) &
1151  !$omp shared(damp_t_MOMY,do_put,nstep)
1152  do j = js, je
1153  do k = ks, ke
1154  damp = - damp_alpha_vely(k,is,j) &
1155  * ( diff(k,is,j) & ! rayleigh damping
1156  - ( diff(k,is,j-1) + diff(k,is,j+1) - diff(k,is,j)*2.0_rp ) &
1157  * 0.25_rp * bnd_smoother_fact ) ! horizontal smoother
1158  momy_t(k,is,j) = momy_tp(k,is,j) & ! tendency from physical step
1159  + damp &
1160  + ( dens_damp(k,is,j) + dens_damp(k,is,j+1) ) * momy(k,is,j) / ( dens(k,is,j) + dens(k,is,j+1) )
1161  if ( do_put ) damp_t_momy(k,is,j) = damp_t_momy(k,is,j) + damp / nstep
1162  enddo
1163  enddo
1164  else
1165  if( spnudge_uv ) then
1166  if( spnudge_uv_divfree ) then
1167 
1168  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1169 !OCL XFILL
1170  do j = js, je
1171  do i = is, ie
1172  do k = ks, ke
1173  diff2(k,i,j) = diff3(k,i,j) * ( dens(k,i,j)+dens(k,i,j+1) ) * 0.5_rp
1174  enddo
1175  enddo
1176  enddo
1177 
1178  else
1179  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1180 !OCL XFILL
1181  do j = js, je
1182  do i = is, ie
1183  do k = ks, ke
1184  diff2(k,i,j) = diff(k,i,j) / ( ( dens(k,i,j)+dens(k,i,j+1) ) * 0.5_rp )
1185  enddo
1186  enddo
1187  enddo
1188 
1189  call dft_g2g(ka,ks,ke,ia,is,ie,ja,js,je,spnudge_uv_lm,spnudge_uv_mm,diff2)
1190 
1191  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1192 !OCL XFILL
1193  do j = js, je
1194  do i = is, ie
1195  do k = ks, ke
1196  diff2(k,i,j) = diff2(k,i,j) * ( dens(k,i,j)+dens(k,i,j+1) ) * 0.5_rp
1197  enddo
1198  enddo
1199  enddo
1200  endif
1201  endif
1202 
1203 !OCL XFILL
1204  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
1205  !$omp private(i,j,k,damp) &
1206  !$omp shared(JS,JE,IS,IE,KS,KE) &
1207  !$omp shared(DENS_damp,DENS,MOMY) &
1208  !$omp shared(DAMP_alpha_VELY,diff,diff2,BND_SMOOTHER_FACT,SPNUDGE_v_alpha,MOMY_tp,MOMY_t) &
1209  !$omp shared(damp_t_MOMY,do_put,nstep)
1210  do j = js, je
1211  do i = is, ie
1212  do k = ks, ke
1213  damp = - damp_alpha_vely(k,i,j) &
1214  * ( diff(k,i,j) & ! rayleigh damping
1215  - ( diff(k,i-1,j) + diff(k,i+1,j) + diff(k,i,j-1) + diff(k,i,j+1) - diff(k,i,j)*4.0_rp ) &
1216  * 0.125_rp * bnd_smoother_fact ) & ! horizontal smoother
1217  - spnudge_v_alpha(k,i,j) * diff2(k,i,j)
1218  momy_t(k,i,j) = momy_tp(k,i,j) & ! tendency from physical step
1219  + damp &
1220  + ( dens_damp(k,i,j) + dens_damp(k,i,j+1) ) * momy(k,i,j) / ( dens(k,i,j) + dens(k,i,j+1) )
1221  if ( do_put ) damp_t_momy(k,i,j) = damp_t_momy(k,i,j) + damp / nstep
1222  enddo
1223  enddo
1224  enddo
1225  end if
1226 !OCL XFILL
1227 
1228  call atmos_dyn_fill_halo( momy_t, 0.0_rp, .false., .true. )
1229  call comm_vars8( momy_t(:,:,:), i_comm_momy_t )
1230 
1231  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1232 !OCL XFILL
1233  do j = js-1, je+2
1234  do i = max(is-1,1), min(ie+2,ia)
1235  do k = ks, ke
1236  diff(k,i,j) = rhot(k,i,j) - damp_pott(k,i,j) * dens(k,i,j)
1237  enddo
1238  enddo
1239  enddo
1240 
1241  call file_history_query( hist_damp(5), do_put )
1242  if ( twod ) then
1243 !OCL XFILL
1244  !$omp parallel do default(none) OMP_SCHEDULE_ &
1245  !$omp private(j,k,damp) &
1246  !$omp shared(IS,JS,JE,KS,KE) &
1247  !$omp shared(DENS_damp,DENS,RHOT) &
1248  !$omp shared(DAMP_alpha_POTT,diff,BND_SMOOTHER_FACT,RHOT_t,RHOT_tp) &
1249  !$omp shared(damp_t_RHOT,do_put,nstep)
1250  do j = js, je
1251  do k = ks, ke
1252  damp = - damp_alpha_pott(k,is,j) &
1253  * ( diff(k,is,j) & ! rayleigh damping
1254  - ( diff(k,is,j-1) + diff(k,is,j+1) - diff(k,is,j)*2.0_rp ) &
1255  * 0.25_rp * bnd_smoother_fact ) ! horizontal smoother
1256  rhot_t(k,is,j) = rhot_tp(k,is,j) & ! tendency from physical step
1257  + damp &
1258  + dens_damp(k,is,j) * rhot(k,is,j) / dens(k,is,j)
1259  if ( do_put ) damp_t_rhot(k,is,j) = damp_t_rhot(k,is,j) + damp / nstep
1260  enddo
1261  enddo
1262  else
1263  if( spnudge_pt ) then
1264 
1265  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1266 !OCL XFILL
1267  do j = js, je
1268  do i = is, ie
1269  do k = ks, ke
1270  diff2(k,i,j) = diff(k,i,j) / dens(k,i,j)
1271  enddo
1272  enddo
1273  enddo
1274 
1275  call dft_g2g(ka,ks,ke,ia,is,ie,ja,js,je,spnudge_pt_lm,spnudge_pt_mm,diff2)
1276 
1277  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1278 !OCL XFILL
1279  do j = js, je
1280  do i = is, ie
1281  do k = ks, ke
1282  diff2(k,i,j) = diff2(k,i,j) * dens(k,i,j)
1283  enddo
1284  enddo
1285  enddo
1286 
1287  else
1288 
1289  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
1290 !OCL XFILL
1291  do j = js, je
1292  do i = is, ie
1293  do k = ks, ke
1294  diff2(k,i,j) = 0
1295  enddo
1296  enddo
1297  enddo
1298 
1299  endif
1300 
1301 !OCL XFILL
1302  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
1303  !$omp private(i,j,k,damp) &
1304  !$omp shared(JS,JE,IS,IE,KS,KE) &
1305  !$omp shared(DENS_damp,DENS,RHOT) &
1306  !$omp shared(DAMP_alpha_POTT,diff,diff2,BND_SMOOTHER_FACT,SPNUDGE_pt_alpha,RHOT_t,RHOT_tp) &
1307  !$omp shared(damp_t_RHOT,do_put,nstep)
1308  do j = js, je
1309  do i = is, ie
1310  do k = ks, ke
1311  damp = - damp_alpha_pott(k,i,j) &
1312  * ( diff(k,i,j) & ! rayleigh damping
1313  - ( diff(k,i-1,j) + diff(k,i+1,j) + diff(k,i,j-1) + diff(k,i,j+1) - diff(k,i,j)*4.0_rp ) &
1314  * 0.125_rp * bnd_smoother_fact ) & ! horizontal smoother
1315  - spnudge_pt_alpha(k,i,j) * diff2(k,i,j)
1316  rhot_t(k,i,j) = rhot_tp(k,i,j) & ! tendency from physical step
1317  + damp &
1318  + dens_damp(k,i,j) * rhot(k,i,j) / dens(k,i,j)
1319  if ( do_put ) damp_t_rhot(k,i,j) = damp_t_rhot(k,i,j) + damp / nstep
1320  enddo
1321  enddo
1322  enddo
1323  end if
1324 !OCL XFILL
1325 
1326  call atmos_dyn_fill_halo( rhot_t, 0.0_rp, .false., .true. )
1327  call comm_vars8( rhot_t(:,:,:), i_comm_rhot_t )
1328 
1329  call comm_wait ( dens_t(:,:,:), i_comm_dens_t, .false. )
1330  call comm_wait ( momz_t(:,:,:), i_comm_momz_t, .false. )
1331  call comm_wait ( momx_t(:,:,:), i_comm_momx_t, .false. )
1332  call comm_wait ( momy_t(:,:,:), i_comm_momy_t, .false. )
1333  call comm_wait ( rhot_t(:,:,:), i_comm_rhot_t, .false. )
1334 
1335  call prof_rapend ("DYN_Large_Tendency", 2)
1336 
1337  call prof_rapstart("DYN_Large_Numfilter", 2)
1338 
1339  !-----< prepare the fluxes of explicit numerical diffusion >-----
1340  if ( nd_coef == 0.0_rp .or. eval_type_numfilter == 'FILTER' ) then
1341 !OCL XFILL
1342  num_diff(:,:,:,:,:) = 0.0_rp
1343  else
1344  call atmos_dyn_fvm_numfilter_flux( num_diff(:,:,:,:,:), & ! [OUT]
1345  dens, momz, momx, momy, rhot, & ! [IN]
1346  cdz, cdx, cdy, fdz, fdx, fdy, twod, dts, & ! [IN]
1347  ref_dens, ref_pott, & ! [IN]
1348  nd_coef, nd_laplacian_num, nd_sfc_fact, nd_use_rs ) ! [IN]
1349  endif
1350 
1351  call prof_rapend ("DYN_Large_Numfilter", 2)
1352 
1353  !-----< calculate the divegence term >-----
1354 
1355  if ( divdmp_coef > 0.0_rp ) then
1356 
1357  call atmos_dyn_divergence( ddiv, & ! (out)
1358  momz, momx, momy, & ! (in)
1359  gsqrt, j13g, j23g, j33g, mapf, & ! (in)
1360  twod, & ! (in)
1361  rcdz, rcdx, rcdy, rfdz, fdz ) ! (in)
1362 
1363  else
1364 
1365 !XFILL
1366  ddiv = 0.0_rp
1367 
1368  end if
1369 
1370  !------------------------------------------------------------------------
1371  ! Start short time integration
1372  !------------------------------------------------------------------------
1373 
1374  call file_history_set_disable( .not. ( llast .and. ( step==nstep ) ) )
1375 
1376  call prof_rapstart("DYN_Short_Tinteg", 2)
1377 
1378  call atmos_dyn_tinteg_short( dens, momz, momx, momy, rhot, prog, & ! (inout)
1379  mflx, tflx, & ! (inout, out)
1380  dens_t, momz_t, momx_t, momy_t, rhot_t, & ! (in)
1381  dpres0, rt2p, corioli, & ! (in)
1382  num_diff, wdamp_coef, divdmp_coef, ddiv, & ! (in)
1383  flag_fct_momentum, flag_fct_t, & ! (in)
1384  flag_fct_along_stream, & ! (in)
1385  cdz, fdz, fdx, fdy, & ! (in)
1386  rcdz, rcdx, rcdy, rfdz, rfdx, rfdy, & ! (in)
1387  phi, gsqrt, j13g, j23g, j33g, mapf, & ! (in)
1388  ref_dens, ref_rhot, & ! (in)
1389  bnd_w, bnd_e, bnd_s, bnd_n, twod, & ! (in)
1390  dts ) ! (in)
1391 
1392  call prof_rapend ("DYN_Short_Tinteg", 2)
1393 
1394 
1395  call file_history_set_disable( .false. )
1396 
1397  if ( nd_coef > 0.0_rp .and. eval_type_numfilter == 'FILTER' ) then
1398  call comm_vars8( dens(:,:,:), i_comm_dens )
1399  call comm_vars8( momz(:,:,:), i_comm_momz )
1400  call comm_vars8( momx(:,:,:), i_comm_momx )
1401  call comm_vars8( momy(:,:,:), i_comm_momy )
1402  call comm_vars8( rhot(:,:,:), i_comm_rhot )
1403  call comm_wait ( dens(:,:,:), i_comm_dens, .false. )
1404  call comm_wait ( momz(:,:,:), i_comm_momz, .false. )
1405  call comm_wait ( momx(:,:,:), i_comm_momx, .false. )
1406  call comm_wait ( momy(:,:,:), i_comm_momy, .false. )
1407  call comm_wait ( rhot(:,:,:), i_comm_rhot, .false. )
1408 
1409  call atmos_dyn_fvm_apply_numfilter( &
1410  num_diff, & ! (out)
1411  dens, momz, momx, momy, rhot, & ! (inout)
1412  cdz, cdx, cdy, fdz, fdx, fdy, & ! (in)
1413  rcdz, rcdx, rcdy, rfdz, rfdx, rfdy, & ! (in)
1414  twod, dts, & ! (in)
1415  gsqrt, mapf, ref_dens, ref_pott, & ! (in)
1416  nd_coef, nd_laplacian_num, nd_sfc_fact, nd_use_rs ) ! (in)
1417  endif
1418 
1419  !$omp parallel do default(none) private(i,j,iv) OMP_SCHEDULE_ collapse(2) &
1420  !$omp shared(JSB,JEB,ISB,IEB,KS,KA,DENS,MOMZ,MOMX,MOMY,RHOT,VA,PROG,KE)
1421  do j = jsb, jeb
1422  do i = isb, ieb
1423  dens( 1:ks-1,i,j) = dens(ks,i,j)
1424  momz( 1:ks-1,i,j) = 0.0_rp
1425  momx( 1:ks-1,i,j) = momx(ks,i,j)
1426  momy( 1:ks-1,i,j) = momy(ks,i,j)
1427  rhot( 1:ks-1,i,j) = rhot(ks,i,j)
1428  do iv = 1, va
1429  prog( 1:ks-1,i,j,iv) = prog(ks,i,j,iv)
1430  end do
1431  dens(ke+1:ka, i,j) = dens(ke,i,j)
1432  momz(ke+1:ka, i,j) = 0.0_rp
1433  momx(ke+1:ka, i,j) = momx(ke,i,j)
1434  momy(ke+1:ka, i,j) = momy(ke,i,j)
1435  rhot(ke+1:ka, i,j) = rhot(ke,i,j)
1436  do iv = 1, va
1437  prog(ke+1:ka, i,j,iv) = prog(ke,i,j,iv)
1438  end do
1439  enddo
1440  enddo
1441 
1442  call comm_vars8( dens(:,:,:), i_comm_dens )
1443  call comm_vars8( momz(:,:,:), i_comm_momz )
1444  call comm_vars8( momx(:,:,:), i_comm_momx )
1445  call comm_vars8( momy(:,:,:), i_comm_momy )
1446  call comm_vars8( rhot(:,:,:), i_comm_rhot )
1447  do iv = 1, va
1448  call comm_vars8( prog(:,:,:,iv), i_comm_prog(iv) )
1449  end do
1450  call comm_wait ( dens(:,:,:), i_comm_dens, .false. )
1451  call comm_wait ( momz(:,:,:), i_comm_momz, .false. )
1452  call comm_wait ( momx(:,:,:), i_comm_momx, .false. )
1453  call comm_wait ( momy(:,:,:), i_comm_momy, .false. )
1454  call comm_wait ( rhot(:,:,:), i_comm_rhot, .false. )
1455  do iv = 1, va
1456  call comm_wait ( prog(:,:,:,iv), i_comm_prog(iv), .false. )
1457  end do
1458 
1459  if ( use_average ) then
1460  !$omp parallel do
1461  do j = jsb, jeb
1462  do i = isb, ieb
1463  do k = ks, ke
1464  dens_av(k,i,j) = dens_av(k,i,j) + dens(k,i,j) / nstep
1465  momz_av(k,i,j) = momz_av(k,i,j) + momz(k,i,j) / nstep
1466  momx_av(k,i,j) = momx_av(k,i,j) + momx(k,i,j) / nstep
1467  momy_av(k,i,j) = momy_av(k,i,j) + momy(k,i,j) / nstep
1468  rhot_av(k,i,j) = rhot_av(k,i,j) + rhot(k,i,j) / nstep
1469  end do
1470  end do
1471  end do
1472  endif
1473 
1474  !$omp parallel
1475  do n = 1, 3
1476  !$omp do
1477  do j = jsb, jeb
1478  do i = isb, ieb
1479  do k = ks, ke
1480  mflx_av(k,i,j,n) = mflx_av(k,i,j,n) + mflx(k,i,j,n)
1481  end do
1482  end do
1483  end do
1484  !$omp end do nowait
1485  end do
1486  !$omp end parallel
1487 
1488  enddo ! dynamical steps
1489 
1490 
1491  !###########################################################################
1492  ! Update Tracers
1493  !###########################################################################
1494 
1495  !$omp parallel
1496  do n = 1, 3
1497 !OCL XFILL
1498  !$omp do
1499  do j = jsb, jeb
1500  do i = isb, ieb
1501  do k = ks, ke
1502  mflx(k,i,j,n) = mflx_av(k,i,j,n) / nstep
1503  end do
1504  end do
1505  end do
1506  !$omp end do nowait
1507  end do
1508  !$omp end parallel
1509 
1510  call comm_vars8( mflx(:,:,:,zdir), i_comm_mflx_z )
1511  if ( .not. twod ) call comm_vars8( mflx(:,:,:,xdir), i_comm_mflx_x )
1512  call comm_vars8( mflx(:,:,:,ydir), i_comm_mflx_y )
1513  call comm_wait ( mflx(:,:,:,zdir), i_comm_mflx_z, .false. )
1514  if ( .not. twod ) call comm_wait ( mflx(:,:,:,xdir), i_comm_mflx_x, .false. )
1515  call comm_wait ( mflx(:,:,:,ydir), i_comm_mflx_y, .false. )
1516 
1517 #ifndef DRY
1518 
1519  !------------------------------------------------------------------------
1520  ! Update each tracer
1521  !------------------------------------------------------------------------
1522 
1523  do iq = 1, qa
1524 
1525  if ( tracer_advc(iq) ) then
1526 
1527  call prof_rapstart("DYN_Large_Numfilter", 2)
1528 
1529  if ( nd_coef_q == 0.0_rp ) then
1530 !OCL XFILL
1531  num_diff_q(:,:,:,:) = 0.0_rp
1532  else
1533  call atmos_dyn_fvm_numfilter_flux_q( num_diff_q(:,:,:,:), & ! [OUT]
1534  dens00, qtrc(:,:,:,iq), iq==i_qv, & ! [IN]
1535  cdz, cdx, cdy, twod, dtl, & ! [IN]
1536  ref_qv, iq, & ! [IN]
1537  nd_coef_q, nd_laplacian_num, nd_sfc_fact, nd_use_rs ) ! [IN]
1538  endif
1539 
1540  call prof_rapend ("DYN_Large_Numfilter", 2)
1541 
1542  call prof_rapstart("DYN_Tracer_Tinteg", 2)
1543 
1544  if ( flag_tracer_split_tend ) then
1545  !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
1546  !$omp shared(KA,IA,JA,iq,QTRC,RHOQ_t,DENS,dtl)
1547  do j = 1, ja
1548  do i = 1, ia
1549  do k = 1, ka
1550  qtrc(k,i,j,iq) = qtrc(k,i,j,iq) &
1551  + rhoq_t(k,i,j,iq) * dtl / dens(k,i,j)
1552  end do
1553  end do
1554  end do
1555  rhoq_tn => zero
1556  else
1557  rhoq_tn => rhoq_t(:,:,:,iq)
1558  end if
1559 
1560  call atmos_dyn_tinteg_tracer( &
1561  qtrc(:,:,:,iq), & ! (inout)
1562  qflx(:,:,:,:), & ! (out)
1563  qtrc0(:,:,:,iq), rhoq_tn, & ! (in)
1564  dens00, dens, & ! (in)
1565  mflx, num_diff_q, & ! (in)
1566  gsqrt, mapf(:,:,:,i_xy), & ! (in)
1567  cdz, rcdz, rcdx, rcdy, & ! (in)
1568  bnd_w, bnd_e, bnd_s, bnd_n, & ! (in)
1569  twod, & ! (in)
1570  dtl, & ! (in)
1571  llast .AND. flag_fct_tracer, & ! (in)
1572  flag_fct_along_stream ) ! (in)
1573 
1574  call prof_rapend ("DYN_Tracer_Tinteg", 2)
1575 
1576  if ( llast ) then
1577  do iv = 1, 3
1578  call file_history_query( hist_qflx(iv,iq), do_put )
1579  if ( do_put .or. monit_lateral_flag(iv) ) then
1580  call multiply_flux_by_metric_xyz( iv, qflx, gsqrt, mapf )
1581  call file_history_put( hist_qflx(iv,iq), qflx(:,:,:,iv) )
1582  end if
1583  end do
1584 
1585  if ( tracer_mass(iq) == 1.0_rp ) then
1586  if ( bnd_w .and. monit_qflx_west > 0 ) then
1587  !$omp parallel do
1588  do j = js, je
1589  do k = ks, ke
1590  qflx_west(k,j) = qflx_west(k,j) + qflx(k,is-1,j,xdir)
1591  end do
1592  end do
1593  end if
1594  if ( bnd_e .and. monit_qflx_east > 0 ) then
1595  !$omp parallel do
1596  do j = js, je
1597  do k = ks, ke
1598  qflx_east(k,j) = qflx_east(k,j) + qflx(k,ie,j,xdir)
1599  end do
1600  end do
1601  end if
1602  if ( bnd_s .and. monit_qflx_south > 0 ) then
1603  !$omp parallel do
1604  do i = is, ie
1605  do k = ks, ke
1606  qflx_south(k,i) = qflx_south(k,i) + qflx(k,i,js-1,ydir)
1607  end do
1608  end do
1609  end if
1610  if ( bnd_n .and. monit_qflx_north > 0 ) then
1611  !$omp parallel do
1612  do i = is, ie
1613  do k = ks, ke
1614  qflx_north(k,i) = qflx_north(k,i) + qflx(k,i,je,ydir)
1615  end do
1616  end do
1617  end if
1618 
1619  end if
1620 
1621  end if
1622 
1623  else
1624 
1625  !$omp parallel do
1626  do j = js, je
1627  do i = is, ie
1628  do k = ks, ke
1629  qtrc(k,i,j,iq) = ( qtrc0(k,i,j,iq) * dens00(k,i,j) &
1630  + rhoq_t(k,i,j,iq) * dtl ) / dens(k,i,j)
1631  end do
1632  end do
1633  end do
1634 
1635  end if
1636 
1637  call comm_vars8( qtrc(:,:,:,iq), i_comm_qtrc(iq) )
1638 
1639  if ( use_average ) then
1640  do j = jsb, jeb
1641  do i = isb, ieb
1642  do k = ks, ke
1643  qtrc_av(k,i,j,iq) = qtrc(k,i,j,iq)
1644  end do
1645  end do
1646  end do
1647  endif
1648 
1649  enddo ! scalar quantities loop
1650 
1651  do iq = 1, qa
1652  call comm_wait ( qtrc(:,:,:,iq), i_comm_qtrc(iq), .false. )
1653  enddo
1654 #endif
1655 
1656  if ( llast ) then
1657  !- history ---------------------------------------
1658  call file_history_put( hist_phys(1), dens_tp )
1659  call file_history_put( hist_phys(2), momz_tp )
1660  call file_history_put( hist_phys(3), momx_tp )
1661  call file_history_put( hist_phys(4), momy_tp )
1662  call file_history_put( hist_phys(5), rhot_tp )
1663 
1664  call file_history_put( hist_damp(1), damp_t_dens )
1665  call file_history_put( hist_damp(2), damp_t_momz )
1666  call file_history_put( hist_damp(3), damp_t_momx )
1667  call file_history_put( hist_damp(4), damp_t_momy )
1668  call file_history_put( hist_damp(5), damp_t_rhot )
1669 
1670  do iv = 1, 3
1671  call file_history_query( hist_mflx(iv), do_put )
1672  if ( do_put .or. monit_lateral_flag(iv) ) then
1673  call multiply_flux_by_metric_xyz( iv, mflx, gsqrt, mapf )
1674  call file_history_put( hist_mflx(iv), mflx(:,:,:,iv) )
1675  end if
1676  end do
1677 
1678  do iv = 1, 3
1679  call file_history_query( hist_tflx(iv), do_put )
1680  if ( do_put ) then
1681  call multiply_flux_by_metric_xyz( iv, tflx, gsqrt, mapf )
1682  call file_history_put( hist_tflx(iv), tflx(:,:,:,iv) )
1683  end if
1684  end do
1685 
1686  !- monitor mass budget ------------------------------------
1687  call monitor_put( monit_damp_mass, damp_t_dens(:,:,:) )
1688  call monitor_put_lateral_flux( monit_mflx_west, bnd_w .and. monit_mflx_west > 0, mflx(:,is-1,:,xdir), zero_x )
1689  call monitor_put_lateral_flux( monit_mflx_east, bnd_e .and. monit_mflx_east > 0, mflx(:,ie,:,xdir), zero_x )
1690  call monitor_put_lateral_flux( monit_mflx_south, bnd_s .and. monit_mflx_south > 0, mflx(:,:,js-1,ydir), zero_y )
1691  call monitor_put_lateral_flux( monit_mflx_north, bnd_n .and. monit_mflx_north > 0, mflx(:,:,je,ydir), zero_y )
1692 
1693  call monitor_put( monit_damp_qtot, dens_tq(:,:,:) )
1694  call monitor_put_lateral_flux( monit_qflx_west, bnd_w, qflx_west, zero_x )
1695  call monitor_put_lateral_flux( monit_qflx_east, bnd_e, qflx_east, zero_x )
1696  call monitor_put_lateral_flux( monit_qflx_south, bnd_s, qflx_south, zero_y )
1697  call monitor_put_lateral_flux( monit_qflx_north, bnd_n, qflx_north, zero_y )
1698  end if
1699 
1700  return
1701  end subroutine atmos_dyn_tstep_large_fvm_heve
1702 
1703  !-- private subroutines --------------------------------------------
1704 
1705  subroutine monitor_put_lateral_flux( MONIT_ID, BND_flag, flx, flx_zero )
1706  use scale_monitor, only: &
1707  monitor_put
1708  implicit none
1709 
1710  integer, intent(in) :: MONIT_ID
1711  logical, intent(in) :: BND_flag
1712  real(RP), target, intent(in) :: flx(:,:)
1713  real(RP), target, intent(in) :: flx_zero(:,:)
1714 
1715  real(RP), pointer :: flx_ptr(:,:)
1716  !------------------------------------------
1717  if ( bnd_flag ) then
1718  flx_ptr => flx
1719  else
1720  flx_ptr => flx_zero
1721  end if
1722  call monitor_put( monit_id, flx_ptr(:,:) )
1723 
1724  return
1725  end subroutine monitor_put_lateral_flux
1726 
1727  subroutine multiply_flux_by_metric_xyz( I_DIR, flx, GSQRT, MAPF )
1728  implicit none
1729  integer, intent(in) :: I_DIR
1730  real(RP), intent(inout) :: flx(KA,IA,JA,3)
1731  real(RP), intent(in) :: GSQRT(KA,IA,JA,7)
1732  real(RP), intent(in) :: MAPF (IA,JA,2,4)
1733 
1734  integer :: i, j, k
1735  !---------------------------------------------------
1736 
1737  if (i_dir == zdir) then
1738  !$omp parallel do
1739  do j = js, je
1740  do i = is, ie
1741  do k = ks-1, ke
1742  flx(k,i,j,zdir) = flx(k,i,j,zdir) * mapf(i,j,1,i_xy) * mapf(i,j,2,i_xy) / gsqrt(k,i,j,i_xyw)
1743  end do
1744  end do
1745  end do
1746  end if
1747 
1748  if (i_dir == xdir) then
1749  !$omp parallel do
1750  do j = js, je
1751  do i = isb, ieb
1752  do k = ks, ke
1753  flx(k,i,j,xdir) = flx(k,i,j,xdir) * mapf(i,j,2,i_uy) / gsqrt(k,i,j,i_uyz)
1754  end do
1755  end do
1756  end do
1757  end if
1758 
1759  if (i_dir == ydir) then
1760  !$omp parallel do
1761  do j = jsb, jeb
1762  do i = is, ie
1763  do k = ks, ke
1764  flx(k,i,j,ydir) = flx(k,i,j,ydir) * mapf(i,j,1,i_xv) / gsqrt(k,i,j,i_xvz)
1765  end do
1766  end do
1767  end do
1768  end if
1769 
1770  return
1771  end subroutine multiply_flux_by_metric_xyz
1772 
1773 end module scale_atmos_dyn_tstep_large_fvm_heve
scale_atmos_grid_cartesc_index::isb
integer, public isb
Definition: scale_atmos_grid_cartesC_index.F90:63
scale_atmos_grid_cartesc_index::i_uy
integer, public i_uy
Definition: scale_atmos_grid_cartesC_index.F90:99
scale_atmos_grid_cartesc_index::ke
integer, public ke
end point of inner domain: z, local
Definition: scale_atmos_grid_cartesC_index.F90:52
scale_atmos_grid_cartesc_index::i_xv
integer, public i_xv
Definition: scale_atmos_grid_cartesC_index.F90:100
scale_prc::prc_abort
subroutine, public prc_abort
Abort Process.
Definition: scale_prc.F90:342
scale_atmos_dyn_tinteg_tracer::atmos_dyn_tinteg_tracer
procedure(tinteg), pointer, public atmos_dyn_tinteg_tracer
Definition: scale_atmos_dyn_tinteg_tracer.F90:74
scale_tracer::qa
integer, public qa
Definition: scale_tracer.F90:34
scale_spnudge::spnudge_uv_lm
integer, public spnudge_uv_lm
Definition: scale_spnudge.F90:26
scale_index
module Index
Definition: scale_index.F90:11
scale_spnudge::spnudge_pt_lm
integer, public spnudge_pt_lm
Definition: scale_spnudge.F90:30
scale_spnudge::spnudge_pt
logical, public spnudge_pt
Definition: scale_spnudge.F90:29
scale_atmos_dyn_fvm_numfilter::atmos_dyn_fvm_numfilter_flux
subroutine, public atmos_dyn_fvm_numfilter_flux(num_diff, DENS, MOMZ, MOMX, MOMY, RHOT, CDZ, CDX, CDY, FDZ, FDX, FDY, TwoD, DT, REF_dens, REF_pott, ND_COEF, ND_LAPLACIAN_NUM, ND_SFC_FACT, ND_USE_RS)
Calculate fluxes with numerical filter for prognostic variables of dynamical core.
Definition: scale_atmos_dyn_fvm_numfilter.F90:390
scale_tracer::tracer_mass
real(rp), dimension(qa_max), public tracer_mass
Definition: scale_tracer.F90:46
scale_const::const_undef2
integer, parameter, public const_undef2
undefined value (INT2)
Definition: scale_const.F90:38
scale_prc_cartesc::prc_has_s
logical, public prc_has_s
Definition: scale_prc_cartesC.F90:50
scale_atmos_dyn_common
module Atmosphere / Dynamics common
Definition: scale_atmos_dyn_common.F90:12
scale_spnudge::spnudge_pt_alpha
real(rp), dimension(:,:,:), allocatable, public spnudge_pt_alpha
Definition: scale_spnudge.F90:39
scale_precision
module PRECISION
Definition: scale_precision.F90:14
scale_atmos_grid_cartesc_index::ka
integer, public ka
Definition: scale_atmos_grid_cartesC_index.F90:47
scale_tracer::tracer_advc
logical, dimension(qa_max), public tracer_advc
Definition: scale_tracer.F90:45
scale_dft::dft_g2g_divfree
subroutine, public dft_g2g_divfree(KA, KS, KE, IA, IS, IE, JA, JS, JE, LM, MM, u, v)
Definition: scale_dft.F90:227
scale_spnudge::spnudge_uv_divfree
logical, public spnudge_uv_divfree
Definition: scale_spnudge.F90:25
scale_file_history::file_history_set_disable
subroutine, public file_history_set_disable(switch)
set switch to turn on/off history
Definition: scale_file_history.F90:1712
scale_const::const_eps
real(rp), public const_eps
small number
Definition: scale_const.F90:33
scale_prof::prof_rapstart
subroutine, public prof_rapstart(rapname_base, level, disable_barrier)
Start raptime.
Definition: scale_prof.F90:159
scale_file_history
module file_history
Definition: scale_file_history.F90:15
scale_spnudge::spnudge_qv_lm
integer, public spnudge_qv_lm
Definition: scale_spnudge.F90:34
scale_prc_cartesc::prc_has_n
logical, public prc_has_n
Definition: scale_prc_cartesC.F90:48
scale_spnudge::spnudge_qv
logical, public spnudge_qv
Definition: scale_spnudge.F90:33
scale_atmos_grid_cartesc_index::jeb
integer, public jeb
Definition: scale_atmos_grid_cartesC_index.F90:66
scale_atmos_dyn_tinteg_short
module Atmosphere / Dynamics Temporal integration
Definition: scale_atmos_dyn_tinteg_short.F90:11
scale_prc_cartesc::prc_has_e
logical, public prc_has_e
Definition: scale_prc_cartesC.F90:49
scale_index::va
integer, public va
Definition: scale_index.F90:35
scale_dft
Definition: scale_dft.F90:3
scale_prc
module PROCESS
Definition: scale_prc.F90:11
scale_spnudge::spnudge_pt_mm
integer, public spnudge_pt_mm
Definition: scale_spnudge.F90:31
scale_atmos_grid_cartesc_index::i_xy
integer, public i_xy
Definition: scale_atmos_grid_cartesC_index.F90:98
scale_atmos_grid_cartesc_index::ie
integer, public ie
end point of inner domain: x, local
Definition: scale_atmos_grid_cartesC_index.F90:54
scale_atmos_grid_cartesc_index::i_uyz
integer, public i_uyz
Definition: scale_atmos_grid_cartesC_index.F90:94
scale_io
module STDIO
Definition: scale_io.F90:10
scale_tracer::k
real(rp), public k
Definition: scale_tracer.F90:44
scale_atmos_grid_cartesc_index
module atmosphere / grid / cartesC index
Definition: scale_atmos_grid_cartesC_index.F90:12
scale_const
module CONSTANT
Definition: scale_const.F90:11
scale_atmos_grid_cartesc_index::ia
integer, public ia
Definition: scale_atmos_grid_cartesC_index.F90:48
scale_comm_cartesc::comm_vars8_init
subroutine, public comm_vars8_init(varname, var, vid)
Register variables.
Definition: scale_comm_cartesC.F90:294
scale_spnudge::spnudge_u_alpha
real(rp), dimension(:,:,:), allocatable, public spnudge_u_alpha
Definition: scale_spnudge.F90:37
scale_debug::check
subroutine, public check(current_line, v)
Undefined value checker.
Definition: scale_debug.F90:56
scale_atmos_grid_cartesc_index::zdir
integer, parameter, public zdir
Definition: scale_atmos_grid_cartesC_index.F90:32
scale_atmos_grid_cartesc_index::i_xvz
integer, public i_xvz
Definition: scale_atmos_grid_cartesC_index.F90:95
scale_atmos_dyn_tstep_large_fvm_heve::atmos_dyn_tstep_large_fvm_heve
subroutine, public atmos_dyn_tstep_large_fvm_heve(DENS, MOMZ, MOMX, MOMY, RHOT, QTRC, PROG, DENS_av, MOMZ_av, MOMX_av, MOMY_av, RHOT_av, QTRC_av, num_diff, num_diff_q, QTRC0, DENS_tp, MOMZ_tp, MOMX_tp, MOMY_tp, RHOT_tp, RHOQ_tp, CORIOLI, CDZ, CDX, CDY, FDZ, FDX, FDY, RCDZ, RCDX, RCDY, RFDZ, RFDX, RFDY, PHI, GSQRT, J13G, J23G, J33G, MAPF, AQ_R, AQ_CV, AQ_CP, AQ_MASS, REF_dens, REF_pott, REF_qv, REF_pres, BND_W, BND_E, BND_S, BND_N, TwoD, ND_COEF, ND_COEF_Q, ND_LAPLACIAN_NUM, ND_SFC_FACT, ND_USE_RS, BND_QA, BND_IQ, BND_SMOOTHER_FACT, DAMP_DENS, DAMP_VELZ, DAMP_VELX, DAMP_VELY, DAMP_POTT, DAMP_QTRC, DAMP_alpha_DENS, DAMP_alpha_VELZ, DAMP_alpha_VELX, DAMP_alpha_VELY, DAMP_alpha_POTT, DAMP_alpha_QTRC, MFLUX_OFFSET_X, MFLUX_OFFSET_Y, wdamp_coef, divdmp_coef, FLAG_TRACER_SPLIT_TEND, FLAG_FCT_MOMENTUM, FLAG_FCT_T, FLAG_FCT_TRACER, FLAG_FCT_ALONG_STREAM, USE_AVERAGE, I_QV, DTLS, DTSS, Llast)
Dynamical Process.
Definition: scale_atmos_dyn_tstep_large_fvm_heve.F90:407
scale_atmos_dyn_common::atmos_dyn_fill_halo
subroutine, public atmos_dyn_fill_halo(var, fill_constval, lateral_halo, top_bottom_halo)
Definition: scale_atmos_dyn_common.F90:129
scale_prc_cartesc
module process / cartesC
Definition: scale_prc_cartesC.F90:11
scale_spnudge::spnudge_qv_mm
integer, public spnudge_qv_mm
Definition: scale_spnudge.F90:35
scale_tracer::tracer_name
character(len=h_short), dimension(qa_max), public tracer_name
Definition: scale_tracer.F90:38
scale_prof
module profiler
Definition: scale_prof.F90:11
scale_atmos_dyn_tstep_large_fvm_heve
module Atmosphere / Dynamics
Definition: scale_atmos_dyn_tstep_large_fvm_heve.F90:12
scale_atmos_grid_cartesc_index::is
integer, public is
start point of inner domain: x, local
Definition: scale_atmos_grid_cartesC_index.F90:53
scale_atmos_dyn_fvm_numfilter
module Atmosphere / Dynamics common
Definition: scale_atmos_dyn_fvm_numfilter.F90:12
scale_spnudge::spnudge_uv_mm
integer, public spnudge_uv_mm
Definition: scale_spnudge.F90:27
scale_atmos_dyn_tinteg_short::atmos_dyn_tinteg_short
procedure(short), pointer, public atmos_dyn_tinteg_short
Definition: scale_atmos_dyn_tinteg_short.F90:109
scale_precision::dp
integer, parameter, public dp
Definition: scale_precision.F90:32
scale_monitor::monitor_reg
subroutine, public monitor_reg(name, desc, unit, itemid, ndims, dim_type, is_tendency)
Search existing item, or matching check between requested and registered item.
Definition: scale_monitor.F90:241
scale_atmos_grid_cartesc_index::ja
integer, public ja
Definition: scale_atmos_grid_cartesC_index.F90:49
scale_atmos_dyn_fvm_numfilter::atmos_dyn_fvm_numfilter_flux_q
subroutine, public atmos_dyn_fvm_numfilter_flux_q(num_diff_q, DENS, QTRC, is_qv, CDZ, CDX, CDY, TwoD, dt, REF_qv, iq, ND_COEF, ND_LAPLACIAN_NUM, ND_SFC_FACT, ND_USE_RS)
Calculate fluxes with numerical filter for tracer variables.
Definition: scale_atmos_dyn_fvm_numfilter.F90:1160
scale_dft::dft_g2g
subroutine, public dft_g2g(KA, KS, KE, IA, IS, IE, JA, JS, JE, LM, MM, f)
Definition: scale_dft.F90:214
scale_tracer
module TRACER
Definition: scale_tracer.F90:12
scale_atmos_dyn_fvm_numfilter::atmos_dyn_fvm_apply_numfilter
subroutine, public atmos_dyn_fvm_apply_numfilter(num_diff, DENS, MOMZ, MOMX, MOMY, RHOT, CDZ, CDX, CDY, FDZ, FDX, FDY, RCDZ, RCDX, RCDY, RFDZ, RFDX, RFDY, TwoD, DT, GSQRT, MAPF, REF_dens, REF_pott, ND_COEF, ND_LAPLACIAN_NUM, ND_SFC_FACT, ND_USE_RS)
Definition: scale_atmos_dyn_fvm_numfilter.F90:1490
scale_atmos_dyn_tstep_large_fvm_heve::monitor_put_lateral_flux
subroutine monitor_put_lateral_flux(MONIT_ID, BND_flag, flx, flx_zero)
Definition: scale_atmos_dyn_tstep_large_fvm_heve.F90:1706
scale_spnudge
Definition: scale_spnudge.F90:2
scale_atmos_grid_cartesc_index::xdir
integer, parameter, public xdir
Definition: scale_atmos_grid_cartesC_index.F90:33
scale_atmos_grid_cartesc_index::ks
integer, public ks
start point of inner domain: z, local
Definition: scale_atmos_grid_cartesC_index.F90:51
scale_atmos_dyn_tstep_large_fvm_heve::multiply_flux_by_metric_xyz
subroutine multiply_flux_by_metric_xyz(I_DIR, flx, GSQRT, MAPF)
Definition: scale_atmos_dyn_tstep_large_fvm_heve.F90:1728
scale_atmos_dyn_tinteg_tracer
module Atmosphere / Dynamics Temporal integration
Definition: scale_atmos_dyn_tinteg_tracer.F90:12
scale_debug
module DEBUG
Definition: scale_debug.F90:11
scale_atmos_dyn_common::atmos_dyn_divergence
subroutine, public atmos_dyn_divergence(DDIV, MOMZ, MOMX, MOMY, GSQRT, J13G, J23G, J33G, MAPF, TwoD, RCDZ, RCDX, RCDY, RFDZ, FDZ)
Definition: scale_atmos_dyn_common.F90:382
scale_comm_cartesc
module COMMUNICATION
Definition: scale_comm_cartesC.F90:11
scale_atmos_grid_cartesc_index::js
integer, public js
start point of inner domain: y, local
Definition: scale_atmos_grid_cartesC_index.F90:55
scale_atmos_grid_cartesc_index::ydir
integer, parameter, public ydir
Definition: scale_atmos_grid_cartesC_index.F90:34
scale_atmos_grid_cartesc_index::ieb
integer, public ieb
Definition: scale_atmos_grid_cartesC_index.F90:64
scale_spnudge::spnudge_qv_alpha
real(rp), dimension(:,:,:), allocatable, public spnudge_qv_alpha
Definition: scale_spnudge.F90:40
scale_file_history::file_history_reg
subroutine, public file_history_reg(name, desc, unit, itemid, standard_name, ndims, dim_type, cell_measures, fill_halo)
Register/Append variable to history file.
Definition: scale_file_history.F90:650
scale_prof::prof_rapend
subroutine, public prof_rapend(rapname_base, level, disable_barrier)
Save raptime.
Definition: scale_prof.F90:217
scale_atmos_dyn_common::atmos_dyn_prep_pres_linearization
subroutine, public atmos_dyn_prep_pres_linearization(DPRES, RT2P, REF_rhot, RHOT, QTRC, REF_pres, AQ_R, AQ_CV, AQ_CP, AQ_MASS)
Definition: scale_atmos_dyn_common.F90:510
scale_const::const_undef
real(rp), public const_undef
Definition: scale_const.F90:41
scale_atmos_grid_cartesc_index::jsb
integer, public jsb
Definition: scale_atmos_grid_cartesC_index.F90:65
scale_prc_cartesc::prc_has_w
logical, public prc_has_w
Definition: scale_prc_cartesC.F90:47
scale_atmos_dyn_tstep_large_fvm_heve::atmos_dyn_tstep_large_fvm_heve_setup
subroutine, public atmos_dyn_tstep_large_fvm_heve_setup(DENS, MOMZ, MOMX, MOMY, RHOT, QTRC, PROG)
Setup.
Definition: scale_atmos_dyn_tstep_large_fvm_heve.F90:123
scale_const::const_ohm
real(rp), public const_ohm
angular velocity of the planet [1/s]
Definition: scale_const.F90:45
scale_io::io_fid_conf
integer, public io_fid_conf
Config file ID.
Definition: scale_io.F90:56
scale_atmos_grid_cartesc_index::je
integer, public je
end point of inner domain: y, local
Definition: scale_atmos_grid_cartesC_index.F90:56
scale_prc_cartesc::prc_twod
logical, public prc_twod
2D experiment
Definition: scale_prc_cartesC.F90:55
scale_spnudge::spnudge_uv
logical, public spnudge_uv
Definition: scale_spnudge.F90:24
scale_spnudge::spnudge_v_alpha
real(rp), dimension(:,:,:), allocatable, public spnudge_v_alpha
Definition: scale_spnudge.F90:38
scale_atmos_grid_cartesc_index::i_xyw
integer, public i_xyw
Definition: scale_atmos_grid_cartesC_index.F90:91
scale_monitor
module MONITOR
Definition: scale_monitor.F90:12
Generated by   doxygen 1.8.17