# 
# Greenwich Mean Sidereal Time
# at 0H GMT for N days 

N=2000;

def jd0(dd,mm,yyyy):
    if mm>2:
        y=yyyy
        m=mm-3
    else:
        y=yyyy-1
        m=mm+9
    jd=1721103.5+int(365.25*y)+int(30.6*m+0.5)+dd
    return(jd)

# IAU formula (1/10 seconds)
def gmst_IAU(dia):
    t=dia/36525.0
    gmst=24110.54841+8640184.812866*t+0.093104*t*t-0.0000062*t*t*t
    gmst=int(gmst*10.0)
    gmst=gmst%(24*3600*10)
    return(gmst)

# initial date
jd_zero=jd0(30,11,2016)

d0=jd_zero-2451545.0

for i in range(N):
    dia=d0+i
    gmst=gmst_IAU(dia)
    print ( " %8d," % gmst),
    if 0==(i+1)%8 : print " \n ",