Changeset 7045


Ignore:
Timestamp:
May 15, 2009, 4:17:16 PM (11 years ago)
Author:
jgriffin
Message:
 
Location:
anuga_work/production/new_south_wales/batemans_bay
Files:
1 added
4 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/new_south_wales/batemans_bay/Arc_asc2raster_GDA94z56.py

    r6988 r7045  
    2222gp.overwriteoutput = 1
    2323
    24 scenario_dir="\\\\nas2\\gemd\\georisk_models\\inundation\\data\\new_south_wales\\batemans_bay_tsunami_scenario_2009\\"
     24scenario_dir="\\\\nas2\\gemd\\georisk_models\\inundation\\data\\new_south_wales\\gosford_tsunami_scenario_2009\\"
    2525output_dir = "anuga\\outputs\\"
    2626
     
    2929##time_dirs = [time_dir1, time_dir2]
    3030
    31 time_dir1 = '20090323_141118_run_final_0_58129_jgriffin'
    32 time_dir2 = '20090323_161308_run_final_0_58115_jgriffin'
    33 time_dir3 = '20090323_161328_run_final_0_58226_jgriffin'
    34 time_dir4 = '20090323_164424_run_final_0_58284_jgriffin'
    35 time_dir5 = '20090323_164755_run_final_0_58286_jgriffin'
     31time_dir1 = '20090514_094907_run_final_0_58242_jgriffin'
     32##time_dir2 = '20090323_161308_run_final_0_58115_jgriffin'
     33##time_dir3 = '20090323_161328_run_final_0_58226_jgriffin'
     34##time_dir4 = '20090323_164424_run_final_0_58284_jgriffin'
     35##time_dir5 = '20090323_164755_run_final_0_58286_jgriffin'
    3636##time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'
    3737
    3838
    39 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5]#, time_dir6]
     39time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5]#, time_dir6]
    4040
    4141for time_dir in time_dirs:
     
    5656   
    5757    #replication dictionary
    58     replicate = (('batemans_bay', ''),
     58    replicate = (('gosford', ''),
    5959                 ('_', ''),('Geordie', 'Geo'),('Sorrento', 'Sor'), ('max','M_'),
    6060                 ('Fremantle', 'Fre'),('Rockingham', 'Roc'),('depth','_dep_'),
     
    6262
    6363    generate_filename = []
    64     input_ascii = glob.glob(folder + '*depth_max.asc')
     64    input_ascii = glob.glob(folder + '*max.asc')
    6565    print time_dir
    6666
  • anuga_work/production/new_south_wales/batemans_bay/export_results_max.py

    r6988 r7045  
    2929##time_dirs = [time_dir1, time_dir2]#, time_dir3]
    3030
    31 time_dir1 = '20090422_162114_run_final_0.0_51436_jgriffin'
     31time_dir1 = '20090514_094907_run_final_0_58242_jgriffin'
    3232##time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt'
    3333##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt'
     
    6363
    6464# one or more key strings from var_equations above
    65 var = ['depth', 'stage', 'speed']#, 'stage', 'speed','depth']
     65var = ['elevation']#, 'stage', 'speed','depth']
    6666
    6767######
     
    9393
    9494            name1 = join(directory, time_dir, project.scenario_name)
    95             name2 = join(directory, time_dir, project.scenario_name)+'_time_6310_0'
    96             name3 = join(directory, time_dir, project.scenario_name)+'_time_12620_0'
    97             name4 = join(directory, time_dir, project.scenario_name)+'_time_18930_0'
    98             name5 = join(directory, time_dir, project.scenario_name)+'_time_25240_0'
    99             name6 = join(directory, time_dir, project.scenario_name)+'_time_31550_0'
    100             name7 = join(directory, time_dir, project.scenario_name)+'_time_37860_0'
    101             name8 = join(directory, time_dir, project.scenario_name)+'_time_44170_0'
    102             name9 = join(directory, time_dir, project.scenario_name)+'_time_50480_0'
    103             name10 = join(directory, time_dir, project.scenario_name)+'_time_56790_0'
     95##            name2 = join(directory, time_dir, project.scenario_name)+'_time_6310_0'
     96##            name3 = join(directory, time_dir, project.scenario_name)+'_time_12620_0'
     97##            name4 = join(directory, time_dir, project.scenario_name)+'_time_18930_0'
     98##            name5 = join(directory, time_dir, project.scenario_name)+'_time_25240_0'
     99##            name6 = join(directory, time_dir, project.scenario_name)+'_time_31550_0'
     100##            name7 = join(directory, time_dir, project.scenario_name)+'_time_37860_0'
     101##            name8 = join(directory, time_dir, project.scenario_name)+'_time_44170_0'
     102##            name9 = join(directory, time_dir, project.scenario_name)+'_time_50480_0'
     103##            name10 = join(directory, time_dir, project.scenario_name)+'_time_56790_0'
    104104
    105105           
    106106
    107             names = [name1, name2, name3,name4,name5,name6,name7,name8,name9,name10]
     107            names = [name1]#, name2, name3,name4,name5,name6,name7,name8,name9,name10]
    108108     
    109109            asc_name = []   
  • anuga_work/production/new_south_wales/batemans_bay/project.py

    r7020 r7045  
    3030tide = 0                # difference between MSL and HAT in metres
    3131zone = 56               # specify zone of model
    32 event_number = 31891    # the event number or the mux file name
     32#event_number = 58242    # Puysegur 1 in 10 000 # the event number or the mux file name
     33event_number = 51436    # New Hebrides 1 in 10 000
     34#event_number = 58349    # Puysegur 1 in 5000
     35#event_number = 58284    # Puysegur 1 in 2000
     36#event_number = 58187    # Puysegur 1 in 1000
     37#event_number = 58113    # Puysegur 1 in 500
     38#event_number = 58025    # Puysegur 1 in 200
     39
    3340alpha = 0.1             # smoothing parameter for mesh
    3441friction=0.01           # manning's friction coefficient
     
    3643finaltime=1000         # final time for simulation
    3744
    38 setup = 'trial'         # This can be one of three values
     45setup = 'final'         # This can be one of three values
    3946                        #    trial - coarsest mesh, fast
    4047                        #    basic - coarse mesh
     
    5259#-------------------------------------------------------------------------------
    5360
    54 output_comment = [setup, tide, event_number, 'lidar']
     61output_comment = [setup, tide, event_number]
    5562
    5663#-------------------------------------------------------------------------------
     
    6269# Format for ascii grids, as produced in ArcGIS + a projection file
    6370ascii_grid_filenames = ['e151_s34_clip2',
    64                         'aus197_topo2raster_umina',
    65                         'aus197_topo2raster_terrigal']   #GA data??
     71                        'aus197_topo2raster_clip',
     72                        'hydro_topo2raster_clip']
     73                        #'aus197_topo2raster_terrigal']   #GA data??
    6674                   
    6775
     
    97105##                   'XYHI341_HSDB_SD100035103_MGA',
    98106##                   'XYHI339_F_HSDB_SD100031669_MGA',
    99                    'AUS197_MGA_AHD_1972.txt',
    100                    'hydro.txt'
     107##                   'AUS197_MGA_AHD_1972.txt',
     108                   'hydro_neg.txt'
    101109                   ]
    102110point_filenames.extend(lidar_filenames)
  • anuga_work/production/new_south_wales/batemans_bay/setup_model.py

    r6970 r7045  
    103103    project.scale_factor = 100
    104104    project.time_thinning = 96
    105     project.step = 240
     105    project.yieldstep = 240
    106106elif project.setup == 'basic':
    107107    print 'basic'
Note: See TracChangeset for help on using the changeset viewer.