Changeset 5413


Ignore:
Timestamp:
Jun 24, 2008, 8:53:59 AM (14 years ago)
Author:
duncan
Message:

Current Hinwood scenario

Location:
anuga_work/development/Hinwood_2008
Files:
1 added
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/Hinwood_2008/plot.py

    r5410 r5413  
    8282    """
    8383    from scenarios import scenarios
    84     outputdir_tag = "_good_tri_area_0.001_A"
     84    outputdir_tag = "_good_tri_area_0.01_A"
    8585    #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!
    8686    # Supported formats: emf, eps, pdf, png, ps, raw, rgba, svg, svgz.
     
    9191        file_exp, file_sim, location_sims, location_exps, outputdir_name, \
    9292                  save_as_list = temp
     93        print "run_data['scenario_id']", run_data['scenario_id']
     94        #location_sims = [location_sims[0]]
     95        #location_exps = [location_exps[0]]
     96        #save_as_list = [save_as_list[0]]
    9397        for loc_sim, loc_exp, save_as, gauge in map(None, location_sims,
    9498                                                    location_exps,
  • anuga_work/development/Hinwood_2008/run_dam.py

    r5410 r5413  
    248248#-------------------------------------------------------------
    249249if __name__ == "__main__":
    250     #slopes = [[-4.5,0.0],[0.0,0.0],[1.285,0.090],[16.1,.960]]
    251     # Note, gauge A has been removed, since it is used as the
    252     # boundary.
    253     run_type = 6
    254 
    255     # T1R5
    256     run_data = {'xleft':[-3.106,0.0],  # Av' of ADV and Gauge A
    257                      'xtoe':[0.0,0.0],
    258                      'xbeach':[1.285,0.090],
    259                      'xright':[16.1,.960],
    260                      'offshore_water_depth':.4,
    261                      'scenario_id':'T1R5',
    262                      'gauge_names':['B','1','2','3','4','5','6','7','8',
    263                                     '9','10','11','12','13','14'],
    264                      'gauge_x':[-0.68, 1.572, 2.572, 3.572, 4.572, 5.572,
    265                                 6.572, 7.572, 8.572, 9.572, 10.572, 11.572,
    266                                 12.572, 13.572, 14.572],
    267                      'gauge_bed_elevation':[-0.400000, -0.293158,
    268                      -0.234473, -0.175788, -0.117104, -0.058419, 0.000266,
    269                      0.058950, 0.117635, 0.176320, 0.235004, 0.293689,
    270                      0.352374, 0.411058, 0.469743]
    271                      }
    272     main( run_data['scenario_id'] + '_boundary.tsm' , run_data,
    273           run_type = run_type,
    274          outputdir_name=run_data['scenario_id'])
    275 
    276     # T1R3
    277     run_data = {'xleft':[-3.106,0.0],  # Av' of ADV and Gauge A
    278                      'xtoe':[0.0,0.0],
    279                      'xbeach':[1.285,0.090],
    280                      'xright':[16.1,.960],
    281                      'offshore_water_depth':.4,
    282                      'scenario_id':'T1R3',
    283                      'gauge_names':['B','1','2','3','4','5','6','7','8',
    284                                     '9','10','11','12','13','14'],
    285                      'gauge_x':[-0.68, 1.572, 2.572, 3.572, 4.572, 5.572,
    286                                 6.572, 7.572, 8.572, 9.572, 10.572, 11.572,
    287                                 12.572, 13.572, 14.572],
    288                      'gauge_bed_elevation':[-0.400000, -0.293158,
    289                      -0.234473, -0.175788, -0.117104, -0.058419, 0.000266,
    290                      0.058950, 0.117635, 0.176320, 0.235004, 0.293689,
    291                      0.352374, 0.411058, 0.469743]
    292                      }
    293     main( run_data['scenario_id'] + '_boundary.tsm'  , run_data,
    294           run_type = run_type,
    295          outputdir_name=run_data['scenario_id'])
    296 
    297     # #T2R7
    298     # xleft is different
    299     run_data = {'xleft':[-4.586,0.0],  # Av' of ADV and Gauge A
    300                      'xtoe':[0.0,0.0],
    301                      'xbeach':[1.285,0.090],
    302                      'xright':[16.1,.960],
    303                      'offshore_water_depth':.4,
    304                      'scenario_id':'T2R7',
    305                      'gauge_names':['B','1','2','3','4','5','6','7','8',
    306                                     '9','10','11','12','13','14'],
    307                      'gauge_x':[-0.68, 1.572, 2.572, 3.572, 4.572, 5.572,
    308                                 6.572, 7.572, 8.572, 9.572, 10.572, 11.572,
    309                                 12.572, 13.572, 14.572],
    310                      'gauge_bed_elevation':[-0.400000, -0.293158,
    311                      -0.234473, -0.175788, -0.117104, -0.058419, 0.000266,
    312                      0.058950, 0.117635, 0.176320, 0.235004, 0.293689,
    313                      0.352374, 0.411058, 0.469743]
    314                      }
    315     main( run_data['scenario_id'] + '_boundary.tsm'  , run_data,
    316           run_type = run_type,
    317           outputdir_name=run_data['scenario_id'])
    318 
    319     # #T2R8
    320     # xleft is different
    321     run_data = {'xleft':[-4.586,0.0],  # Av' of ADV and Gauge A
    322                      'xtoe':[0.0,0.0],
    323                      'xbeach':[1.285,0.090],
    324                      'xright':[16.1,.960],
    325                      'offshore_water_depth':.4,
    326                      'scenario_id':'T2R8',
    327                      'gauge_names':['B','1','2','3','4','5','6','7','8',
    328                                     '9','10','11','12','13','14'],
    329                      'gauge_x':[-0.68, 1.572, 2.572, 3.572, 4.572, 5.572,
    330                                 6.572, 7.572, 8.572, 9.572, 10.572, 11.572,
    331                                 12.572, 13.572, 14.572],
    332                      'gauge_bed_elevation':[-0.400000, -0.293158,
    333                      -0.234473, -0.175788, -0.117104, -0.058419, 0.000266,
    334                      0.058950, 0.117635, 0.176320, 0.235004, 0.293689,
    335                      0.352374, 0.411058, 0.469743]
    336                      }
    337     main( run_data['scenario_id'] + '_boundary.tsm'  , run_data,
    338           run_type = run_type,
    339           outputdir_name=run_data['scenario_id'])
    340 
    341     # #T3R29
    342     # xleft is different
    343     run_data = {'xleft':[-3.875,0.0],  # Av' of ADV and Gauge A
    344                 'xtoe':[0.0,0.0],
    345                 'xbeach':[1.285,0.090],
    346                 'xright':[16.1,.440],
    347                 'offshore_water_depth':.336,
    348                 'scenario_id':'T3R29',
    349                 'gauge_names':['1','2','3','4','5','6','7','8',
    350                                '9','10','11','12','13','14','B'],
    351                 'gauge_x':[1.572, 2.572, 3.572, 4.572, 5.572,
    352                            6.572, 7.572, 8.572, 9.572, 10.572, 11.572,
    353                            12.572, 13.572, 14.572, -0.325],
    354                 'gauge_bed_elevation':[-0.237263, -0.213789, -0.190315,
    355                                        -0.166841, -0.143368, -0.119894,
    356                                        -0.096420, -0.072946, -0.049472,
    357                                        -0.025998, -0.002524, 0.020949,
    358                                        0.044423, 0.067897, -0.336000]
    359                      }
    360     main( run_data['scenario_id'] + '_boundary.tsm'  , run_data,
    361           run_type = run_type,
    362           outputdir_name=run_data['scenario_id'])
    363 
    364     # #T3R28
    365     run_data = {'xleft':[-3.875,0.0],  # Av' of ADV and Gauge A
    366                 'xtoe':[0.0,0.0],
    367                 'xbeach':[1.285,0.090],
    368                 'xright':[16.1,.440],
    369                 'offshore_water_depth':.336,
    370                 'scenario_id':'T3R28',
    371                 'gauge_names':['1','2','3','4','5','6','7','8',
    372                                '9','10','11','12','13','14','B'],
    373                 'gauge_x':[1.572, 2.572, 3.572, 4.572, 5.572,
    374                            6.572, 7.572, 8.572, 9.572, 10.572, 11.572,
    375                            12.572, 13.572, 14.572, -0.325],
    376                 'gauge_bed_elevation':[-0.237263, -0.213789, -0.190315,
    377                                        -0.166841, -0.143368, -0.119894,
    378                                        -0.096420, -0.072946, -0.049472,
    379                                        -0.025998, -0.002524, 0.020949,
    380                                        0.044423, 0.067897, -0.336000]
    381                      }
    382     main( run_data['scenario_id'] + '_boundary.tsm'  , run_data,
    383           run_type = run_type,
    384           outputdir_name=run_data['scenario_id'])
    385 
    386     # #T4R31
    387     # xleft is different
    388     run_data = {'xleft':[-2.43,0.0],  # Av' of ADV and Gauge A
    389                 'xtoe':[0.0,0.0],
    390                 'xbeach':[1.285,0.090],
    391                 'xright':[16.1,.440],
    392                 'offshore_water_depth':.336,
    393                 'scenario_id':'T4R31',
    394                 'gauge_names':['1','2','3','4','5','6','7','8',
    395                                '9','10','11','12','13','14','B'],
    396                 'gauge_x':[1.572, 2.572, 3.572, 4.572, 5.572,
    397                            6.572, 7.572, 8.572, 9.572, 10.572, 11.572,
    398                            12.572, 13.572, 14.572, -0.325],
    399                 'gauge_bed_elevation':[-0.237263, -0.213789, -0.190315,
    400                                        -0.166841, -0.143368, -0.119894,
    401                                        -0.096420, -0.072946, -0.049472,
    402                                        -0.025998, -0.002524, 0.020949,
    403                                        0.044423, 0.067897, -0.336000]
    404                      }
    405     main( run_data['scenario_id'] + '_boundary.tsm'  , run_data,
    406           run_type = run_type,
    407           outputdir_name=run_data['scenario_id'])
    408    
    409     # #T4R32
    410     run_data = {'xleft':[-2.43,0.0],  # Av' of ADV and Gauge A
    411                 'xtoe':[0.0,0.0],
    412                 'xbeach':[1.285,0.090],
    413                 'xright':[16.1,.440],
    414                 'offshore_water_depth':.336,
    415                 'scenario_id':'T4R32',
    416                 'gauge_names':['1','2','3','4','5','6','7','8',
    417                                '9','10','11','12','13','14','B'],
    418                 'gauge_x':[1.572, 2.572, 3.572, 4.572, 5.572,
    419                            6.572, 7.572, 8.572, 9.572, 10.572, 11.572,
    420                            12.572, 13.572, 14.572, -0.325],
    421                 'gauge_bed_elevation':[-0.237263, -0.213789, -0.190315,
    422                                        -0.166841, -0.143368, -0.119894,
    423                                        -0.096420, -0.072946, -0.049472,
    424                                        -0.025998, -0.002524, 0.020949,
    425                                        0.044423, 0.067897, -0.336000]
    426                      }
    427     main( run_data['scenario_id'] + '_boundary.tsm'  , run_data,
    428           run_type = run_type,
    429           outputdir_name=run_data['scenario_id'])
     250   
     251    from scenarios import scenarios
     252   
     253    run_type = 5
     254    for run_data in [scenarios[0]]:
     255        main( run_data['scenario_id'] + '_boundary.tsm'  , run_data,
     256              run_type = run_type,
     257              outputdir_name=run_data['scenario_id'])
     258   
Note: See TracChangeset for help on using the changeset viewer.