wt = range(1880,2011)
print "yeardata = "+repr(wt)+";"

t = open( "RadF.txt" ).readlines()
cols = [4,7,1,9,10,2,3,5,6,8]
print "// "+repr(wt)
print "gisshead = ['Solar','Volcanoes','GHG(mixed)','Aerosol(dir)','Aerosol(ind)','Ozone','H2O(strato)','LandUse','SnowAlbedo','BlackCarbon'];"
wt = [[round(float(t[i+4].split()[j]),4) for i in range(131)] for j in cols]
print "gissdata = ["
for i in range(len(wt)):
  print "  "+repr(wt[i])+","
print "];"

t = open( "RCP85_MIDYEAR_RADFORCING.DAT" ).readlines()
cols = [3,2,5,41,48,49,50,51,52,53]
print "// "+repr(wt)
print "cmiphead = ['Solar','Volcanoes','GHG(mixed)','Aerosol(dir)','Aerosol(ind)','Ozone(strato)','Ozone(tropo)','H2O(strato)','LandUse','BlackCarbon(snow)'];"
wt = [[round(float(t[i+175].split()[j]),4) for i in range(131)] for j in cols]
print "cmipdata = ["
for i in range(len(wt)):
  print "  "+repr(wt[i])+","
print "];"

#t = open( "otto.txt" ).readlines()
#print "ottohead = ['Solar(NULL)','Volcanoes(NULL)','All'];"
#w = [round(float(t[i].split()[1]),4) for i in range(131)]
#print "ottodata = ["
#for i in range(2):
#  print "  "+repr([0.0 for i in range(131)])+","
#print "  "+repr(w)+","
#print "];"

t = open( "Temp.txt" ).readlines()
wt = [float(t[i].split()[1]) for i in range(131)]
print "tempdata = "+repr(wt)+";"

#print "invarhead = ['+ENSO (MEI+4m)','-AMO (Kaplan v2) 121m smooth','-AMO (Kaplan v2): Warning - see notes','-AMO (Kaplan v2+3m)','-AMO (Kaplan v2+6m)','-AMO (Trenberth/ERSST)','-AMO (Trenberth/ERSST+3m)','-AMO (Trenberth/ERSST+6m)','-AMO (Trenberth/ERSST+12m)','-PDO (Kaplan)'];"
print "invarhead = ['+ENSO (MEI+4m)'];"
print "invardata = ["
t = open( "mei.lag4" ).readlines()
wt = [float(t[i].split()[1]) for i in range(131)]
print "  "+repr(wt)+","
#t = open( "amo-ma.txt" ).readlines()
#wt = [float(t[i].split()[1]) for i in range(131)]
#print "  "+repr(wt)+","
#t = open( "amo.txt" ).readlines()
#wt = [float(t[i].split()[1]) for i in range(131)]
#print "  "+repr(wt)+","
#t = open( "amo-lag3.txt" ).readlines()
#wt = [float(t[i].split()[1]) for i in range(131)]
#print "  "+repr(wt)+","
#t = open( "amo-lag6.txt" ).readlines()
#wt = [float(t[i].split()[1]) for i in range(131)]
#print "  "+repr(wt)+","
#t = open( "amo2.txt" ).readlines()
#wt = [float(t[i].split()[1]) for i in range(131)]
#print "  "+repr(wt)+","
#t = open( "amo2-lag3.txt" ).readlines()
#wt = [float(t[i].split()[1]) for i in range(131)]
#print "  "+repr(wt)+","
#t = open( "amo2-lag6.txt" ).readlines()
#wt = [float(t[i].split()[1]) for i in range(131)]
#print "  "+repr(wt)+","
#t = open( "amo2-lag12.txt" ).readlines()
#wt = [float(t[i].split()[1]) for i in range(131)]
#print "  "+repr(wt)+","
#t = open( "pdo.txt" ).readlines()
#wt = [float(t[i].split()[1]) for i in range(131)]
#print "  "+repr(wt)+","
print "];"