apply_coverage_factors Subroutine

    subroutine apply_coverage_factors
        !! Apply coverage factors to volumes
        !! author: J. Morris, CCFE, Culham Science Centre
        !! Apply coverage factors to volumes
        !
        ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

        use build_variables, only: blareaob, blarea, blareaib, shareaib, shareaob, &
          sharea
        use fwbs_variables, only: fdiv, fhcd, volblkto, volblkti, volblkt, fvolsi, &
          fvolso, volshld, vdewin, fvoldw
        use physics_variables, only: idivrt

        implicit none

        ! Apply blanket coverage factors
        if (idivrt == 2) then
          ! double null configuration
          blareaob = blarea*(1.0D0-2.0D0*fdiv-fhcd) - blareaib
        else
          ! single null configuration
          blareaob = blarea*(1.0D0-fdiv-fhcd) - blareaib
        end if

        blarea = blareaib + blareaob

        volblkto = volblkt*(1.0D0-fdiv-fhcd) - volblkti
        volblkt = volblkti + volblkto

        ! Apply shield coverage factors
        shareaib = fvolsi*shareaib
        shareaob = fvolso*shareaob
        sharea = shareaib + shareaob

        volshldi = fvolsi*volshldi
        volshldo = fvolso*volshldo
        volshld = volshldi + volshldo

        ! Apply vacuum vessel coverage factor
        ! moved from dshaped_* and elliptical_* to keep coverage factor
        ! changes in the same location.
        vdewin = fvoldw*vdewin

    end subroutine apply_coverage_factors