Changeset 8286


Ignore:
Timestamp:
Dec 15, 2011, 10:51:25 PM (8 years ago)
Author:
steve
Message:

working version of sww_merge

Location:
trunk/anuga_core/source
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/anuga_core/source/anuga/file/sww.py

    r8281 r8286  
    703703        #outfile.createDimension('numbers_in_range', 2)
    704704
    705         print 'store parallel data'
     705        #print 'store parallel data'
    706706        outfile.number_of_global_triangles = number_of_global_triangles
    707707        outfile.number_of_global_nodes = number_of_global_nodes
     
    712712        outfile.createVariable('tri_full_flag', netcdf_int, ('number_of_volumes',))
    713713
    714         print tri_l2g.shape
    715         print tri_l2g
    716         print outfile.variables['tri_l2g'].shape
     714        #print tri_l2g.shape
     715        #print tri_l2g
     716        #print outfile.variables['tri_l2g'].shape
    717717
    718718        outfile.variables['tri_l2g'][:] = tri_l2g.astype(num.int32)
    719719
    720         print node_l2g.shape
    721         print node_l2g
    722         print outfile.variables['node_l2g'].shape
     720        #print node_l2g.shape
     721        #print node_l2g
     722        #print outfile.variables['node_l2g'].shape
    723723
    724724        outfile.variables['node_l2g'][:] = node_l2g.astype(num.int32)
    725725
    726         print tri_full_flag.shape
    727         print tri_full_flag
    728         print outfile.variables['tri_full_flag'].shape
     726        #print tri_full_flag.shape
     727        #print tri_full_flag
     728        #print outfile.variables['tri_full_flag'].shape
    729729
    730730        outfile.variables['tri_full_flag'][:] = tri_full_flag.astype(num.int32)
  • trunk/anuga_core/source/anuga_parallel/run_parallel_sw_merimbula.py

    r8284 r8286  
    4242#--------------------------------------------------------------------------
    4343
    44 #mesh_filename = "merimbula_10785_1.tsh" ; x0 = 756000.0 ; x1 = 756500.0
     44mesh_filename = "merimbula_10785_1.tsh" ; x0 = 756000.0 ; x1 = 756500.0
    4545#mesh_filename = "merimbula_43200.tsh"   ; x0 = 756000.0 ; x1 = 756500.0
    4646#mesh_filename = "test-100.tsh" ; x0 = 0.25 ; x1 = 0.5
    47 mesh_filename = "test-20.tsh" ; x0 = 250.0 ; x1 = 350.0
     47#mesh_filename = "test-20.tsh" ; x0 = 250.0 ; x1 = 350.0
    4848yieldstep = 20
    49 finaltime = 40
     49finaltime = 1000
    5050verbose = True
    5151
     
    8585    domain = create_domain_from_file(mesh_filename)
    8686    domain.set_quantity('stage', Set_Stage(x0, x1, 2.0))
    87     domain.set_quantity('elevation', Set_Elevation(500.0))
     87    #domain.set_quantity('elevation', Set_Elevation(500.0))
    8888
    89     print domain.statistics()
    90     print domain.get_extent()
    91     print domain.get_extent(absolute=True)
    92     print domain.geo_reference
     89    #print domain.statistics()
     90    #print domain.get_extent()
     91    #print domain.get_extent(absolute=True)
     92    #print domain.geo_reference
    9393else:
    9494    domain = None
     
    109109domain.set_name('merimbula')
    110110
    111 for p in range(numprocs):
    112     if myid == p:
    113         print 'P%d'%p
    114         print domain.get_extent()
    115         print domain.get_extent(absolute=True)
    116         print domain.geo_reference
    117         print domain.s2p_map
    118         print domain.p2s_map
    119         print domain.tri_l2g
    120         print domain.node_l2g
    121     else:
    122         pass
    123    
    124     barrier()
     111#for p in range(numprocs):
     112#    if myid == p:
     113#        print 'P%d'%p
     114#        print domain.get_extent()
     115#        print domain.get_extent(absolute=True)
     116#        print domain.geo_reference
     117#        print domain.s2p_map
     118#        print domain.p2s_map
     119#        print domain.tri_l2g
     120#        print domain.node_l2g
     121#    else:
     122#        pass
     123#
     124#    barrier()
    125125
    126126
Note: See TracChangeset for help on using the changeset viewer.