lat_c = root_2000.variables['lat_centre'][:]
lon_c = root_2000.variables['lon_centre'][:]
size = len(lat_c)*len(lon_c)
m = Basemap(projection='cyl',llcrnrlat=lat_e[0],urcrnrlat=lat_e[-1],llcrnrlon=lon_e[0],urcrnrlon=lon_e[-1],resolution='c')
print 'b'
(diurnal_mag_2000,diurnal_ph_2000,seasonal_mag_2000,seasonal_ph_2000,mean_2000,p1_2000,p5_2000,p25_2000,p50_2000,p75_2000,p95_2000,p99_2000,diurnal_ave_waveform_2000,seasonal_waveform_2000,full_ave_waveform_2000,pc_var_daily_2000,pc_var_seasonal_2000,pc_var_full_2000,pc_var_noise_2000,total_var_2000,
diurnal_mag_spring_2000,diurnal_ph_spring_2000,mean_spring_2000,p1_spring_2000,p5_spring_2000,p25_spring_2000,p50_spring_2000,p75_spring_2000,p95_spring_2000,p99_spring_2000,diurnal_waveform_spring_2000,
diurnal_mag_summer_2000,diurnal_ph_summer_2000,mean_summer_2000,p1_summer_2000,p5_summer_2000,p25_summer_2000,p50_summer_2000,p75_summer_2000,p95_summer_2000,p99_summer_2000,diurnal_waveform_summer_2000,
diurnal_mag_autumn_2000,diurnal_ph_autumn_2000,mean_autumn_2000,p1_autumn_2000,p5_autumn_2000,p25_autumn_2000,p50_autumn_2000,p75_autumn_2000,p95_autumn_2000,p99_autumn_2000,diurnal_waveform_autumn_2000,
diurnal_mag_winter_2000,diurnal_ph_winter_2000,mean_winter_2000,p1_winter_2000,p5_winter_2000,p25_winter_2000,p50_winter_2000,p75_winter_2000,p95_winter_2000,p99_winter_2000,diurnal_waveform_winter_2000,
seasonal_mag_day_2000,seasonal_ph_day_2000,mean_day_2000,p1_day_2000,p5_day_2000,p25_day_2000,p50_day_2000,p75_day_2000,p95_day_2000,p99_day_2000,seasonal_waveform_day_2000,
seasonal_mag_night_2000,seasonal_ph_night_2000,mean_night_2000,p1_night_2000,p5_night_2000,p25_night_2000,p50_night_2000,p75_night_2000,p95_night_2000,p99_night_2000,seasonal_waveform_night_2000,
daily_h3_mag_2000,daily_h2_mag_2000,daily_h1_mag_2000,daily_mag_2000,annual_h3_mag_2000,annual_h2_mag_2000,annual_h1_mag_2000,annual_mag_2000) = modules.get_periodic_specific('na',f2000)
if model != 'GISSE2R':
(diurnal_mag_2100,diurnal_ph_2100,seasonal_mag_2100,seasonal_ph_2100,mean_2100,p1_2100,p5_2100,p25_2100,p50_2100,p75_2100,p95_2100,p99_2100,diurnal_ave_waveform_2100,seasonal_waveform_2100,full_ave_waveform_2100,pc_var_daily_2100,pc_var_seasonal_2100,pc_var_full_2100,pc_var_noise_2100,total_var_2100,
diurnal_mag_spring_2100,diurnal_ph_spring_2100,mean_spring_2100,p1_spring_2100,p5_spring_2100,p25_spring_2100,p50_spring_2100,p75_spring_2100,p95_spring_2100,p99_spring_2100,diurnal_waveform_spring_2100,
diurnal_mag_summer_2100,diurnal_ph_summer_2100,mean_summer_2100,p1_summer_2100,p5_summer_2100,p25_summer_2100,p50_summer_2100,p75_summer_2100,p95_summer_2100,p99_summer_2100,diurnal_waveform_summer_2100,
diurnal_mag_autumn_2100,diurnal_ph_autumn_2100,mean_autumn_2100,p1_autumn_2100,p5_autumn_2100,p25_autumn_2100,p50_autumn_2100,p75_autumn_2100,p95_autumn_2100,p99_autumn_2100,diurnal_waveform_autumn_2100,
diurnal_mag_winter_2100,diurnal_ph_winter_2100,mean_winter_2100,p1_winter_2100,p5_winter_2100,p25_winter_2100,p50_winter_2100,p75_winter_2100,p95_winter_2100,p99_winter_2100,diurnal_waveform_winter_2100,
seasonal_mag_day_2100,seasonal_ph_day_2100,mean_day_2100,p1_day_2100,p5_day_2100,p25_day_2100,p50_day_2100,p75_day_2100,p95_day_2100,p99_day_2100,seasonal_waveform_day_2100,
seasonal_mag_night_2100,seasonal_ph_night_2100,mean_night_2100,p1_night_2100,p5_night_2100,p25_night_2100,p50_night_2100,p75_night_2100,p95_night_2100,p99_night_2100,seasonal_waveform_night_2100,
daily_h3_mag_2100,daily_h2_mag_2100,daily_h1_mag_2100,daily_mag_2100,annual_h3_mag_2100,annual_h2_mag_2100,annual_h1_mag_2100,annual_mag_2100) = modules.get_periodic_specific('na',f2100)
print 'c'
if model != 'CMAM':
(diurnal_mag_2100e,diurnal_ph_2100e,seasonal_mag_2100e,seasonal_ph_2100e,mean_2100e,p1_2100e,p5_2100e,p25_2100e,p50_2100e,p75_2100e,p95_2100e,p99_2100e,diurnal_ave_waveform_2100e,seasonal_waveform_2100e,full_ave_waveform_2100e,pc_var_daily_2100e,pc_var_seasonal_2100e,pc_var_full_2100e,pc_var_noise_2100e,total_var_2100e,
p99_summer_2000, diurnal_waveform_summer_2000, diurnal_mag_autumn_2000,
diurnal_ph_autumn_2000, mean_autumn_2000, p1_autumn_2000, p5_autumn_2000,
p25_autumn_2000, p50_autumn_2000, p75_autumn_2000, p95_autumn_2000,
p99_autumn_2000, diurnal_waveform_autumn_2000, diurnal_mag_winter_2000,
diurnal_ph_winter_2000, mean_winter_2000, p1_winter_2000, p5_winter_2000,
p25_winter_2000, p50_winter_2000, p75_winter_2000, p95_winter_2000,
p99_winter_2000, diurnal_waveform_winter_2000, seasonal_mag_day_2000,
seasonal_ph_day_2000, mean_day_2000, p1_day_2000, p5_day_2000,
p25_day_2000, p50_day_2000, p75_day_2000, p95_day_2000, p99_day_2000,
seasonal_waveform_day_2000, seasonal_mag_night_2000,
seasonal_ph_night_2000, mean_night_2000, p1_night_2000, p5_night_2000,
p25_night_2000, p50_night_2000, p75_night_2000, p95_night_2000,
p99_night_2000, seasonal_waveform_night_2000, daily_h3_mag_2000,
daily_h2_mag_2000, daily_h1_mag_2000, daily_mag_2000, annual_h3_mag_2000,
annual_h2_mag_2000, annual_h1_mag_2000,
annual_mag_2000) = modules.get_periodic_specific('na', f2000)
if model != 'GISSE2R':
(diurnal_mag_2100, diurnal_ph_2100, seasonal_mag_2100,
seasonal_ph_2100, mean_2100, p1_2100, p5_2100, p25_2100, p50_2100,
p75_2100, p95_2100, p99_2100, diurnal_ave_waveform_2100,
seasonal_waveform_2100, full_ave_waveform_2100, pc_var_daily_2100,
pc_var_seasonal_2100, pc_var_full_2100, pc_var_noise_2100,
total_var_2100, diurnal_mag_spring_2100, diurnal_ph_spring_2100,
mean_spring_2100, p1_spring_2100, p5_spring_2100, p25_spring_2100,
p50_spring_2100, p75_spring_2100, p95_spring_2100, p99_spring_2100,
diurnal_waveform_spring_2100, diurnal_mag_summer_2100,
diurnal_ph_summer_2100, mean_summer_2100, p1_summer_2100,
p5_summer_2100, p25_summer_2100, p50_summer_2100, p75_summer_2100,
p95_summer_2100, p99_summer_2100, diurnal_waveform_summer_2100,
diurnal_mag_autumn_2100, diurnal_ph_autumn_2100, mean_autumn_2100,
,obs_seasonal_mag_night,obs_seasonal_phase,obs_seasonal_phase_day,obs_seasonal_phase_night,obs_seasonal_phase_min\
,obs_seasonal_phase_min_day,obs_seasonal_phase_min_night,obs_daily_ff,obs_daily_ff_spring,obs_daily_ff_summer,obs_daily_ff_autumn\
,obs_daily_ff_winter,obs_seasonal_ff,obs_seasonal_ff_day,obs_seasonal_ff_night,obs_daily_waveform,obs_daily_waveform_spring\
,obs_daily_waveform_summer,obs_daily_waveform_autumn,obs_daily_waveform_winter,obs_seasonal_waveform,obs_seasonal_waveform_day\
,obs_seasonal_waveform_night,obs_full_waveform,obs_periodic_variance_daily,obs_periodic_variance_daily_spring,obs_periodic_variance_daily_summer\
,obs_periodic_variance_daily_autumn,obs_periodic_variance_daily_winter,obs_periodic_variance_seasonal,obs_periodic_variance_seasonal_day\
,obs_periodic_variance_seasonal_night,obs_periodic_variance_all,model_daily_mag,model_daily_mag_spring,model_daily_mag_summer,model_daily_mag_autumn\
,model_daily_mag_winter,model_daily_phase,model_daily_phase_spring,model_daily_phase_spring,model_daily_phase_summer\
,model_daily_phase_autumn,model_daily_phase_winter,model_daily_phase_min,model_daily_phase_min_spring,model_daily_phase_min_summer\
,model_daily_phase_min_autumn,model_daily_phase_min_winter,model_ave,model_ave_spring,model_ave_summer,model_ave_autumn,model_ave_winter\
,model_ave_day,model_ave_night,model_seasonal_mag,model_seasonal_mag_day,model_seasonal_phase,model_seasonal_phase_day,model_seasonal_phase_min\
,model_seasonal_phase_min_day,model_daily_ff,model_daily_ff_spring,model_daily_ff_summer,model_daily_ff_autumn,model_daily_ff_winter,model_seasonal_ff\
,model_seasonal_ff_day,model_seasonal_ff_night,model_daily_waveform,model_daily_waveform_spring,model_daily_waveform_summer,model_daily_waveform_autumn\
,model_daily_waveform_winter,model_seasonal_waveform,model_seasonal_waveform_day,model_seasonal_waveform_night,model_full_waveform,model_periodic_variance_daily\
,model_periodic_variance_daily_spring,model_periodic_variance_daily_summer,model_periodic_variance_daily_autumn,model_periodic_variance_daily_winter,model_periodic_variance_seasonal\
,model_periodic_variance_seasonal_day,model_periodic_variance_seasonal_night,model_periodic_variance_all = modules.get_periodic_specific(periodic_obs_fname,periodic_model_fname,obs_refs)
#------------------------------------------
tags = modules.get_tags(obs_refs)
loc_colors = {
'ANT': 'magenta',
'ARC': 'purple',
'AF': 'orange',
'AS': 'blue',
'EU': 'red',
'OC': 'black',
'O': 'yellow',
'NA': 'green',
'SA': 'pink'
}
model_p99_summer, model_diurnal_waveform_summer, model_diurnal_mag_autumn,
model_diurnal_ph_autumn, model_mean_autumn, model_p1_autumn, model_p5_autumn,
model_p25_autumn, model_p50_autumn, model_p75_autumn, model_p95_autumn,
model_p99_autumn, model_diurnal_waveform_autumn, model_diurnal_mag_winter,
model_diurnal_ph_winter, model_mean_winter, model_p1_winter, model_p5_winter,
model_p25_winter, model_p50_winter, model_p75_winter, model_p95_winter,
model_p99_winter, model_diurnal_waveform_winter, model_seasonal_mag_day,
model_seasonal_ph_day, model_mean_day, model_p1_day, model_p5_day,
model_p25_day, model_p50_day, model_p75_day, model_p95_day, model_p99_day,
model_seasonal_waveform_day, model_seasonal_mag_night,
model_seasonal_ph_night, model_mean_night, model_p1_night, model_p5_night,
model_p25_night, model_p50_night, model_p75_night, model_p95_night,
model_p99_night, model_seasonal_waveform_night, model_daily_h3_mag,
model_daily_h2_mag, model_daily_h1_mag, model_daily_mag, model_annual_h3_mag,
model_annual_h2_mag, model_annual_h1_mag,
model_annual_mag) = modules.get_periodic_specific(obs_periodic_fname,
model_periodic_fname)
#group data
#xy_params
xy_params = [
obs_diurnal_mag, obs_diurnal_ph, obs_seasonal_mag, obs_seasonal_ph,
obs_mean, obs_p1, obs_p5, obs_p25, obs_p50, obs_p75, obs_p95, obs_p99,
obs_pc_var_daily, obs_pc_var_seasonal, obs_pc_var_full,
obs_diurnal_mag_spring, obs_diurnal_ph_spring, obs_mean_spring,
obs_p1_spring, obs_p5_spring, obs_p25_spring, obs_p50_spring,
obs_p75_spring, obs_p95_spring, obs_p99_spring, obs_diurnal_mag_summer,
obs_diurnal_ph_summer, obs_mean_summer, obs_p1_summer, obs_p5_summer,
obs_p25_summer, obs_p50_summer, obs_p75_summer, obs_p95_summer,
obs_p99_summer, obs_diurnal_mag_autumn, obs_diurnal_ph_autumn,
obs_mean_autumn, obs_p1_autumn, obs_p5_autumn, obs_p25_autumn,
(obs_diurnal_mag,obs_diurnal_ph,obs_seasonal_mag,obs_seasonal_ph,obs_mean,obs_p1,obs_p5,obs_p25,obs_p50,obs_p75,obs_p95,obs_p99,obs_diurnal_ave_waveform,obs_diurnal_ave_waveform_extended,obs_diurnal_season_waveform_extended,obs_seasonal_waveform,obs_seasonal_waveform_extended,obs_full_ave_waveform,obs_full_season_waveform,obs_pc_var_daily,obs_pc_var_seasonal,obs_pc_var_full, obs_diurnal_mag_spring,obs_diurnal_ph_spring,obs_mean_spring,obs_p1_spring,obs_p5_spring,obs_p25_spring,obs_p50_spring,obs_p75_spring,obs_p95_spring,obs_p99_spring,obs_diurnal_waveform_spring, obs_diurnal_mag_summer,obs_diurnal_ph_summer,obs_mean_summer,obs_p1_summer,obs_p5_summer,obs_p25_summer,obs_p50_summer,obs_p75_summer,obs_p95_summer,obs_p99_summer,obs_diurnal_waveform_summer, obs_diurnal_mag_autumn,obs_diurnal_ph_autumn,obs_mean_autumn,obs_p1_autumn,obs_p5_autumn,obs_p25_autumn,obs_p50_autumn,obs_p75_autumn,obs_p95_autumn,obs_p99_autumn,obs_diurnal_waveform_autumn, obs_diurnal_mag_winter,obs_diurnal_ph_winter,obs_mean_winter,obs_p1_winter,obs_p5_winter,obs_p25_winter,obs_p50_winter,obs_p75_winter,obs_p95_winter,obs_p99_winter,obs_diurnal_waveform_winter, obs_seasonal_mag_day,obs_seasonal_ph_day,obs_mean_day,obs_p1_day,obs_p5_day,obs_p25_day,obs_p50_day,obs_p75_day,obs_p95_day,obs_p99_day,obs_seasonal_waveform_day, obs_seasonal_mag_night,obs_seasonal_ph_night,obs_mean_night,obs_p1_night,obs_p5_night,obs_p25_night,obs_p50_night,obs_p75_night,obs_p95_night,obs_p99_night,obs_seasonal_waveform_night, obs_daily_h3_mag,obs_daily_h2_mag,obs_daily_h1_mag,obs_daily_mag,obs_annual_h3_mag,obs_annual_h2_mag,obs_annual_h1_mag,obs_annual_mag, model_diurnal_mag,model_diurnal_ph,model_seasonal_mag,model_seasonal_ph,model_mean,model_p1,model_p5,model_p25,model_p50,model_p75,model_p95,model_p99,model_diurnal_ave_waveform,model_diurnal_ave_waveform_extended,model_diurnal_season_waveform_extended,model_seasonal_waveform,model_seasonal_waveform_extended,model_full_ave_waveform,model_full_season_waveform,model_pc_var_daily,model_pc_var_seasonal,model_pc_var_full, model_diurnal_mag_spring,model_diurnal_ph_spring,model_mean_spring,model_p1_spring,model_p5_spring,model_p25_spring,model_p50_spring,model_p75_spring,model_p95_spring,model_p99_spring,model_diurnal_waveform_spring, model_diurnal_mag_summer,model_diurnal_ph_summer,model_mean_summer,model_p1_summer,model_p5_summer,model_p25_summer,model_p50_summer,model_p75_summer,model_p95_summer,model_p99_summer,model_diurnal_waveform_summer, model_diurnal_mag_autumn,model_diurnal_ph_autumn,model_mean_autumn,model_p1_autumn,model_p5_autumn,model_p25_autumn,model_p50_autumn,model_p75_autumn,model_p95_autumn,model_p99_autumn,model_diurnal_waveform_autumn, model_diurnal_mag_winter,model_diurnal_ph_winter,model_mean_winter,model_p1_winter,model_p5_winter,model_p25_winter,model_p50_winter,model_p75_winter,model_p95_winter,model_p99_winter,model_diurnal_waveform_winter, model_seasonal_mag_day,model_seasonal_ph_day,model_mean_day,model_p1_day,model_p5_day,model_p25_day,model_p50_day,model_p75_day,model_p95_day,model_p99_day,model_seasonal_waveform_day, model_seasonal_mag_night,model_seasonal_ph_night,model_mean_night,model_p1_night,model_p5_night,model_p25_night,model_p50_night,model_p75_night,model_p95_night,model_p99_night,model_seasonal_waveform_night, model_daily_h3_mag,model_daily_h2_mag,model_daily_h1_mag,model_daily_mag,model_annual_h3_mag,model_annual_h2_mag,model_annual_h1_mag,model_annual_mag) = modules.get_periodic_specific(obs_periodic_fname,model_periodic_fname) #group data #xy_params xy_params = [obs_diurnal_mag,obs_diurnal_ph,obs_seasonal_mag,obs_seasonal_ph,obs_mean,obs_p1,obs_p5,obs_p25,obs_p50,obs_p75,obs_p95,obs_p99,obs_pc_var_daily,obs_pc_var_seasonal,obs_pc_var_full, obs_diurnal_mag_spring,obs_diurnal_ph_spring,obs_mean_spring,obs_p1_spring,obs_p5_spring,obs_p25_spring,obs_p50_spring,obs_p75_spring,obs_p95_spring,obs_p99_spring, obs_diurnal_mag_summer,obs_diurnal_ph_summer,obs_mean_summer,obs_p1_summer,obs_p5_summer,obs_p25_summer,obs_p50_summer,obs_p75_summer,obs_p95_summer,obs_p99_summer, obs_diurnal_mag_autumn,obs_diurnal_ph_autumn,obs_mean_autumn,obs_p1_autumn,obs_p5_autumn,obs_p25_autumn,obs_p50_autumn,obs_p75_autumn,obs_p95_autumn,obs_p99_autumn, obs_diurnal_mag_winter,obs_diurnal_ph_winter,obs_mean_winter,obs_p1_winter,obs_p5_winter,obs_p25_winter,obs_p50_winter,obs_p75_winter,obs_p95_winter,obs_p99_winter, obs_seasonal_mag_day,obs_seasonal_ph_day,obs_mean_day,obs_p1_day,obs_p5_day,obs_p25_day,obs_p50_day,obs_p75_day,obs_p95_day,obs_p99_day, obs_seasonal_mag_night,obs_seasonal_ph_night,obs_mean_night,obs_p1_night,obs_p5_night,obs_p25_night,obs_p50_night,obs_p75_night,obs_p95_night,obs_p99_night, obs_daily_h3_mag,obs_daily_h2_mag,obs_daily_h1_mag,obs_daily_mag,obs_annual_h3_mag,obs_annual_h2_mag,obs_annual_h1_mag,obs_annual_mag, model_diurnal_mag,model_diurnal_ph,model_seasonal_mag,model_seasonal_ph,model_mean,model_p1,model_p5,model_p25,model_p50,model_p75,model_p95,model_p99,model_pc_var_daily,model_pc_var_seasonal,model_pc_var_full, model_diurnal_mag_spring,model_diurnal_ph_spring,model_mean_spring,model_p1_spring,model_p5_spring,model_p25_spring,model_p50_spring,model_p75_spring,model_p95_spring,model_p99_spring, model_diurnal_mag_summer,model_diurnal_ph_summer,model_mean_summer,model_p1_summer,model_p5_summer,model_p25_summer,model_p50_summer,model_p75_summer,model_p95_summer,model_p99_summer,