Changeset 2076


Ignore:
Timestamp:
Nov 25, 2005, 2:47:32 PM (18 years ago)
Author:
duncan
Message:

update

Location:
production/gippsland_2005
Files:
1 added
3 edited

Legend:

Unmodified
Added
Removed
  • production/gippsland_2005/create_mesh.py

    r2059 r2076  
    1010from coordinate_transforms.geo_reference import Geo_reference
    1111
    12 def create_mesh(max_area):
     12def create_mesh(max_area, mesh_file=None):
    1313    #Lower left corner
    1414    xllcorner = 570066
     
    5353
    5454    m.generateMesh('q28.0z', maxArea = max_area)
    55     #m.export_mesh_file('lakes_small_100_test.tsh')
    56     return m, len(m.getTriangulation())
     55    if mesh_file is None:   
     56        return m, len(m.getTriangulation())
     57    else:
     58        m.export_mesh_file(mesh_file)
     59        return len(m.getTriangulation())
     60
    5761
    5862#-------------------------------------------------------------
  • production/gippsland_2005/project.py

    r2058 r2076  
    1010#Making assumptions about the location of scenario data
    1111# Assumes the INUNDATIONHOME environ variable exists
     12
    1213scenario_dir_name = 'gippsland_storm_surge_scenario_2005'
    13 
    14 basename = 'test'
    15 
     14basename = 'lakes_100'
    1615boundary_basename = 'csiro_bondary'
    1716
     
    3231vcur_dir =  boundarydir + 'vcur_expanded'
    3332
    34 meshname = meshdir + basename
    35 demname = datadir + basename
    36 boundaryname = boundarydir + boundary_basename
    37 outputname = outputdir + basename  #Used by post processing
     33meshname = meshdir + basename + '.tsh'
     34mesh_elevname = meshdir + basename + '_elev' + '.tsh'
     35demname = datadir + basename + '.dem'
     36#demname = datadir + 'lakes_100'
     37pointname = datadir + basename + '.pts'
     38boundaryname = boundarydir + boundary_basename + '.sww'
     39outputname = outputdir + basename  + '.sww' #Used by post processing
    3840
    3941gauge_filename = 'all_bld_ind.csv'
  • production/gippsland_2005/run_small_gippsland.py

    r2063 r2076  
    4444# NOTE: This function is dependent on a lot of files, and I haven't listed
    4545# them, since there is so many and I don't think they will be changing.
     46
     47print "project.boundaryname + '.sww'",project.boundaryname + '.sww'
    4648cache(asc_csiro2sww,
    4749      (project.bath_dir,
     
    5860       'zscale': 1,                 #Enhance tsunami
    5961       'fail_on_NaN': False},
    60        #evaluate = True,
     62       #evaluate = True,                 
    6163      verbose = True)
    62 
    6364
    6465#Create the mesh ..
    6566
    6667#cache this, with dependancy on create_mesh.py
    67 mesh, triagle_count = create_mesh(10000000)
    68 
     68mesh, triagle_count = create_mesh(100000)
     69#mesh..export_mesh_file(meshname + '.tsh')
    6970
    7071#meshname = project.meshname + '_%d.msh' %resolution
     
    7475#Setup domain
    7576
    76 domain = cache(pmesh_instance_to_domain_instance, (mesh, Domain),
    77                #dependencies = [meshname],                     
    78                verbose = True)               
     77domain = pmesh_instance_to_domain_instance(mesh, Domain)
     78
     79# this was crashing outside of python
     80#domain = cache(pmesh_instance_to_domain_instance, (mesh, Domain),
     81#               #dependencies = [meshname],     
     82#               evaluate = True,                     
     83#               verbose = True)               
    7984
    8085
     
    116121Br = Reflective_boundary(domain)
    117122Bt = Transmissive_boundary(domain)
    118 Bd = Dirichlet_boundary([tide,0,0])
     123Bd = Dirichlet_boundary([project.tide,0,0])
    119124Bw = Time_boundary(domain=domain,
    120125                   f=lambda t: [(60<t<660)*4, 0, 0])
Note: See TracChangeset for help on using the changeset viewer.