Changeset 4422


Ignore:
Timestamp:
Apr 24, 2007, 12:03:30 PM (17 years ago)
Author:
nick
Message:

good version of dampier and update the export_results.py

Location:
anuga_work/production/dampier_2006
Files:
4 edited

Legend:

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

    r4347 r4422  
    55from os import sep
    66
    7 time_dir = '20061127_055714_run' #
     7time_dir = '20070418_035124_run' #
     8cellsize = 25
     9timestep = None
    810directory = project.output_dir
    911#name = directory + time_dir + sep + project.scenario_name
    10 name = directory+sep+project.scenario_name
     12name = directory+time_dir+sep+project.scenario_name
    1113
    12 is_parallel = True
     14is_parallel = False
    1315if is_parallel == True: nodes = 4
    1416print 'output dir:', name
    1517
    16 which_var = 1
     18which_var = 2
    1719if which_var == 0:  # Stage
    1820    outname = name + '_stage'
     
    3537    quantityname = 'elevation'  #Elevation
    3638
     39
     40
    3741if is_parallel == True:
     42#    print 'is_parallel',is_parallel
    3843    for i in range(0,nodes):
    3944        namei = name + '_P%d_%d' %(i,nodes)
     
    4247        sww2dem(namei, basename_out = outnamei,
    4348                    quantity = quantityname,
    44                     timestep = 1,
    45                     cellsize = 100,     
    46                     #easting_min = project.e_min_area,
    47                     #easting_max = project.e_max_area,
    48                     #northing_min = project.n_min_area,
    49                     #northing_max = project.n_max_area,       
     49                    timestep = timestep,
     50                    cellsize = cellsize,     
     51                    easting_min = project.e_min_area,
     52                    easting_max = project.e_max_area,
     53                    northing_min = project.n_min_area,
     54                    northing_max = project.n_max_area,       
    5055                    reduction = max,
    5156                    verbose = True,
     
    5560    sww2dem(name, basename_out = outname,
    5661                quantity = quantityname,
    57                 cellsize = 20,     
     62                timestep = timestep,
     63                cellsize = cellsize,     
    5864                easting_min = project.e_min_area,
    5965                easting_max = project.e_max_area,
  • anuga_work/production/dampier_2006/get_timeseries.py

    r4373 r4422  
    2424# sww file created from URS boundary
    2525production_dirs = {#'boundaries': '10000 yr wave height'
    26                    '20070413_033455_run': 'URS boundary normal'
     26                   '20070418_035124_run': 'URS boundary normal'
    2727                   }
    2828
    2929   
    30 is_parallel = True
     30is_parallel = False
    3131if is_parallel == True:
    3232    nodes = 8
     
    8989       
    9090    texname, elev_output = sww2timeseries(swwfiles,
    91                                           project.gauges_dir_name_simple,
     91                                          project.gauges_dir_name,
     92#                                          project.gauges_dir_name_simple,
    9293                                          production_dirs,
    9394                                          report = False,
  • anuga_work/production/dampier_2006/project.py

    r4401 r4422  
    2121
    2222#time stuff
    23 time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
     23time = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
     24#time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
    2425gtime = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
    2526build_time = time+'_build'
  • anuga_work/production/dampier_2006/run_dampier.py

    r4401 r4422  
    100100#domain = cache(Domain, (meshes_dir_name), {'use_cache':True, 'verbose':True}, verbose=True)
    101101#above don't work
    102 domain = Domain(meshes_dir_name, use_cache=True, verbose=True)
     102domain = Domain(meshes_dir_name, use_cache=False, verbose=True)
    103103print 'domain id', id(domain)
    104104print 'myhash', myhash(domain)     
     
    106106print domain.statistics()
    107107
    108 boundaries_dir_name=project.boundaries_dir_name
    109 
    110 print 'starting to create boundary conditions'
    111 
    112 from anuga.shallow_water.data_manager import urs2sww
     108#boundaries_dir_name=project.boundaries_dir_name
     109
     110#print 'starting to create boundary conditions'
     111
     112#from anuga.shallow_water.data_manager import urs2sww
    113113
    114114
     
    122122    from polygon import *
    123123    #following sets the stage/water to be offcoast only
    124     IC = Polygon_function( [(project.poly_mainland, 0.)], default = tide)
    125 #    IC = Polygon_function( [(project.poly_mainland, 0.)], default = 200)
     124#    IC = Polygon_function( [(project.poly_mainland, 0.)], default = tide)
     125    IC = Polygon_function( [(project.poly_mainland, 0.)], default = 200)
    126126    domain.set_quantity('stage', IC)
    127127    domain.set_quantity('friction', 0.01)
     
    159159domain.set_maximum_allowed_speed(0.1) # Allow a little runoff (0.1 is OK)
    160160print 'domain id', id(domain)
    161 #domain.beta_h = 0
     161domain.beta_h = 0
    162162#domain.limit2007 = 1
    163163
     
    170170
    171171Bf = Field_boundary(project.boundaries_dir_namea + '.sww',
    172                     domain, time_thinning=12, mean_stage=tide,
     172                    domain, time_thinning=1, mean_stage=tide,
    173173                    use_cache=True, verbose=True)
    174174
Note: See TracChangeset for help on using the changeset viewer.