setgradientboundarycondition_DGModel2D_t Subroutine

  subroutine setgradientboundarycondition_DGModel2D_t(this)
    !! Apply registered boundary conditions for the solution gradient.
    !! Each boundary condition method loops over its own
    !! boundary faces.
    implicit none
    class(DGModel2D_t),intent(inout) :: this
    ! Local
    type(BoundaryCondition),pointer :: bc
    procedure(SELF_bcMethod),pointer :: apply_bc

    bc => this%parabolicBCs%head
    do while(associated(bc))
      apply_bc => bc%bcMethod
      call apply_bc(bc,this)
      bc => bc%next
    enddo

  endsubroutine setgradientboundarycondition_DGModel2D_t