Changeset 4055


Ignore:
Timestamp:
Dec 5, 2006, 2:41:16 PM (17 years ago)
Author:
duncan
Message:

friction scripts; update

Location:
anuga_work/development
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/friction_UA_flume_2006/elevation.xya

    r4047 r4055  
    1 elevation                               
    2 0       ,       0       ,       0
    3 0.1     ,       0       ,       0
    4 0.2     ,       0       ,       0
    5 0.3     ,       0       ,       0
    6 0.4     ,       0       ,       0
    7 0.5     ,       0       ,       0
    8 0.6     ,       0       ,       0
    9 0.7     ,       0       ,       0
    10 0.8     ,       0       ,       0
    11 0.9     ,       0       ,       0
    12 1       ,       0       ,       0
    13 1.1     ,       0       ,       0
    14 1.2     ,       0       ,       0
    15 1.3     ,       0       ,       0
    16 1.4     ,       0       ,       0
    17 1.5     ,       0       ,       0
    18 1.6     ,       0       ,       0
    19 1.7     ,       0       ,       0
    20 1.8     ,       0       ,       0
    21 1.9     ,       0       ,       0
    22 2       ,       0       ,       0
    23 2.1     ,       0       ,       0
    24 2.2     ,       0       ,       0
    25 2.3     ,       0       ,       0
    26 2.4     ,       0       ,       0
    27 2.5     ,       0       ,       0
    28 2.6     ,       0       ,       0
    29 2.7     ,       0       ,       0
    30 2.8     ,       0       ,       0
    31 2.9     ,       0       ,       0
    32 3       ,       0       ,       0
    33 3.1     ,       0       ,       0
    34 3.2     ,       0       ,       0
    35 3.3     ,       0       ,       0
    36 3.4     ,       0       ,       0
    37 3.5     ,       0       ,       0
    38 3.6     ,       0       ,       0
    39 3.7     ,       0       ,       0
    40 3.8     ,       0       ,       0
    41 3.9     ,       0       ,       0
    42 4       ,       0       ,       0
    43 4.1     ,       0       ,       0
    44 4.2     ,       0       ,       0
    45 4.3     ,       0       ,       0
    46 4.4     ,       0       ,       0
    47 4.5     ,       0       ,       0
    48 4.6     ,       0       ,       0
    49 4.7     ,       0       ,       0
    50 4.8     ,       0       ,       0
    51 4.9     ,       0       ,       0
    52 5       ,       0       ,       0
    53 5.1     ,       0       ,       0.01
    54 5.2     ,       0       ,       0.02
    55 5.3     ,       0       ,       0.03
    56 5.4     ,       0       ,       0.04
    57 5.5     ,       0       ,       0.05
    58 5.6     ,       0       ,       0.06
    59 5.7     ,       0       ,       0.07
    60 5.8     ,       0       ,       0.08
    61 5.9     ,       0       ,       0.09
    62 6       ,       0       ,       0.1
    63 6.1     ,       0       ,       0.11
    64 6.2     ,       0       ,       0.12
    65 6.3     ,       0       ,       0.13
    66 6.4     ,       0       ,       0.14
    67 6.5     ,       0       ,       0.15
    68 6.6     ,       0       ,       0.16
    69 6.7     ,       0       ,       0.17
    70 6.8     ,       0       ,       0.18
    71 6.9     ,       0       ,       0.19
    72 7       ,       0       ,       0.2
    73 7.1     ,       0       ,       0.21
    74 7.2     ,       0       ,       0.22
    75 7.3     ,       0       ,       0.23
    76 7.4     ,       0       ,       0.24
    77 7.5     ,       0       ,       0.25
    78 7.6     ,       0       ,       0.26
    79 7.7     ,       0       ,       0.27
    80 7.8     ,       0       ,       0.28
    81 7.9     ,       0       ,       0.29
    82 8       ,       0       ,       0.3
    83 8.1     ,       0       ,       0.31
    84 8.2     ,       0       ,       0.32
    85 8.3     ,       0       ,       0.33
    86 8.4     ,       0       ,       0.34
    87 8.5     ,       0       ,       0.35
    88 8.6     ,       0       ,       0.36
    89 8.7     ,       0       ,       0.37
    90 8.8     ,       0       ,       0.38
    91 8.9     ,       0       ,       0.39
    92 9       ,       0       ,       0.4
    93 9.1     ,       0       ,       0.41
    94 9.2     ,       0       ,       0.42
    95 9.3     ,       0       ,       0.43
    96 9.4     ,       0       ,       0.44
    97 9.5     ,       0       ,       0.45
    98 9.6     ,       0       ,       0.46
    99 9.7     ,       0       ,       0.47
    100 9.8     ,       0       ,       0.48
    101 9.9     ,       0       ,       0.49
    102 10      ,       0       ,       0.5
    103 10.1    ,       0       ,       0.51
    104 10.2    ,       0       ,       0.52
    105 10.3    ,       0       ,       0.53
    106 10.4    ,       0       ,       0.54
    107 10.5    ,       0       ,       0.55
    108 10.6    ,       0       ,       0.56
    109 10.7    ,       0       ,       0.57
    110 10.8    ,       0       ,       0.58
    111 10.9    ,       0       ,       0.59
    112 11      ,       0       ,       0.6
    113 11.1    ,       0       ,       0.61
    114 11.2    ,       0       ,       0.62
    115 11.3    ,       0       ,       0.63
    116 11.4    ,       0       ,       0.64
    117 11.5    ,       0       ,       0.65
    118 11.6    ,       0       ,       0.66
    119 11.7    ,       0       ,       0.67
    120 11.8    ,       0       ,       0.68
    121 11.9    ,       0       ,       0.69
    122 12      ,       0       ,       0.7
    123 12.1    ,       0       ,       0.71
    124 12.2    ,       0       ,       0.72
    125 12.3    ,       0       ,       0.73
    126 12.4    ,       0       ,       0.74
    127 12.5    ,       0       ,       0.75
    128 12.6    ,       0       ,       0.76
    129 12.7    ,       0       ,       0.77
    130 12.8    ,       0       ,       0.78
    131 12.9    ,       0       ,       0.79
    132 13      ,       0       ,       0.8
     1elevation,,
     20,0,0
     30.1,0,0
     40.2,0,0
     50.3,0,0
     60.4,0,0
     70.5,0,0
     80.6,0,0
     90.7,0,0
     100.8,0,0
     110.9,0,0
     121,0,0
     131.1,0,0
     141.2,0,0
     151.3,0,0
     161.4,0,0
     171.5,0,0
     181.6,0,0
     191.7,0,0
     201.8,0,0
     211.9,0,0
     222,0,0
     232.1,0,0
     242.2,0,0
     252.3,0,0
     262.4,0,0
     272.5,0,0
     282.6,0,0
     292.7,0,0
     302.8,0,0
     312.9,0,0
     323,0,0
     333.1,0,0
     343.2,0,0
     353.3,0,0
     363.4,0,0
     373.5,0,0
     383.6,0,0
     393.7,0,0
     403.8,0,0
     413.9,0,0
     424,0,0
     434.1,0,0.01
     444.2,0,0.02
     454.3,0,0.03
     464.4,0,0.04
     474.5,0,0.05
     484.6,0,0.06
     494.7,0,0.07
     504.8,0,0.08
     514.9,0,0.09
     525,0,0.1
     535.1,0,0.11
     545.2,0,0.12
     555.3,0,0.13
     565.4,0,0.14
     575.5,0,0.15
     585.6,0,0.16
     595.7,0,0.17
     605.8,0,0.18
     615.9,0,0.19
     626,0,0.2
     636.1,0,0.21
     646.2,0,0.22
     656.3,0,0.23
     666.4,0,0.24
     676.5,0,0.25
     686.6,0,0.26
     696.7,0,0.27
     706.8,0,0.28
     716.9,0,0.29
     727,0,0.3
     737.1,0,0.31
     747.2,0,0.32
     757.3,0,0.33
     767.4,0,0.34
     777.5,0,0.35
     787.6,0,0.36
     797.7,0,0.37
     807.8,0,0.38
     817.9,0,0.39
     828,0,0.4
     838.1,0,0.41
     848.2,0,0.42
     858.3,0,0.43
     868.4,0,0.44
     878.5,0,0.45
     888.6,0,0.46
     898.7,0,0.47
     908.8,0,0.48
     918.9,0,0.49
     929,0,0.5
     939.1,0,0.51
     949.2,0,0.52
     959.3,0,0.53
     969.4,0,0.54
     979.5,0,0.55
     989.6,0,0.56
     999.7,0,0.57
     1009.8,0,0.58
     1019.9,0,0.59
     10210,0,0.6
     10310.1,0,0.61
     10410.2,0,0.62
     10510.3,0,0.63
     10610.4,0,0.64
     10710.5,0,0.65
     10810.6,0,0.66
     10910.7,0,0.67
     11010.8,0,0.68
     11110.9,0,0.69
     11211,0,0.7
     11311.1,0,0.71
     11411.2,0,0.72
     11511.3,0,0.73
     11611.4,0,0.74
     11711.5,0,0.75
     11811.6,0,0.76
     11911.7,0,0.77
     12011.8,0,0.78
     12111.9,0,0.79
     12212,0,0.8
     12312.1,0,0.81
  • anuga_work/development/friction_UA_flume_2006/parallel_run_dam.py

    r4051 r4055  
    6363        pro_instance = project.Project(['data','flumes','dam_2006'],
    6464                                        outputdir_name=outputdir_name)
     65        print "The output dir is", pro_instance.outputdir
     66        copy_code_files(pro_instance.outputdir,__file__,
     67                        dirname(project.__file__) \
     68                        + sep + project.__name__+'.py')
     69        copy (pro_instance.codedir + 'run_dam.py',
     70              pro_instance.outputdir + 'run_dam.py')
     71        copy (pro_instance.codedir + 'create_mesh.py',
     72              pro_instance.outputdir + 'create_mesh.py')
    6573    else:
    6674        # wait for a bit..I don't know if I need this
     
    8189
    8290    # creates copy of code in output dir
    83     print "The output dir is", pro_instance.outputdir
    84     copy_code_files(pro_instance.outputdir,__file__,
    85                     dirname(project.__file__) \
    86                     + sep + project.__name__+'.py')
    87     copy (pro_instance.codedir + 'run_dam.py',
    88           pro_instance.outputdir + 'run_dam' + str(friction)+ '.py')
    89     copy (pro_instance.codedir + 'create_mesh.py',
    90           pro_instance.outputdir + 'create_mesh.py')
    9191    if is_trial_run is False:
    92         start_screen_catcher(pro_instance.outputdir, int(friction*100))
     92        start_screen_catcher(pro_instance.outputdir, rank, pypar.size())
    9393
    9494    print 'USER:    ', pro_instance.user
     
    119119    domain.set_datadir(pro_instance.outputdir)
    120120    domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
    121     domain.set_minimum_storable_height(0.01)
     121    domain.set_minimum_storable_height(0.001)
    122122    #domain.set_store_vertices_uniquely(True)  # for writting to sww
    123123
     
    140140
    141141    Br = Reflective_boundary(domain)
    142     Bd = Dirichlet_boundary([0,0,0])  # to drain the water out.
    143     domain.set_boundary( {'wall': Br, 'edge': Bd} )
     142    #Bd = Dirichlet_boundary([0,0,0])  # to drain the water out.
     143    domain.set_boundary( {'wall': Br, 'edge': Br} )
    144144
    145145    #-------------------------------------------------------------------------
     
    187187    size = pypar.size() #Number of MPI processes
    188188    node = pypar.get_processor_name()
    189     frictions = [0.0, 0.01, 0.0125]
     189    frictions = [0.0, 0.005, 0.0075, 0.01, 0.0125, 0.014, 0.015,
     190                 0.016, 0.017, 0.018,
     191                 0.019, 0.02,
     192                 0.0225, 0.025, 0.03, 0.035, 0.04]
    190193    for i, friction in enumerate(frictions):
    191194        if i%size == rank:
    192195            print "I am processor %d of %d on node %s" %(rank, size, node)
    193             #main(friction, is_trial_run = True, outputdir_name='friction_UA')
     196            main(friction, is_trial_run = False, outputdir_name='friction_UA_min_sww_fix')
  • anuga_work/development/friction_UA_flume_2006/run_dam.py

    r4051 r4055  
    4040        outputdir_name += '_test'
    4141        yieldstep = 1
    42         finaltime = 31
     42        finaltime = 15
    4343    else:
    4444        yieldstep = 0.01
    45         finaltime = 31
     45        finaltime = 15
    4646           
    4747    pro_instance = project.Project(['data','flumes','dam_2006'],
     
    152152#-------------------------------------------------------------
    153153if __name__ == "__main__":
    154      frictions = [0.04]
     154     frictions = [0.00,0.02]
    155155     for i, friction in enumerate(frictions):
    156          main(friction, is_trial_run = True, outputdir_name='friction_UA_2')
     156         main(friction, is_trial_run = True, outputdir_name='friction_UA')
    157157
  • anuga_work/development/friction_dam_2006/calc_norm.py

    r4047 r4055  
    130130#-------------------------------------------------------------
    131131if __name__ == "__main__":
    132     calc_norm('fromD_0.2_Slope0.05_x0.2-0.6.csv',
     132    calc_norm('fromD_0.2_Slope0.05_x0.2-0.6_clean.csv',
    133133              is_trial_run = False,
    134134              outputdir_name='friction_set_A')
  • anuga_work/development/friction_dam_2006/run_dam.py

    r4047 r4055  
    152152#-------------------------------------------------------------
    153153if __name__ == "__main__":
    154     for friction in [ 0.0275]:
     154    for friction in [ 0.05]:
    155155        main(friction, is_trial_run = False, outputdir_name='friction_set_A')
    156156   
Note: See TracChangeset for help on using the changeset viewer.