setCpTarget.f90

Go to the documentation of this file.

subroutine setcptargets(cptarget, npts, sps_in)
    use constants
    use blockpointers
    use flowvarrefstate
    use inputtimespectral
    use communication
    use inputphysics
    use utils, only: setpointers
    implicit none
    !
    !      Arguments.
    !
    integer(kind=intType), intent(in) :: npts, sps_in
    real(kind=realtype), intent(in) :: cptarget(npts)
    !
    !      Local variables.
    !
    integer(kind=intType) :: mm, nn, i, j, ii, sps
    !integer(kind=intType) :: iBeg, iEnd, jBeg, jEnd
 
    sps = sps_in
 
    ii = 0
    domains: do nn = 1, ndom
        call setpointers(nn, 1_inttype, sps)
 
        ! Loop over the number of viscous boundary subfaces of this block.
        ! According to preprocessing/viscSubfaceInfo, visc bocos are numbered
        ! before other bocos. Therefore, mm_nViscBocos == mm_nBocos
        bocos: do mm = 1, nbocos
            if (bctype(mm) == eulerwall .or. bctype(mm) == nswalladiabatic .or. &
                bctype(mm) == nswallisothermal) then
                do j = bcdata(mm)%jnBeg, bcdata(mm)%jnEnd
                    do i = bcdata(mm)%inBeg, bcdata(mm)%inEnd
                        ii = ii + 1
                        bcdata(mm)%CpTarget(i, j) = cptarget(ii)
                    end do
                end do
            end if
 
        end do bocos
    end do domains
 
end subroutine setcptargets
 
subroutine getcptargets(cptarget, npts, sps)
 
    use constants
    use blockpointers, only: ndom, nbocos, bcdata, bctype
!  use flowVarRefState, only : TRef
    use utils, only: setpointers
    implicit none
 
    ! Input Variables
    integer(kind=intType), intent(in) :: npts, sps
    real(kind=realtype), intent(out) :: cptarget(npts)
 
    ! Local Variables
    integer(kind=intType) :: mm, nn, i, j, ii
    integer(kind=intType) :: iBeg, iEnd, jBeg, jEnd
 
    ii = 0
    domains: do nn = 1, ndom
        call setpointers(nn, 1_inttype, sps)
        ! Loop over the number of viscous boundary subfaces of this block.
        bocos: do mm = 1, nbocos
            wall: if (bctype(mm) == eulerwall .or. bctype(mm) == nswalladiabatic .or. &
                      bctype(mm) == nswallisothermal) then
                do j = bcdata(mm)%jnBeg, bcdata(mm)%jnEnd
                    do i = bcdata(mm)%inBeg, bcdata(mm)%inEnd
                        ii = ii + 1
                        cptarget(ii) = bcdata(mm)%CpTarget(i, j)
                    end do
                end do
            end if wall
        end do bocos
    end do domains
end subroutine getcptargets