Changeset 6596


Ignore:
Timestamp:
Mar 24, 2009, 2:14:11 PM (15 years ago)
Author:
jgriffin
Message:

changed event number

Location:
anuga_work/production/australia_ph2/sydney
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/sydney/Arc_asc2raster_GDA94z56.py

    r6529 r6596  
    3030
    3131time_dir1 = '20090311_145740_run_final_0_58152_2956_Tb_jgriffin'
    32 time_dir2 = '20090311_145828_run_final_0_58152_2956_Bf_jgriffin'
    33 time_dir3 = '20090313_154712_run_final_0_58222_2956_Bf_jgriffin'
    34 time_dir4 = '20090313_154810_run_final_0_58326_2956_Bf_jgriffin'
    35 time_dir5 = '20090316_142137_run_final_0_58326_2956_Bf_jgriffin'
     32##time_dir2 = '20090311_145828_run_final_0_58152_2956_Bf_jgriffin'
     33##time_dir3 = '20090313_154712_run_final_0_58222_2956_Bf_jgriffin'
     34##time_dir4 = '20090313_154810_run_final_0_58326_2956_Bf_jgriffin'
     35##time_dir5 = '20090316_142137_run_final_0_58326_2956_Bf_jgriffin'
    3636##time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'
    3737
    3838
    39 time_dirs = [time_dir3, time_dir4 , time_dir3, time_dir4, time_dir5]#, time_dir6]
     39time_dirs = [time_dir1]#, time_dir4 , time_dir3, time_dir4, time_dir5]#, time_dir6]
    4040
    4141for time_dir in time_dirs:
     
    4848    #ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_conditions_ocean1"
    4949   
    50     print 'Process: Create File GDB'
    51     gp.CreateFileGDB_management(folder, "raster")
     50##    print 'Process: Create File GDB'
     51##    gp.CreateFileGDB_management(folder, "raster")
    5252
    5353    gp.Workspace = raster_gbd
  • anuga_work/production/australia_ph2/sydney/export_results_max.py

    r6529 r6596  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090316_142137_run_final_0_58326_2956_Bf_jgriffin'
     24time_dir1 = '20090311_145740_run_final_0_58152_2956_Tb_jgriffin'
    2525##time_dir2 = '20090311_152929_run_final_0_70385_1203_Tb_jgriffin'
    2626##time_dir3 = '20090311_153009_run_final_0_70454_1203_Bf_jgriffin'
     
    2828time_dirs = [time_dir1]#,time_dir2,time_dir3,time_dir4] 
    2929
    30 cellsize = 250
     30cellsize = 20
    3131##cellsize = 5
    3232
    33 #timestep = None    # None means no timestep!
    34 timestep = 0
     33timestep = None    # None means no timestep!
     34#timestep = 0
    3535
    3636######
     
    5252
    5353# one or more key strings from var_equations above
    54 var = ['elevation','stage']
     54var = ['stage']
    5555
    5656######
     
    6969            northing_min = None
    7070            northing_max = None
     71##            easting_min = 290000
     72##            easting_max = 316000
     73##            northing_min = 6117000
     74##            northing_max = 6170000
    7175        else:
    7276            try:
     
    8286
    8387            name1 = join(directory, time_dir, project.scenario_name)
    84 ##            name2 = join(directory, time_dir, project.scenario_name)+'_time_66360_0'
     88            name2 = join(directory, time_dir, project.scenario_name)+'_time_66360_0'
    8589##            name3 = join(directory, time_dir, project.scenario_name)+'_time_26640_0'
    8690##            name4 = join(directory, time_dir, project.scenario_name)+'_time_39960_0'
  • anuga_work/production/australia_ph2/sydney/get_timeseries.py

    r6176 r6596  
    1212"""
    1313
    14 from os import sep
     14from os import sep, rename, listdir, system
    1515from anuga.abstract_2d_finite_volumes.util import sww2csv_gauges,csv2timeseries_graphs
    1616import project
    1717
    18 timestamp1='20080815_103442_run_final_0.0_polyline_alpha0.1_kvanputt'
    19 timestamp2='20080815_103336_run_final_0.6_polyline_alpha0.1_kvanputt'
     18directory = project.output_folder
    2019
    21 timestamps = [timestamp1, timestamp2]
    22 for timestamp in timestamps:
     20##time_dir1 = '20081009_184721_run_final_0.6_68693_alpha0.1_kvanputt'+sep+'gauge'
     21##time_dir2 = '20081031_120316_run_final_0_68693_alpha0.1_kvanputt'+sep+'gauge'
     22##time_dir3 = '20081031_120515_run_final_0_27255_alpha0.1_kvanputt'+sep+'gauge'
     23##time_dir4 = '20081031_120642_run_final_0_27283_alpha0.1_kvanputt'+sep+'gauge'
     24##time_dir5 = '20081031_120753_run_final_0.6_27283_alpha0.1_kvanputt'+sep+'gauge'
     25##time_dir6 = '20081111_113832_run_final_0.6_27255_alpha0.1_kvanputt'+sep+'gauge'
     26##
     27##time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    2328
    24     filename1=project.output_dir+timestamp+sep+project.scenario_name+'.sww'
    25     filename2=project.output_dir+timestamp+sep+'sww2'+sep+project.scenario_name+'_time_39900_0.sww'
     29##time_dir1 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt'
     30time_dir6 = '20090311_145740_run_final_0_58152_2956_Tb_jgriffin'
    2631
    27     filenames = [filename1, filename2]
    28     for filename in filenames:
     32time_dirs = [time_dir6]
    2933
    30         gauges = [project.gauges_dir_name, project.gauges_dir_name2]
    31         for gauge in gauges:
    32            
    33             sww2csv_gauges(filename,gauge,
    34                             quantities = ['stage','speed','depth','elevation'],
    35                             verbose=True)
    36                        
     34for time_dir in time_dirs:
     35    name = directory+sep+time_dir+sep+project.scenario_name
     36    gauge = project.gauges
     37    sww2csv_gauges(name+'.sww',gauge,
     38                   quantities = ['stage','speed','depth','elevation'],
     39                   verbose=True)
     40 
    3741
    38    
    39 
    40 
  • anuga_work/production/australia_ph2/sydney/project.py

    r6529 r6596  
    2929tide = 0                # difference between MSL and HAT
    3030zone = 56              # specify zone of model
    31 event_number = 58326    # the event number or the mux file name
     31event_number = 58152    # the event number or the mux file name
    3232alpha = 0.1             # smoothing parameter for mesh
    3333friction=0.01           # manning's friction coefficient
     
    3737# index is only used when wave = Tb
    3838index = 2956            # index from the PTHA
    39 wave = 'Bf'             # Bf (sts wave) Tb (index wave)
     39wave = 'Tb'             # Bf (sts wave) Tb (index wave)
    4040
    4141setup = 'final'         # This can be one of three values
Note: See TracChangeset for help on using the changeset viewer.