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])
示例#4
0
#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])