Changeset 6713 for anuga_work/development/anuga_1d/test_discontinuous_z.py
 Timestamp:
 Apr 3, 2009, 12:34:47 PM (14 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

anuga_work/development/anuga_1d/test_discontinuous_z.py
r6698 r6713 85 85 domain.set_quantity('elevation',elevation_arbitrary2) 86 86 domain.set_boundary({'exterior':Reflective_boundary(domain)}) 87 X=domain.vertices 88 C=domain.centroids 87 X=domain.vertices.flat 88 C=domain.centroids.flat 89 89 90 90 yieldstep = 1.0 … … 148 148 domain.write_time() 149 149 print "integral", domain.quantities['height'].get_integral() 150 StageQ=domain.quantities['stage'].vertex_values 150 StageQ=domain.quantities['stage'].vertex_values.flat 151 151 StageUp=domain.quantities['stage'].explicit_update 152 MomentumQ=domain.quantities['xmomentum'].vertex_values 152 MomentumQ=domain.quantities['xmomentum'].vertex_values.flat 153 153 MomentumUp=domain.quantities['xmomentum'].explicit_update 154 154 155 ElevationQ=domain.quantities['elevation'].vertex_values 156 VelocityQ=domain.quantities['velocity'].vertex_values 155 ElevationQ=domain.quantities['elevation'].vertex_values.flat 156 VelocityQ=domain.quantities['velocity'].vertex_values.flat 157 157 158 158 print "momentum", MomentumQ … … 161 161 #StageUp=domain.quantities['stage'].explicit_update 162 162 print "MomUpdates", MomentumUp 163 163 164 hold(False) 164 165 plot1 = subplot(311) 166 #plot(X,StageQ) 165 167 plot(X,StageQ, X,ElevationQ) 166 168 plot1.set_ylim([1.0,8.0]) … … 178 180 plot3 = subplot(313) 179 181 plot(X,VelocityQ) 180 #plot 2.set_ylim([5.0,30.0])182 #plot3.set_ylim([5.0,30.0]) 181 183 legend( ('Numerical Solution', 'for velocity'), 'upper right', shadow=False) 182 184 xlabel('Position')
Note: See TracChangeset
for help on using the changeset viewer.