subroutine SAVETQBANKS

      implicit none
      
#include "skparm.h"
#include "skhead.h"
#include "sktq.h"
#include "zbs.h"

      common/tqtempbank/
     +     ldat,ndata_save,ibuf_save,ibufz_save,
     +     ndataa_save, ibufa_save, ibufza_save

      integer ldat
      integer ndata_save, ibuf_save(10)
      integer ibufz_save(3,MAXPM*30)
      integer ndataa_save, ibufa_save(10)
      integer ibufza_save(3,MAXPMA*30)

      if(sk_geometry.ge.4) then
C  Store TQREAL and TQAREAL
         CALL KZGET1('TQREAL',0,ndata_save,ibuf_save)
         CALL KZLDAT('TQREAL',ldat)
         CALL LBFCOPY(izbs(ldat+ndata_save+1),ibufz_save,3*ibuf_save(1))
         CALL KZGET1('TQAREAL',0,ndataa_save,ibufa_save)
         CALL KZLDAT('TQAREAL',ldat)
         CALL LBFCOPY(izbs(ldat+ndataa_save+1),ibufza_save,3*ibufa_save(1))

      else
         print *,'SAVETQBANKS: Only SK4 is supported!'
      endif

      return
      end

      
      subroutine RESTORETQBANKS

      implicit none
      
#include "skparm.h"
#include "skhead.h"
#include "sktq.h"
#include "zbs.h"

      common/tqtempbank/
     +     ldat,ndata_save,ibuf_save,ibufz_save,
     +     ndataa_save, ibufa_save, ibufza_save

      integer ldat
      integer ndata_save, ibuf_save(10)
      integer ibufz_save(3,MAXPM*30)
      integer ndataa_save, ibufa_save(10)
      integer ibufza_save(3,MAXPMA*30)

      integer   LUN
      parameter(LUN=10) ! no meaning

      if(sk_geometry.ge.4) then
C  Restore TQREAL and TQAREAL
         CALL KZREP1('TQREAL',0,ndata_save,ibuf_save)
         CALL KZLDAT('TQREAL',ldat)
         CALL LBFCOPY(ibufz_save,izbs(ldat+ndata_save+1),3*ibuf_save(1))
         CALL KZREP1('TQAREAL',0,ndataa_save,ibufa_save)
         CALL KZLDAT('TQAREAL',ldat)
         CALL LBFCOPY(ibufza_save,izbs(ldat+ndataa_save+1),3*ibufa_save(1))

         call skread(-LUN)

      else
         print *,'RESTORETQBANKS: Only SK4 is supported!'
      endif

      return
      end