Changeset 7606


Ignore:
Timestamp:
Jan 12, 2010, 1:30:41 PM (8 years ago)
Author:
jgriffin
Message:

Added extra events for probabilistic run

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

Legend:

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

    r7601 r7606  
    1515from anuga.shallow_water.data_manager import urs2sts
    1616
    17 import numpy as num
     17import Numeric as num
    1818
    1919
  • anuga_work/production/new_south_wales/batemans_bay/export_results_max.py

    r7577 r7606  
    2323
    2424#Specify output directories
    25 time_dir1 = '20090601_172248_run_final_0.0_58284__250m_jgriffin'
     25time_dir1 = '20090528_153025_run_final_0.0_58346_jgriffin'
    2626##time_dir2 = '20090529_143458_run_final_0.0_58346_jgriffin'
    2727
     
    3030# sww filename extensions ie. if batemans_bay_time_37860_0.sww, input into list 37860
    3131# make sure numbers are in sequential order
    32 times = []
     32times = [37860]
    3333
    3434#Modify the cellsize value to set the size of the raster you require
    3535#Take into account mesh size when aplying this paramater
    36 cellsize = 250   #dependent on data resolution in area of interest.
     36cellsize = 20   #dependent on data resolution in area of interest.
    3737
    3838#Now set the timestep at which you want the raster generated.
     
    5757                 'depth':     'stage-elevation',
    5858                 'speed':     '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6)',
     59                 '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)',
    5960                 'elevation': 'elevation' }
    6061
    6162# one or more key strings from var_equations above
    6263#var = ['depth', 'speed','stage']
    63 var = ['stage','speed','elevation']
     64var = ['energy']
    6465######
    6566# Start script, running through variables, area, folder, sww file (determine by times)
  • anuga_work/production/new_south_wales/batemans_bay/project.py

    r7600 r7606  
    5353else:   
    5454    event_number = 58284
    55 #event_number_list = [31583,31602,31714,31913] #7871,#1.5-2.0m
    56 ##event_number_list = [51268,51375,51436,51452] #31997, #1.5-2.0m
    57 ##event_number_list =  [58286,58318,58337,58350] #58176,#1.5-2.0m
    58 
    59 ##event_number_list = [31831,31897,31961,31973]#31708,
    60 ##event_number_list = [51270,51328,51415,51460]#31981,
    61 ##event_number_list = [58274,58331,58348,58357,58360]#51470,
     55##event_number_list = [31602,31714,31913] #7871,31583,#1.5-2.0m
     56##event_number_list = [51375,51436,51452] #31997,51268, #1.5-2.0m
     57##event_number_list =  [58318,58337,58350] #58176,58286,#1.5-2.0m
     58##event_number_list = [31897,31961,31973]#31708,31831,#2.1-3.0m
     59
     60##event_number_list = [51328,51415,51460]#31981,51270,#2.1-3.0m
     61##event_number_list = [58331,58348,58357,58360]#51470,58274,#2.1-3.0m
     62   
     63##event_number_list = [2009,2045,2060,5716] #1964,#0.3-1.0m
     64##event_number_list = [5890,6017,6052,6143] #5867,#0.3-1.0m
     65   
     66##event_number_list = [6262,6360,6358,6464] #6157,#0.3-1.0m
     67##event_number_list = [6512,7741,7771,7813] #6475,#0.3-1.0m
     68##event_number_list = [15442,15491,15505,16964]#7865,#0.3-1.0m
     69##event_number_list = [16990,17006,17020,20766]#16968,#0.3-1.0m
     70
     71##event_number_list = [20993,22691,22784,22847] #20902,#0.3-1.0m
     72##event_number_list = [24300,24356,24382,28082] #24247,#0.3-1.0m
     73##event_number_list = [30207,30379,30532,30541] #28084,#0.3-1.0m
     74##event_number_list = [31012,31173,31177,31361] #30858,#0.3-1.0m
     75##event_number_list = [31536,31593,31964,31803]#31441,#0.3-1.0m
     76
     77# Re-run those that didn't work first time
     78event_number_list = [15442,5890,6017,6052]
    6279
    6380
  • anuga_work/production/new_south_wales/batemans_bay/run_multiple_events.py

    r7592 r7606  
    99for event_number in project.event_number_list:
    1010    s = 'python2.5 run_model.py %d' % event_number
    11     print s
     11    #print s
    1212    os.system(s)
Note: See TracChangeset for help on using the changeset viewer.