Changeset 3115


Ignore:
Timestamp:
Jun 8, 2006, 2:10:15 PM (18 years ago)
Author:
ole
Message:

Got parallel wollongong going and flagged two issues with FIXME

File:
1 edited

Legend:

Unmodified
Added
Removed
  • production/wollongong_2006/run_flagstaff.py

    r3105 r3115  
    5454print 'I am processor %d of %d on node %s' %(myid, numprocs, processor_name)
    5555
    56 bounding_rectangle = zeros(4, Float) # Buffer for results
    57 
    5856
    5957#------------------------------------------------------------------------------
     
    9492                                    interior_regions=project.interior_regions)
    9593   
    96     ## Add buildings
    97     #dict = importUngenerateFile(project.buildings_filename)
    98     #Segment.set_default_tag('wall') # This should bind to a Reflective boundary
    99     #mesh.addVertsSegs(dict)
    100 
    101  
    10294    # Add buildings
    10395    # This should bind to a Reflective boundary
     
    117109    print domain.statistics()
    118110
    119     domain.set_name(project.basename)
     111    #domain.set_name(project.basename)
    120112    domain.set_datadir(project.outputdir)
    121     domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
    122 
    123     bounding_rectangle[:] = domain.get_extent(absolute=True)
     113    #domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
     114    domain.set_quantities_to_be_stored(None)
     115
     116    domain.set_quantity('elevation',
     117                        filename=project.demname + '.pts',
     118                        use_cache=True,
     119                        verbose=True)   
     120
    124121
    125122    # Subdivide the mesh
     
    161158
    162159
    163 pypar.broadcast(bounding_rectangle, 0)
    164 
    165160# Build the domain for this processor
    166161domain = Parallel_Domain(points, vertices, boundary,
     
    168163                         ghost_recv_dict = ghost_recv_dict)
    169164
     165# FIXME (Ole): Name currently has to be set here to get the processor number
     166# right. It would be easy to build into Parallel_Domain
     167domain.set_name(project.basename)
    170168
    171169#------------------------------------------------------------------------------
     
    175173domain.set_quantity('stage', project.initial_sealevel)
    176174domain.set_quantity('friction', 0.03)
    177 domain.set_quantity('elevation',
    178                     filename=project.demname + '.pts',
    179                     use_cache=True,
    180                     verbose=True)
     175
     176#
     177# FIXME (Ole): This one segfaults which is bad, because set_quantity is
     178# time consuming and should be done here rather than on processor 0
     179#
     180#domain.set_quantity('elevation',
     181#                    filename=project.demname + '.pts',
     182#                    use_cache=True,
     183#                    verbose=True)
    181184
    182185
     
    196199                     'ghost': None})
    197200
    198                            
     201
     202
     203print 'P%d: Ready to evolve. Value of store is %s' %(myid, str(domain.store))
     204
     205
    199206#------------------------------------------------------------------------------
    200207# Evolve system through time
Note: See TracChangeset for help on using the changeset viewer.