yinterpolate Subroutine

public subroutine yinterpolate(input, output, ks, kf)

Uses

  • proc~~yinterpolate~~UsesGraph proc~yinterpolate modinlet::yinterpolate module~modglobal modglobal proc~yinterpolate->module~modglobal

Arguments

Type IntentOptional Attributes Name
real, intent(in), dimension(jbdum:jedum,ks:kf,1:nstore) :: input
real, intent(inout), dimension(jb :je ,ks:kf,1:nstore) :: output
integer, intent(in) :: ks
integer, intent(in) :: kf

Called by

proc~~yinterpolate~~CalledByGraph proc~yinterpolate modinlet::yinterpolate proc~readinletfile modinlet::readinletfile proc~readinletfile->proc~yinterpolate proc~inletgen modinlet::inletgen proc~inletgen->proc~readinletfile proc~inletgennotemp modinlet::inletgennotemp proc~inletgennotemp->proc~readinletfile proc~readinitfiles modstartup::readinitfiles proc~readinitfiles->proc~readinletfile proc~boundary modboundary::boundary proc~readinitfiles->proc~boundary proc~boundary->proc~inletgen proc~boundary->proc~inletgennotemp proc~startup modstartup::startup proc~startup->proc~readinitfiles program~dalesurban DALESURBAN program~dalesurban->proc~boundary program~dalesurban->proc~startup

Contents

Source Code

yinterpolate

Source Code

  subroutine yinterpolate(input,output,ks,kf)
  use modglobal, only : jb,je, nstore
  integer, intent(in) :: ks
  integer, intent(in) :: kf
  real, dimension(jbdum:jedum,ks:kf,1:nstore), intent(in)  :: input
  real, dimension(jb   :je   ,ks:kf,1:nstore), intent(inout) :: output
  integer j

    do j=jb,je
!      if (np==0 .and. yloclowf(j)==)
      output(j,:,:) = input(yloclowf(j),:,:) + (input(ylocupf(j),:,:) - input(yloclowf(j),:,:)) / (yfdum(ylocupf(j)) - yfdum(yloclowf(j))) * (yf(j) - yfdum(yloclowf(j)))
    end do
  end subroutine yinterpolate