Changeset 5609


Ignore:
Timestamp:
Aug 5, 2008, 1:47:42 PM (11 years ago)
Author:
kristy
Message:

updated run file with optimal alpha

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/busselton/run_busselton.py

    r5578 r5609  
    4141from anuga.fit_interpolate.benchmark_least_squares import mem_usage
    4242from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
     43from anuga.geospatial_data.geospatial_data import find_optimal_smoothing_parameter
    4344
    4445# Application specific imports
     
    114115    #barrier()
    115116   
     117         covariance_value,alpha = find_optimal_smoothing_parameter (data_file= kwargs['elevation_file'],
     118                                alpha_list=[0.001, 0.01, 0.1, 0.15, 0.2, 0.25, 0.3, 0.4, 0.5],
     119                                mesh_file = project.meshes_dir_name+'.msh')
     120        print 'optimal alpha', covariance_value,alpha       
    116121
    117122    #-------------------------------------------------------------------------
     
    181186    print 'domain id', id(domain)
    182187   
    183     print'set_boundary'
    184 
    185188    boundary_urs_out=project.boundaries_dir_name
    186189   
    187190    print 'Available boundary tags', domain.get_boundary_tags()
    188     Bf = File_boundary(boundary_urs_out+'.sts',
    189                    domain, time_thinning=1,
     191    Bf = File_boundary(boundary_urs_out+'.sts', # Change from file_boundary
     192                   domain, mean_stage= project.tide,
     193                   time_thinning=1,
    190194                   use_cache=True,
    191                    verbose = True,
    192                    boundary_polygon=bounding_polygon)
     195                   verbose = True)
     196                   #boundary_polygon=bounding_polygon)
    193197   
    194198    Br = Reflective_boundary(domain)
    195199    Bd = Dirichlet_boundary([kwargs['tide'],0,0])
    196     #Bw = Dirichlet_boundary([kwargs['tide']+10.0,0,0]) # To be deleted for FESA runs
    197 
     200
     201    print dir(Bf)
    198202    print 'finished reading boundary file'
    199203
     
    215219        domain.write_time()
    216220        domain.write_boundary_statistics(tags = 'ocean')
    217 
    218 ##        # To be deleted for FESA runs
    219 ##        if allclose(t, 240):
    220 ##            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bw})
    221 ##
    222 ##        if allclose(t, 1440):
    223 ##            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bd})
    224 
    225221
    226222
     
    255251    kwargs['output_dir']=project.output_run_time_dir
    256252    kwargs['bathy_file']=project.combined_dir_name+'.txt'
    257 #    kwargs['bathy_file']=project.combined_small_dir_name + '.pts'
    258  
    259     kwargs['file_name']=project.home+'detail.csv'
     253  kwargs['file_name']=project.home+'detail.csv'
    260254    kwargs['aa_scenario_name']=project.scenario_name
    261255    kwargs['ab_time']=project.time
Note: See TracChangeset for help on using the changeset viewer.