machine_path_prefix = '/p/user_pub/e3sm/e3sm_diags_data/' #param.reference_data_path = os.path.join(machine_path_prefix, 'obs_for_e3sm_diags/time-series/') #param.test_data_path = os.path.join(machine_path_prefix, 'test_model_data_for_acme_diags/time-series/CESM1-CAM5_cmip/') param.reference_data_path = '/p/user_pub/e3sm/zhang40/e3sm_cmip6_xmls/v1_water_cycle/amip/' param.test_data_path = '/p/user_pub/e3sm/zhang40/e3sm_cmip6_xmls/v1_water_cycle/amip/r1i1p1f1' param.test_name = 'e3sm_v1_amip_r1' #For compy #prefix = '/compyfs/www/zhan429/doc_examples/' #For cori #prefix = '/global/project/projectdirs/acme/www/zhang40/' #For acme1 #prefix = '/var/www/acme/acme-diags/zhang40/tests/' #param.results_dir = os.path.join(prefix, 'area_mean_amip_r1-3') param.multiprocessing = True param.num_workers = 25 # We're passing in this new object as well, in # addition to the CoreParameter object. ts_param = AreaMeanTimeSeriesParameter() ts_param.ref_names = ['r2i1p1f1', 'r3i1p1f1'] #This setting plot model data only ts_param.start_yr = '1870' ts_param.end_yr = '2014' runner.sets_to_run = ['area_mean_time_series'] runner.run_diags([param, ts_param])
#For cori #machine_path_prefix = '/global/project/projectdirs/acme/acme_diags' param.test_data_path = os.path.join( machine_path_prefix, 'test_model_data_for_acme_diags/time-series/E3SM_v1_historical_r1/') #param.reference_data_path = os.path.join(machine_path_prefix, 'obs_for_e3sm_diags/time-series/') param.reference_data_path = os.path.join( machine_path_prefix, 'test_model_data_for_acme_diags/time-series/') param.test_name = 'e3sm_v1' param.rune_type = 'model-vs-model' prefix = '/compyfs/www/zhan429/examples/' param.results_dir = os.path.join(prefix, 'area_mean_with_models') param.multiprocessing = True param.num_workers = 40 # We're passing in this new object as well, in # addtion to the CoreParameter object. ts_param = AreaMeanTimeSeriesParameter() ts_param.ref_names = [ 'E3SM_v1_historical_r2', 'E3SM_v1_historical_r3', 'E3SM_v1_historical_r4', 'E3SM_v1_historical_r5' ] #This setting plot model data only ts_param.start_yr = '1850' ts_param.end_yr = '2014' runner.sets_to_run = ['area_mean_time_series'] runner.run_diags([param, ts_param])
param.test_start_yr = '2002' param.test_end_yr = '2004' # if the time range is invalid for reference datasets, only model results will be plotted param.ref_timeseries_input = True param.ref_start_yr = '2002' param.ref_end_yr = '2004' prefix = '/compyfs/www/zhan429/run_with_api/' param.results_dir = os.path.join(prefix, 'area_mean') attrs_core = vars(param) #print (', '.join("%s: %s" % item for item in attrs_core.items())) print('-------------------------------') ts_param = AreaMeanTimeSeriesParameter() ts_param.ref_names = ['none'] #param.multiprocessing = True #param.num_workers = 40 # We're passing in this new object as well, in # addtion to the CoreParameter object. runner.sets_to_run = ['area_mean_time_series'] print('*******************************') #attrs = vars(ts_param) #print (', '.join("%s: %s" % item for item in attrs.items())) #runner.run_diags([param])
#machine_path_prefix = '/global/project/projectdirs/acme/acme_diags' #for acme1 machine_path_prefix = '/p/user_pub/e3sm/e3sm_diags_data/' #param.reference_data_path = os.path.join(machine_path_prefix, 'obs_for_e3sm_diags/time-series/') #param.test_data_path = os.path.join(machine_path_prefix, 'test_model_data_for_acme_diags/time-series/CESM1-CAM5_cmip/') param.reference_data_path = '/p/user_pub/e3sm/zhang40/e3sm_cmip6_xmls/v1_water_cycle/abrupt-4xCO2/' param.test_data_path = '/p/user_pub/e3sm/zhang40/e3sm_cmip6_xmls/v1_water_cycle/piControl/r1i1p1f1' param.test_name = 'e3sm_v1_piControl' #For compy #prefix = '/compyfs/www/zhan429/doc_examples/' #For cori #prefix = '/global/project/projectdirs/acme/www/zhang40/' #For acme1 #prefix = '/var/www/acme/acme-diags/zhang40/tests/' #param.results_dir = os.path.join(prefix, 'area_mean_picontrol_r1') param.multiprocessing = True param.num_workers = 25 # We're passing in this new object as well, in # addition to the CoreParameter object. ts_param = AreaMeanTimeSeriesParameter() ts_param.ref_names = [] #This setting plot model data only ts_param.start_yr = '0001' ts_param.end_yr = '0500' runner.sets_to_run = ['area_mean_time_series'] runner.run_diags([param, ts_param])