Changeset 5092


Ignore:
Timestamp:
Feb 29, 2008, 10:06:48 AM (17 years ago)
Author:
duncan
Message:

updated simulations

Location:
anuga_work/development/Hinwood_2008
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/Hinwood_2008/boundary.csv

    r5076 r5092  
    1 time sec,stage m,"X velocity, m/s","Y velocity, m/s"
    2 0,0.2,0.23,0.01
    3 0.5,0.205,0.23,0.01
    4 1,0.21,0.23,0.01
    5 1.5,0.215,0.23,0.01
    6 2,0.22,0.23,0.01
    7 2.5,0.225,0.23,0.01
    8 3,0.23,0.23,0.01
    9 3.5,0.235,0.23,0.01
    10 4,0.24,0.23,0.01
    11 4.5,0.245,0.23,0.01
    12 5,0.25,0.23,0.01
    13 5.5,0.255,0.23,0.01
    14 6,0.26,0.23,0.01
    15 6.5,0.265,0.23,0.01
    16 7,0.27,0.23,0.01
    17 7.5,0.27,0.23,0.01
    18 8,0.265,0.23,0.01
    19 8.5,0.26,0.23,0.01
    20 9,0.255,0.23,0.01
    21 9.5,0.25,0.23,0.01
    22 10,0.245,0.23,0.01
    23 10.5,0.24,0.23,0.01
    24 11,0.235,0.23,0.01
    25 11.5,0.23,0.23,0.01
    26 12,0.225,0.23,0.01
    27 12.5,0.22,0.23,0.01
    28 13,0.215,0.23,0.01
    29 13.5,0.21,0.23,0.01
    30 14,0.205,0.23,0.01
    31 14.5,0.2,0.23,0.01
    32 15,0.195,0.23,0.01
     1time sec,stage m,X velocity m/s,Y velocity m/s
     20,0.4,0,0.01
     30.1,0.401996668,0.005990005,0.01
     40.2,0.403973387,0.01192016,0.01
     50.3,0.405910404,0.017731212,0.01
     60.4,0.407788367,0.023365101,0.01
     70.5,0.409588511,0.028765532,0.01
     80.6,0.411292849,0.033878548,0.01
     90.7,0.412884354,0.038653061,0.01
     100.8,0.414347122,0.043041365,0.01
     110.9,0.415666538,0.046999615,0.01
     121,0.41682942,0.050488259,0.01
     131.1,0.417824147,0.053472442,0.01
     141.2,0.418640782,0.055922345,0.01
     151.3,0.419271164,0.057813491,0.01
     161.4,0.419708995,0.059126984,0.01
     171.5,0.4199499,0.059849699,0.01
     181.6,0.419991472,0.059974416,0.01
     191.7,0.419833296,0.059499889,0.01
     201.8,0.419476953,0.058430858,0.01
     211.9,0.418926002,0.056778005,0.01
     222,0.418185949,0.054557846,0.01
     232.1,0.417264187,0.051792562,0.01
     242.2,0.416169928,0.048509784,0.01
     252.3,0.414914104,0.044742313,0.01
     262.4,0.413509264,0.040527791,0.01
     272.5,0.411969443,0.035908329,0.01
     282.6,0.410310027,0.030930082,0.01
     292.7,0.408547598,0.025642793,0.01
     302.8,0.406699763,0.020099289,0.01
     312.9,0.404784987,0.01435496,0.01
     323,0.4028224,0.0084672,0.01
     333.1,0.400831613,0.00249484,0.01
     343.2,0.398832517,-0.003502449,0.01
     353.3,0.396845086,-0.009464742,0.01
     363.4,0.394889178,-0.015332466,0.01
     373.5,0.392984335,-0.021046994,0.01
     383.6,0.391149591,-0.026551227,0.01
     393.7,0.389403277,-0.031790168,0.01
     403.8,0.387762842,-0.036711473,0.01
     413.9,0.386244677,-0.04126597,0.01
     424,0.38486395,-0.04540815,0.01
     434.1,0.383634458,-0.049096627,0.01
     444.2,0.382568485,-0.052294546,0.01
     454.3,0.381676681,-0.054969956,0.01
     464.4,0.380967959,-0.057096124,0.01
     474.5,0.380449398,-0.058651807,0.01
     484.6,0.38012618,-0.05962146,0.01
     494.7,0.380001535,-0.059995395,0.01
     504.8,0.380076708,-0.059769877,0.01
     514.9,0.380350948,-0.058947157,0.01
     525,0.380821515,-0.057535456,0.01
     535.1,0.381483706,-0.055548881,0.01
     545.2,0.382330907,-0.053007279,0.01
     555.3,0.383354651,-0.049936047,0.01
     565.4,0.38454471,-0.046365869,0.01
     575.5,0.385889193,-0.04233242,0.01
     585.6,0.387374667,-0.037875998,0.01
     595.7,0.388986289,-0.033041133,0.01
     605.8,0.390707956,-0.027876131,0.01
     615.9,0.392522467,-0.0224326,0.01
     626,0.39441169,-0.01676493,0.01
     636.1,0.39635675,-0.01092975,0.01
     646.2,0.398338212,-0.004985364,0.01
     656.3,0.400336278,0.001008834,0.01
     666.4,0.402330984,0.006992952,0.01
     676.5,0.4043024,0.012907199,0.01
     686.6,0.406230827,0.018692482,0.01
     696.7,0.408096998,0.024290995,0.01
     706.8,0.409882267,0.029646801,0.01
     716.9,0.411568795,0.034706386,0.01
     727,0.413139732,0.039419196,0.01
     737.1,0.414579381,0.043738142,0.01
     747.2,0.415873357,0.047620072,0.01
     757.3,0.417008732,0.051026197,0.01
     767.4,0.417974162,0.053922486,0.01
     777.5,0.41876,0.056279999,0.01
     787.6,0.419358393,0.05807518,0.01
     797.7,0.419763365,0.059290094,0.01
     807.8,0.419970867,0.059912601,0.01
     817.9,0.419978827,0.059936481,0.01
     828,0.419787165,0.059361495,0.01
     838.1,0.419397796,0.058193389,0.01
     848.2,0.418814611,0.056443833,0.01
     858.3,0.418043437,0.05413031,0.01
     868.4,0.417091978,0.051275934,0.01
     878.5,0.415969742,0.047909227,0.01
     888.6,0.414687942,0.044063826,0.01
     898.7,0.413259385,0.039778154,0.01
     908.8,0.411698344,0.035095032,0.01
     918.9,0.410020417,0.030061251,0.01
     929,0.40824237,0.024727109,0.01
     939.1,0.406381967,0.019145902,0.01
     949.2,0.404457798,0.013373395,0.01
     959.3,0.402489088,0.007467265,0.01
     969.4,0.400495509,0.001486526,0.01
     979.5,0.398496978,-0.004509067,0.01
     989.6,0.396513464,-0.010459607,0.01
     999.7,0.394564787,-0.016305638,0.01
     1009.8,0.392670417,-0.021988748,0.01
     1019.9,0.390849282,-0.027452154,0.01
     10210,0.389119578,-0.032641267,0.01
     10310.1,0.387498587,-0.037504239,0.01
     10410.2,0.386002506,-0.041992481,0.01
     10510.3,0.384646284,-0.046061149,0.01
     10610.4,0.383443471,-0.049669588,0.01
     10710.5,0.382406085,-0.052781746,0.01
     10810.6,0.381544492,-0.055366525,0.01
     10910.7,0.3808673,-0.057398101,0.01
     11010.8,0.380381275,-0.058856174,0.01
     11110.9,0.380091275,-0.059726175,0.01
     11211,0.380000196,-0.059999412,0.01
     11311.1,0.380108948,-0.059673155,0.01
     11411.2,0.380416445,-0.058750664,0.01
     11511.3,0.380919615,-0.057241155,0.01
     11611.4,0.381613429,-0.055159712,0.01
     11711.5,0.382490957,-0.05252713,0.01
     11811.6,0.383543428,-0.049369716,0.01
     11911.7,0.384760328,-0.045719015,0.01
     12011.8,0.386129498,-0.041611505,0.01
     12111.9,0.387637258,-0.037088227,0.01
     12212,0.389268542,-0.032194375,0.01
     12312.1,0.391007051,-0.026978848,0.01
     12412.2,0.392835414,-0.021493757,0.01
     12512.3,0.394735364,-0.015793907,0.01
     12612.4,0.396687916,-0.009936251,0.01
     12712.5,0.398673562,-0.003979314,0.01
     12812.6,0.400672461,0.002017383,0.01
     12912.7,0.402664641,0.007993922,0.01
     13012.8,0.404630197,0.01389059,0.01
     13112.9,0.406549489,0.019648466,0.01
     13213,0.408403341,0.025210022,0.01
     13313.1,0.410173229,0.030519688,0.01
     13413.2,0.41184147,0.035524411,0.01
     13513.3,0.413391395,0.040174186,0.01
     13613.4,0.414807518,0.044422553,0.01
     13713.5,0.416075689,0.048227066,0.01
     13813.6,0.417183236,0.051549709,0.01
     13913.7,0.418119095,0.054357285,0.01
     14013.8,0.418873913,0.05662174,0.01
     14113.9,0.41944015,0.05832045,0.01
     14214,0.419812147,0.059436441,0.01
     14314.1,0.419986188,0.059958563,0.01
     14414.2,0.419960533,0.059881599,0.01
     14514.3,0.419735439,0.059206318,0.01
     14614.4,0.419313156,0.057939467,0.01
     14714.5,0.418697901,0.056093703,0.01
     14814.6,0.417895823,0.05368747,0.01
     14914.7,0.416914937,0.05074481,0.01
     15014.8,0.415765041,0.047295124,0.01
     15114.9,0.414457627,0.043372881,0.01
     15215,0.413005757,0.03901727,0.01
  • anuga_work/development/Hinwood_2008/create_mesh.py

    r5076 r5092  
    55from anuga.pmesh.mesh import *
    66from anuga.coordinate_transforms.geo_reference import Geo_reference
     7
     8xslope = 5.0 # Distance between the boundary ande the start of the slope
    79
    810
     
    1517    """
    1618    #Basic geometry
     19    global xslope
    1720   
    1821    xright  = 19.0
    1922    ybottom = 0
    20     ytop    = 0.45
     23    ytop    = 1.00
    2124    xleft = 0.0
    22     xslope = 4.0
     25    ###xslope = slope
    2326
    2427    #Outline
  • anuga_work/development/Hinwood_2008/run_dam.py

    r5078 r5092  
    4242def elevation_function(x,y):
    4343    from Numeric import zeros, size, Float
    44     slope = 4 ## Bit of a magic Number
    45    
     44
     45    xslope = create_mesh.xslope  #4 ## Bit of a magic Number
     46    print "xslope",xslope
    4647    z = zeros(size(x), Float)
    4748    for i in range(len(x)):
    48         if x[i] < slope:
     49        if x[i] < xslope:
    4950            z[i] = 0.0  #WARNING: the code in prepare_time_boundary
    5051                        # that calc's momentum assumes this is 0.0
    5152        else:
    52             z[i] = (x[i]-slope)*0.1
     53            z[i] = (x[i]-xslope)*(1./16.)
    5354    return z
    5455
     
    124125    #-------------------------------------------------------------------------
    125126
    126     domain.set_quantity('stage', 0.06)
     127    domain.set_quantity('stage', 0.4)
    127128    domain.set_quantity('friction', friction)
    128129    domain.set_quantity('elevation', elevation_function)
     
    151152    print 'finished'
    152153
    153     points = [[2.8,0.225],  #-1.8m from SWL
    154               [5.1,0.225],  #0.5m from SWL
    155               [6.6,0.225],  #2m from SWL
    156               [6.95,0.255], #2.35m from SWL
    157               [7.6,0.255],  #3m from SWL
    158               [8.2,0.255],  #3.5m from SWL
    159               [9.2,0.255]  #4.5m from SWL
     154    points = [[2.8,0.5],  #-1.8m from SWL
     155              [5.1,0.5],  #0.5m from SWL
     156              [6.6,0.5],  #2m from SWL
     157              [6.95,0.5], #2.35m from SWL
     158              [7.6,0.5],  #3m from SWL
     159              [8.2,0.5],  #3.5m from SWL
     160              [9.2,0.5]  #4.5m from SWL
    160161              ]
    161162
     
    177178if __name__ == "__main__":
    178179    main( is_trial_run = True,
    179          outputdir_name='Hinwood_draft')
     180         outputdir_name='Hinwood_low_stage_low_velocity_draft')
Note: See TracChangeset for help on using the changeset viewer.