def METHOD_NAME():
    input={}
    # grid
    input['x']= [0.0, 1.0, 4.0] 
    # stationary flow
    input['u'] = [1.0, 1.0, 1.0, 1.0, 1.0]
    # cross sectional area
    input['a'] = [1.0, 1.0, 1.0, 1.0, 1.0]
    # initial concentrations
    input['c'] = [0.1, 0.2, 0.3, 0.4, 0.5]
    # simulation timespan
    input['refdate'] = '01 dec 2000'
    #unit is always seconds
    input['unit'] = 'seconds'
    input['time'] = [0.0, 1.0, 10.0]
    # sources mass/m^3/s
    input['source_locations'] = [2]
    input['source_labels'] = ['default']
    input['source_values']= {}
    input['source_values']['default'] = [5.0]
    # boundaries
    input['bound_labels']=['left', 'right']
    input['bound_locations']=[0, -1]
    input['bound_values']={}
    input['bound_values']['left']=[-1000.0, 0.01, 0.02, 0.03]
    input['bound_values']['right']=[0.0]
    #output (index based and 0 based)
    input['output_file'] = 'default.output'
    input['output_locations'] = [1, 2]
    input['output_labels']=['defaultOutput1', 'defaultOutput2']
    return input