*-- Author : S.Burke / J.V. Morris SUBROUTINE FKSAVE(IFLAG,JFIRST,JLAST) ********************************************************************** * * * Save/restore smoothed vector (+ residuals) from JFIRST to * * JLAST (inclusive) * * * * Save if IFLAG is >= 0; restore otherwise * * * * Note that there is NO error checking. * * * * Also note that this does not save the entire status (i.e. * * projected and filtered vectors). * * * ********************************************************************** *KEEP,FKNPL. * * Per-track values can go in H1WORK; note that LTRUE and LFIRST must * be set at least per event. * * This is about 36k words long; the remaining common blocks are * about 3.6k in total. Some of this could be in /H1WORK/, but the * blocks would have to be reorganised. * * /FKPROJ/ * /FKFILT/ * /FKSMTH/ * /FKINT/ * /FKRSID/ * /FKTRUE/ * /FKDBG/ *KEEP,FKFLAG. *KEEP,FKSMTH. *KEEP,FKRSID. *KEEP,FKTRUE. *KEND. ********************************************************************** ********************************************************************** CALL UCOPY(SSMT(1,JMIN),SSAVE,10*JDIFF) CALL UCOPY(CSMT(1,1,JMIN),CSAVE,50*JDIFF) CALL UCOPY(RSMT(1,JMIN),SRSAVE,4*JDIFF) CALL UCOPY(CRSMT(1,1,JMIN),CRSAVE,8*JDIFF) CALL UCOPY(CHISMT(JMIN),CHSAVE,2*JDIFF) CALL UCOPY(RTRUE(1,JMIN),RTSAVE,10*JDIFF) CALL UCOPY(CHITRU(JMIN),CTSAVE,2*JDIFF) CALL UCOPY(SSAVE,SSMT(1,JMIN),10*JDIFF) CALL UCOPY(CSAVE,CSMT(1,1,JMIN),50*JDIFF) CALL UCOPY(SRSAVE,RSMT(1,JMIN),4*JDIFF) CALL UCOPY(CRSAVE,CRSMT(1,1,JMIN),8*JDIFF) CALL UCOPY(CHSAVE,CHISMT(JMIN),2*JDIFF) CALL UCOPY(RTSAVE,RTRUE(1,JMIN),10*JDIFF) CALL UCOPY(CTSAVE,CHITRU(JMIN),2*JDIFF) *