Changeset 7608


Ignore:
Timestamp:
Jan 15, 2010, 3:37:10 PM (15 years ago)
Author:
jgriffin
Message:

Added bed shear stress calculation to export_results_max.py
Added more probabilistic events to project.py

Location:
anuga_work/production/new_south_wales/batemans_bay
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • anuga_work/production/new_south_wales/batemans_bay/Arc_asc2raster_GDA94z56.py

    r7577 r7608  
    2525output_dir = "anuga\\outputs\\"
    2626
    27 time_dir1 = '20090601_172248_run_final_0.0_58284__250m_jgriffin'
     27time_dir1 = '20091223_101813_run_final_0.0_6157_jgriffin'
    2828
    2929time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5]#, time_dir6, time_dir7, time_dir8, time_dir9]
    … …  
    3838    #ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_conditions_ocean1"
    3939   
    40 ##    print 'Process: Create File GDB'
    41 ##    gp.CreateFileGDB_management(folder, "raster")
     40    print 'Process: Create File GDB'
     41    gp.CreateFileGDB_management(folder, "raster")
    4242
    4343    gp.Workspace = raster_gbd
    … …  
    5252
    5353    generate_filename = []
    54     input_ascii = glob.glob(folder + '*.asc')
     54    input_ascii = glob.glob(folder + '*max.asc')
    5555    print time_dir
    5656
    … …  
    7878##        output_extract = output_DEM + 'E'
    7979##        print 'Output Extract ',output_extract
    80 ##   
     80####   
    8181##        print 'Process: Extract by Mask'
    8282##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
    83 ##
     83####
    8484##       

  • anuga_work/production/new_south_wales/batemans_bay/export_results_max.py

    r7606 r7608  
    2222directory = project.output_folder
    2323
     24time_dirs =['20090422_162054_run_final_0.0_58368_jgriffin',
     25            '20090422_162114_run_final_0.0_51436_jgriffin',
     26            '20090520_144937_run_final_0.0_58272_jgriffin',
     27            '20090520_145156_run_final_1.0_58272_jgriffin',
     28            '20090520_145453_run_final_0.0_58115_jgriffin',
     29            '20090520_145510_run_final_0.0_58129_jgriffin',
     30            '20090520_145608_run_final_0.0_58226_jgriffin',
     31            '20090520_145616_run_final_0.0_58284_jgriffin',
     32            '20090520_145700_run_final_0.0_58286_jgriffin',
     33            '20090525_093040_run_final_1.0_51445_jgriffin',
     34            '20090528_143200_run_final_1.0_51204_jgriffin',
     35            '20090528_143323_run_final_1.0_58346_jgriffin',
     36            '20090528_153025_run_final_0.0_58346_jgriffin',
     37            '20090528_153109_run_final_0.0_51204_jgriffin',
     38            '20090529_143332_run_final_0.0_51378_jgriffin',
     39            '20090529_143419_run_final_0.0_51077_jgriffin',
     40            '20090529_143442_run_final_0.0_51347_jgriffin',
     41            '20090529_143458_run_final_0.0_51292_jgriffin',
     42            '20090529_143527_run_final_0.0_51424_jgriffin']
     43           
     44
    2445#Specify output directories
    25 time_dir1 = '20090528_153025_run_final_0.0_58346_jgriffin'
     46##time_dirs =[
     47##['20091216_122203_run_final_0.0_58176_jgriffin',
     48##'20091216_122218_run_final_0.0_31997_jgriffin',
     49##'20091216_122348_run_final_0.0_7871_jgriffin',
     50##'20091216_124238_run_final_0.0_31708_jgriffin',
     51##'20091216_124254_run_final_0.0_31981_jgriffin',
     52##'20091216_124340_run_final_0.0_51470_jgriffin',
     53##'20091221_110506_run_final_0.0_58274_jgriffin',
     54##'20091221_110525_run_final_0.0_51270_jgriffin',
     55##'20091221_110609_run_final_0.0_58286_jgriffin']
     56##'20091221_110617_run_final_0.0_51268_jgriffin',
     57##'20091221_110624_run_final_0.0_31831_jgriffin',
     58##'20091221_110632_run_final_0.0_31583_jgriffin',
     59##'20091223_100248_run_final_0.0_5867_jgriffin',
     60##'20091223_100422_run_final_0.0_1964_jgriffin',
     61##'20091223_101813_run_final_0.0_6157_jgriffin',
     62##'20091223_101905_run_final_0.0_6475_jgriffin',
     63##'20091223_102012_run_final_0.0_7865_jgriffin',
     64##'20091223_102123_run_final_0.0_16968_jgriffin',
     65##'20091223_102158_run_final_0.0_20902_jgriffin',
     66##'20091223_102326_run_final_0.0_24247_jgriffin']
     67##'20091223_102727_run_final_0.0_28084_jgriffin',
     68##'20091223_102851_run_final_0.0_30858_jgriffin',
     69##'20091223_102951_run_final_0.0_31441_jgriffin',
     70##'20100105_103540_run_final_0.0_31602_jgriffin',
     71##'20100105_103553_run_final_0.0_51375_jgriffin',
     72##'20100105_103701_run_final_0.0_31897_jgriffin',
     73##'20100105_103701_run_final_0.0_58318_jgriffin',
     74##'20100105_134917_run_final_0.0_51328_jgriffin',
     75##'20100105_134929_run_final_0.0_58331_jgriffin']
     76##'20100105_135108_run_final_0.0_6262_jgriffin',
     77##'20100105_135137_run_final_0.0_6512_jgriffin',
     78##'20100105_135246_run_final_0.0_16990_jgriffin',
     79##'20100105_135253_run_final_0.0_20993_jgriffin',
     80##'20100105_135319_run_final_0.0_24300_jgriffin',
     81##'20100105_135401_run_final_0.0_30207_jgriffin',
     82##'20100105_135419_run_final_0.0_31012_jgriffin',
     83##'20100105_135440_run_final_0.0_31536_jgriffin',
     84##'20100105_150548_run_final_0.0_15491_jgriffin',
     85##'20100105_170027_run_final_0.0_6143_jgriffin']
     86
     87##time_dirs =['20100109_002128_run_final_0.0_31593_jgriffin',
     88##            '20100109_042200_run_final_0.0_17006_jgriffin',
     89##            '20100109_101855_run_final_0.0_31714_jgriffin',
     90##            '20100109_103842_run_final_0.0_6360_jgriffin',
     91##            '20100109_105331_run_final_0.0_7741_jgriffin',
     92##            '20100109_111451_run_final_0.0_58337_jgriffin',
     93##            '20100109_124401_run_final_0.0_15505_jgriffin',
     94##            '20100109_132407_run_final_0.0_30379_jgriffin',
     95##            '20100109_142643_run_final_0.0_31173_jgriffin',
     96##            '20100109_145548_run_final_0.0_24356_jgriffin',
     97##            '20100109_150711_run_final_0.0_51436_jgriffin',
     98##            '20100109_162600_run_final_0.0_22691_jgriffin']
     99           
     100
     101
    26102##time_dir2 = '20090529_143458_run_final_0.0_58346_jgriffin'
    27103
    28 time_dirs = [time_dir1]#, time_dir2]   
     104#time_dirs = [time_dir1]#, time_dir2]   
    29105
    30106# sww filename extensions ie. if batemans_bay_time_37860_0.sww, input into list 37860
    … …  
    46122######
    47123
    48 #area = ['', '']  # strings must match keys in var_equations below
    49 area = ['All']      # 'All' means no special areas - the whole thing
     124area = ['AOI']  # strings must match keys in var_equations below
     125##area = ['All']      # 'All' means no special areas - the whole thing
    50126
    51127######
    … …  
    53129# This would not normally change.
    54130######
     131# Added bed shear stress quantity
     132# Note that mannings n (friction value) is taken as 0.01, as in the model run
     133# density of water is 1000
     134
     135
     136# Quantities
    55137var_equations = {'stage':     'stage',
    56138                 'momentum':  '(xmomentum**2 + ymomentum**2)**0.5',
    … …  
    58140                 'speed':     '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6)',
    59141                 'energy':  '(((xmomentum/(stage-elevation+1.e-6))**2 + (ymomentum/(stage-elevation+1.e-6))**2)*0.5*1000*(stage-elevation+1.e-6))+(9.81*stage*1000)',
     142                 'bed_shear_stress': ('(((1/depth**(7./3.))*1000*9.81*0.01**2*(xmomentum/(stage-elevation+1.e-6))*((xmomentum/(stage-elevation+1.e-6))**2+(ymomentum/(stage-elevation+1.e-6))**2)**0.5)**2+'
     143                                      '((1/depth**(7./3.))*1000*9.81*0.01**2*(ymomentum/(stage-elevation+1.e-6))*((xmomentum/(stage-elevation+1.e-6))**2+(ymomentum/(stage-elevation+1.e-6))**2)**0.5)**2)**0.5'),
    60144                 'elevation': 'elevation' }
    61145
    62146# one or more key strings from var_equations above
    63147#var = ['depth', 'speed','stage']
    64 var = ['energy']
     148var = ['bed_shear_stress']
    65149######
    66150# Start script, running through variables, area, folder, sww file (determine by times)

  • anuga_work/production/new_south_wales/batemans_bay/project.py

    r7606 r7608  
    5252    event_number = int(sys.argv[1])
    5353else:   
    54     event_number = 58284
     54    event_number = 74057#60194#74552#
    5555##event_number_list = [31602,31714,31913] #7871,31583,#1.5-2.0m
    5656##event_number_list = [51375,51436,51452] #31997,51268, #1.5-2.0m
    … …  
    7676
    7777# Re-run those that didn't work first time
    78 event_number_list = [15442,5890,6017,6052]
     78event_number_list = [5890,6017,6052]#15442,
     79
     80# more events start 12 Jan 2010 #0.3-1.0m
     81##event_number_list = [31852,31940,31990,35179]
     82##event_number_list = [35377,35451,35542,35557]
     83##event_number_list = [35598,35652,41034,41217]
     84##event_number_list = [41281,41340,41395,41415]
     85##event_number_list = [45841,46039,46221,46297]
     86##event_number_list = [46439,46536,46533,46642]
     87##event_number_list = [46681,49097,49118,49123]
     88##event_number_list = [50621,50797,50873]
     89##event_number_list = [49591,50998,51116,51193]
     90##event_number_list = [51224,51299,51313,51314]
     91##event_number_list = [51319,51336,51351,51389]
     92##event_number_list = [51401,51421,51437,51455]
     93##event_number_list = [53713,53789,53823,53858]
     94##event_number_list = [53872,53901,53910,53917]
     95##event_number_list = [53948,53971,53997,57150]
     96##event_number_list = [57262,57277,57385,57424]
     97##event_number_list = [57478,57560,56566,57573]
     98##event_number_list = [57608,57610,57942,57987]
     99##event_number_list = [58022,58068,58100,58105]
     100##event_number_list = [58146,58172,58180,58201]
     101##event_number_list = [58216,58245,58287,58294]
     102##event_number_list = [58321,60110,60143,60165]
     103##event_number_list = [60194,60209,60232,60248]
     104##event_number_list = [61850,61882,61912,61916]
     105##event_number_list = [61931,61944,61955,62674]
     106##event_number_list = [62746,62759,62780,62792]
     107##event_number_list = [62821,62848,62861,62882]
     108##event_number_list = [62904,63731,63736,72010]
     109##event_number_list = [72022,72036,73960,74015]
     110##event_number_list = [74057,74078,74079,74123]
     111##event_number_list = [74132,74146,74552]
     112
     113# Re-run those that didn't work first time
     114##event_number_list = [60209,74057,74146]
     115##event_number_list = [53971,58068]
     116##event_number_list = [58022,53948]#74552
     117##event_number_list = [60194,62904]
     118##event_number_list = [72022,74057]
    79119
    80120
    … …  
    184224
    185225# ASCII export grid for Batemans Bay
    186 ##xminBatemans_Bay =
    187 ##xmaxBatemans_Bay =
    188 ##yminBatemans_Bay =
    189 ##ymaxBatemans_Bay =
     226xminAOI = 239690
     227xmaxAOI = 255206
     228yminAOI = 6036317
     229ymaxAOI = 6050529
    190230
    191231
Note: See TracChangeset for help on using the changeset viewer.