Changeset 5387


Ignore:
Timestamp:
Jun 2, 2008, 8:59:09 AM (16 years ago)
Author:
kristy
Message:
 
Location:
anuga_work/production
Files:
4 edited

Legend:

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

    r5381 r5387  
    9090G2 = Geospatial_data(file_name = coast_in_dir_name + '.txt',verbose=True)
    9191G3 = Geospatial_data(file_name = coast_in_dir_name1 + '.txt',verbose=True)
    92 #G3 = Geospatial_data(file_name = island_dir_name + '.pts',verbose=True)
    9392G_off = Geospatial_data(file_name = offshore_in_dir_name + '.txt',verbose=True)
    9493
     
    108107
    109108print'project.combined_dir_name + .txt',project.combined_dir_name + '.txt'
    110 #G_all=Geospatial_data(file_name = project.combined_dir_name + '.txt')
    111 print'split'
    112 G_all_1 = G.split(.10)
    113 print'export 1'
    114 G_all_1.export_points_file(project.combined_dir_name+'_small' + '.txt')
    115 #print'export 2'
    116 #G_all_2.export_points_file(project.combined_dir_name+'_other1' + '.xya')
    117109
    118 '''
    119110#-------------------------------------------------------------------------
    120111# Convert URS to SWW file for boundary conditions
    121112#-------------------------------------------------------------------------
    122113print 'starting to create boundary conditions'
     114from anuga.shallow_water.data_manager import urs2sww, urs_ungridded2sww
     115
    123116boundaries_in_dir_name = project.boundaries_in_dir_name
     117print 'boundaries_in_dir_name',project.boundaries_in_dir_name
    124118
    125 from anuga.shallow_water.data_manager import urs2sww
    126 
    127 print 'minlat=project.north_boundary, maxlat=project.south_boundary',project.north_boundary, project.south_boundary
    128 print 'minlon= project.west_boundary, maxlon=project.east_boundary',project.west_boundary, project.east_boundary
    129119
    130120#import sys; sys.exit()
    131121
    132 #if access(project.boundaries_dir,F_OK) == 0:
    133 #    mkdir (project.boundaries_dir)
     122urs_ungridded2sww(project.boundaries_in_dir_name, project.boundaries_in_dir_name,
     123                  verbose=True, mint=4000, maxt=80000, zscale=1)
    134124
    135 from caching import cache
    136 cache(urs2sww,
    137       (boundaries_in_dir_name,
    138        project.boundaries_dir_name1),
    139       {'verbose': True,
    140        'minlat': project.south_boundary,
    141        'maxlat': project.north_boundary,
    142        'minlon': project.west_boundary,
    143        'maxlon': project.east_boundary,
    144 #       'minlat': project.south,
    145 #       'maxlat': project.north,
    146 #       'minlon': project.west,
    147 #       'maxlon': project.east,
    148        'mint': 0, 'maxt': 40000,
    149 #       'origin': domain.geo_reference.get_origin(),
    150        'mean_stage': project.tide,
    151 #       'zscale': 1,                 #Enhance tsunami
    152        'fail_on_NaN': False},
    153        verbose = True,
    154        )
    155 #       dependencies = source_dir + project.boundary_basename + '.sww')
    156 
    157 '''
    158 
    159 
    160 
    161 
    162 
    163 
    164 
  • anuga_work/production/busselton/project.py

    r5381 r5387  
    4242finaltime=10000
    4343export_cellsize=50
    44 setup='trial'
    45 source='other'
     44setup='final'
     45source='exmouth'
    4646
    4747
     
    108108tide_dir = anuga_dir+'tide_data'+sep
    109109
    110 if source =='dampier':
    111     boundaries_name = 'broome_3854_17042007' #Dampier gun
    112     boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'dampier'+sep+'1_10000'+sep
    113 
    114 if source=='onslow':
    115     boundaries_name = 'broome_3859_16052007' #onslow_hedland_broome gun
    116     boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'onslow_hedland_broome'+sep+'1_10000'+sep
    117    
    118110if source=='exmouth':
    119     boundaries_name = 'broome_3103_18052007' #exmouth gun
     111    boundaries_name = 'busselton_3103_30052008' #exmouth gun
    120112    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
    121113
    122 if source=='other':
    123     boundaries_name = 'test' #exmouth gun
     114if source=='test':
     115    boundaries_name = 'other' #exmouth gun
    124116    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
    125117
     
    190182trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
    191183
    192 poly_mainland=read_polygon(topographies_in_dir+'initial_condition_south.csv')
     184poly_mainland=read_polygon(polygons_dir+'initial_condition.csv')
    193185
    194186print 'min number triangles', trigs_min
  • anuga_work/production/perth/build_perth.py

    r5386 r5387  
    141141
    142142urs_ungridded2sww(project.boundaries_in_dir_name, project.boundaries_in_dir_name,
    143                   verbose=True, mint=4000, maxt=10000, zscale=1)
     143                  verbose=True, mint=4000, maxt=80000, zscale=1)
    144144
    145145
  • anuga_work/production/perth/project.py

    r5386 r5387  
    130130
    131131#boundaries locations
    132 boundaries_in_dir_name = boundaries_in_dir + scenario_name
     132boundaries_in_dir_name = boundaries_in_dir + boundaries_name
    133133boundaries_dir = anuga_dir+'boundaries'+sep
    134134boundaries_dir_name = boundaries_dir + scenario_name
Note: See TracChangeset for help on using the changeset viewer.