Changeset 3636


Ignore:
Timestamp:
Sep 20, 2006, 4:22:56 PM (18 years ago)
Author:
ole
Message:

Work on Okushiri validation (towards recreating 18 Aug 2005 results)

Location:
anuga_validation/okushiri_2005
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_validation/okushiri_2005/create_mesh.py

    r3535 r3636  
    113113   
    114114#
    115 
    116     #Fine resolution. Needs alpha >= 0.001   
     115   
    117116    base_resolution = 1
    118117
     
    122121    mid = m.addRegionEN(point_mbottom[0]+.1, ybottom+.1)
    123122    mid.setMaxArea(0.0005*base_resolution)
    124    
     123 
    125124
    126125    inner = m.addRegionEN(island_3[0]+.1, island_3[1]+.1)
    127     inner.setMaxArea(0.00007*base_resolution)
     126    inner.setMaxArea(0.0003*base_resolution)
    128127
    129128
    130129    gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
    131     gulleys.setMaxArea(0.00002*base_resolution)
     130    gulleys.setMaxArea(0.00005*base_resolution)
    132131
    133132
     133    # From r 1709 11 August 2005
     134    #base_resolution = 100
     135    #
     136    #ocean = m.addRegionEN(xleft+.1, ybottom+.1)
     137    #ocean.setMaxArea(0.01*base_resolution)
     138    #
     139    #mid = m.addRegionEN(point_mbottom[0]+.1, ybottom+.1)
     140    #mid.setMaxArea(0.001*base_resolution)
    134141
     142    #inner = m.addRegionEN(island_3[0]+.1, island_3[1]+.1)
     143    #inner.setMaxArea(0.0001*base_resolution)
     144
     145
     146    #gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
     147    #gulleys.setMaxArea(0.00001*base_resolution)       
     148   
     149   
    135150    m.generateMesh('pzq28.0za1000000a')
    136151
  • anuga_validation/okushiri_2005/extract_timeseries.py

    r3514 r3636  
    44import sys
    55from os import sep
    6 from caching import cache
     6from anuga.caching import cache
    77
    88sys.path.append('..'+sep+'..'+sep)
     
    5151
    5252
    53 from anuga.pyvolution.util import file_function
     53from anuga.abstract_2d_finite_volumes.util import file_function
    5454from anuga.utilities.numerical_tools import ensure_numeric
    5555gauge_values = [ensure_numeric(input_stage),
Note: See TracChangeset for help on using the changeset viewer.