subroutine where(x,y,fact) common/init/ilun,factor,imode,ipenst call plot(x,y,4) fact=factor return end