Changeset 1708
- Timestamp:
- Aug 11, 2005, 3:15:54 PM (18 years ago)
- Location:
- inundation/ga/storm_surge/validation/LWRU2
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
inundation/ga/storm_surge/validation/LWRU2/create_mesh.py
r1707 r1708 79 79 #Localised refined area for gulleys 80 80 dict = {} 81 dict['points'] = [ [4.5, 1.8], 82 [5.0, 1.8], 83 [5.0, 2.5], 84 [4.5, 2.5] ] 85 81 xl = 4.8 82 xr = 5.3 83 yb = 1.6 84 yt = 2.3 85 p0 = [xl, yb] 86 p1 = [xr, yb] 87 p2 = [xr, yt] 88 p3 = [xl, yt] 89 90 dict['points'] = [p0, p1, p2, p3] 91 92 86 93 dict['segments'] = [[0,1], [1,2], [2,3], [3,0]] 87 94 dict['segment_tags'] = ['', '', '', ''] … … 89 96 90 97 91 base_resolution = 1 98 base_resolution = 100 92 99 93 100 ocean = m.addRegionEN(xleft+1, ybottom+1) … … 101 108 inner.setMaxArea(0.0001*base_resolution) 102 109 103 gulleys = m.addRegionEN(dict['points'][0]+0.1, dict['points'][1]+0.1) 110 111 gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1) 104 112 gulleys.setMaxArea(0.00001*base_resolution) 105 113 -
inundation/ga/storm_surge/validation/LWRU2/lwru2.py
r1707 r1708 81 81 82 82 try: 83 #raise Exception83 raise Exception 84 84 fid = open(picklefile) 85 85 print 'Read pickled domain' … … 169 169 for t in domain.evolve(yieldstep = 0.05, finaltime = 22.5): 170 170 domain.write_time() 171 Stage_cv = domain.quantities['stage'].centroid_values172 print Stage_cv[0]173 print 'Function value',function(t)171 #Stage_cv = domain.quantities['stage'].centroid_values 172 #print Stage_cv[0] 173 #print 'Function value',function(t) 174 174 175 175 print 'That took %.2f seconds' %(time.time()-t0)
Note: See TracChangeset
for help on using the changeset viewer.