FE-Project
Loading...
Searching...
No Matches
Functions/Subroutines
scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu Module Reference

module FElib / Fluid dyn solver / Atmosphere / Global nonhydrostatic model / HEVE More...

Functions/Subroutines

subroutine, public atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_init (mesh)
subroutine, public atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_final ()
subroutine, public atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_cal_tend_shallow_atm (dens_dt, momx_dt, momy_dt, momz_dt, rhot_dt, ddens_, momx_, momy_, momz_, drhot_, dpres_, dens_hyd, pres_hyd, pres_hyd_ref, therm_hyd, coriolis, rtot, cvtot, cptot, dphyddx, dphyddy, element3d_operation, dx, dy, dz, sx, sy, sz, lift, lmesh, elem, lmesh2d, elem2d)
subroutine, public atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_cal_tend_deep_atm (dens_dt, momx_dt, momy_dt, momz_dt, rhot_dt, ddens_, momx_, momy_, momz_, drhot_, dpres_, dens_hyd, pres_hyd, pres_hyd_ref, therm_hyd, coriolis, rtot, cvtot, cptot, dphyddx, dphyddy, element3d_operation, dx, dy, dz, sx, sy, sz, lift, lmesh, elem, lmesh2d, elem2d)

Detailed Description

module FElib / Fluid dyn solver / Atmosphere / Global nonhydrostatic model / HEVE

Description
HEVE DGM scheme for Global Atmospheric Dynamical process. The governing equations is a fully compressible nonhydrostatic equations, which consist of mass, momentum, and thermodynamics (density * potential temperature conservation) equations.
Author
Yuta Kawai, Xuanzhengbo Ren, and Team SCALE

Function/Subroutine Documentation

◆ atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_init()

subroutine, public scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu::atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_init ( class(meshbase3d), intent(in) mesh)

Definition at line 75 of file scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu.F90.

76 implicit none
77 class(MeshBase3D), intent(in) :: mesh
78 !--------------------------------------------
79
80 call atm_dyn_dgm_nonhydro3d_common_init( mesh )
81
82 return

References scale_atm_dyn_dgm_nonhydro3d_common::atm_dyn_dgm_nonhydro3d_common_init().

◆ atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_final()

subroutine, public scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu::atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_final

Definition at line 86 of file scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu.F90.

87 implicit none
88 !--------------------------------------------
89
90 call atm_dyn_dgm_nonhydro3d_common_final()
91 return

References scale_atm_dyn_dgm_nonhydro3d_common::atm_dyn_dgm_nonhydro3d_common_final().

◆ atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_cal_tend_shallow_atm()

subroutine, public scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu::atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_cal_tend_shallow_atm ( real(rp), dimension(elem%np,lmesh%nea), intent(out) dens_dt,
real(rp), dimension(elem%np,lmesh%nea), intent(out) momx_dt,
real(rp), dimension(elem%np,lmesh%nea), intent(out) momy_dt,
real(rp), dimension(elem%np,lmesh%nea), intent(out) momz_dt,
real(rp), dimension(elem%np,lmesh%nea), intent(out) rhot_dt,
real(rp), dimension(elem%np,lmesh%nea), intent(in) ddens_,
real(rp), dimension(elem%np,lmesh%nea), intent(in) momx_,
real(rp), dimension(elem%np,lmesh%nea), intent(in) momy_,
real(rp), dimension(elem%np,lmesh%nea), intent(in) momz_,
real(rp), dimension(elem%np,lmesh%nea), intent(in) drhot_,
real(rp), dimension(elem%np,lmesh%nea), intent(in) dpres_,
real(rp), dimension(elem%np,lmesh%nea), intent(in) dens_hyd,
real(rp), dimension(elem%np,lmesh%nea), intent(in) pres_hyd,
real(rp), dimension(elem%np,lmesh%nea), intent(in) pres_hyd_ref,
real(rp), dimension(elem%np,lmesh%nea), intent(in) therm_hyd,
real(rp), dimension(elem2d%np,lmesh2d%nea), intent(in) coriolis,
real(rp), dimension (elem%np,lmesh%nea), intent(in) rtot,
real(rp), dimension(elem%np,lmesh%nea), intent(in) cvtot,
real(rp), dimension(elem%np,lmesh%nea), intent(in) cptot,
real(rp), dimension(elem%np,lmesh%nea), intent(in) dphyddx,
real(rp), dimension(elem%np,lmesh%nea), intent(in) dphyddy,
class(elementoperationbase3d), intent(in) element3d_operation,
type(sparsemat), intent(in) dx,
type(sparsemat), intent(in) dy,
type(sparsemat), intent(in) dz,
type(sparsemat), intent(in) sx,
type(sparsemat), intent(in) sy,
type(sparsemat), intent(in) sz,
type(sparsemat), intent(in) lift,
class(localmesh3d), intent(in) lmesh,
class(elementbase3d), intent(in) elem,
class(localmesh2d), intent(in) lmesh2d,
class(elementbase2d), intent(in) elem2d )

Definition at line 97 of file scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu.F90.

104
105 use scale_atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_gpu, only: &
106 get_ebnd_flux => atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_get_generalhvc_gpu
107 use scale_element_operation_gpu_driver, only: &
108 elementoperationgpudriver
109 implicit none
110
111 class(LocalMesh3D), intent(in) :: lmesh
112 class(ElementBase3D), intent(in) :: elem
113 class(LocalMesh2D), intent(in) :: lmesh2D
114 class(ElementBase2D), intent(in) :: elem2D
115 class(ElementOperationBase3D), intent(in) :: element3D_operation
116 type(SparseMat), intent(in) :: Dx, Dy, Dz, Sx, Sy, Sz, Lift
117 real(RP), intent(out) :: DENS_dt(elem%Np,lmesh%NeA)
118 real(RP), intent(out) :: MOMX_dt(elem%Np,lmesh%NeA)
119 real(RP), intent(out) :: MOMY_dt(elem%Np,lmesh%NeA)
120 real(RP), intent(out) :: MOMZ_dt(elem%Np,lmesh%NeA)
121 real(RP), intent(out) :: RHOT_dt(elem%Np,lmesh%NeA)
122 real(RP), intent(in) :: DDENS_(elem%Np,lmesh%NeA)
123 real(RP), intent(in) :: MOMX_(elem%Np,lmesh%NeA)
124 real(RP), intent(in) :: MOMY_(elem%Np,lmesh%NeA)
125 real(RP), intent(in) :: MOMZ_(elem%Np,lmesh%NeA)
126 real(RP), intent(in) :: DRHOT_(elem%Np,lmesh%NeA)
127 real(RP), intent(in) :: DPRES_(elem%Np,lmesh%NeA)
128 real(RP), intent(in) :: DENS_hyd(elem%Np,lmesh%NeA)
129 real(RP), intent(in) :: PRES_hyd(elem%Np,lmesh%NeA)
130 real(RP), intent(in) :: PRES_hyd_ref(elem%Np,lmesh%NeA)
131 real(RP), intent(in) :: THERM_hyd(elem%Np,lmesh%NeA)
132 real(RP), intent(in) :: CORIOLIS(elem2D%Np,lmesh2D%NeA)
133 real(RP), intent(in) :: Rtot (elem%Np,lmesh%NeA)
134 real(RP), intent(in) :: CVtot(elem%Np,lmesh%NeA)
135 real(RP), intent(in) :: CPtot(elem%Np,lmesh%NeA)
136 real(RP), intent(in) :: DPhydDx(elem%Np,lmesh%NeA)
137 real(RP), intent(in) :: DPhydDy(elem%Np,lmesh%NeA)
138
139 type(ElementOperationGPUDriver) :: element3D_operation_driver
140 real(RP) :: tend_tmp(elem%Nnode_h1D**2,5,elem%Nnode_v,lmesh%Ne)
141 real(RP) :: Flux2D(elem%Nnode_h1D**2,5,elem%Nnode_v,lmesh%Ne,2)
142 real(RP) :: FluxZ_store(elem%Nnode_h1D**2,5,elem%Nnode_v,lmesh%Ne)
143 real(RP) :: del_flux(elem%NfpTot,PRGVAR_NUM,lmesh%Ne)
144 real(RP) :: drho(elem%Np,lmesh%Ne)
145
146 integer :: IndexH2Dto3D(elem%Np)
147
148 integer :: NeS, NeE, Nnode_h1D, Nnode_v
149 !------------------------------------------------------------------------
150
151 indexh2dto3d(:) = elem%IndexH2Dto3D(:)
152 call element3d_operation_driver%Init( element3d_operation )
153 nes = lmesh%NeS
154 nee = lmesh%NeE
155 nnode_h1d = elem%Nnode_h1D
156 nnode_v = elem%Nnode_v
157
158 !$acc data present(DDENS_,MOMX_,MOMY_,MOMZ_,DRHOT_,DPRES_, &
159 !$acc DENS_hyd,PRES_hyd,THERM_hyd, &
160 !$acc CORIOLIS,Rtot,CVtot,CPtot,DPhydDx,DPhydDy, &
161 !$acc DENS_dt,MOMX_dt,MOMY_dt,MOMZ_dt,RHOT_dt, lmesh,elem) &
162 !$acc create(del_flux,Flux2D,FluxZ_store,tend_tmp,drho) copyin(IndexH2Dto3D)
163
164 call prof_rapstart('cal_dyn_tend_bndflux', 3)
165 call get_ebnd_flux( &
166 del_flux, & ! (out)
167 ddens_, momx_, momy_, momz_, drhot_, dpres_, & ! (in)
168 dens_hyd, pres_hyd, therm_hyd, rtot, cvtot, cptot, & ! (in)
169 lmesh%Gsqrt, lmesh%GIJ(:,:,1,1), lmesh%GIJ(:,:,1,2), lmesh%GIJ(:,:,2,2), & ! (in)
170 lmesh%GsqrtH, lmesh%gam, lmesh%GI3(:,:,1), lmesh%GI3(:,:,2), & ! (in)
171 lmesh%normal_fn(:,:,1), lmesh%normal_fn(:,:,2), lmesh%normal_fn(:,:,3), & ! (in)
172 lmesh%vmapM, lmesh%vmapP, elem%IndexH2Dto3D_bnd, & ! (in)
173 lmesh, elem, lmesh2d, elem2d ) ! (in)
174 !$acc wait(1)
175 call prof_rapend('cal_dyn_tend_bndflux', 3)
176
177 !-----
178 call prof_rapstart('cal_dyn_tend_interior', 3)
179 call cal_tend_interior_gpu( &
180 dens_dt, momx_dt, momy_dt, momz_dt, rhot_dt, & ! (out)
181 ddens_, momx_, momy_, momz_, drhot_, dpres_, del_flux, & ! (in)
182 dens_hyd, therm_hyd, dphyddx, dphyddy, & ! (in)
183 lmesh%Gsqrt,lmesh%GsqrtH,lmesh%GI3(:,:,1),lmesh%GI3(:,:,2), & ! (in)
184 lmesh%Escale(:,:,1,1), lmesh%Escale(:,:,2,2), lmesh%Escale(:,:,3,3), & ! (in)
185 lmesh2d%pos_en(:,:,1), lmesh2d%pos_en(:,:,2), & ! (in)
186 flux2d, fluxz_store, tend_tmp, drho, & ! (in)
187 element3d_operation_driver, lmesh%EMap3Dto2D, indexh2dto3d, & ! (in)
188 lmesh,elem, lmesh%NeS, lmesh%NeE, lmesh%NeA, lmesh%Ne2DA, elem%Nnode_h1D, elem%Nnode_v ) ! (in) ! (in)
189 !$acc wait(1)
190 call prof_rapend('cal_dyn_tend_interior', 3)
191
192 !$acc end data
193
194 call element3d_operation_driver%Final()
195 return
scale_atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_gpu
module FElib / Fluid dyn solver / Atmosphere / Nonhydrostatic model / HEVE / Numflux
Definition scale_atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_gpu.F90:11
scale_atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_gpu::atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_get_generalhvc_gpu
subroutine, public atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_get_generalhvc_gpu(del_flux, ddens_, momx_, momy_, momz_, drhot_, dpres, dens_hyd, pres_hyd, therm_hyd, rtot, cvtot, cptot, gsqrt, g11, g12, g22, gsqrth, gam, g13, g23, nx, ny, nz, vmapm, vmapp, im2dto3d, lmesh, elem, lmesh2d, elem2d)
Definition scale_atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_gpu.F90:260
scale_element_operation_gpu_driver
module FElib / Element / Driver for operation with 3D tensor product elements using GPU
Definition scale_element_operation_gpu_driver.F90:16
scale_element_operation_gpu_driver::elementoperationgpudriver
Driver for element operation with 3D tensor product elements using GPU.
Definition scale_element_operation_gpu_driver.F90:39

References scale_atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_gpu::atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_get_generalhvc_gpu(), scale_element_operation_gpu_driver::elementoperationgpu_div_kplane(), scale_element_operation_gpu_driver::elementoperationgpu_divvar5(), scale_element_operation_gpu_driver::elementoperationgpu_divvar5_kplane2(), scale_element_operation_gpu_driver::elementoperationgpu_divvar5_z_lift(), and scale_element_operation_gpu_driver::elementoperationgpu_vfilterpm1().

◆ atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_cal_tend_deep_atm()

subroutine, public scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu::atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_cal_tend_deep_atm ( real(rp), dimension(elem%np,lmesh%nea), intent(out) dens_dt,
real(rp), dimension(elem%np,lmesh%nea), intent(out) momx_dt,
real(rp), dimension(elem%np,lmesh%nea), intent(out) momy_dt,
real(rp), dimension(elem%np,lmesh%nea), intent(out) momz_dt,
real(rp), dimension(elem%np,lmesh%nea), intent(out) rhot_dt,
real(rp), dimension(elem%np,lmesh%nea), intent(in) ddens_,
real(rp), dimension(elem%np,lmesh%nea), intent(in) momx_,
real(rp), dimension(elem%np,lmesh%nea), intent(in) momy_,
real(rp), dimension(elem%np,lmesh%nea), intent(in) momz_,
real(rp), dimension(elem%np,lmesh%nea), intent(in) drhot_,
real(rp), dimension(elem%np,lmesh%nea), intent(in) dpres_,
real(rp), dimension(elem%np,lmesh%nea), intent(in) dens_hyd,
real(rp), dimension(elem%np,lmesh%nea), intent(in) pres_hyd,
real(rp), dimension(elem%np,lmesh%nea), intent(in) pres_hyd_ref,
real(rp), dimension(elem%np,lmesh%nea), intent(in) therm_hyd,
real(rp), dimension(elem2d%np,lmesh2d%nea), intent(in) coriolis,
real(rp), dimension (elem%np,lmesh%nea), intent(in) rtot,
real(rp), dimension(elem%np,lmesh%nea), intent(in) cvtot,
real(rp), dimension(elem%np,lmesh%nea), intent(in) cptot,
real(rp), dimension(elem%np,lmesh%nea), intent(in) dphyddx,
real(rp), dimension(elem%np,lmesh%nea), intent(in) dphyddy,
class(elementoperationbase3d), intent(in) element3d_operation,
type(sparsemat), intent(in) dx,
type(sparsemat), intent(in) dy,
type(sparsemat), intent(in) dz,
type(sparsemat), intent(in) sx,
type(sparsemat), intent(in) sy,
type(sparsemat), intent(in) sz,
type(sparsemat), intent(in) lift,
class(localmesh3d), intent(in) lmesh,
class(elementbase3d), intent(in) elem,
class(localmesh2d), intent(in) lmesh2d,
class(elementbase2d), intent(in) elem2d )

Definition at line 399 of file scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu.F90.

406
407 use scale_atm_dyn_dgm_nonhydro3d_rhot_heve_numflux, only: &
408 get_ebnd_flux => atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_get_generalhvc_asis
409 use scale_const, only: &
410 ohm => const_ohm
411 implicit none
412
413 class(LocalMesh3D), intent(in) :: lmesh
414 class(ElementBase3D), intent(in) :: elem
415 class(LocalMesh2D), intent(in) :: lmesh2D
416 class(ElementBase2D), intent(in) :: elem2D
417 class(ElementOperationBase3D), intent(in) :: element3D_operation
418 type(SparseMat), intent(in) :: Dx, Dy, Dz, Sx, Sy, Sz, Lift
419 real(RP), intent(out) :: DENS_dt(elem%Np,lmesh%NeA)
420 real(RP), intent(out) :: MOMX_dt(elem%Np,lmesh%NeA)
421 real(RP), intent(out) :: MOMY_dt(elem%Np,lmesh%NeA)
422 real(RP), intent(out) :: MOMZ_dt(elem%Np,lmesh%NeA)
423 real(RP), intent(out) :: RHOT_dt(elem%Np,lmesh%NeA)
424 real(RP), intent(in) :: DDENS_(elem%Np,lmesh%NeA)
425 real(RP), intent(in) :: MOMX_(elem%Np,lmesh%NeA)
426 real(RP), intent(in) :: MOMY_(elem%Np,lmesh%NeA)
427 real(RP), intent(in) :: MOMZ_(elem%Np,lmesh%NeA)
428 real(RP), intent(in) :: DRHOT_(elem%Np,lmesh%NeA)
429 real(RP), intent(in) :: DPRES_(elem%Np,lmesh%NeA)
430 real(RP), intent(in) :: DENS_hyd(elem%Np,lmesh%NeA)
431 real(RP), intent(in) :: PRES_hyd(elem%Np,lmesh%NeA)
432 real(RP), intent(in) :: PRES_hyd_ref(elem%Np,lmesh%NeA)
433 real(RP), intent(in) :: THERM_hyd(elem%Np,lmesh%NeA)
434 real(RP), intent(in) :: CORIOLIS(elem2D%Np,lmesh2D%NeA)
435 real(RP), intent(in) :: Rtot (elem%Np,lmesh%NeA)
436 real(RP), intent(in) :: CVtot(elem%Np,lmesh%NeA)
437 real(RP), intent(in) :: CPtot(elem%Np,lmesh%NeA)
438 real(RP), intent(in) :: DPhydDx(elem%Np,lmesh%NeA)
439 real(RP), intent(in) :: DPhydDy(elem%Np,lmesh%NeA)
440
441 real(RP) :: Fx(elem%Np), Fy(elem%Np), Fz(elem%Np), LiftDelFlx(elem%Np)
442 real(RP) :: DPRES_hyd(elem%Np), GradPhyd_x(elem%Np), GradPhyd_y(elem%Np)
443 real(RP) :: del_flux(elem%NfpTot,lmesh%Ne,PRGVAR_NUM)
444 real(RP) :: del_flux_hyd(elem%NfpTot,lmesh%Ne,2)
445 real(RP) :: RHOT_(elem%Np)
446 real(RP) :: rdens_(elem%Np), u_(elem%Np), v_(elem%Np), w_(elem%Np), wt_(elem%np), drho(elem%Np)
447
448 real(RP) :: G11(elem%Np), G12(elem%Np), G22(elem%Np)
449 real(RP) :: GsqrtV(elem%Np), RGsqrtV(elem%Np), Rgam2(elem%Np)
450 real(RP) :: X2D(elem2D%Np,lmesh2D%Ne), Y2D(elem2D%Np,lmesh2D%Ne)
451 real(RP) :: X(elem%Np), Y(elem%Np), twoOVdel2(elem%Np)
452 real(RP) :: OM1(elem%Np), OM2(elem%Np), OM3(elem%Np), DEL(elem%Np), R(elem%Np)
453 logical :: is_panel1to4
454 real(RP) :: s
455
456 integer :: ke, ke2d
457 integer :: p
458
459 real(RP) :: rgamm
460 real(RP) :: rP0
461 real(RP) :: P0ovR
462 !------------------------------------------------------------------------
463
464 call prof_rapstart('cal_dyn_tend_bndflux', 3)
465 call get_ebnd_flux( &
466 del_flux, del_flux_hyd, & ! (out)
467 ddens_, momx_, momy_, momz_, drhot_, dpres_, dens_hyd, pres_hyd, & ! (in)
468 rtot, cvtot, cptot, & ! (in)
469 lmesh%Gsqrt, lmesh%GIJ(:,:,1,1), lmesh%GIJ(:,:,1,2), lmesh%GIJ(:,:,2,2), & ! (in)
470 lmesh%GsqrtH, lmesh%gam, lmesh%GI3(:,:,1), lmesh%GI3(:,:,2), & ! (in)
471 lmesh%normal_fn(:,:,1), lmesh%normal_fn(:,:,2), lmesh%normal_fn(:,:,3), & ! (in)
472 lmesh%vmapM, lmesh%vmapP, elem%IndexH2Dto3D_bnd, & ! (in)
473 lmesh, elem, lmesh2d, elem2d ) ! (in)
474 call prof_rapend('cal_dyn_tend_bndflux', 3)
475
476 !-----
477 call prof_rapstart('cal_dyn_tend_interior', 3)
478 rgamm = cvdry / cpdry
479 rp0 = 1.0_rp / pres00
480 p0ovr = pres00 / rdry
481
482 s = 2.0_rp * ohm
483 is_panel1to4 = .true.
484 if ( lmesh%panelID == 5 ) then
485 is_panel1to4 = .false.
486 else if ( lmesh%panelID == 6 ) then
487 is_panel1to4 = .false.
488 s = - s
489 end if
490
491 !$omp parallel private( &
492 !$omp RHOT_, rdens_, u_, v_, w_, wt_, &
493 !$omp Fx, Fy, Fz, LiftDelFlx, &
494 !$omp drho, DPRES_hyd, GradPhyd_x, GradPhyd_y, &
495 !$omp G11, G12, G22, Rgam2, GsqrtV, RGsqrtV, &
496 !$omp X, Y, twoOVdel2, &
497 !$omp OM1, OM2, OM3, DEL, R, ke, ke2D )
498
499 !$omp do
500 do ke2d = lmesh2d%NeS, lmesh2d%NeE
501 x2d(:,ke2d) = tan(lmesh2d%pos_en(:,ke2d,1))
502 y2d(:,ke2d) = tan(lmesh2d%pos_en(:,ke2d,2))
503 end do
504
505 !$omp do
506 do ke = lmesh%NeS, lmesh%NeE
507 !--
508 ke2d = lmesh%EMap3Dto2D(ke)
509 rgam2(:) = 1.0_rp / lmesh%gam(:,ke)**2
510 g11(:) = lmesh%GIJ(elem%IndexH2Dto3D,ke2d,1,1) * rgam2(:)
511 g12(:) = lmesh%GIJ(elem%IndexH2Dto3D,ke2d,1,2) * rgam2(:)
512 g22(:) = lmesh%GIJ(elem%IndexH2Dto3D,ke2d,2,2) * rgam2(:)
513 gsqrtv(:) = lmesh%Gsqrt(:,ke) * rgam2(:) / lmesh%GsqrtH(elem%IndexH2Dto3D,ke2d)
514 rgsqrtv(:) = 1.0_rp / gsqrtv(:)
515
516 !--
517 rhot_(:) = p0ovr * ( pres_hyd(:,ke) * rp0 )**rgamm + drhot_(:,ke)
518 ! DPRES_(:) = PRES00 * ( Rtot(:,ke) * rP0 * RHOT_(:) )**( CPtot(:,ke) / CVtot(:,ke) ) &
519 ! - PRES_hyd(:,ke)
520
521 rdens_(:) = 1.0_rp / ( ddens_(:,ke) + dens_hyd(:,ke) )
522 u_(:) = momx_(:,ke) * rdens_(:)
523 v_(:) = momy_(:,ke) * rdens_(:)
524 w_(:) = momz_(:,ke) * rdens_(:)
525 wt_(:) = w_(:) * rgsqrtv(:) + lmesh%GI3(:,ke,1) * u_(:) + lmesh%GI3(:,ke,2) * v_(:)
526
527 x(:) = x2d(elem%IndexH2Dto3D,ke2d)
528 y(:) = y2d(elem%IndexH2Dto3D,ke2d)
529 del(:) = sqrt( 1.0_rp + x(:)**2 + y(:)**2 )
530 twoovdel2(:) = 2.0_rp / ( 1.0_rp + x(:)**2 + y(:)**2 )
531
532 r(:) = rplanet * lmesh%gam(:,ke)
533
534 !- pnl=1~4: OM1: 0, OM2: s del / (r (1+Y^2)) , s OM3 : Y /del
535 !- pnl=5,6: OM1: - s X del/(r (1+X^2)), OM2: - s Y del/(r(1+Y^2)), OM3 : s/del
536 if ( is_panel1to4 ) then
537 om1(:) = 0.0_rp
538 om2(:) = s * del(:) / ( r(:) * ( 1.0_rp + y(:)**2 ) )
539 om3(:) = s * y(:) / del(:)
540 else
541 om1(:) = - s * x(:) * del(:) / ( r(:) * ( 1.0_rp + x(:)**2 ) )
542 om2(:) = - s * y(:) * del(:) / ( r(:) * ( 1.0_rp + y(:)**2 ) )
543 om3(:) = s / del(:)
544 end if
545
546 drho(:) = matmul(intrpmat_vpordm1, ddens_(:,ke))
547
548 !-- Gradient hydrostatic pressure
549
550 dpres_hyd(:) = pres_hyd(:,ke) - pres_hyd_ref(:,ke)
551
552 call sparsemat_matmul(dx, gsqrtv(:) * dpres_hyd(:), fx)
553 call sparsemat_matmul(dz, gsqrtv(:) * lmesh%GI3(:,ke,1) * dpres_hyd(:), fz)
554 call sparsemat_matmul(lift, lmesh%Fscale(:,ke) * del_flux_hyd(:,ke,1), liftdelflx)
555 gradphyd_x(:) = lmesh%Escale(:,ke,1,1) * fx(:) &
556 + lmesh%Escale(:,ke,3,3) * fz(:) &
557 + liftdelflx(:)
558
559 call sparsemat_matmul(dy, gsqrtv(:) * dpres_hyd(:), fy)
560 call sparsemat_matmul(dz, gsqrtv(:) * lmesh%GI3(:,ke,2) * dpres_hyd(:), fz)
561 call sparsemat_matmul(lift, lmesh%Fscale(:,ke) * del_flux_hyd(:,ke,2), liftdelflx)
562 gradphyd_y(:) = lmesh%Escale(:,ke,2,2) * fy(:) &
563 + lmesh%Escale(:,ke,3,3) * fz(:) &
564 + liftdelflx(:)
565
566 !-- DENS
567 call sparsemat_matmul(dx, lmesh%Gsqrt(:,ke) * momx_(:,ke), fx)
568 call sparsemat_matmul(dy, lmesh%Gsqrt(:,ke) * momy_(:,ke), fy)
569 call sparsemat_matmul(dz, lmesh%Gsqrt(:,ke) * ( ddens_(:,ke) + dens_hyd(:,ke) ) * wt_(:), fz)
570 call sparsemat_matmul(lift, lmesh%Fscale(:,ke) * del_flux(:,ke,dens_vid), liftdelflx)
571
572 dens_dt(:,ke) = - ( &
573 lmesh%Escale(:,ke,1,1) * fx(:) &
574 + lmesh%Escale(:,ke,2,2) * fy(:) &
575 + lmesh%Escale(:,ke,3,3) * fz(:) &
576 + liftdelflx(:) ) / lmesh%Gsqrt(:,ke)
577
578 !-- MOMX
579 call sparsemat_matmul(dx, lmesh%Gsqrt(:,ke) * ( u_(:) * momx_(:,ke) + g11(:) * dpres_(:,ke) ), fx)
580 call sparsemat_matmul(dy, lmesh%Gsqrt(:,ke) * ( v_(:) * momx_(:,ke) + g12(:) * dpres_(:,ke) ), fy)
581 call sparsemat_matmul(dz, lmesh%Gsqrt(:,ke) * ( wt_(:) * momx_(:,ke) &
582 + ( lmesh%GI3(:,ke,1) * g11(:) + lmesh%GI3(:,ke,2) * g12(:) ) * dpres_(:,ke) ), fz)
583 call sparsemat_matmul(lift, lmesh%Fscale(:,ke) * del_flux(:,ke,momx_vid), liftdelflx)
584
585 momx_dt(:,ke) = &
586 - ( lmesh%Escale(:,ke,1,1) * fx(:) &
587 + lmesh%Escale(:,ke,2,2) * fy(:) &
588 + lmesh%Escale(:,ke,3,3) * fz(:) &
589 + liftdelflx(:) ) / lmesh%Gsqrt(:,ke) &
590 - twoovdel2(:) * y(:) * & !-> metric terms
591 ( x(:) * y(:) * u_(:) - ( 1.0_rp + y(:)**2 ) * v_(:) ) * momx_(:,ke) & !
592 - 2.0_rp * u_(:) * momz_(:,ke) / r(:) & !<-
593 - ( g11(:) * gradphyd_x(:) + g12(:) * gradphyd_y(:) ) * rgsqrtv(:) & !-> gradient hydrostatic pressure
594 - lmesh%Gsqrt(:,ke) * ( g11(:) * ( om2(:) * momz_(:,ke) - om3(:) * momy_(:,ke) ) & !-> Coriolis term
595 - g12(:) * ( om1(:) * momz_(:,ke) - om3(:) * momx_(:,ke) ) )
596
597 !-- MOMY
598 call sparsemat_matmul(dx, lmesh%Gsqrt(:,ke) * ( u_(:) * momy_(:,ke) + g12(:) * dpres_(:,ke) ), fx)
599 call sparsemat_matmul(dy, lmesh%Gsqrt(:,ke) * ( v_(:) * momy_(:,ke) + g22(:) * dpres_(:,ke) ), fy)
600 call sparsemat_matmul(dz, lmesh%Gsqrt(:,ke) * ( wt_(:) * momy_(:,ke) &
601 + ( lmesh%GI3(:,ke,1) * g12(:) + lmesh%GI3(:,ke,2) * g22(:) ) * dpres_(:,ke) ), fz)
602 call sparsemat_matmul(lift, lmesh%Fscale(:,ke) * del_flux(:,ke,momy_vid), liftdelflx)
603
604 momy_dt(:,ke) = &
605 - ( lmesh%Escale(:,ke,1,1) * fx(:) &
606 + lmesh%Escale(:,ke,2,2) * fy(:) &
607 + lmesh%Escale(:,ke,3,3) * fz(:) &
608 + liftdelflx(:) ) / lmesh%Gsqrt(:,ke) &
609 - twoovdel2(:) * x(:) * & !-> metric terms
610 ( - (1.0_rp + x(:)**2) * u_(:) + x(:) * y(:) * v_(:) ) * momy_(:,ke) & !
611 - 2.0_rp * v_(:) * momz_(:,ke) / r(:) & !<-
612 - ( g12(:) * gradphyd_x(:) + g22(:) * gradphyd_y(:) ) * rgsqrtv(:) & !-> gradient hydrostatic pressure
613 - lmesh%Gsqrt(:,ke) * ( g12(:) * ( om2(:) * momz_(:,ke) - om3(:) * momy_(:,ke) ) & !-> Coriolis term
614 - g22(:) * ( om1(:) * momz_(:,ke) - om3(:) * momx_(:,ke) ) )
615
616 !-- MOMZ
617 call sparsemat_matmul(dx, lmesh%Gsqrt(:,ke) * u_(:) * momz_(:,ke), fx)
618 call sparsemat_matmul(dy, lmesh%Gsqrt(:,ke) * v_(:) * momz_(:,ke), fy)
619 call sparsemat_matmul(dz, lmesh%Gsqrt(:,ke) * ( wt_(:) * momz_(:,ke) + rgsqrtv(:) * dpres_(:,ke) ), fz)
620 call sparsemat_matmul(lift, lmesh%Fscale(:,ke) * del_flux(:,ke,momz_vid), liftdelflx)
621
622 momz_dt(:,ke) = &
623 - ( lmesh%Escale(:,ke,1,1) * fx(:) &
624 + lmesh%Escale(:,ke,2,2) * fy(:) &
625 + lmesh%Escale(:,ke,3,3) * fz(:) &
626 + liftdelflx(:) ) / lmesh%Gsqrt(:,ke) &
627 - 0.25_rp * r(:) * twoovdel2(:)**2 * ( 1.0_rp * x(:)**2 ) * ( 1.0_rp * y(:)**2 ) & !-> metric terms
628 * ( - ( 1.0_rp + x(:)**2 ) * momx_(:,ke) * u_(:) & !
629 + 2.0_rp * x(:) * y(:) * momx_(:,ke) * v_(:) & !
630 - ( 1.0_rp + y(:)**2 ) * momy_(:,ke) * v_(:) ) & !<-
631 + 2.0_rp * dpres_(:,ke) / r(:) & !-> metric term with gradient of pressure deviaition
632 - lmesh%Gsqrt(:,ke) * ( om1(:) * momy_(:,ke) - om2(:) * momx_(:,ke) ) & !-> Coriolis term
633 - grav * rgam2(:) * drho(:) !-> buoyancy term
634
635 !-- RHOT
636 call sparsemat_matmul(dx, lmesh%Gsqrt(:,ke) * u_(:) * rhot_(:), fx)
637 call sparsemat_matmul(dy, lmesh%Gsqrt(:,ke) * v_(:) * rhot_(:), fy)
638 call sparsemat_matmul(dz, lmesh%Gsqrt(:,ke) * wt_(:) * rhot_(:), fz)
639 call sparsemat_matmul(lift, lmesh%Fscale(:,ke) * del_flux(:,ke,rhot_vid), liftdelflx)
640
641 rhot_dt(:,ke) = &
642 - ( lmesh%Escale(:,ke,1,1) * fx(:) &
643 + lmesh%Escale(:,ke,2,2) * fy(:) &
644 + lmesh%Escale(:,ke,3,3) * fz(:) &
645 + liftdelflx(:) ) / lmesh%Gsqrt(:,ke)
646 end do
647 !$omp end do
648 !$omp end parallel
649 call prof_rapend('cal_dyn_tend_interior', 3)
650
651 return
scale_atm_dyn_dgm_nonhydro3d_rhot_heve_numflux
module FElib / Fluid dyn solver / Atmosphere / Nonhydrostatic model / HEVE / Numflux
Definition scale_atm_dyn_dgm_nonhydro3d_rhot_heve_numflux.F90:11
scale_atm_dyn_dgm_nonhydro3d_rhot_heve_numflux::atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_get_generalhvc_asis
subroutine, public atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_get_generalhvc_asis(del_flux, del_flux_hyd, ddens_, momx_, momy_, momz_, drhot_, dpres_, dens_hyd, pres_hyd, rtot, cvtot, cptot, gsqrt, g11, g12, g22, gsqrth, gam, g13, g23, nx, ny, nz, vmapm, vmapp, im2dto3d, lmesh, elem, lmesh2d, elem2d)
Definition scale_atm_dyn_dgm_nonhydro3d_rhot_heve_numflux.F90:1357

References scale_atm_dyn_dgm_nonhydro3d_rhot_heve_numflux::atm_dyn_dgm_nonhydro3d_rhot_heve_numflux_get_generalhvc_asis(), and scale_atm_dyn_dgm_nonhydro3d_common::intrpmat_vpordm1.

Generated by doxygen 1.14.0