Changeset 1708


Ignore:
Timestamp:
Aug 11, 2005, 3:15:54 PM (18 years ago)
Author:
ole
Message:

Added refined region around gulleys to mesh

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  
    7979    #Localised refined area for gulleys
    8080    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       
    8693    dict['segments'] = [[0,1], [1,2], [2,3], [3,0]]
    8794    dict['segment_tags'] = ['', '', '', '']
     
    8996
    9097   
    91     base_resolution = 1
     98    base_resolution = 100
    9299
    93100    ocean = m.addRegionEN(xleft+1, ybottom+1)
     
    101108    inner.setMaxArea(0.0001*base_resolution)
    102109
    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)
    104112    gulleys.setMaxArea(0.00001*base_resolution)   
    105113   
  • inundation/ga/storm_surge/validation/LWRU2/lwru2.py

    r1707 r1708  
    8181
    8282try:
    83     #raise Exception
     83    raise Exception
    8484    fid = open(picklefile)
    8585    print 'Read pickled domain'
     
    169169for t in domain.evolve(yieldstep = 0.05, finaltime = 22.5):
    170170    domain.write_time()
    171     Stage_cv = domain.quantities['stage'].centroid_values
    172     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)
    174174
    175175print 'That took %.2f seconds' %(time.time()-t0)
Note: See TracChangeset for help on using the changeset viewer.