Changeset 3423


Ignore:
Timestamp:
Jul 25, 2006, 4:57:26 PM (19 years ago)
Author:
duncan
Message:

sync scenarios

Location:
development/dam_2006
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • development/dam_2006/project.py

    r3421 r3423  
    3838#community_scenario = gaugedir + 'community_onslow.csv'
    3939#gaugetimeseries = gaugedir + 'onslow'
    40 depth_filename =
     40depth_filename = outputtimedir + "depth.csv"
     41velocity_filename = outputtimedir + "velocity.csv"
    4142
    4243
  • development/dam_2006/run_dam.py

    r3421 r3423  
    2727from pyvolution.util import Screen_Catcher
    2828from pyvolution.region import Set_region
     29from fit_interpolate.interpolate import interpolate_sww2csv
     30import create_mesh
    2931
    3032# Application specific imports
    3133import project                 # Definition of file names and polygons
    32 import create_mesh
    3334
    3435def main():
     
    4344    print'output dir', project.outputtimedir
    4445
     46    #FIXME this isn't working
    4547    #normal screen output is stored in
    4648    screen_output_name = project.outputtimedir + "screen_output.txt"
     
    5355    create_mesh.generate(project.mesh_filename,
    5456                         0.75,
    55                          is_course=True) # this creates the mesh
    56                          #is_course=False) # this creates the mesh
     57                         #is_course=True) # this creates the mesh
     58                         is_course=False) # this creates the mesh
    5759
    5860    head,tail = path.split(project.mesh_filename)
     
    106108    t0 = time.time()
    107109
    108     for t in domain.evolve(yieldstep = 0.1, finaltime = 15):
     110    for t in domain.evolve(yieldstep = 0.1, finaltime = 25):
    109111        domain.write_time()
    110112   
     
    122124    # Calc gauge info
    123125    #-------------------------------------------------------------------------
    124     #interpolate_sww2csv(project.basename +".sww",
    125     #                    points,
    126     #                    "depth.csv",
    127     #                    "velocity.csv")
     126    interpolate_sww2csv(project.outputtimedir + project.basename +".sww",
     127                        points,
     128                        project.depth_filename,
     129                        project.velocity_filename)
    128130
    129131#-------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.