Changeset 2001
- Timestamp:
- Nov 4, 2005, 1:15:36 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
inundation/analytical solutions/Analytical_solution_circular_hydraulic_jump.py
r1995 r2001 80 80 if distance[j] <= r: 81 81 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]) 84 83 return z 84 85 85 domain.set_quantity('elevation', bed_z) 86 86 … … 123 123 124 124 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 125 129 126 130 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) 136 134 137 135 w = domain.quantities['stage'] … … 141 139 import time 142 140 143 f = open(" output.txt","w")141 f = open("circular_hydraulic_jump_true.txt","w") 144 142 145 143 t0 = time.time() 146 for t in domain.evolve(yieldstep = .01, finaltime = 0.1):144 for t in domain.evolve(yieldstep = 1.0, finaltime = 1000.): 147 145 domain.write_time() 148 146 … … 150 148 radial_momentum = domain.create_quantity_from_expression(exp) 151 149 152 print 'outer stage ', w.get_values(location='centroids',153 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]]) 157 155 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 163 160 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)) 168 163 169 164 170 165 f.close() 171 166 172 173 #vxmom.update_quantity('xmomentum')174 #vymom.update_quantity('ymomentum')175 176 #Old version177 #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 185 167 print 'That took %.2f seconds' %(time.time()-t0)
Note: See TracChangeset
for help on using the changeset viewer.