INTEGER maxcube
PARAMETER( maxcube = 1 )
IF ( nipass.GT.maxcube ) STOP 'maxcube needs to be = 3'
passkey = ipass + (k-1) *maxcube
& + (igadkey-1)*maxcube*Nr