- if (diff.gt.0) then
- do i=1,diff
- ihpb(nhpb+i)=ihpb(nss+i)
- jhpb(nhpb+i)=jhpb(nss+i)
- forcon(nhpb+i)=forcon(nss+i)
- dhpb(nhpb+i)=dhpb(nss+i)
- enddo
- else if (diff.lt.0) then
- do i=diff,-1
- ihpb(nss+i)=ihpb(nhpb+i)
- jhpb(nss+i)=jhpb(nhpb+i)
- forcon(nss+i)=forcon(nhpb+i)
- dhpb(nss+i)=dhpb(nhpb+i)
- enddo
- endif
-
- nhpb=nhpb+diff