! This FUNCTION is used to return the acceleration of incident waves ! of point IWaveKin in the xi- (KDirection=1), yi- (KDirection=2), or ! zi- (KDirection=3) direction, respectively, at time ZTime to the ! calling program.-! This FUNCTION