SetCoriolis_betaplane_LinearShallowWater2D_t Subroutine

  subroutine SetCoriolis_betaplane_LinearShallowWater2D_t(this,f0,beta)
    implicit none
    class(LinearShallowWater2D_t),intent(inout) :: this
    real(prec),intent(in) :: f0
    real(prec),intent(in) :: beta
    ! Local
    integer :: iel
    integer :: i
    integer :: j
    real(prec) :: y

    do concurrent(i=1:this%solution%N+1,j=1:this%solution%N+1, &
                  iel=1:this%mesh%nElem)
      y = this%geometry%x%interior(i,j,iel,1,2)
      this%fCori%interior(i,j,iel,1) = f0+beta*y
    enddo
    call this%fCori%UpdateDevice()

  endsubroutine SetCoriolis_betaplane_LinearShallowWater2D_t