Changeset 2001


Ignore:
Timestamp:
Nov 4, 2005, 1:15:36 PM (18 years ago)
Author:
chris
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • inundation/analytical solutions/Analytical_solution_circular_hydraulic_jump.py

    r1995 r2001  
    8080            if distance[j] <= r:
    8181                if distance[j+1] > r:
    82                     z[i] = bed[0]
    83                     #bed[j] + (bed[j+1] - bed[j])/(distance[j+1] - distance[j])*(r - distance[j])
     82                    z[i] = bed[j] + (bed[j+1] - bed[j])/(distance[j+1] - distance[j])*(r - distance[j])
    8483    return z
     84
    8585domain.set_quantity('elevation', bed_z)
    8686
     
    123123
    124124
     125# domain.initialise_visualiser()
     126# domain.visualiser.coloring['stage'] = True
     127# domain.visualiser.scale_z['stage'] = 2.0
     128# domain.visualiser.scale_z['elevation'] = 0.05
    125129
    126130
    127 domain.initialise_visualiser()
    128     #domain.visualiser.coloring['stage'] = True
    129 domain.visualiser.scale_z['stage'] = 2.0
    130 domain.visualiser.scale_z['elevation'] = 0.05
    131 
    132 
    133 from realtime_visualisation_new import Visualiser
    134 #vxmom = Visualiser(domain,title='xmomentum',scale_z=10.0)
    135 #vymom = Visualiser(domain,title='ymomentum',scale_z=10.0)
     131# from realtime_visualisation_new import Visualiser
     132# vxmom = Visualiser(domain,title='xmomentum',scale_z=10.0)
     133# vymom = Visualiser(domain,title='ymomentum',scale_z=10.0)
    136134
    137135w = domain.quantities['stage']
     
    141139import time
    142140
    143 f = open("output.txt","w")
     141f = open("circular_hydraulic_jump_true.txt","w")
    144142
    145143t0 = time.time()
    146 for t in domain.evolve(yieldstep = .01, finaltime = 0.1):
     144for t in domain.evolve(yieldstep = 1.0, finaltime = 1000.):
    147145    domain.write_time()
    148146
     
    150148    radial_momentum = domain.create_quantity_from_expression(exp)
    151149
    152     print 'outer stage      ', w.get_values(location='centroids',
    153                                             indices=[typical_outer[0]])
    154     print '      radial mom ', \
    155           radial_momentum.get_values(location='centroids',
    156                                      indices=[typical_outer[0]])
     150    outer_stage = w.get_values(location='centroids', indices=[typical_outer[0]])
     151    outer_radial_mom = radial_momentum.get_values(location='centroids', indices=[typical_outer[0]])
     152   
     153    inner_stage = w.get_values(location='centroids', indices=[typical_inner[0]])
     154    inner_radial_mom = radial_momentum.get_values(location='centroids', indices=[typical_inner[0]])
    157155
    158     print 'inner stage      ', w.get_values(location='centroids',
    159                                             indices=[typical_inner[0]])
    160     print '      radial mom ', \
    161           radial_momentum.get_values(location='centroids',
    162                                      indices=[typical_inner[0]])
     156    print inner_stage
     157    print inner_radial_mom
     158    print outer_stage
     159    print outer_radial_mom
    163160
    164 
    165     f.write('time = %25.15e wall clock time %g \n' % (domain.time, time.time()))
    166     f.write('%g \n' % w.get_values(location='centroids',
    167                                    indices=[typical_outer[0]])[0])
     161#    f.write('time = %25.15e wall clock time %g \n' % (domain.time, time.time()))
     162    f.write('%10.3f %25.15e %25.15e %25.15e %25.15e \n' % (domain.time, inner_stage, inner_radial_mom, outer_stage, outer_radial_mom))
    168163
    169164
    170165f.close()
    171166
    172 
    173     #vxmom.update_quantity('xmomentum')
    174     #vymom.update_quantity('ymomentum')
    175 
    176     #Old version
    177     #print 'outer stage      ',domain.quantities['stage'].get_values(location='centroids',indices=[typical_outer[0]])
    178     #print '      radial mom ',sqrt(domain.quantities['xmomentum'].get_values(location='centroids',indices=[typical_outer[0]])[0]**2 +
    179     #                               domain.quantities['ymomentum'].get_values(location='centroids',indices=[typical_outer[0]])[0]**2)
    180     #print 'inner stage      ',domain.quantities['stage'].get_values(location='centroids',indices=[typical_inner[0]])
    181     #print '      radial mom ',sqrt(domain.quantities['xmomentum'].get_values(location='centroids',indices=[typical_inner[0]])[0]**2 +
    182     #                               domain.quantities['ymomentum'].get_values(location='centroids',indices=[typical_inner[0]])[0]**2)
    183 
    184 
    185167print 'That took %.2f seconds' %(time.time()-t0)
Note: See TracChangeset for help on using the changeset viewer.