Changeset 5626


Ignore:
Timestamp:
Aug 7, 2008, 3:49:36 PM (15 years ago)
Author:
kristy
Message:

Updated

Location:
anuga_work/production
Files:
6 edited

Legend:

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

    r5624 r5626  
    7474
    7575scenario_name=project.scenario_name
    76 order_filename=os.path.join(project.boundaries_dir, 'thinned_bound_order.txt')
     76order_filename=os.path.join(project.boundaries_dir, 'thinned_bounding_order.txt')
    7777
    7878print 'reading', order_filename
  • anuga_work/production/busselton/project.py

    r5607 r5626  
    160160#gauges
    161161gauge_name = 'busselton.csv'
     162gauge_name2 = 'thinned_MGA50+1.csv'
     163
    162164gauges_dir = home+state+sep+scenario+sep+'anuga'+sep+'gauges'+sep
    163165gauges_dir_name = gauges_dir + gauge_name
     166gauges_dir_name2 = gauges_dir + gauge_name2
    164167
    165168buildings_filename = gauges_dir + 'Busselton_res_Project.csv'
  • anuga_work/production/busselton/run_busselton.py

    r5609 r5626  
    115115    #barrier()
    116116   
    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       
     117##        covariance_value,alpha = find_optimal_smoothing_parameter (data_file = kwargs['bathy_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       
    121121
    122122    #-------------------------------------------------------------------------
     
    189189   
    190190    print 'Available boundary tags', domain.get_boundary_tags()
    191     Bf = File_boundary(boundary_urs_out+'.sts', # Change from file_boundary
    192                    domain, mean_stage= project.tide,
     191    Bf = Field_boundary(boundary_urs_out+'.sts', # Change from file_boundary
     192                   domain, mean_stage=project.tide,
    193193                   time_thinning=1,
    194194                   use_cache=True,
    195                    verbose = True)
    196                    #boundary_polygon=bounding_polygon)
     195                   verbose = True,
     196                   boundary_polygon=bounding_polygon)
    197197   
    198198    Br = Reflective_boundary(domain)
     
    251251    kwargs['output_dir']=project.output_run_time_dir
    252252    kwargs['bathy_file']=project.combined_dir_name+'.txt'
    253   kwargs['file_name']=project.home+'detail.csv'
     253    kwargs['file_name']=project.home+'detail.csv'
    254254    kwargs['aa_scenario_name']=project.scenario_name
    255255    kwargs['ab_time']=project.time
  • anuga_work/production/perth/export_results.py

    r5592 r5626  
    1313#time_dir = '20080725_173911_run_final_0.6_polyline_alpha0.1_kvanputt'
    1414#time_dir = '20080728_112519_run_final_0.6_polyline_alpha0.2_kvanputt'
    15 time_dir = '20080728_154123_run_final_0.6_polyline_alpha0.15_kvanputt'
     15time_dir = '20080807_134943_run_final_0.0_polyline_alpha0.1_kvanputt'
    1616
    1717cellsize = 25
  • anuga_work/production/perth/project.py

    r5607 r5626  
    157157#gauges
    158158gauge_name = 'perth.csv'
    159 gauge_name2 = 'thinned_MGA50+50.csv'
     159gauge_name2 = 'thinned_MGA50+1.csv'
    160160
    161161gauges_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'gauges'+sep
     
    211211res_sorrento_gauge = 500*res_factor
    212212
     213#Polygon designed to incorporate Dredge Area from Fremantle to
     214#Rockingham the steep incline was making the mesh go to 0
     215poly_dredge = read_polygon(polygons_dir+'DredgeArea.csv')
     216res_dredge = 1000*res_factor
    213217
    214218#assert zone == refzone
     
    218222                     ,[poly_rockingham,res_rockingham],[poly_geordie_bay,res_geordie_bay]
    219223                     ,[poly_sorrento_gauge,res_sorrento_gauge],[poly_rottnest_in, res_rottnest_in]
    220                      ,[poly_rottnest_ex, res_rottnest_ex], [poly_garden, res_garden]]
     224                     ,[poly_rottnest_ex, res_rottnest_ex], [poly_garden, res_garden]
     225                     ,[poly_dredge, res_dredge]]
    221226
    222227   
  • anuga_work/production/perth/run_perth.py

    r5617 r5626  
    210210                   time_thinning=1,
    211211                   use_cache=True,
    212                    verbose = True)
    213                   # boundary_polygon=bounding_polygon)
     212                   verbose = True,
     213                  boundary_polygon=bounding_polygon)
    214214
    215215    Br = Reflective_boundary(domain)
Note: See TracChangeset for help on using the changeset viewer.