Changeset 5092
- Timestamp:
- Feb 29, 2008, 10:06:48 AM (16 years ago)
- Location:
- anuga_work/development/Hinwood_2008
- Files:
-
- 3 edited
-
boundary.csv (modified) (1 diff)
-
create_mesh.py (modified) (2 diffs)
-
run_dam.py (modified) (4 diffs)
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 1 time sec,stage m,X velocity m/s,Y velocity m/s 2 0,0.4,0,0.01 3 0.1,0.401996668,0.005990005,0.01 4 0.2,0.403973387,0.01192016,0.01 5 0.3,0.405910404,0.017731212,0.01 6 0.4,0.407788367,0.023365101,0.01 7 0.5,0.409588511,0.028765532,0.01 8 0.6,0.411292849,0.033878548,0.01 9 0.7,0.412884354,0.038653061,0.01 10 0.8,0.414347122,0.043041365,0.01 11 0.9,0.415666538,0.046999615,0.01 12 1,0.41682942,0.050488259,0.01 13 1.1,0.417824147,0.053472442,0.01 14 1.2,0.418640782,0.055922345,0.01 15 1.3,0.419271164,0.057813491,0.01 16 1.4,0.419708995,0.059126984,0.01 17 1.5,0.4199499,0.059849699,0.01 18 1.6,0.419991472,0.059974416,0.01 19 1.7,0.419833296,0.059499889,0.01 20 1.8,0.419476953,0.058430858,0.01 21 1.9,0.418926002,0.056778005,0.01 22 2,0.418185949,0.054557846,0.01 23 2.1,0.417264187,0.051792562,0.01 24 2.2,0.416169928,0.048509784,0.01 25 2.3,0.414914104,0.044742313,0.01 26 2.4,0.413509264,0.040527791,0.01 27 2.5,0.411969443,0.035908329,0.01 28 2.6,0.410310027,0.030930082,0.01 29 2.7,0.408547598,0.025642793,0.01 30 2.8,0.406699763,0.020099289,0.01 31 2.9,0.404784987,0.01435496,0.01 32 3,0.4028224,0.0084672,0.01 33 3.1,0.400831613,0.00249484,0.01 34 3.2,0.398832517,-0.003502449,0.01 35 3.3,0.396845086,-0.009464742,0.01 36 3.4,0.394889178,-0.015332466,0.01 37 3.5,0.392984335,-0.021046994,0.01 38 3.6,0.391149591,-0.026551227,0.01 39 3.7,0.389403277,-0.031790168,0.01 40 3.8,0.387762842,-0.036711473,0.01 41 3.9,0.386244677,-0.04126597,0.01 42 4,0.38486395,-0.04540815,0.01 43 4.1,0.383634458,-0.049096627,0.01 44 4.2,0.382568485,-0.052294546,0.01 45 4.3,0.381676681,-0.054969956,0.01 46 4.4,0.380967959,-0.057096124,0.01 47 4.5,0.380449398,-0.058651807,0.01 48 4.6,0.38012618,-0.05962146,0.01 49 4.7,0.380001535,-0.059995395,0.01 50 4.8,0.380076708,-0.059769877,0.01 51 4.9,0.380350948,-0.058947157,0.01 52 5,0.380821515,-0.057535456,0.01 53 5.1,0.381483706,-0.055548881,0.01 54 5.2,0.382330907,-0.053007279,0.01 55 5.3,0.383354651,-0.049936047,0.01 56 5.4,0.38454471,-0.046365869,0.01 57 5.5,0.385889193,-0.04233242,0.01 58 5.6,0.387374667,-0.037875998,0.01 59 5.7,0.388986289,-0.033041133,0.01 60 5.8,0.390707956,-0.027876131,0.01 61 5.9,0.392522467,-0.0224326,0.01 62 6,0.39441169,-0.01676493,0.01 63 6.1,0.39635675,-0.01092975,0.01 64 6.2,0.398338212,-0.004985364,0.01 65 6.3,0.400336278,0.001008834,0.01 66 6.4,0.402330984,0.006992952,0.01 67 6.5,0.4043024,0.012907199,0.01 68 6.6,0.406230827,0.018692482,0.01 69 6.7,0.408096998,0.024290995,0.01 70 6.8,0.409882267,0.029646801,0.01 71 6.9,0.411568795,0.034706386,0.01 72 7,0.413139732,0.039419196,0.01 73 7.1,0.414579381,0.043738142,0.01 74 7.2,0.415873357,0.047620072,0.01 75 7.3,0.417008732,0.051026197,0.01 76 7.4,0.417974162,0.053922486,0.01 77 7.5,0.41876,0.056279999,0.01 78 7.6,0.419358393,0.05807518,0.01 79 7.7,0.419763365,0.059290094,0.01 80 7.8,0.419970867,0.059912601,0.01 81 7.9,0.419978827,0.059936481,0.01 82 8,0.419787165,0.059361495,0.01 83 8.1,0.419397796,0.058193389,0.01 84 8.2,0.418814611,0.056443833,0.01 85 8.3,0.418043437,0.05413031,0.01 86 8.4,0.417091978,0.051275934,0.01 87 8.5,0.415969742,0.047909227,0.01 88 8.6,0.414687942,0.044063826,0.01 89 8.7,0.413259385,0.039778154,0.01 90 8.8,0.411698344,0.035095032,0.01 91 8.9,0.410020417,0.030061251,0.01 92 9,0.40824237,0.024727109,0.01 93 9.1,0.406381967,0.019145902,0.01 94 9.2,0.404457798,0.013373395,0.01 95 9.3,0.402489088,0.007467265,0.01 96 9.4,0.400495509,0.001486526,0.01 97 9.5,0.398496978,-0.004509067,0.01 98 9.6,0.396513464,-0.010459607,0.01 99 9.7,0.394564787,-0.016305638,0.01 100 9.8,0.392670417,-0.021988748,0.01 101 9.9,0.390849282,-0.027452154,0.01 102 10,0.389119578,-0.032641267,0.01 103 10.1,0.387498587,-0.037504239,0.01 104 10.2,0.386002506,-0.041992481,0.01 105 10.3,0.384646284,-0.046061149,0.01 106 10.4,0.383443471,-0.049669588,0.01 107 10.5,0.382406085,-0.052781746,0.01 108 10.6,0.381544492,-0.055366525,0.01 109 10.7,0.3808673,-0.057398101,0.01 110 10.8,0.380381275,-0.058856174,0.01 111 10.9,0.380091275,-0.059726175,0.01 112 11,0.380000196,-0.059999412,0.01 113 11.1,0.380108948,-0.059673155,0.01 114 11.2,0.380416445,-0.058750664,0.01 115 11.3,0.380919615,-0.057241155,0.01 116 11.4,0.381613429,-0.055159712,0.01 117 11.5,0.382490957,-0.05252713,0.01 118 11.6,0.383543428,-0.049369716,0.01 119 11.7,0.384760328,-0.045719015,0.01 120 11.8,0.386129498,-0.041611505,0.01 121 11.9,0.387637258,-0.037088227,0.01 122 12,0.389268542,-0.032194375,0.01 123 12.1,0.391007051,-0.026978848,0.01 124 12.2,0.392835414,-0.021493757,0.01 125 12.3,0.394735364,-0.015793907,0.01 126 12.4,0.396687916,-0.009936251,0.01 127 12.5,0.398673562,-0.003979314,0.01 128 12.6,0.400672461,0.002017383,0.01 129 12.7,0.402664641,0.007993922,0.01 130 12.8,0.404630197,0.01389059,0.01 131 12.9,0.406549489,0.019648466,0.01 132 13,0.408403341,0.025210022,0.01 133 13.1,0.410173229,0.030519688,0.01 134 13.2,0.41184147,0.035524411,0.01 135 13.3,0.413391395,0.040174186,0.01 136 13.4,0.414807518,0.044422553,0.01 137 13.5,0.416075689,0.048227066,0.01 138 13.6,0.417183236,0.051549709,0.01 139 13.7,0.418119095,0.054357285,0.01 140 13.8,0.418873913,0.05662174,0.01 141 13.9,0.41944015,0.05832045,0.01 142 14,0.419812147,0.059436441,0.01 143 14.1,0.419986188,0.059958563,0.01 144 14.2,0.419960533,0.059881599,0.01 145 14.3,0.419735439,0.059206318,0.01 146 14.4,0.419313156,0.057939467,0.01 147 14.5,0.418697901,0.056093703,0.01 148 14.6,0.417895823,0.05368747,0.01 149 14.7,0.416914937,0.05074481,0.01 150 14.8,0.415765041,0.047295124,0.01 151 14.9,0.414457627,0.043372881,0.01 152 15,0.413005757,0.03901727,0.01 -
anuga_work/development/Hinwood_2008/create_mesh.py
r5076 r5092 5 5 from anuga.pmesh.mesh import * 6 6 from anuga.coordinate_transforms.geo_reference import Geo_reference 7 8 xslope = 5.0 # Distance between the boundary ande the start of the slope 7 9 8 10 … … 15 17 """ 16 18 #Basic geometry 19 global xslope 17 20 18 21 xright = 19.0 19 22 ybottom = 0 20 ytop = 0.4523 ytop = 1.00 21 24 xleft = 0.0 22 xslope = 4.025 ###xslope = slope 23 26 24 27 #Outline -
anuga_work/development/Hinwood_2008/run_dam.py
r5078 r5092 42 42 def elevation_function(x,y): 43 43 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 46 47 z = zeros(size(x), Float) 47 48 for i in range(len(x)): 48 if x[i] < slope:49 if x[i] < xslope: 49 50 z[i] = 0.0 #WARNING: the code in prepare_time_boundary 50 51 # that calc's momentum assumes this is 0.0 51 52 else: 52 z[i] = (x[i]- slope)*0.153 z[i] = (x[i]-xslope)*(1./16.) 53 54 return z 54 55 … … 124 125 #------------------------------------------------------------------------- 125 126 126 domain.set_quantity('stage', 0. 06)127 domain.set_quantity('stage', 0.4) 127 128 domain.set_quantity('friction', friction) 128 129 domain.set_quantity('elevation', elevation_function) … … 151 152 print 'finished' 152 153 153 points = [[2.8,0. 225], #-1.8m from SWL154 [5.1,0. 225], #0.5m from SWL155 [6.6,0. 225], #2m from SWL156 [6.95,0. 255], #2.35m from SWL157 [7.6,0. 255], #3m from SWL158 [8.2,0. 255], #3.5m from SWL159 [9.2,0. 255] #4.5m from SWL154 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 160 161 ] 161 162 … … 177 178 if __name__ == "__main__": 178 179 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.
