Ignore:
Timestamp:
Oct 14, 2005, 5:22:01 PM (19 years ago)
Author:
ole
Message:

Suggested using the new expression functionality for printing the radial momentum.
I have checked the calculated results are the same.
The old code is still there but has been commented out.

File:
1 edited

Legend:

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

    r1876 r1929  
    1717from math import pi, sqrt
    1818from mesh_factory import strang_mesh
     19
    1920
    2021#---------
     
    134135#vymom = Visualiser(domain,title='ymomentum',scale_z=10.0)
    135136
     137w = domain.quantities['stage']
    136138
    137139#----------
     
    142144for t in domain.evolve(yieldstep = .01, finaltime = 10):
    143145    domain.write_time()
     146
     147    exp = '(xmomentum**2 + ymomentum**2)**0.5'
     148    radial_momentum = domain.create_quantity_from_expression(exp)
     149
     150    print 'outer stage      ', w.get_values(location='centroids',
     151                                            indices=[typical_outer[0]])
     152    print '      radial mom ', \
     153          radial_momentum.get_values(location='centroids',
     154                                     indices=[typical_outer[0]])
     155
     156    print 'inner stage      ', w.get_values(location='centroids',
     157                                            indices=[typical_inner[0]])
     158    print '      radial mom ', \
     159          radial_momentum.get_values(location='centroids',
     160                                     indices=[typical_inner[0]])   
     161   
     162
     163   
     164
     165   
     166   
    144167    #vxmom.update_quantity('xmomentum')
    145168    #vymom.update_quantity('ymomentum')
    146     print 'outer stage      ',domain.quantities['stage'].get_values(location='centroids',indices=[typical_outer[0]])
    147     print '      radial mom ',sqrt(domain.quantities['xmomentum'].get_values(location='centroids',indices=[typical_outer[0]])[0]**2 +
    148                                    domain.quantities['ymomentum'].get_values(location='centroids',indices=[typical_outer[0]])[0]**2)
    149     print 'inner stage      ',domain.quantities['stage'].get_values(location='centroids',indices=[typical_inner[0]])
    150     print '      radial mom ',sqrt(domain.quantities['xmomentum'].get_values(location='centroids',indices=[typical_inner[0]])[0]**2 +
    151                                    domain.quantities['ymomentum'].get_values(location='centroids',indices=[typical_inner[0]])[0]**2)
     169
     170    #Old version
     171    #print 'outer stage      ',domain.quantities['stage'].get_values(location='centroids',indices=[typical_outer[0]])
     172    #print '      radial mom ',sqrt(domain.quantities['xmomentum'].get_values(location='centroids',indices=[typical_outer[0]])[0]**2 +
     173    #                               domain.quantities['ymomentum'].get_values(location='centroids',indices=[typical_outer[0]])[0]**2)
     174    #print 'inner stage      ',domain.quantities['stage'].get_values(location='centroids',indices=[typical_inner[0]])
     175    #print '      radial mom ',sqrt(domain.quantities['xmomentum'].get_values(location='centroids',indices=[typical_inner[0]])[0]**2 +
     176    #                               domain.quantities['ymomentum'].get_values(location='centroids',indices=[typical_inner[0]])[0]**2)
    152177
    153178
Note: See TracChangeset for help on using the changeset viewer.