rfmtftoc.f90

Go to the documentation of this file.

! Copyright (C) 2016 J. K. Dewhurst, S. Sharma and E. K. U. Gross.
! This file is distributed under the terms of the GNU General Public License.
! See the file COPYING for license details.

pure subroutine rfmtftoc(nrc,nrci,rfmt,rfcmt)
use modmain
implicit none
! arguments
integer, intent(in) :: nrc,nrci
real(8), intent(in) :: rfmt(*)
real(8), intent(out) :: rfcmt(*)
! local variables
integer irc,i,j,k
i=1
j=1
k=lmmaxi*lradstp
do irc=1,nrci
  rfcmt(i:i+lmmaxi-1)=rfmt(j:j+lmmaxi-1)
  i=i+lmmaxi
  j=j+k
end do
j=j+(lradstp-1)*(lmmaxo-lmmaxi)
k=lmmaxo*lradstp
do irc=nrci+1,nrc
  rfcmt(i:i+lmmaxo-1)=rfmt(j:j+lmmaxo-1)
  i=i+lmmaxo
  j=j+k
end do
end subroutine