Python calc_subregion_area_mean_and_stdの例

コード例 #1

ファイル: run_RCMES.py プロジェクト: maheshguru/climate

    model_datasets.append(dsp.ensemble(model_datasets))
    model_names.append('ENS')
""" Step 6: Generate subregion average and standard deviation """
if config['use_subregions']:
    # sort the subregion by region names and make a list
    subregions = sorted(config['subregions'].items(),
                        key=operator.itemgetter(0))

    # number of subregions
    nsubregion = len(subregions)

    print('Calculating spatial averages and standard deviations of ',
          str(nsubregion), ' subregions')

    obs_subregion_mean, obs_subregion_std, subregion_array = (
        utils.calc_subregion_area_mean_and_std([obs_dataset], subregions))
    model_subregion_mean, model_subregion_std, subregion_array = (
        utils.calc_subregion_area_mean_and_std(model_datasets, subregions))
""" Step 7: Write a netCDF file """
workdir = config['workdir']
if workdir[-1] != '/':
    workdir = workdir + '/'
print 'Writing a netcdf file: ', workdir + config['output_netcdf_filename']
if not os.path.exists(workdir):
    os.system("mkdir -p " + workdir)

if config['use_subregions']:
    dsp.write_netcdf_multiple_datasets_with_subregions(
        obs_dataset,
        obs_name,
        model_datasets,

コード例 #2

ファイル: run_RCMES.py プロジェクト: feihugis/climate

print 'Generating multi-model ensemble'
if len(model_datasets) >= 2.:
    model_datasets.append(dsp.ensemble(model_datasets))
    model_names.append('ENS')

""" Step 6: Generate subregion average and standard deviation """
if config['use_subregions']:
    # sort the subregion by region names and make a list
    subregions= sorted(config['subregions'].items(),key=operator.itemgetter(0))

    # number of subregions
    nsubregion = len(subregions)

    print 'Calculating spatial averages and standard deviations of ',str(nsubregion),' subregions'

    ref_subregion_mean, ref_subregion_std, subregion_array = utils.calc_subregion_area_mean_and_std([ref_dataset], subregions) 
    model_subregion_mean, model_subregion_std, subregion_array = utils.calc_subregion_area_mean_and_std(model_datasets, subregions) 

""" Step 7: Write a netCDF file """
workdir = config['workdir']
if workdir[-1] != '/':
    workdir = workdir+'/'
print 'Writing a netcdf file: ',workdir+config['output_netcdf_filename']
if not os.path.exists(workdir):
    os.system("mkdir "+workdir)

if config['use_subregions']:
    dsp.write_netcdf_multiple_datasets_with_subregions(ref_dataset, ref_name, model_datasets, model_names,
                                                       path=workdir+config['output_netcdf_filename'],
                                                       subregions=subregions, subregion_array = subregion_array, 
                                                       ref_subregion_mean=ref_subregion_mean, ref_subregion_std=ref_subregion_std,

コード例 #3

if len(model_datasets) >= 2.:
    model_datasets.append(dsp.ensemble(model_datasets))
    model_names.append('ENS')
""" Step 6: Generate subregion average and standard deviation """
if config['use_subregions']:
    # sort the subregion by region names and make a list
    subregions = sorted(config['subregions'].items(),
                        key=operator.itemgetter(0))

    # number of subregions
    nsubregion = len(subregions)

    print 'Calculating spatial averages and standard deviations of ', str(
        nsubregion), ' subregions'

    ref_subregion_mean, ref_subregion_std, subregion_array = utils.calc_subregion_area_mean_and_std(
        [ref_dataset], subregions)
    model_subregion_mean, model_subregion_std, subregion_array = utils.calc_subregion_area_mean_and_std(
        model_datasets, subregions)
""" Step 7: Write a netCDF file """
workdir = config['workdir']
if workdir[-1] != '/':
    workdir = workdir + '/'
print 'Writing a netcdf file: ', workdir + config['output_netcdf_filename']
if not os.path.exists(workdir):
    os.system("mkdir " + workdir)

if config['use_subregions']:
    dsp.write_netcdf_multiple_datasets_with_subregions(
        ref_dataset,
        ref_name,
        model_datasets,

コード例 #4

ファイル: run_RCMES.py プロジェクト: MichaelArthurAnderson/climate

    target_datasets.append(dsp.ensemble(target_datasets))
    target_names.append('ENS')

""" Step 5: Generate subregion average and standard deviation """
if config['use_subregions']:
    # sort the subregion by region names and make a list
    subregions= sorted(config['subregions'].items(),key=operator.itemgetter(0))

    # number of subregions
    nsubregion = len(subregions)

    print('Calculating spatial averages and standard deviations of {} subregions'
          .format(nsubregion))

    reference_subregion_mean, reference_subregion_std, subregion_array = (
        utils.calc_subregion_area_mean_and_std([reference_dataset], subregions))
    target_subregion_mean, target_subregion_std, subregion_array = (
        utils.calc_subregion_area_mean_and_std(target_datasets, subregions))

""" Step 6: Write a netCDF file """
workdir = config['workdir']
if workdir[-1] != '/':
    workdir = workdir+'/'
print('Writing a netcdf file: ',workdir+config['output_netcdf_filename'])
if not os.path.exists(workdir):
    os.system("mkdir -p "+workdir)

if config['use_subregions']:
    dsp.write_netcdf_multiple_datasets_with_subregions(
        reference_dataset, reference_name, target_datasets, target_names,
        path=workdir+config['output_netcdf_filename'],

コード例 #5

ファイル: run_RCMES.py プロジェクト: sgoggins/climate

    target_names.append('ENS')
""" Step 5: Generate subregion average and standard deviation """
if config['use_subregions']:
    # sort the subregion by region names and make a list
    subregions = sorted(config['subregions'].items(),
                        key=operator.itemgetter(0))

    # number of subregions
    nsubregion = len(subregions)

    print(
        'Calculating spatial averages and standard deviations of {} subregions'
        .format(nsubregion))

    reference_subregion_mean, reference_subregion_std, subregion_array = (
        utils.calc_subregion_area_mean_and_std([reference_dataset],
                                               subregions))
    target_subregion_mean, target_subregion_std, subregion_array = (
        utils.calc_subregion_area_mean_and_std(target_datasets, subregions))
""" Step 6: Write a netCDF file """
workdir = config['workdir']
if workdir[-1] != '/':
    workdir = workdir + '/'
print('Writing a netcdf file: ', workdir + config['output_netcdf_filename'])
if not os.path.exists(workdir):
    os.system("mkdir -p " + workdir)

if config['use_subregions']:
    dsp.write_netcdf_multiple_datasets_with_subregions(
        reference_dataset,
        reference_name,
        target_datasets,