copyright (C) 2001 MSC-RPN COMM %%%MC2%%% ***s/r inv_posiz -- *subroutine inv_posiz (p,h,href,hint,ktop,kbot,ng,nk,nkref) 5,1 implicit none * integer ng,nk,nkref real p(ng,nk,2),href(ng,nkref),hint(ng,nk), $ h(ng,0:nkref+2),ktop(ng,nk),kbot(ng,nk) * integer i,k real dum,w1(ng,nk,2),w2(ng,0:nkref+2) * *---------------------------------------------------------------------- * do k=1,nkref/2 do i=1,ng dum = href (i,k) href(i,k) = href (i,nkref-k+1) href(i,nkref-k+1) = dum end do end do * call posiz3
( p,h,href,hint,ktop,kbot,ng,nk,nkref ) * do k=1,nkref/2 do i=1,ng dum = href (i,k) href(i,k) = href (i,nkref-k+1) href(i,nkref-k+1) = dum end do end do * *---------------------------------------------------------------------- * return end