Compare the VER computed by the TIE-GCM model with MSIS/IRI nightglow

20 Mar 2009 (f10.7a = 70, Dst = 0, Kp = 1.3)

fig = figure(figsize=(8,3))

################################# RED #############################

mx = max(A_gcm_red.max(),A_mine_red.max())
mn = 1e-4*mx

subplot(2,2,1)
pcolormesh(T,ALT,A_gcm_red, cmap='plasma', norm=matplotlib.colors.LogNorm(vmin=mn,vmax=mx))
cb = colorbar(aspect=10)
cb.set_label('VER [ph/cm^3/s]',rotation=-90, labelpad=10)
axis('tight')
title('Scott Red')
ylabel('Altitude [km]')
gca().set_xticklabels([])

subplot(2,2,3)
pcolormesh(T,ALT,A_mine_red, cmap='plasma', norm=matplotlib.colors.LogNorm(vmin=mn,vmax=mx))
cb = colorbar(aspect=10)
cb.set_label('VER [ph/cm^3/s]',rotation=-90, labelpad=10)
axis('tight')
title('Brian Red')
xlabel('Time [sec]')
ylabel('Altitude [km]')

################################# GREEN #############################

mx = max(A_gcm_green.max(),A_mine_green.max())
mn = 1e-4*mx

subplot(2,2,2)
pcolormesh(T,ALT,A_gcm_green, cmap='viridis', norm=matplotlib.colors.LogNorm(vmin=mn,vmax=mx))
cb = colorbar(aspect=10)
cb.set_label('VER [ph/cm^3/s]',rotation=-90, labelpad=10)
axis('tight')
title('Scott Green')
ylabel('Altitude [km]')
gca().set_xticklabels([])

subplot(2,2,4)
pcolormesh(T,ALT,A_mine_green, cmap='viridis', norm=matplotlib.colors.LogNorm(vmin=mn,vmax=mx))
cb = colorbar(aspect=10)
cb.set_label('VER [ph/cm^3/s]',rotation=-90, labelpad=10)
axis('tight')
title('Brian Green')
xlabel('Time [sec]')
ylabel('Altitude [km]')

tight_layout()

Integrate vertically to get column brightness [Rayleigh]

DALT = zeros((M,N))
DALT[:-1,:] = diff(ALT,axis=0)*1e3 # convert to m
DALT[-1,:] = DALT[-2,:]

figure(figsize=(8,3))
################################# RED #######################################
subplot(1,2,1)
B_gcm_red =  sum(A_gcm_red*1e6*DALT, axis=0)*1e-10 # R
B_mine_red = sum(A_mine_red*1e6*DALT, axis=0)*1e-10 # R
plot(tvec, B_gcm_red,  'r-',  label = 'Scott')
plot(tvec, B_mine_red, 'r--', label = 'Brian',dashes=(2,2))
legend(loc='best')
xlabel('Time [sec]')
ylabel('Vertical Column Brightness [R]')
gca().set_yscale('log')
grid()

################################# GREEN #######################################
subplot(1,2,2)
B_gcm_green =  sum(A_gcm_green*1e6*DALT, axis=0)*1e-10 # R
B_mine_green = sum(A_mine_green*1e6*DALT, axis=0)*1e-10 # R
plot(tvec, B_gcm_green,  'g-',  label = 'Scott')
plot(tvec, B_mine_green, 'g--', label = 'Brian',dashes=(2,2))
legend(loc='best')
xlabel('Time [sec]')
ylabel('Vertical Column Brightness [R]')
gca().set_yscale('log')
grid()