Teq210_25 = Tatm_sac * em210_25
Teq210_50 = Tatm_sac * em210_50
Teq210_75 = Tatm_sac * em210_75
clf()
xlim(0,13)
plot(month, Teq210_25, 'ro-')
plot(month, Teq210_50, 'go-')
plot(month, Teq210_75, 'bo-')
#### Concordia
atm_DomeC_summer = {'Name':'Atm Dome C ', 'T':200, 'e':0.025, 'trans':1}
atm_DomeC_winter = {'Name':'Atm Dome C ', 'T':260, 'e':0.06, 'trans':1}
val_DomeC_summer = loading.give_NET(atm=atm_DomeC_summer, verbose=True)
net_150_DomeC_summer = val_DomeC_summer[3][0]
net_220_DomeC_summer = val_DomeC_summer[3][1]
val_DomeC_winter = loading.give_NET(atm=atm_DomeC_winter, verbose=True)
net_150_DomeC_winter = val_DomeC_winter[3][0]
net_220_DomeC_winter = val_DomeC_winter[3][1]
#### Chorillos
net_150_chorillos_25 = np.zeros(len(month))
net_150_chorillos_50 = np.zeros(len(month))
net_150_chorillos_75 = np.zeros(len(month))
net_220_chorillos_25 = np.zeros(len(month))
net_220_chorillos_50 = np.zeros(len(month))
net_220_chorillos_75 = np.zeros(len(month))
xlabel('Month')
ylabel(r'Atmospheric Brightness Temperature (K)')
plot(month, em220_ch * TAvchorillo, 'b', lw=3, label ='Chorillos - 220 GHz')
plot(month, em150_ch * TAvchorillo, 'r', lw=3, label ='Chorillos - 150 GHz')
plot(month, em220_chaj * TAvchorillo, 'b--', lw=3, label ='Chajnantor (1998 ext from 857 GHz) - 220 GHz')
plot(month, em150_chaj * TAvchorillo, 'r--', lw=3, label ='Chajnantor (1998 ext from 857 GHz) - 150 GHz')
ylim(0,110)
legend()
grid()
#savefig('atm_brightnesscho-chaj.png')
from Sites import loading
theatm = {'Name':'Dome C ', 'T':Tdomec[0], 'e':em150_dc[0], 'trans':1}
vals = loading.give_NET(atm=theatm, verbose=True, bol_efficiency=0.8, freqGHz=150.)
theatm = {'Name':'Dome C ', 'T':Tdomec[0], 'e':em220_dc[0], 'trans':1}
vals = loading.give_NET(atm=theatm, verbose=True, bol_efficiency=0.8, freqGHz=220.)
fraction = 0.76
net150_dc = np.zeros(12)
net220_dc = np.zeros(12)
net150_ch = np.zeros(12)
net220_ch = np.zeros(12)
margin = 0.8
for i in np.arange(12):
theatm = {'Name':'Atm Chorillos ', 'T':TAvchorillo[i], 'e':em150_ch[i], 'trans':1}
