Type | Intent | Optional | 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 |
subroutine yinterpolateh(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, jj
do j=jb,je
! if (np==0 .and. yloclowf(j)==)
output(j,:,:) = input(yloclowh(j),:,:) + (input(ylocuph(j),:,:) - input(yloclowh(j),:,:)) / (yhdum(ylocuph(j)) - yhdum(yloclowh(j))) * (yh(j) - yhdum(yloclowh(j)))
end do
end subroutine yinterpolateh