run_zonal_mean_2d_strat_model_model.py