- Timestamp:
- Oct 14, 2005, 5:22:01 PM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
inundation/analytical solutions/Analytical_solution_circular_hydraulic_jump.py
r1876 r1929 17 17 from math import pi, sqrt 18 18 from mesh_factory import strang_mesh 19 19 20 20 21 #--------- … … 134 135 #vymom = Visualiser(domain,title='ymomentum',scale_z=10.0) 135 136 137 w = domain.quantities['stage'] 136 138 137 139 #---------- … … 142 144 for t in domain.evolve(yieldstep = .01, finaltime = 10): 143 145 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 144 167 #vxmom.update_quantity('xmomentum') 145 168 #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) 152 177 153 178
Note: See TracChangeset
for help on using the changeset viewer.