Changeset 6672


Ignore:
Timestamp:
Mar 31, 2009, 5:39:21 PM (15 years ago)
Author:
myall
Message:

running and exporting scripts

Location:
anuga_work/production/australia_ph2
Files:
12 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/adelaide/project.py

    r6617 r6672  
    2929tide = 0                # difference between MSL and HAT
    3030event_number = 64322    # the event number
    31 # events: 7867 (altiplano), 27346 (java), 46591 (midAmerica),
    32 # 64322 (sandwich), 65368 (shetland),
    33 # other events: 51466 (newHebrides), 57537 (peru), 58334 (puysegur), 68791 (sumatra)
    34 # 11100 (andaman),
     31# events: 27346 (java 2.1e-6), 58334 (puysegur 8.8e-7), 64322 (sandwich 8.2e-6),
     32# other events: 51466 (newHebrides 7.2e-7), 57537 (peru 4.7e-7), 46591 (midAmerica 1.3e-6),
     33# 68791 (sumatra 7.5e-7), 11100 (andaman 4.9e-7), 65368 (shetland 2.5e-6), 7867 (altiplano 1.0e-6)
    3534alpha = 0.1             # smoothing parameter for mesh
    3635friction=0.01           # manning's friction coefficient
     
    4241wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    4342
    44 setup = 'trial'         # This can be one of three values
     43setup = 'final'         # This can be one of three values
    4544                        #    trial - coarsest mesh, fast
    4645                        #    basic - coarse mesh
  • anuga_work/production/australia_ph2/albany/project.py

    r6617 r6672  
    2727# One or all can be changed each time the run_model script is executed
    2828tide = 0                # difference between MSL and HAT
     29central_meridian = 117.0 # central meridian for projection (optional)
    2930event_number = 64344    # the event number or the mux file name
    3031# selected events: 27319 (java), 64344 (sandwich), 46697 (midamerica),
     
    4748                        #    final - fine mesh, slowest
    4849
    49 internal_polygon = True
     50internal_polygon = False
    5051
    5152#-------------------------------------------------------------------------------
  • anuga_work/production/australia_ph2/broome/project.py

    r6617 r6672  
    3636# index is only used when wave = Tb
    3737index = 1693            # index from the PTHA - Y2000 5.4259m
    38 wave = 'Bf'             # Bf (sts wave) Tb (index wave)
     38wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    3939
    4040
  • anuga_work/production/australia_ph2/carnarvon/export_results_max.py

    r6412 r6672  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090224_102606_run_trial_0_58362_kvanputt'
    25 ##time_dir2 = ''
     24time_dir1 = '20090324_150124_run_final_0_27285_109_Bf_mhingee'
     25time_dir2 = '20090324_161452_run_final_0_27285_109_Tb_mhingee'
    2626##time_dir3 = ''
    2727##time_dir4 = ''
     
    2929##time_dir6 = ''
    3030##
    31 time_dirs = [time_dir1] #time_dir2]# , time_dir4, time_dir5, time_dir6]
     31time_dirs = [time_dir1, time_dir2] # , time_dir4, time_dir5, time_dir6]
    3232 
    3333
     
    3535##cellsize = 5
    3636
    37 #timestep = None    # None means no timestep!
    38 timestep = 0
     37timestep = None    # None means no timestep!
     38##timestep = 0
    3939
    4040######
     
    5656
    5757# one or more key strings from var_equations above
    58 var = ['elevation']
     58var = ['stage','speed']
    5959
    6060######
     
    8686
    8787            name1 = join(directory, time_dir, project.scenario_name)
    88             #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0'
    89             #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0'
    90 
    91             names = [name1] #, name2]#, name3]
     88            name2 = join(directory, time_dir, project.scenario_name)+'_time_8760_0'
     89            name3 = join(directory, time_dir, project.scenario_name)+'_time_17520_0'
     90            name4 = join(directory, time_dir, project.scenario_name)+'_time_26280_0'
     91            name5 = join(directory, time_dir, project.scenario_name)+'_time_35040_0'
     92            name6 = join(directory, time_dir, project.scenario_name)+'_time_43800_0'
     93            name7 = join(directory, time_dir, project.scenario_name)+'_time_52560_0'
     94            name8 = join(directory, time_dir, project.scenario_name)+'_time_61320_0'
     95            name9 = join(directory, time_dir, project.scenario_name)+'_time_70080_0'
     96            name10 = join(directory, time_dir, project.scenario_name)+'_time_78840_0'
     97           
     98            names = [name1, name2, name3, name4, name5, name6, name7, name8, name9, name10]
    9299     
    93            # asc_name = []   
     100            asc_name = []   
    94101               
    95102            for name in names:
     
    112119                            format = 'asc')
    113120
    114 ##                asc_name.append(outname + '.asc')
    115 ##   
    116 ##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
    117 ##
    118 ##            print 'max asc outname ', maxasc_outname
    119 ##            print 'asc_name ', str(asc_name)
    120 ##           
    121 ##            MaxAsc(maxasc_outname, asc_name)
     121                asc_name.append(outname + '.asc')
     122   
     123            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     124
     125            print 'max asc outname ', maxasc_outname
     126            print 'asc_name ', str(asc_name)
     127           
     128            MaxAsc(maxasc_outname, asc_name)
  • anuga_work/production/australia_ph2/carnarvon/project.py

    r6617 r6672  
    4444                        #    final - fine mesh, slowest
    4545
    46 internal_polygon = False
     46internal_polygon = True
    4747
    4848#-------------------------------------------------------------------------------
  • anuga_work/production/australia_ph2/ceduna/project.py

    r6617 r6672  
    3838friction=0.01           # manning's friction coefficient
    3939starttime=0             # start time for simulation
    40 finaltime=1000          # final time for simulation
     40finaltime=80000          # final time for simulation
    4141
    4242# index is only used when wave = Tb
  • anuga_work/production/australia_ph2/darwin/export_results_max.py

    r6532 r6672  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090316_120226_run_final_0_17898_1833_Bf_mhingee'
    25 ##time_dir2 = '20090316_120344_run_final_0_64977_1833_Bf_mhingee'
    26 ##time_dir3 = '20090316_120457_run_final_0_70906_1833_Bf_mhingee'
    27 ##time_dir4 = '20090316_120558_run_final_0_71261_1833_Bf_mhingee'
     24time_dir1 = '20090327_103428_run_final_0_70906_1833_Bf_None_mhingee'
     25##time_dir2 = ''
     26##time_dir3 = ''
     27##time_dir4 = ''
    2828##time_dir5 = ''
    2929##time_dir6 = ''
  • anuga_work/production/australia_ph2/darwin/project.py

    r6617 r6672  
    2828tide = 0                # difference between MSL and HAT
    2929central_meridian = 129.0 # for projection (necessary for nonstandard)
    30 event_number = 70906    # the event number or the mux file name
     30event_number = 17898    # the event number or the mux file name
    3131# contributing events: 17898 (Flores), 64977 (Seram), 70906 (Sumba_normal),
    3232# 71261 (Timor_OTF)
  • anuga_work/production/australia_ph2/hobart/project.py

    r6476 r6672  
    2727# One or all can be changed each time the run_model script is executed
    2828tide = 0                # difference between MSL and HAT
    29 zone = 55               # specify zone of model
     29central_meridian = 147.0 # Central meridian for projection (optional)
    3030
    3131alpha = 0.1             # smoothing parameter for mesh
     
    3535
    3636# index is only used when wave = Tb
    37 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     37wave = 'Bf'             # Bf (sts wave) Tb (index wave)
     38
    3839# for melbourne (west) side:
    39 event_number = 58349    # the event number or the mux file name
     40##event_number = 58349    # the event number or the mux file name
    4041# contributing sources:
    4142# 51328 (newhebrides), 58349 (puysegur), 64535 (sandwich)
    4243# other possible: 31966 (kermandec), 65350 (shetland)
    43 index = 1983            # index from the PTHA - west side, close to melbourne; 1983 has 1 in 2000 wave height 0.534m
     44##index = 1983            # index from the PTHA - west side, close to melbourne; 1983 has 1 in 2000 wave height 0.534m
    4445
    4546# for east, Hobart, side - point opposite hobart
    46 ##event_number = 51202    # the event number or the mux file name
    47 # contributing events: 64434 (sandwich), 58117 (puysegur), 51202 (new hebrides)
    48 # Others: 31899 (kermandec), 46657 (midamerica), 57586 (peru),
    49 ##index = 2473            # index from PTHA - east, opposite Hobart; 1.112m for 1 in 2000
     47event_number = 64461    # the event number or the mux file name
     48# contributing events:  51450 (newHebrides 7.3e-6), 58078 (puysegur 5.2e-6), 64461 (sandwich 8.2e-6)
     49# Others: 31816 (kermadec 4.7e-7), 46651 (midamerica 3.5e-7),
     50#
     51index = 2479            # index from PTHA - Y2000 1.0977m
     52
     53# index = 2460   # for Bichen
    5054
    5155setup = 'final'         # This can be one of three values
     
    5357                        #    basic - coarse mesh
    5458                        #    final - fine mesh, slowest
     59
     60internal_polygon = False
    5561
    5662#-------------------------------------------------------------------------------
     
    6571#-------------------------------------------------------------------------------
    6672
    67 output_comment = [setup, tide, event_number, index, wave]
     73if internal_polygon:
     74    internal_poly_comment = 'internal'
     75else:
     76    internal_poly_comment = None
     77
     78output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
    6879
    6980#-------------------------------------------------------------------------------
     
    7788
    7889# Format for point is x,y,elevation (with header)
    79 point_filenames = ['mel_grid250m.txt',
    80                    'bight_250m.txt'] # 250m grid 2005
     90point_filenames = ['grid_250m.txt'] # 250m grid 2005
    8191
    8292### Add csv header list to all files in point_filenames
     
    94104# Used in run_model.py
    95105# Format for points easting,northing (no header)                   
    96 interior_regions_data = [['tasmania.csv', 1000000]]
     106interior_regions_data = []
     107
     108# add an internal polygon to force different mesh generation
     109# used to test for discretisation error when building elevation
     110# make sure file is in same folder as interior regions and bouding polygon;
     111# format is same (2 column .csv; easting, northing; no header)
     112if internal_polygon:
     113    interior_regions_data.append(['internal_polygon.csv',
     114                                  bounding_polygon_maxarea])
    97115
    98116# LAND - used to set the initial stage/water to be offcoast only
  • anuga_work/production/australia_ph2/hobart/run_model.py

    r6511 r6672  
    156156                        verbose=True)
    157157    domain.set_boundary({'back': Br,
    158                      'side': Bd,
     158                     'side': Bt,
    159159                     'ocean': Bf})
    160160   
     
    163163
    164164    domain.set_boundary({'back': Br,
    165                          'side': Bd,
     165                         'side': Bt,
    166166                         'ocean': Tb})
    167167else:
  • anuga_work/production/australia_ph2/strahan/export_results_max.py

    r6619 r6672  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090302_104336_run_final_0_27255_mhingee'
    25 ##time_dir2 = ''
    26 ##time_dir3 = ''
    27 ##time_dir4 = ''
    28 ##time_dir5 = ''
    29 ##time_dir6 = ''
    30 ##
    31 time_dirs = [time_dir1] #time_dir2]# , time_dir4, time_dir5, time_dir6]
     24time_dir1 = '20090327_143655_run_final_0_64214_2044_Tb_None_mhingee'
     25time_dir2 = '20090327_143836_run_final_0_58337_2044_Tb_None_mhingee'
     26time_dir3 = '20090330_095007_run_final_0_68779_2044_Tb_None_mhingee'
     27time_dir4 = '20090330_095719_run_final_0_64214_2044_Bf_None_mhingee'
     28time_dir5 = '20090330_095845_run_final_0_58337_2044_Bf_None_mhingee'
     29time_dir6 = '20090330_143221_run_final_0_68779_2044_Bf_None_mhingee'
     30
     31time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    3232 
    3333
     
    3636
    3737#timestep = None    # None means no timestep!
    38 timestep = 0
     38timestep = None
    3939
    4040######
     
    5656
    5757# one or more key strings from var_equations above
    58 var = ['elevation','stage']
     58var = ['stage','speed']
    5959
    6060######
     
    9191            names = [name1] #, name2]#, name3]
    9292     
    93            # asc_name = []   
     93            asc_name = []   
    9494               
    9595            for name in names:
     
    112112                            format = 'asc')
    113113
    114 ##                asc_name.append(outname + '.asc')
    115 ##   
    116 ##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
    117 ##
    118 ##            print 'max asc outname ', maxasc_outname
    119 ##            print 'asc_name ', str(asc_name)
    120 ##           
    121 ##            MaxAsc(maxasc_outname, asc_name)
     114                asc_name.append(outname + '.asc')
     115   
     116            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     117
     118            print 'max asc outname ', maxasc_outname
     119            print 'asc_name ', str(asc_name)
     120           
     121            MaxAsc(maxasc_outname, asc_name)
  • anuga_work/production/australia_ph2/strahan/project.py

    r6619 r6672  
    2828tide = 0                # difference between MSL and HAT
    2929central_meridian = 147.0 # Central meridian for projection (optional)
    30 event_number =     # the event number or the mux file name
    31 # main events: 64214 (sandwich 1.4e-5), 58337 (puysegur 3e-6),
     30event_number = 68779    # the event number or the mux file name
     31# main events: 64214 (sandwich 1.4e-5), 58337 (puysegur 3e-6), 68779 (sumatra 7.53-7)
    3232# others: 27295 (java 3e-7), 31817 (kermadec 5e-7), 46559 (midAmerica 4.2e-7)
    3333# 51460 (newHebrides 7e-7), 57477 (peru 5.4e-7), 65322 (shetland 3.6e-7),
    34 # 68779 (sumatra 7.53-7)
     34
    3535alpha = 0.1             # smoothing parameter for mesh
    3636friction=0.01           # manning's friction coefficient
     
    4040# index is only used when wave = Tb
    4141index = 2044            # index from the PTHA - Y2000 is 0.60374m
    42 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     42wave = 'Bf'             # Bf (sts wave) Tb (index wave)
    4343
    4444setup = 'final'         # This can be one of three values
     
    7777
    7878# Format for point is x,y,elevation (with header)
    79 point_filenames = ['grid250m_pts.txt',
    80                    'coast_250.txt'] # 250m grid 2005
     79point_filenames = ['grid_250m.txt'] # 250m grid 2005
    8180
    8281### Add csv header list to all files in point_filenames
     
    125124# Thinned ordering file from Hazard Map (geographic)
    126125# Format is index,latitude,longitude (with header)
    127 urs_order_filename = 'urs_order_austtg.csv'
     126urs_order_filename = 'urs_order.csv'
    128127
    129128# Landward bounding points
Note: See TracChangeset for help on using the changeset viewer.