subroutine oldwkday(ij,iy,dayawk)
        implicit none  
        integer, intent(in):: ij,iy   !iy at least 4 digits
        integer, intent(out):: dayawk
!	returns  day of week (mon=1, sun=7) for supplied day of yr and year 
	logical leap
	integer centdy,leapad,ijday(12)
	data ijday /0,31,59,90,120,151,181,212,243,273,304,334/
	leapad=(iy-1901)/4 ! add day after each leap year
	leap=(mod(iy,4).eq.0).and.(mod(iy,100).ne.0) ! leap: wrong in 2100
	if(mod(iy,400).eq.0) leap=.true.
	centdy=(iy-1900)*365+leapad+ij  !days since 1900
	dayawk=mod(centdy-1,7)+1 !
	return
	end