scale_atm_dyn_dgm_globalsw Module Reference

module FElib / Fluid dyn solver / Atmosphere / Global shallow water More...

Functions/Subroutines
subroutine, public	atm_dyn_dgm_globalsw_init (mesh)
subroutine, public	atm_dyn_dgm_globalsw_final ()
subroutine, public	atm_dyn_dgm_globalsw_cal_tend (h_dt, u_dt, v_dt, h_, u_, v_, hs_, u1_, u2_, coriolis, dx, dy, sx, sy, lift, lmesh, elem)

Detailed Description

module FElib / Fluid dyn solver / Atmosphere / Global shallow water

Description

DGM scheme for Atmospheric dynamical process.

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

atm_dyn_dgm_globalsw_init()

subroutine, public scale_atm_dyn_dgm_globalsw::atm_dyn_dgm_globalsw_init

(

class(meshbase2d), intent(in)

mesh

)

Definition at line 63 of file scale_atm_dyn_dgm_globalsw.F90.

    implicit none
    class(MeshBase2D), intent(in) :: mesh
    !--------------------------------------------
    return

atm_dyn_dgm_globalsw_final()

subroutine, public scale_atm_dyn_dgm_globalsw::atm_dyn_dgm_globalsw_final

Definition at line 71 of file scale_atm_dyn_dgm_globalsw.F90.

    implicit none
    !--------------------------------------------
    return

atm_dyn_dgm_globalsw_cal_tend()

subroutine, public scale_atm_dyn_dgm_globalsw::atm_dyn_dgm_globalsw_cal_tend	(	real(rp), dimension(elem%np,lmesh%nea), intent(out)	h_dt,
		real(rp), dimension(elem%np,lmesh%nea), intent(out)	u_dt,
		real(rp), dimension(elem%np,lmesh%nea), intent(out)	v_dt,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	h_,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	u_,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	v_,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	hs_,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	u1_,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	u2_,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	coriolis,
		type(sparsemat), intent(in)	dx,
		type(sparsemat), intent(in)	dy,
		type(sparsemat), intent(in)	sx,
		type(sparsemat), intent(in)	sy,
		type(sparsemat), intent(in)	lift,
		class(localmesh2d), intent(in)	lmesh,
		class(elementbase2d), intent(in)	elem )

Definition at line 80 of file scale_atm_dyn_dgm_globalsw.F90.

 
    implicit none
 
    class(LocalMesh2D), intent(in) :: lmesh
    class(ElementBase2D), intent(in) :: elem
    type(SparseMat), intent(in) :: Dx, Dy, Sx, Sy, Lift
    real(RP), intent(out) :: h_dt(elem%Np,lmesh%NeA)
    real(RP), intent(out) :: U_dt(elem%Np,lmesh%NeA)
    real(RP), intent(out) :: V_dt(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: h_(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: U_(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: V_(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: hs_(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: u1_(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: u2_(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: CORIOLIS(elem%Np,lmesh%NeA)
 
    real(RP) :: Fx(elem%Np), Fy(elem%Np), LiftDelFlx(elem%Np)
    real(RP) :: del_flux(elem%NfpTot,lmesh%Ne,PROG_VARS_NUM)
    real(RP) :: del_flux_aux(elem%NfpTot,lmesh%Ne,1)
    real(RP) :: u1_dt(elem%Np), u2_dt(elem%Np)
    real(RP) :: VOR(elem%Np), E(elem%Np)
    integer :: ke       
    !------------------------------------------------------------------------
 
    call prof_rapstart('cal_dyn_tend_bndflux', 3)
    call cal_del_flux_dyn( del_flux, del_flux_aux,         & ! (out)
      h_, u_, v_, hs_, u1_, u2_,                           & ! (in)
      lmesh%Gsqrt, lmesh%GIJ(:,:,1,1), lmesh%GIJ(:,:,2,2), & ! (in)
      lmesh%normal_fn(:,:,1), lmesh%normal_fn(:,:,2),      & ! (in)
      lmesh%vmapM, lmesh%vmapP,                            & ! (in)
      lmesh, elem )                                          ! (in)
    call prof_rapend('cal_dyn_tend_bndflux', 3)
 
    !-----
    call prof_rapstart('cal_dyn_tend_interior', 3)
 
    !$omp parallel do private( &
    !$omp E, VOR,                          &
    !$omp u1_dt, u2_dt, Fx, Fy, LiftDelFlx )
    do ke = lmesh%NeS, lmesh%NeE
      !--
      call sparsemat_matmul(dx, u2_(:,ke), fx)
      call sparsemat_matmul(dy, u1_(:,ke), fy)
      call sparsemat_matmul(lift, lmesh%Fscale(:,ke) * del_flux_aux(:,ke,1), liftdelflx)
 
      vor(:) = ( &
          lmesh%Escale(:,ke,1,1) * fx(:) &
        - lmesh%Escale(:,ke,2,2) * fy(:) &
        + liftdelflx(:) ) / lmesh%Gsqrt(:,ke)
      
      e(:) = grav * ( hs_(:,ke) + h_(:,ke) )                        &
           + 0.5_rp * ( u1_(:,ke) * u_(:,ke) + u2_(:,ke) * v_(:,ke) )
      
      !-- h
      call sparsemat_matmul(dx, lmesh%Gsqrt(:,ke) * h_(:,ke) * u_(:,ke), fx)
      call sparsemat_matmul(dy, lmesh%Gsqrt(:,ke) * h_(:,ke) * v_(:,ke), fy)
      call sparsemat_matmul(lift, lmesh%Fscale(:,ke) * del_flux(:,ke,vars_h_id), liftdelflx)
 
      h_dt(:,ke) = - ( &
              lmesh%Escale(:,ke,1,1) * fx(:) &
            + lmesh%Escale(:,ke,2,2) * fy(:) &
            + liftdelflx(:) ) / lmesh%Gsqrt(:,ke)
      
      !-- u1
      call sparsemat_matmul(dx, e(:), fx)
      call sparsemat_matmul(lift, lmesh%Fscale(:,ke) * del_flux(:,ke,vars_u1_id), liftdelflx)
 
      u1_dt(:) = - (  &
            lmesh%Escale(:,ke,1,1) * fx(:)                                &
          + liftdelflx(:)                                                 &
          - lmesh%Gsqrt(:,ke) * ( coriolis(:,ke) + vor(:) ) * v_(:,ke)    )
      
      !-- u2
      call sparsemat_matmul(dy, e(:), fy)
      call sparsemat_matmul(lift, lmesh%Fscale(:,ke) * del_flux(:,ke,vars_u2_id), liftdelflx)
 
      u2_dt(:) = - (  &
            lmesh%Escale(:,ke,2,2) * fy(:)                                &
          + liftdelflx(:)                                                 &
          + lmesh%Gsqrt(:,ke) * ( coriolis(:,ke) + vor(:) ) * u_(:,ke)    )
    
      !--
      u_dt(:,ke) = lmesh%GIJ(:,ke,1,1) * u1_dt(:) + lmesh%GIJ(:,ke,1,2) * u2_dt(:) 
      v_dt(:,ke) = lmesh%GIJ(:,ke,2,1) * u1_dt(:) + lmesh%GIJ(:,ke,2,2) * u2_dt(:) 
    end do
 
    call prof_rapend('cal_dyn_tend_interior', 3)
 
    return