Changeset 7225


Ignore:
Timestamp:
Jun 19, 2009, 12:38:12 PM (15 years ago)
Author:
jgriffin
Message:
 
Location:
anuga_work/production/new_south_wales
Files:
5 added
10 edited

Legend:

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

    r7076 r7225  
    2525output_dir = "anuga\\outputs\\"
    2626
    27 ##time_dir1 = '20081209_155431_run_final_0_27255_250m_none_lfountai'
    28 ##time_dir2 = '20081210_100528_run_final_0_68693_250m_none_kvanputt'
    29 ##time_dirs = [time_dir1, time_dir2]
    30 
    31 time_dir1 = '20090520_144937_run_final_0.0_58272_jgriffin'
    32 time_dir2 = '20090520_145453_run_final_0.0_58115_jgriffin'
    33 time_dir3 = '20090520_145156_run_final_1.0_58272_jgriffin'
    34 time_dir4 = '20090520_145208_run_final_0.0_51445_jgriffin'
    35 time_dir5 = '20090520_145510_run_final_0.0_58129_jgriffin'
    36 time_dir6 = '20090520_145608_run_final_0.0_58226_jgriffin'
    37 time_dir7 = '20090520_145616_run_final_0.0_58284_jgriffin'
    38 time_dir8 = '20090520_145700_run_final_0.0_58286_jgriffin'
    39 time_dir9 = '20090525_093040_run_final_1.0_51445_jgriffin'
     27time_dir5 = '20090529_143332_run_final_0.0_51378_jgriffin'
     28time_dir1 = '20090529_143419_run_final_0.0_51077_jgriffin'
     29time_dir2 = '20090529_143442_run_final_0.0_51347_jgriffin'
     30time_dir3 = '20090529_143458_run_final_0.0_51292_jgriffin'
     31time_dir4 = '20090529_143527_run_final_0.0_51424_jgriffin'
    4032
    4133
    42 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6, time_dir7, time_dir8, time_dir9]
     34time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5]#, time_dir6, time_dir7, time_dir8, time_dir9]
    4335
    4436for time_dir in time_dirs:
  • anuga_work/production/new_south_wales/batemans_bay/export_results_max.py

    r7139 r7225  
    2929##time_dirs = [time_dir1, time_dir2]#, time_dir3]
    3030
    31 time_dir1 = '20090528_153109_run_final_0.0_51204_jgriffin'
    32 ##time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt'
     31time_dir1 = '20090529_143527_run_final_0.0_51424_jgriffin'
     32##time_dir2 = '20090529_143458_run_final_0.0_51292_jgriffin'
    3333##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt'
    3434##time_dir4 = '20081211_162433_run_final_0.6_68693_alpha0.1_kvanputt'
     
    3636##time_dir6 = '20081211_162744_run_final_0_27283_alpha0.1_kvanputt'
    3737##
    38 time_dirs = [time_dir1] #5, time_dir4]# , time_dir4, time_dir5, time_dir6]
     38time_dirs = [time_dir1]#,time_dir2] #5, time_dir4]# , time_dir4, time_dir5, time_dir6]
    3939 
    4040
  • anuga_work/production/new_south_wales/batemans_bay/get_timeseries.py

    r7132 r7225  
    1818directory = project.output_folder
    1919
    20 time_dir1 = '20090520_144937_run_final_0.0_58272_jgriffin'
    21 time_dir2 = '20090520_145453_run_final_0.0_58115_jgriffin'
    22 time_dir3 = '20090520_145156_run_final_1.0_58272_jgriffin'
    23 time_dir4 = '20090520_145208_run_final_0.0_51445_jgriffin'
    24 time_dir5 = '20090520_145510_run_final_0.0_58129_jgriffin'
    25 time_dir6 = '20090520_145608_run_final_0.0_58226_jgriffin'
    26 time_dir7 = '20090520_145616_run_final_0.0_58284_jgriffin'
    27 time_dir8 = '20090520_145700_run_final_0.0_58286_jgriffin'
    28 time_dir9 = '20090525_093040_run_final_1.0_51445_jgriffin'
     20time_dir1 = '20090529_143419_run_final_0.0_51077_jgriffin'
     21time_dir2 = '20090529_143442_run_final_0.0_51347_jgriffin'
     22time_dir3 = '20090529_143458_run_final_0.0_51292_jgriffin'
     23time_dir4 = '20090529_143527_run_final_0.0_51424_jgriffin'
     24time_dir5 = '20090529_143332_run_final_0.0_51378_jgriffin'
     25
     26#time_dir1 = '20090528_143200_run_final_1.0_51204_jgriffin'
     27#time_dir2 = '20090528_143323_run_final_1.0_58346_jgriffin'
     28#time_dir3 = '20090528_153025_run_final_0.0_58346_jgriffin'
     29#time_dir4 = '20090528_153109_run_final_0.0_51204_jgriffin'
     30##time_dir5 = '20090520_145510_run_final_0.0_58129_jgriffin'
     31##time_dir6 = '20090520_145608_run_final_0.0_58226_jgriffin'
     32##time_dir7 = '20090520_145616_run_final_0.0_58284_jgriffin'
     33##time_dir8 = '20090520_145700_run_final_0.0_58286_jgriffin'
     34##time_dir9 = '20090525_093040_run_final_1.0_51445_jgriffin'
    2935
    3036
    31 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6, time_dir7, time_dir8, time_dir9]
     37time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5]#, time_dir6, time_dir7, time_dir8, time_dir9]
    3238
    3339
  • anuga_work/production/new_south_wales/batemans_bay/project.py

    r7139 r7225  
    4343##event_number = 51204     #1 in 10000 yr New Hebrides
    4444
    45 ######event_number = 58368    #1 in 100 000 yr Puysegur!!!
    46 ######event_number = 51436    #1 in 100 000 yr New Hebrides!!!
    47 ######event_number = 58272      #1 in 10000 yr Puysegur????
    48 ######event_number = 51445     #1 in 10000 yr New Hebrides????
     45######event_number = 58368    #1 in 100 000 yr Puysegur
     46######event_number = 51436    #1 in 100 000 yr New Hebrides
     47######event_number = 58272      #1 in ~15000 yr Puysegur
     48######event_number = 51445     #1 in ~15000 yr New Hebrides
    4949
    5050alpha = 0.1             # smoothing parameter for mesh
     
    9696                   
    9797# Format for point is x,y,elevation (with header)
    98 point_filenames = ['SD100031996_jgriffin_clip.csv',
     98point_filenames = ['SD100031996_clip.csv',
    9999                   'tomaga_offshore_AHD_MGA_1997.csv',
    100100                   'Batemans_BBHD_MGA_1995.csv',
     
    136136# Format latitude,longitude etc (geographic)
    137137##building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     138
     139# AREA OF IMAGES - Extent of each image to find out highest runup
     140# Header - easting,northing,id,value
     141# Used in get_runup.py
     142images_filename = 'images.csv'
    138143
    139144# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
     
    272277##building_exposure = join(gauges_folder, building_exposure_filename)
    273278
     279# The absolute pathname for the image file
     280# Used for get_runup.py
     281if images_filename:
     282    images = join(polygons_folder, images_filename)
     283
    274284# full path to where MUX files (or meta-files) live
    275285mux_input = join(event_folder, mux_input_filename)
  • anuga_work/production/new_south_wales/gosford/Arc_asc2raster_GDA94z56.py

    r7069 r7225  
    2929##time_dirs = [time_dir1, time_dir2]
    3030
    31 time_dir1 = '20090521_114544_run_final_0.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'
    36 ##time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'
     31time_dir1 = '20090615_160527_run_final_1.0_51436_jgriffin'
     32time_dir2 = '20090615_160610_run_final_0.0_51436_jgriffin'
     33time_dir3 = '20090615_160900_run_final_0.0_58242_jgriffin'
     34time_dir4 = '20090615_160939_run_final_1.0_58242_jgriffin'
     35time_dir5 = '20090615_161039_run_final_0.0_58284_jgriffin'
     36time_dir6 = '20090615_161145_run_final_0.0_58349_jgriffin'
     37time_dir7 = '20090615_161318_run_final_0.0_58187_jgriffin'
     38time_dir8 = '20090615_161339_run_final_0.0_58113_jgriffin'
     39time_dir9 = '20090615_161420_run_final_0.0_58025_jgriffin'
     40##time_dir10 = '20090527_103636_run_final_0.0_58242_jgriffin'
     41##time_dir11 = '20090527_103718_run_final_1.0_58242_jgriffin'
     42##time_dir12 = '20090527_103748_run_final_1.0_51436_jgriffin'
     43##time_dir13 = '20090527_103817_run_final_0.0_51436_jgriffin'
    3744
    38 
    39 time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5]#, time_dir6]
     45time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6, time_dir7, time_dir8, time_dir9]#, time_dir10, time_dir11, time_dir12, time_dir13]
    4046
    4147for time_dir in time_dirs:
  • anuga_work/production/new_south_wales/gosford/build_elevation.py

    r7069 r7225  
    5555        convert_dem_from_ascii2netcdf(absolute_filename,
    5656                                      basename_out=absolute_filename,
    57                                       use_cache=False,
     57                                      use_cache=True,
    5858                                      verbose=True)
    5959        dem2pts(absolute_filename, use_cache=True, verbose=True)
  • anuga_work/production/new_south_wales/gosford/export_results_max.py

    r7132 r7225  
    2929##time_dirs = [time_dir1, time_dir2]#, time_dir3]
    3030
    31 time_dir1 = '20090526_091618_run_final_0.0_58025_jgriffin'
    32 ##time_dir2 = '20090526_091528_run_final_0.0_58113_jgriffin'
     31time_dir1 = '20090615_160610_run_final_0.0_51436_jgriffin'
     32##time_dir2 = '20090615_161318_run_final_0.0_58187_jgriffin'
    3333##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt'
    3434##time_dir4 = '20081211_162433_run_final_0.6_68693_alpha0.1_kvanputt'
     
    3636##time_dir6 = '20081211_162744_run_final_0_27283_alpha0.1_kvanputt'
    3737##
    38 time_dirs = [time_dir1]# time_dir2] #5, time_dir4]# , time_dir4, time_dir5, time_dir6]
     38time_dirs = [time_dir1]#, time_dir2] #5, time_dir4]# , time_dir4, time_dir5, time_dir6]
    3939 
    4040
     
    6363
    6464# one or more key strings from var_equations above
    65 var = ['depth', 'speed','stage']
    66 
     65#var = ['depth', 'speed','stage']
     66var = ['depth', 'speed', 'stage']
    6767######
    6868# Start running the various conversions we require.
  • anuga_work/production/new_south_wales/gosford/get_timeseries.py

    r7069 r7225  
    1818directory = project.output_folder
    1919
    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]
     20##time_dir1 = '20090615_160527_run_final_1.0_51436_jgriffin'
     21##time_dir2 = '20090615_160610_run_final_0.0_51436_jgriffin'
     22##time_dir3 = '20090615_160900_run_final_0.0_58242_jgriffin'
     23##time_dir4 = '20090615_160939_run_final_1.0_58242_jgriffin'
     24##time_dir5 = '20090615_161039_run_final_0.0_58284_jgriffin'
     25##time_dir6 = '20090615_161145_run_final_0.0_58349_jgriffin'
     26##time_dir7 = '20090615_161318_run_final_0.0_58187_jgriffin'
     27##time_dir8 = '20090615_161339_run_final_0.0_58113_jgriffin'
     28time_dir9 = '20090615_161420_run_final_0.0_58025_jgriffin'
    2829
    29 ##time_dir1 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt'
    30 time_dir6 = '20090306_165253_run_final_0_58368_jgriffin'
    3130
    32 time_dirs = [time_dir6]
     31time_dirs = [time_dir9]
     32##time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6, time_dir7, time_dir8, time_dir9]
     33##time_dirs = [time_dir7, time_dir8, time_dir9, time_dir10, time_dir11, time_dir12, time_dir13]
     34
    3335
    3436for time_dir in time_dirs:
     
    3638    gauge = project.gauges
    3739    sww2csv_gauges(name+'.sww',gauge,
    38                    quantities = ['stage','speed','depth','elevation'],
     40                   quantities = ['stage','speed','depth'],
    3941                   verbose=True)
    4042 
  • anuga_work/production/new_south_wales/gosford/project.py

    r7132 r7225  
    99from os.path import join, exists
    1010import csv
     11#import run_model_later
    1112
    1213
     
    2526#-------------------------------------------------------------------------------
    2627
    27 
    28 # Model specific parameters.
    29 # One or all can be changed each time the run_model script is executed
    30 tide = 0.0                # difference between MSL and HAT in metres
     28tide = 0.0   
     29             # difference between MSL and HAT in metres
    3130zone = 56               # specify zone of model
    3231#event_number = 58242    # Puysegur 1 in 10 000 # the event number or the mux file name
    33 event_number = 51436    # New Hebrides 1 in 10 000
     32#event_number = 51436    # New Hebrides 1 in 10 000
    3433#event_number = 58349    # Puysegur 1 in 5000
    3534#event_number = 58284    # Puysegur 1 in 2000
    3635#event_number = 58187    # Puysegur 1 in 1000
    3736#event_number = 58113    # Puysegur 1 in 500
    38 #event_number = 58025    # Puysegur 1 in 200
    39 
     37event_number = 58025    # Puysegur 1 in 200
    4038
    4139alpha = 0.1             # smoothing parameter for mesh
     
    144142##building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
    145143
     144# AREA OF IMAGES - Extent of each image to find out highest runup
     145# Header - easting,northing,id,value
     146# Used in get_runup.py
     147images_filename = 'images.csv'
     148
    146149# BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
    147150# NOTE: when files are put together the points must be in sequence
     
    279282##building_exposure = join(gauges_folder, building_exposure_filename)
    280283
     284# The absolute pathname for the image file
     285# Used for get_runup.py
     286if images_filename:
     287    images = join(polygons_folder, images_filename)
     288
    281289# full path to where MUX files (or meta-files) live
    282290mux_input = join(event_folder, mux_input_filename)
  • anuga_work/production/new_south_wales/gosford/run_model.py

    r7069 r7225  
    176176    print domain.timestepping_statistics()
    177177    print domain.boundary_statistics(tags='ocean')
     178    print domain.volumetric_balance_statistics()
    178179
    179180print 'Simulation took %.2f seconds' % (time.time()-t0)
Note: See TracChangeset for help on using the changeset viewer.