Changeset 7132


Ignore:
Timestamp:
May 29, 2009, 8:52:31 AM (11 years ago)
Author:
jgriffin
Message:
 
Location:
anuga_work/production/new_south_wales
Files:
5 edited

Legend:

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

    r6597 r7132  
    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]
     20time_dir1 = '20090520_144937_run_final_0.0_58272_jgriffin'
     21time_dir2 = '20090520_145453_run_final_0.0_58115_jgriffin'
     22time_dir3 = '20090520_145156_run_final_1.0_58272_jgriffin'
     23time_dir4 = '20090520_145208_run_final_0.0_51445_jgriffin'
     24time_dir5 = '20090520_145510_run_final_0.0_58129_jgriffin'
     25time_dir6 = '20090520_145608_run_final_0.0_58226_jgriffin'
     26time_dir7 = '20090520_145616_run_final_0.0_58284_jgriffin'
     27time_dir8 = '20090520_145700_run_final_0.0_58286_jgriffin'
     28time_dir9 = '20090525_093040_run_final_1.0_51445_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_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6, time_dir7, time_dir8, time_dir9]
     32
    3333
    3434for time_dir in time_dirs:
     
    3636    gauge = project.gauges
    3737    sww2csv_gauges(name+'.sww',gauge,
    38                    quantities = ['stage','speed','depth','elevation'],
     38                   quantities = ['stage','speed','depth'],
    3939                   verbose=True)
    4040 
  • anuga_work/production/new_south_wales/batemans_bay/project.py

    r7070 r7132  
    2626# Model specific parameters.
    2727# One or all can be changed each time the run_model script is executed
    28 tide = 1.0                # difference between MSL and HAT (1.0)
     28tide = 0.0                # difference between MSL and HAT (1.0)
    2929
    3030# the event number or the mux file name
     
    3434##event_number = 58284    #1 in 2000 yr
    3535##event_number = 58286    #1 in 5000 yr
    36 ##event_number = 58272      #1 in 10000 yr Puysegur
    37 event_number = 51445     #1 in 10000 yr New Hebrides
     36####event_number = 58272      #1 in 10000 yr Puysegur????
     37####event_number = 51445     #1 in 10000 yr New Hebrides????
     38##event_number = 58346      #1 in 10000 yr Puysegur
     39event_number = 51204     #1 in 10000 yr New Hebrides
    3840##event_number = 58368    #1 in 100 000 yr Puysegur!!!
    3941##event_number = 51436    #1 in 100 000 yr New Hebrides!!!
  • anuga_work/production/new_south_wales/gosford/export_results_max.py

    r7069 r7132  
    2929##time_dirs = [time_dir1, time_dir2]#, time_dir3]
    3030
    31 time_dir1 = '20090521_120108_run_final_1.0_58242_jgriffin'
    32 ##time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt'
     31time_dir1 = '20090526_091618_run_final_0.0_58025_jgriffin'
     32##time_dir2 = '20090526_091528_run_final_0.0_58113_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/gosford/prepare_data_hydro.py

    r7069 r7132  
    1515
    1616#folder = '//nas/gemd/georisk_models/inundation/data/new_south_wales/gosford_tsunami_scenario_2009/original_reference_data/GosfordLidar/Ground/'
    17 folder = '//nas/gemd/georisk_models/inundation/data/new_south_wales/gosford_tsunami_scenario_2009/anuga/topographies/original/'
     17folder = '//nas/gemd/georisk_models/inundation/data/new_south_wales/gosford_tsunami_scenario_2009/elevation/'
    1818
    1919
    2020#input_grd = glob.glob(folder + '*.grd')
    21 input_grd = [folder+ 'hydro.txt']
     21input_grd = [folder+ 'hydro_clip.csv']
    2222
    2323for file_in in input_grd:
     
    4747        fields = line.split(',')
    4848
    49         X[i] = float(fields[0])  #easting
    50         Y[i] = float(fields[1])  #northing
    51         Z[i] = -1*abs(float(fields[2]))  #elevation
     49        X[i] = float(fields[2])  #easting
     50        Y[i] = float(fields[3])  #northing
     51        Z[i] = -1*abs(float(fields[4]))  #elevation
    5252   
    5353
  • anuga_work/production/new_south_wales/gosford/project.py

    r7069 r7132  
    2828# Model specific parameters.
    2929# One or all can be changed each time the run_model script is executed
    30 tide = 1.0                # difference between MSL and HAT in metres
     30tide = 0.0                # difference between MSL and HAT in metres
    3131zone = 56               # specify zone of model
    32 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 = 58242    # Puysegur 1 in 10 000 # the event number or the mux file name
     33event_number = 51436    # New Hebrides 1 in 10 000
    3434#event_number = 58349    # Puysegur 1 in 5000
    3535#event_number = 58284    # Puysegur 1 in 2000
    3636#event_number = 58187    # Puysegur 1 in 1000
    3737#event_number = 58113    # Puysegur 1 in 500
    38 ##event_number = 58025    # Puysegur 1 in 200
     38#event_number = 58025    # Puysegur 1 in 200
    3939
    4040
     
    7070# Format for ascii grids, as produced in ArcGIS + a projection file
    7171ascii_grid_filenames = ['e151_s34_clip2',
    72                         'aus197_topo2raster_clip',
    73                         'hydro_topo2raster_clip']
     72                        'aus197_topo2raster_clip2',
     73                        'hydro_topo2raster_clip2']
    7474                        #'aus197_topo2raster_terrigal']   #GA data??
    7575                   
     
    8282topographies_folder = join(anuga_folder, 'topographies/original')
    8383
     84#This reads a file containing lidar filenames and appends to a list
    8485file_list = join(topographies_folder, 'file_list.csv')
    85 ###file_list2 = join(topographies_folder, 'file_list2.csv')
    8686fid = csv.reader(open(file_list))
    87 ###fid2 = csv.reader(open(file_list2))
    8887lidar_filenames = []
    8988for row in fid:
    9089    lidar_filenames.append(','.join(row))
    91 ##for row in fid2:
    92 ##    lidar_filenames.append(','.join(row))
    93 ###print lidar_filenames
    94 #lidar_filenames = list(lidar_filenames)
    95 #print lidar_filenames
     90
    9691
    9792# Format for point is x,y,elevation (with header)
     
    107102##                   'XYHI339_F_HSDB_SD100031669_MGA',
    108103##                   'AUS197_MGA_AHD_1972.txt',
    109                    'hydro_neg.txt'
     104                   'hydro_clip_neg.txt'
    110105                   ]
    111106point_filenames.extend(lidar_filenames)
     
    128123# Format for points easting,northing (no header)
    129124
    130 interior_regions_data =[['aoi_umina.csv', 500],
     125interior_regions_data =[['aoi_umina_large.csv', 500],
    131126                        ['aoi_terrigal.csv', 500],
    132                         ['aos_umina.csv', 2500],
     127                        ['aos_umina_large.csv', 2500],
    133128                        ['aos_terrigal.csv', 2500]]
    134129
Note: See TracChangeset for help on using the changeset viewer.