Changeset 8286
- Timestamp:
- Dec 15, 2011, 10:51:25 PM (11 years ago)
- Location:
- trunk/anuga_core/source
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/anuga_core/source/anuga/file/sww.py
r8281 r8286 703 703 #outfile.createDimension('numbers_in_range', 2) 704 704 705 print 'store parallel data'705 #print 'store parallel data' 706 706 outfile.number_of_global_triangles = number_of_global_triangles 707 707 outfile.number_of_global_nodes = number_of_global_nodes … … 712 712 outfile.createVariable('tri_full_flag', netcdf_int, ('number_of_volumes',)) 713 713 714 print tri_l2g.shape715 print tri_l2g716 print outfile.variables['tri_l2g'].shape714 #print tri_l2g.shape 715 #print tri_l2g 716 #print outfile.variables['tri_l2g'].shape 717 717 718 718 outfile.variables['tri_l2g'][:] = tri_l2g.astype(num.int32) 719 719 720 print node_l2g.shape721 print node_l2g722 print outfile.variables['node_l2g'].shape720 #print node_l2g.shape 721 #print node_l2g 722 #print outfile.variables['node_l2g'].shape 723 723 724 724 outfile.variables['node_l2g'][:] = node_l2g.astype(num.int32) 725 725 726 print tri_full_flag.shape727 print tri_full_flag728 print outfile.variables['tri_full_flag'].shape726 #print tri_full_flag.shape 727 #print tri_full_flag 728 #print outfile.variables['tri_full_flag'].shape 729 729 730 730 outfile.variables['tri_full_flag'][:] = tri_full_flag.astype(num.int32) -
trunk/anuga_core/source/anuga_parallel/run_parallel_sw_merimbula.py
r8284 r8286 42 42 #-------------------------------------------------------------------------- 43 43 44 #mesh_filename = "merimbula_10785_1.tsh" ; x0 = 756000.0 ; x1 = 756500.044 mesh_filename = "merimbula_10785_1.tsh" ; x0 = 756000.0 ; x1 = 756500.0 45 45 #mesh_filename = "merimbula_43200.tsh" ; x0 = 756000.0 ; x1 = 756500.0 46 46 #mesh_filename = "test-100.tsh" ; x0 = 0.25 ; x1 = 0.5 47 mesh_filename = "test-20.tsh" ; x0 = 250.0 ; x1 = 350.047 #mesh_filename = "test-20.tsh" ; x0 = 250.0 ; x1 = 350.0 48 48 yieldstep = 20 49 finaltime = 4049 finaltime = 1000 50 50 verbose = True 51 51 … … 85 85 domain = create_domain_from_file(mesh_filename) 86 86 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)) 88 88 89 print domain.statistics()90 print domain.get_extent()91 print domain.get_extent(absolute=True)92 print domain.geo_reference89 #print domain.statistics() 90 #print domain.get_extent() 91 #print domain.get_extent(absolute=True) 92 #print domain.geo_reference 93 93 else: 94 94 domain = None … … 109 109 domain.set_name('merimbula') 110 110 111 for p in range(numprocs):112 if myid == p:113 print 'P%d'%p114 print domain.get_extent()115 print domain.get_extent(absolute=True)116 print domain.geo_reference117 print domain.s2p_map118 print domain.p2s_map119 print domain.tri_l2g120 print domain.node_l2g121 else:122 pass123 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() 125 125 126 126
Note: See TracChangeset
for help on using the changeset viewer.