Changeset 4212


Ignore:
Timestamp:
Feb 5, 2007, 11:28:49 AM (17 years ago)
Author:
nick
Message:

update to dampier

Location:
anuga_work/production/dampier_2006
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/dampier_2006/build_dampier.py

    r4196 r4212  
    6161coast_in_dir_name = project.coast_in_dir_name
    6262offshore_in_dir_name = project.offshore_in_dir_name
    63 #offshore1_in_dir_name = project.offshore1_in_dir_name
    64 #offshore2_in_dir_name = project.offshore2_in_dir_name
     63offshore1_in_dir_name = project.offshore1_in_dir_name
     64offshore2_in_dir_name = project.offshore2_in_dir_name
    6565
    6666onshore_dir_name = project.onshore_dir_name
    6767coast_dir_name = project.coast_dir_name
    6868offshore_dir_name = project.offshore_dir_name
    69 #offshore1_dir_name = project.offshore1_dir_name
    70 #offshore2_dir_name = project.offshore2_dir_name
    71 
     69offshore1_dir_name = project.offshore1_dir_name
     70offshore2_dir_name = project.offshore2_dir_name
     71'''
    7272# creates DEM from asc data
    7373print "creates DEMs from asc data"
     
    7575                              basename_out=onshore_dir_name,
    7676                              use_cache=True, verbose=True)
    77 '''
     77
    7878convert_dem_from_ascii2netcdf(offshore1_in_dir_name,
    7979                              basename_out=offshore1_dir_name,
     
    8282                              basename_out=offshore2_dir_name,
    8383                              use_cache=True, verbose=True)
    84 '''
     84
    8585#creates pts file for onshore DEM
    8686print "creates pts file for onshore DEM"
    8787dem2pts(onshore_dir_name, use_cache=True, verbose=True)
    88 '''
     88
    8989#creates pts file for clipped DEMs
    9090dem2pts(offshore1_dir_name,
     
    9494        basename_out=offshore2_dir_name,
    9595        use_cache=True, verbose=True)
    96 '''
     96
    9797print'create Geospatial onshore objects from topographies'
    9898#print'create Geospatial coastal objects from topographies'
     
    101101#print'create Geospatial offshore objects from topographies'
    102102G = Geospatial_data(file_name = onshore_dir_name + '.pts') +\
    103      Geospatial_data(file_name = coast_in_dir_name + '.txt') +\
    104      Geospatial_data(file_name = offshore_in_dir_name + '.txt')
     103    Geospatial_data(file_name = coast_in_dir_name + '.txt') +\
     104    Geospatial_data(file_name = offshore_in_dir_name + '.txt') +\
     105    Geospatial_data(file_name = offshore1_dir_name + '.pts') +\
     106    Geospatial_data(file_name = offshore2_dir_name + '.pts')
     107     
    105108
    106109print'add all geospatial objects'
     
    124127G_small.export_points_file(project.combined_dir_name + '_small' + '.pts')
    125128G_other.export_points_file(project.combined_dir_name + '_other' + '.pts')
     129'''
     130print 'start reading:',project.combined_dir_name + '.pts'
     131G = Geospatial_data(file_name = project.combined_dir_name + '.pts')
     132print 'start export',project.combined_dir_name + '.txt'
     133G.export_points_file(project.combined_dir_name + '.txt')
    126134
    127135
     136
     137
     138
  • anuga_work/production/dampier_2006/project.py

    r4196 r4212  
    2727print 'gtime: ', gtime
    2828
    29 tide = 0.6
     29tide = 2.4
    3030
    3131#Making assumptions about the location of scenario data
     
    4040coast_name = 'coastline_edited_w_DEM'
    4141offshore_name = 'clipped_bathy'
    42 #offshore1_name = 'elev_501'
    43 #offshore2_name = 'inferrec_e'
     42offshore1_name = 'elev_501'
     43offshore2_name = 'inferrec_e'
    4444
    4545#final topo name
     
    5555coast_in_dir_name = topographies_in_dir + coast_name
    5656offshore_in_dir_name = topographies_in_dir + offshore_name
    57 #offshore1_in_dir_name = topographies_in_dir + offshore1_name
    58 #offshore2_in_dir_name = topographies_in_dir + offshore2_name
     57offshore1_in_dir_name = topographies_in_dir + offshore1_name
     58offshore2_in_dir_name = topographies_in_dir + offshore2_name
    5959
    6060onshore_dir_name = topographies_dir + onshore_name
    6161coast_dir_name = topographies_dir + coast_name
    6262offshore_dir_name = topographies_dir + offshore_name
    63 #offshore1_dir_name = topographies_dir + offshore1_name
    64 #offshore2_dir_name = topographies_dir + offshore2_name
     63offshore1_dir_name = topographies_dir + offshore1_name
     64offshore2_dir_name = topographies_dir + offshore2_name
    6565
    6666#final topo files
  • anuga_work/production/dampier_2006/run_dampier.py

    r4196 r4212  
    145145
    146146    domain.set_quantity('elevation',
    147                     filename = project.combined_dir_name + '.pts',
     147#                    filename = project.combined_dir_name + '.pts',
     148# MUST USE TXT FILES FOR CACHING TO WORK!
     149                    filename = project.combined_dir_name + '.txt',
    148150                    use_cache = True,
    149151                    verbose = True,
     
    179181print 'Reading Boundary file'
    180182Bf = File_boundary(boundaries_dir_name + '.sww',
    181                   domain, time_thinning=10, use_cache=True, verbose=True)
     183                  domain, time_thinning=1, use_cache=True, verbose=True)
    182184
    183185print 'finished reading boundary file'
     
    201203t0 = time.time()
    202204
    203 for t in domain.evolve(yieldstep = 120, finaltime = 34000):
     205for t in domain.evolve(yieldstep = 120, finaltime = 9000):
    204206    domain.write_time()
    205207    domain.write_boundary_statistics(tags = 'ocean')
     
    209211#    domain.write_boundary_statistics(tags = 'ocean')
    210212     
    211 #for t in domain.evolve(yieldstep = 60, finaltime = 28800
    212 #                       ,skip_initial_step = True):
    213 #    domain.write_time()
    214 #    domain.write_boundary_statistics(tags = 'ocean')   
    215 
    216 #for t in domain.evolve(yieldstep = 120, finaltime = 34800
    217 #                       ,skip_initial_step = True):
    218 #    domain.write_time()
    219 #    domain.write_boundary_statistics(tags = 'ocean')   
     213for t in domain.evolve(yieldstep = 60, finaltime = 28800
     214                       ,skip_initial_step = True):
     215    domain.write_time()
     216    domain.write_boundary_statistics(tags = 'ocean')   
     217
     218for t in domain.evolve(yieldstep = 120, finaltime = 34800
     219                       ,skip_initial_step = True):
     220    domain.write_time()
     221    domain.write_boundary_statistics(tags = 'ocean')   
    220222
    221223x, y = domain.get_maximum_inundation_location()
Note: See TracChangeset for help on using the changeset viewer.