import requests
import re
import numpy as np

airmasses = [1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.1,2.2,2.3]
pwv_values = [15.0]
for airmass in airmasses:
    for pwv_value in pwv_values:
        url = 'http://www.eso.org/observing/etc/bin/simu/skycalc'
        payload = {'INS.NAME':'SKYCALC','INS.MODE':'swspectr','POSTFILE.FLAG':0,'COORD.DAY':8,'COORD.MONTH':5,'COORD.YEAR':2017,'COORD.UT.HOUR':0,'COORD.UT.MIN':0,'COORD.UT.SEC':0,'COORD.RA':60,'COORD.DEC':0,'SKYMODEL.TARGET.ALT':90,'SKYMODEL.TARGET.AIRMASS': airmass,'SKYMODEL.SEASON':2,'SKYMODEL.TIME':2,'SKYMODEL.PWV.MODE':'pwv','SKYMODEL.PWV':pwv_value, 'SKYMODEL.MSOLFLUX':'130.00','SKYMODEL.INCL.MOON':'N', 'SKYMODEL.MOON.SUN.SEP':90,'SKYMODEL.MOON.TARGET.SEP':45,'SKYMODEL.MOON.ALT':45,'SKYMODEL.MOON.EARTH.DIST':1,
                'SKYMODEL.INCL.STARLIGHT':'N','SKYMODEL.INCL.ZODIACAL':'N','SKYMODEL.ECL.LON':135,'SKYMODEL.ECL.LAT':90,'SKYMODEL.INCL.MOLEC.EMIS.LOWER.ATM':'Y',
                'SKYMODEL.INCL.MOLEC.EMIS.UPPER.ATM':'Y','SKYMODEL.INCL.AIRGLOW':'Y','SKYMODEL.INCL.THERMAL':'N','SKYMODEL.THERMAL.T1':0,'SKYMODEL.THERMAL.E1':0,
                'SKYMODEL.THERMAL.T2':0,'SKYMODEL.THERMAL.E2':0,'SKYMODEL.THERMAL.T3':0,'SKYMODEL.THERMAL.E3':0,'SKYMODEL.WAVELENGTH.MIN':305,'SKYMODEL.WAVELENGTH.MAX':2510,
                'SKYMODEL.WAVELENGTH.GRID.MODE':'fixed_spectral_resolution','SKYMODEL.WAVELENGTH.DELTA':0.1,'SKYMODEL.WAVELENGTH.RESOLUTION':20000,'SKYMODEL.LSF.KERNEL.TYPE':'none',
                'SKYMODEL.LSF.GAUSS.FWHM':20,'SKYMODEL.LSF.BOXCAR.WIDTH':10,'SKYCALC.RAD.PLOT.FLAG':1,'SKYCALC.DETAILEDRAD.PLOT.FLAG':1,'SKYCALC.TRANS.PLOT.FLAG':1,
                'SKYCALC.DETAILEDTRANS.PLOT.FLAG':1,'SKYCALC.MAG.FLAG':1,'SKYCALC.LSF.PLOT.FLAG':1,}

        r = requests.post(url, data=payload)


        #import json
        #r = requests.post(url, data=json.dumps(payload))

        # Response, status etc
        result = re.search(r'/observing/(.*)/skytable.fits', r.text)
        url = result.group(0)
        filename = "paranal_sky_model_%03d_%02d.fits" % (10*pwv_value,10*airmass)
        print 'wget -O %s -c http://www.eso.org%s' % (filename,url)