import os
from e3sm_diags.parameter.core_parameter import CoreParameter
from e3sm_diags.run import runner


param = CoreParameter()

param.run_type = 'model_vs_model'
param.reference_data_path = '/pscratch/sd/x/xsongsio/DATA/optimization/O2_diag/clim_rgr/'
param.reference_name = '20230804.v3alpha02.zm_opt_O2.1iteration.F2010.chrysalis'
param.test_data_path = '/pscratch/sd/x/xsongsio/DATA/optimization/O3_diag/clim_rgr/'
param.test_name = '20230804.v3alpha02.zm_opt_O3.1iteration.F2010.chrysalis'
#param.seasons = ["ANN"]   #all seasons ["ANN","DJF", "MAM", "JJA", "SON"] will run,if comment out"

prefix = '/global/cfs/cdirs/e3sm/www/xsongsio/v3alpha02.chrysalis_O3-O2/'
param.results_dir = os.path.join(prefix, '1yr')
param.test_title = ' '
param.reference_title = ' '
param.diff_title = '1iteration_chrysalis_O3 - 1iteration_chrysalis_O2'
# Use the following if running in parallel:
param.multiprocessing = True
param.num_workers = 32

# Use below to run all core sets of diags:
runner.sets_to_run = ['lat_lon','zonal_mean_xy', 'zonal_mean_2d', 'polar', 'cosp_histogram', 'meridional_mean_2d']
#runner.sets_to_run = ['lat_lon','zonal_mean_2d', 'polar', 'cosp_histogram', 'meridional_mean_2d']
# Use below to run lat_lon map only:
#runner.sets_to_run = ['lat_lon']
runner.run_diags([param])