Changeset 7907


Ignore:
Timestamp:
Aug 2, 2010, 2:11:03 PM (9 years ago)
Author:
steve
Message:

Problem with plotting vertices. Had to use flat structure

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/anuga_work/development/sudi/sw_1d/shock_detector_shv/dam_h_2sides_RUN.py

    r7906 r7907  
    5757C=domain.centroids
    5858import time
    59 yieldstep=finaltime=1.0
     59yieldstep=finaltime=5.0
    6060t0=time.time()
    6161i=1
    6262
    63 while finaltime < 1.01:
     63while finaltime < 5.01:
    6464    for t in domain.evolve(yieldstep=yieldstep, finaltime=finaltime):
    6565        domain.write_time()
     
    7878    #print "u_error %.10f" %(u_error)
    7979    #print 'That took %.2f seconds' %(time.time()-t0)
    80     X = domain.vertices
    81     StageQ = domain.quantities['stage'].vertex_values
    82     XmomQ = domain.quantities['xmomentum'].vertex_values
    83     VelQ = domain.quantities['velocity'].vertex_values
    84     BedQ = domain.quantities['elevation'].vertex_values
     80    X = domain.vertices.flat
     81    StageQ = domain.quantities['stage'].vertex_values.flat
     82    XmomQ = domain.quantities['xmomentum'].vertex_values.flat
     83    VelQ = domain.quantities['velocity'].vertex_values.flat
     84    BedQ = domain.quantities['elevation'].vertex_values.flat
    8585    #h, uh, u = analytical_sol(X.flat, domain.time)
    8686
Note: See TracChangeset for help on using the changeset viewer.