#!/usr/bin/env python3

ANG = 1.e-10

ucoord = 0.0
vcoord = 0.0
vis2data = 0.2
vis2err = 0.01
eff_wave = 6625*ANG
eff_band = 100*ANG

f = open("Vis.dat", "w")
f.write("# jd ucoord vcoord eff_wave eff_band vis2data vis2err\n")

t1 = 2.5+0.25
t2 = t1
dt = 0.01
u1 = 0.0
u2 = 330.0
du = 1.0
v1 = u1
v2 = u2
dv = du
eps = 1.0e-8

t = t1
while t < t2+eps:
    u = u1
    while u < u2:
        u += du
        v = v1
        while v < v2:
            v += dv
            f.write("%.8f  %.8e  %.8e  %.8e  %.8e  %.8f  %.8f\n" % (t, u, v, eff_wave, eff_band, vis2data, vis2err))
    t += dt

f.close()


