run_zonal_mean_2d_strat.py