Changeset 4581


Ignore:
Timestamp:
Jul 3, 2007, 4:28:38 PM (16 years ago)
Author:
sexton
Message:

minor updates for model interrogation

Location:
anuga_work/production
Files:
16 edited

Legend:

Unmodified
Added
Removed

  • anuga_work/production/broome_2006/export_results.py

    r4504 r4581  
    55from os import sep
    66
    7 #time_dir = '20070423_040235_run' #
    8 #time_dir = '20070509_015227_run_final_4.9_nbartzis' # HAT 1 in 10000 yr wave height for Dampier
    9 time_dir = '20070518_060438_run_final_0_nbartzis' # MSL 1 in 10000 yr wave height for Dampier
     7time_dir = '20070615_063022_run_final_4.9_exmouth_nbartzis' # HAT exmouth
     8#time_dir = '20070615_062730_run_final_0_exmouth_nbartzis' # MSL exmouth
     9#time_dir = '20070618_064014_run_final_4.9_onslow_nbartzis' # HAT onslow
     10#time_dir = '20070622_004008_run_final_0_onslow_nbartzis' # MSL onslow
     11#time_dir = '20070615_063137_run_final_4.9_dampier_nbartzis' # HAT dampier
     12#time_dir = '20070621_002922_run_final_0_dampier_nbartzis' # MSL dampier
     13
     14
    1015cellsize = 25
    1116timestep = None
    … …  
    1419name = directory+time_dir+sep+project_urs.scenario_name
    1520
    16 is_parallel = True
     21#is_parallel = True
     22is_parallel = False
    1723if is_parallel == True: nodes = 4
    1824print 'output dir:', name
    1925
    20 var = [2,3,4]
     26var = [3]
    2127
    2228for which_var in var:
    … …  
    6470                    timestep = timestep,
    6571                    cellsize = cellsize,     
    66                     #easting_min = project.e_min_area,
    67                     #easting_max = project.e_max_area,
    68                     #northing_min = project.n_min_area,
    69                     #northing_max = project.n_max_area,       
     72                    easting_min = project_urs.eastingmin,
     73                    easting_max = project_urs.eastingmax,
     74                    northing_min = project_urs.northingmin,
     75                    northing_max = project_urs.northingmax,       
    7076                    reduction = max,
    7177                    verbose = True,

  • anuga_work/production/broome_2006/get_timeseries.py

    r4580 r4581  
    1818
    1919# bulli event
    20 production_dirs = {'20070518_060438_run_final_0_nbartzis': ''}
     20production_dirs = {'20070702_050939_run_final_4.9_exmouth_nbartzis': ''}
    2121
    22 swwfiles = {}   
    23 for label_id in production_dirs.keys():
     22# 20070702_050939_run_final_4.9_exmouth_nbartzis HAT exmouth
     23# 20070615_062730_run_final_0_exmouth_nbartzis # MSL exmouth
     24# 20070618_064014_run_final_4.9_onslow_nbartzis # HAT onslow
     25# 20070622_004008_run_final_0_onslow_nbartzis # MSL onslow
     26# 20070615_063137_run_final_4.9_dampier_nbartzis # HAT dampier
     27# 20070621_002922_run_final_0_dampier_nbartzis # MSL dampier
     28
     29is_parallel = True
     30
     31if is_parallel is True:
     32    no_nodes = 4
    2433   
    25     file_loc = project.outputdir + label_id + sep
    26     swwfile = file_loc + project.basename + '.sww'
    27     swwfiles[swwfile] = label_id
     34if is_parallel is False:
     35   
     36    swwfiles = {}   
     37    for label_id in production_dirs.keys():
    2838       
    29     texname, elev_output = sww2timeseries(swwfiles,
    30                                           project.beach_gauges,
    31                                           production_dirs,
    32                                           report = False,
    33                                           plot_quantity = ['stage', 'speed'],
    34                                           generate_fig = True,
    35                                           surface = False,
    36                                           time_min = None,
    37                                           time_max = None,
    38                                           time_unit = 'hours',
    39                                           title_on = False,
    40                                           use_cache = True,
    41                                           verbose = True)
     39        file_loc = project.output_dir + label_id + sep
     40        swwfile = file_loc + project.scenario_name + '.sww'
     41        swwfiles[swwfile] = label_id
     42           
     43        texname, elev_output = sww2timeseries(swwfiles,
     44                                              project.beach_gauges,
     45                                              production_dirs,
     46                                              report = False,
     47                                              plot_quantity = ['stage', 'speed'],
     48                                              generate_fig = False,
     49                                              surface = False,
     50                                              time_min = None,
     51                                              time_max = None,
     52                                              time_unit = 'hours',
     53                                              title_on = False,
     54                                              use_cache = True,
     55                                              verbose = True)
     56
     57if is_parallel is True:                                     
     58
     59    for count in range(no_nodes):
     60        swwfiles = {}
     61        for label_id in production_dirs.keys():
     62
     63            file_loc = project.output_dir + label_id + sep
     64            extra = '_P%i_%i' %(count,no_nodes)
     65            swwfile = file_loc + project.scenario_name + extra + '.sww'
     66            swwfiles[swwfile] = label_id
     67
     68            texname, elev_output = sww2timeseries(swwfiles,
     69                                                  project.beach_gauges,
     70                                                  production_dirs,
     71                                                  report = False,
     72                                                  plot_quantity = ['stage', 'speed'],
     73                                                  generate_fig = False,
     74                                                  surface = False,
     75                                                  time_min = None,
     76                                                  time_max = None,
     77                                                  time_unit = 'hours',
     78                                                  title_on = False,
     79                                                  use_cache = True,
     80                                                  verbose = True)
    4281                                     

  • anuga_work/production/broome_2006/project.py

    r4580 r4581  
    114114gauge_name = 'broome_gauges.csv'
    115115gauges_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'gauges'+sep
     116beach_gauges = gauges_dir + 'beach_gauges.csv'
    116117gauges_dir_name = gauges_dir + gauge_name
    117118
    … …  
    153154res_3 = 2000*res_factor
    154155
    155 poly_4 = read_polygon(polygons_dir+'Broome_inner_town_pts.csv')
     156#poly_4 = read_polygon(polygons_dir+'Broome_inner_town_pts.csv')
    156157#res_4 = 5000
    157 res_4 = 500*res_factor
     158#res_4 = 500*res_factor
    158159#assert zone == refzone
    159160
    160 interior_regions = [[poly_0,res_0],[poly_1,res_1],[poly_2,res_2]
    161                      ,[poly_3,res_3],[poly_4,res_4]]
     161#interior_regions = [[poly_0,res_0],[poly_1,res_1],[poly_2,res_2]
     162                     #,[poly_3,res_3],[poly_4,res_4]]
    162163
    163 trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
     164#trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
    164165
    165 print 'min number triangles', trigs_min
     166#print 'min number triangles', trigs_min
    166167
    167 poly_mainland = read_polygon(polygons_dir+'Initial_Condition.csv')
     168#poly_mainland = read_polygon(polygons_dir+'Initial_Condition.csv')
    168169
    169170###################################################################
    … …  
    172173
    173174# exporting asc grid
    174 e_min_area = 412000.0
    175 e_max_area = 423000.0
    176 n_min_area = 8007000.0
    177 n_max_area = 8022000.0
     175e_min_area = 412900.0
     176e_max_area = 422700.0
     177n_min_area = 8008555.0
     178n_max_area = 8021350.0
    178179
    179180

  • anuga_work/production/broome_2006/run_building_inundation.py

    r4504 r4581  
    2323
    2424# Inputs
    25 #timestampdir = '20070509_015227_run_final_4.9_nbartzis' # HAT
    26 timestampdir = '20070518_060438_run_final_0_nbartzis' # MSL
    27 file_loc = project_urs.output_dir + timestampdir + sep
     25#time_dir = '20070615_063022_run_final_4.9_exmouth_nbartzis' # HAT exmouth
     26#time_dir = '20070615_062730_run_final_0_exmouth_nbartzis' # MSL exmouth
     27#time_dir = '20070618_064014_run_final_4.9_onslow_nbartzis' # HAT onslow
     28#time_dir = '20070622_004008_run_final_0_onslow_nbartzis' # MSL onslow
     29time_dir = '20070615_063137_run_final_4.9_dampier_nbartzis' # HAT dampier
     30#time_dir = '20070621_002922_run_final_0_dampier_nbartzis' # MSL dampier
     31
     32file_loc = project_urs.output_dir + time_dir + sep
    2833swwfile = file_loc + project_urs.scenario_name #+ '.sww'
    2934buildings_filename = project_urs.buildings_filename

  • anuga_work/production/dampier_2006/export_results.py

    r4498 r4581  
    55from os import sep
    66
    7 time_dir = '20070426_045818_run' # HAT 1 in Dampier 10000 yr - CIPMA scenario
    8 #time_dir = '20070419_065050_run' # HAT 1 in Dampier 10000 yr
    9 #time_dir = '20070419_065018_run' # MSL 1 in Dampier 10000 yr
    10 #time_dir = '' # HAT 1 in Broome 10000 yr
    11 #time_dir = '' # MSL 1 in Broome 10000 yr
    12 #time_dir = '' # HAT 1 in Pt Hedland 10000 yr
    13 #time_dir = '' # MSL 1 in Pt Hedland 10000 yr
    14 #time_dir = '' # HAT 1 in Onslow 10000 yr
    15 #time_dir = '' # MSL 1 in Onslow 10000 yr
    16 #time_dir = '' # HAT 1 in Exmouth 10000 yr
    17 #time_dir = '' # MSL 1 in Exmouth 10000 yr
    18 cellsize = 20
     7#time_dir = '20070426_045818_run' # HAT 1 in Dampier 10000 yr - CIPMA scenario
     8#time_dir = '20070419_065050_run_final_2.4_dampier_nbartzis' # HAT 1 in Dampier 10000 yr
     9#time_dir = '20070419_065018_run_final_0.0_dampier_nbartzis' # MSL 1 in Dampier 10000 yr
     10#time_dir = '20070618_063824_run_final_2.4_onslow_nbartzis' # HAT 1 in Onslow 10000 yr
     11#time_dir = '20070618_050232_run_final_0.0_onslow_nbartzis' # MSL 1 in Onslow 10000 yr
     12#time_dir = '20070620_044413_run_final_2.4_exmouth_nbartzis' # HAT 1 in Exmouth 10000 yr
     13time_dir = '20070620_044352_run_final_0.0_exmouth_nbartzis' # MSL 1 in Exmouth 10000 yr
     14#time_dir = '20070618_063824_run_final_2.4_onslow_nbartzis' # HAT onslow
     15#time_dir = '20070618_050232_run_final_0.0_onslow_nbartzis' # MSL onslow
     16cellsize = 25
    1917timestep = None
    2018directory = project.output_dir
    … …  
    2725
    2826#var = [2,3] # depth and speed
    29 var = [3]
     27var = [0,2,3]
    3028
    3129for which_var in var:
    3230    if which_var == 0:  # Stage
    33         outname = name + '_stage'
     31        outname = name + '_stage_dampier'
    3432        quantityname = 'stage'
    3533
    … …  
    3937
    4038    if which_var == 2:  # Depth
    41         outname = name + '_depth_facility'
     39        outname = name + '_depth_dampier'
    4240        quantityname = 'stage-elevation' 
    4341
    4442    if which_var == 3:  # Speed
    45         outname = name + '_speed_facility'
     43        outname = name + '_speed_dampier'
    4644        quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
    4745
    4846    if which_var == 4:  # Elevation
    49         outname = name + '_elevation'
     47        outname = name + '_elevation_dampier'
    5048        quantityname = 'elevation'  #Elevation
    5149
    … …  
    7573                    timestep = timestep,
    7674                    cellsize = cellsize,     
    77                     #easting_min = project.e_min_area,
    78                     #easting_max = project.e_max_area,
    79                     #northing_min = project.n_min_area,
    80                     #northing_max = project.n_max_area,       
     75                    easting_min = project.e_min_area_d,
     76                    easting_max = project.e_max_area_d,
     77                    northing_min = project.n_min_area_d,
     78                    northing_max = project.n_max_area_d,       
    8179                    reduction = max,
    8280                    verbose = True,

  • anuga_work/production/dampier_2006/project.py

    r4504 r4581  
    117117#gauges
    118118gauge_name = 'dampier_gauges_up2.csv' #'dampier.csv'
     119
    119120gauge_name_simple = 'dampier_gauges_simple.csv'
    120121gauge_name_test = 'dampier_gauges_up_test.csv' #'dampier.csv'
    … …  
    123124gauges_dir_name_simple = gauges_dir + gauge_name_simple
    124125gauges_dir_name_test = gauges_dir + gauge_name_test
     126beach_gauges = gauges_dir + 'beach_gauges.csv'
    125127
    126128community_filename = gauges_dir + 'CHINS_v2.csv'

  • anuga_work/production/dampier_2006/report/dampier_2006_report_CIPMA.tex

    r4477 r4581  
    104104\input{anuga_setup}
    105105
     106\clearpage
     107
    106108   \section{Metadata}
    107109     \label{sec:metadata}
    108110     \input{metadata}
    109      
    110 \section{ANUGA modelling parameters}
    111 \label{sec:anugasetup}
    112 \input{anuga_setup}
    113111
    114112\clearpage

  • anuga_work/production/dampier_2006/run_building_inundation.py

    r4504 r4581  
    2424# Inputs
    2525#timestampdir = '20070419_065050_run' # HAT 1 in 10000 yr Dampier event
    26 timestampdir = '20070419_065018_run' # MSL 1 in 10000 yr Dampier event
     26#timestampdir = '20070419_065018_run' # MSL 1 in 10000 yr Dampier event
     27#timestampdir = '20070618_063824_run_final_2.4_onslow_nbartzis' # HAT onslow
     28#timestampdir = '20070618_050232_run_final_0.0_onslow_nbartzis' # MSL onslow
     29#timestampdir = '20070620_044352_run_final_0.0_exmouth_nbartzis' # MSL exmouth
     30timestampdir = '20070620_044413_run_final_2.4_exmouth_nbartzis' # HAT exmouth
    2731file_loc = project.output_dir + timestampdir + sep
    2832swwfilebasename = file_loc + project.scenario_name #+ '.sww'

  • anuga_work/production/exmouth_2006/get_timeseries.py

    r4580 r4581  
    1818
    1919# bulli event
    20 production_dirs = {'20070528_061822_run_final_1.4_dampier_nbartzis': ''}
     20production_dirs = {'20070619_042043_run_final_1.4_onslow_nbartzis': ''}
     21# 20070618_035101_run_final_1.4_dampier_nbartzis # HAT dampier
     22# 20070618_035114_run_final_0_dampier_nbartzis # MSL dampier
     23# 20070619_042542_run_final_1.4_exmouth_nbartzis # HAT exmouth parallel 4
     24# 20070619_042417_run_final_0_exmouth_nbartzis # MSL exmouth parallel 4
     25# 20070619_042043_run_final_1.4_onslow_nbartzis # HAT onslow
     26# 20070619_042140_run_final_0_onslow_nbartzis # MSL onslow parallel 2
    2127
    22 swwfiles = {}   
    23 for label_id in production_dirs.keys():
     28is_parallel = True
     29
     30if is_parallel is True:
     31    no_nodes = 4
    2432   
    25     file_loc = project.outputdir + label_id + sep
    26     swwfile = file_loc + project.basename + '.sww'
    27     swwfiles[swwfile] = label_id
     33if is_parallel is False:
     34   
     35    swwfiles = {}   
     36    for label_id in production_dirs.keys():
    2837       
    29     texname, elev_output = sww2timeseries(swwfiles,
    30                                           project.beach_gauges,
    31                                           production_dirs,
    32                                           report = False,
    33                                           plot_quantity = ['stage', 'speed'],
    34                                           generate_fig = True,
    35                                           surface = False,
    36                                           time_min = None,
    37                                           time_max = None,
    38                                           time_unit = 'hours',
    39                                           title_on = False,
    40                                           use_cache = True,
    41                                           verbose = True)
     38        file_loc = project.output_dir + label_id + sep
     39        swwfile = file_loc + project.scenario_name + '.sww'
     40        swwfiles[swwfile] = label_id
     41           
     42        texname, elev_output = sww2timeseries(swwfiles,
     43                                              project.beach_gauges,
     44                                              production_dirs,
     45                                              report = False,
     46                                              plot_quantity = ['stage', 'speed'],
     47                                              generate_fig = False,
     48                                              surface = False,
     49                                              time_min = None,
     50                                              time_max = None,
     51                                              time_unit = 'hours',
     52                                              title_on = False,
     53                                              use_cache = True,
     54                                              verbose = True)
     55
     56if is_parallel is True:                                     
     57
     58    for count in range(no_nodes):
     59        swwfiles = {}
     60        for label_id in production_dirs.keys():
     61
     62            file_loc = project.output_dir + label_id + sep
     63            extra = '_P%i_%i' %(count,no_nodes)
     64            swwfile = file_loc + project.scenario_name + extra + '.sww'
     65            swwfiles[swwfile] = label_id
     66
     67            texname, elev_output = sww2timeseries(swwfiles,
     68                                                  project.beach_gauges,
     69                                                  production_dirs,
     70                                                  report = False,
     71                                                  plot_quantity = ['stage', 'speed'],
     72                                                  generate_fig = False,
     73                                                  surface = False,
     74                                                  time_min = None,
     75                                                  time_max = None,
     76                                                  time_unit = 'hours',
     77                                                  title_on = False,
     78                                                  use_cache = True,
     79                                                  verbose = True)
    4280                                     

  • anuga_work/production/exmouth_2006/project.py

    r4580 r4581  
    140140beach_gauges = 'beach_gauges.csv'
    141141gauge_name = 'exmouth_gauges.csv'
     142
    142143gauges_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'gauges'+sep
     144beach_gauges = gauges_dir + 'beach_gauges.csv'
    143145gauges_dir_name = gauges_dir + gauge_name
    144146

  • anuga_work/production/exmouth_2006/run_building_inundation.py

    r4504 r4581  
    2323
    2424# Inputs
    25 #timestampdir = '' # HAT
    26 #timestampdir = '' # MSL
     25#timestamp_dir = '20070615_024233_run_final_1.4_dampier_nbartzis' # HAT dampier
     26#timestamp_dir = '20070615_024322_run_final_0_dampier_nbartzis' # MSL dampier
     27#timestamp_dir = '20070604_015217_run_final_1.4_exmouth_nbartzis' # HAT exmouth parallel 4
     28#timestamp_dir = '20070615_023718_run_final_0_exmouth_nbartzis' # MSL exmouth parallel 2
     29#timestamp_dir = '20070615_023954_run_final_1.4_onslow_nbartzis' # HAT onslow
     30#timestamp_dir = '20070615_023908_run_final_0_onslow_nbartzis' # MSL onslow parallel 2
     31
    2732file_loc = project.output_dir + timestampdir + sep
    2833swwfile = file_loc + project.scenario_name + '.sww'

  • anuga_work/production/onslow_2006/export_results.py

    r4521 r4581  
    55from os import sep
    66
    7 time_dir = '20070507_215717_run_final_1.5_nbartzis' # HAT 1 in Dampier 10000 yr
    8 #time_dir = '20070517_061654_run_final_0.0_nbartzis' # MSL 1 in Dampier 10000 yr
     7#time_dir = '20070613_034754_run_final_1.5_dampier_nbartzis' # HAT 1 in Dampier 10000 yr
     8#time_dir = '20070518_021050_run_final_0.0_nbartzis' # MSL 1 in Dampier 10000 yr
    99#time_dir = '' # HAT 1 in Broome 10000 yr
    1010#time_dir = '' # MSL 1 in Broome 10000 yr
    1111#time_dir = '' # HAT 1 in Pt Hedland 10000 yr
    1212#time_dir = '' # MSL 1 in Pt Hedland 10000 yr
    13 #time_dir = '' # HAT 1 in Onslow 10000 yr
    14 #time_dir = '' # MSL 1 in Onslow 10000 yr
    15 #time_dir = '' # HAT 1 in Exmouth 10000 yr
    16 #time_dir = '' # MSL 1 in Exmouth 10000 yr
     13#time_dir = '20070530_082920_run_final_1.5_onslow_nbartzis' # HAT 1 in Onslow 10000 yr
     14#time_dir = '20070531_002753_run_final_0.0_onslow_nbartzis' # MSL 1 in Onslow 10000 yr
     15time_dir = '20070608_062811_run_final_1.5_exmouth_nbartzis' # HAT 1 in Exmouth 10000 yr
     16#time_dir = '20070605_234905_run_final_0.0_exmouth_nbartzis' # MSL 1 in Exmouth 10000 yr
    1717#time_dir = '20060704_063005' # 2006 simulation HAT Mw9
    1818cellsize = 25
    … …  
    2424#name = directory+time_dir+sep+'source'
    2525
    26 #is_parallel = False
    27 is_parallel = True
     26from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
     27
     28elevgrid = directory+time_dir+sep+'onslow_elevation_town'
     29
     30convert_dem_from_ascii2netcdf(elevgrid, use_cache=True, verbose=True)
     31
     32dem2pts(elevgrid,run_up_elev=True,run_up_loc=True)
     33
     34"""
     35is_parallel = False
     36#is_parallel = True
    2837if is_parallel == True: nodes = 4
    2938print 'output dir:', name
    3039
    31 var = [2,3,4] # depth and speed
    32 #var = [3] # depth and speed
     40#var = [2,3,4] # depth and speed
     41#var = [2] # depth
     42var = [0]
    3343
    3444for which_var in var:
    … …  
    4656
    4757    if which_var == 3:  # Speed
    48         outname = name + '_speed_150'
     58        outname = name + '_speed'
    4959        quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
    5060
    … …  
    7686                    timestep = timestep,
    7787                    cellsize = cellsize,     
    78                     #easting_min = project.e_min_area,
    79                     #easting_max = project.e_max_area,
    80                     #northing_min = project.n_min_area,
    81                     #northing_max = project.n_max_area,       
     88                    easting_min = project.e_min_area,
     89                    easting_max = project.e_max_area,
     90                    northing_min = project.n_min_area,
     91                    northing_max = project.n_max_area,       
    8292                    reduction = max,
    8393                    verbose = True,
    8494                    format = 'asc')
     95
     96"""

  • anuga_work/production/onslow_2006/project.py

    r4580 r4581  
    7575#for MOST
    7676gauge_filename = gaugedir + 'gauge_location_onslow.csv'
    77 gauge_filename_nsw = gaugedir + 'gauges_nsw.csv'
    7877beach_gauges = gaugedir + 'beach_gauges.csv'
     78boundary_gauges = gaugedir + 'boundary_gauges.csv'
    7979td_gauges = gaugedir + 'gauges_td.csv'
     80db_point = gaugedir + 'db_points.csv'
    8081gauge_filename_bindi = gaugedir + 'gauge_location_bindi.csv'
    8182gauges50 = gaugedir + '50_gauges.xya'
    … …  
    140141n_min_area = 7601550
    141142n_max_area = 7620000
     143
     144# for test grid for ITEM ppt
     145#e_min_area = 207600
     146#e_max_area = 320000
     147#n_min_area = 7600350
     148#n_max_area = 7691200
    142149
    143150# region to export to make elevation map: JS 22/9/06 - NOTE, this won't

  • anuga_work/production/pt_hedland_2006/export_results.py

    r4504 r4581  
    66
    77#time_dir = '20070518_020117_run_final_3.6_nbartzis' # HAT 1 in Dampier 10000 yr
    8 time_dir = '20070520_225029_run_final_0.0_nbartzis' # MSL 1 in Dampier 10000 yr
    9 #time_dir = '' # HAT 1 in Broome 10000 yr
    10 #time_dir = '' # MSL 1 in Broome 10000 yr
    11 #time_dir = '' # HAT 1 in Pt Hedland 10000 yr
    12 #time_dir = '' # MSL 1 in Pt Hedland 10000 yr
    13 #time_dir = '' # HAT 1 in Onslow 10000 yr
    14 #time_dir = '' # MSL 1 in Onslow 10000 yr
    15 #time_dir = '' # HAT 1 in Exmouth 10000 yr
    16 #time_dir = '' # MSL 1 in Exmouth 10000 yr
     8#time_dir = '20070520_225029_run_final_0.0_nbartzis' # MSL 1 in Dampier 10000 yr
     9time_dir = '20070613_061053_run_final_3.6_onslow_nbartzis' # HAT 1 in Onslow 10000 yr
     10#time_dir = '20070613_061125_run_final_0.0_onslow_nbartzis' # MSL 1 in Onslow 10000 yr
     11#time_dir = '20070613_044204_run_final_3.6_exmouth_nbartzis' # HAT 1 in Exmouth 10000 yr
     12#time_dir = '20070613_225854_run_final_0.0_exmouth_nbartzis' # MSL 1 in Exmouth 10000 yr
    1713
    1814cellsize = 25
    … …  
    2218name = directory+time_dir+sep+project_urs.scenario_name
    2319
    24 is_parallel = True
     20is_parallel = False
    2521if is_parallel == True: nodes = 4
    2622print 'output dir:', name
    2723
    28 var = [2,3,4] # depth and speed
     24var = [0,2,3,4] # depth and speed
    2925
    3026for which_var in var:
    … …  
    7268                    timestep = timestep,
    7369                    cellsize = cellsize,     
    74                     easting_min = project.e_min_area,
    75                     easting_max = project.e_max_area,
    76                     northing_min = project.n_min_area,
    77                     northing_max = project.n_max_area,       
     70                    easting_min = project_urs.e_min_area,
     71                    easting_max = project_urs.e_max_area,
     72                    northing_min = project_urs.n_min_area,
     73                    northing_max = project_urs.n_max_area,       
    7874                    reduction = max,
    7975                    verbose = True,

  • anuga_work/production/pt_hedland_2006/project_urs.py

    r4544 r4581  
    131131
    132132#gauges
    133 gauge_name = 'gauge_location_port_hedland.csv' 
     133gauge_name = 'gauge_location_port_hedland.csv'
    134134gauge_name_test = 'gauge_checking_test.csv'
    135135gauges_dir = anuga_dir+'gauges'+sep
    136136gauges_dir_name = gauges_dir + gauge_name
    137137gauges_dir_name_test = gauges_dir + gauge_name_test
     138beach_gauges = gauges_dir + 'beach_gauges.csv'
    138139
    139140tide_dir = anuga_dir+'tide_data'+sep

  • anuga_work/production/pt_hedland_2006/run_building_inundation.py

    r4504 r4581  
    2323
    2424# Inputs
    25 #timestampdir = '20070518_020117_run_final_3.6_nbartzis' # HAT
    26 timestampdir = '20070520_225029_run_final_0.0_nbartzis' # MSL
    27 file_loc = project_urs.output_dir + timestampdir + sep
     25
     26#timestamp_dir = '20070518_020117_run_final_3.6_nbartzis' # HAT 1 in Dampier 10000 yr
     27#timestamp_dir = '20070520_225029_run_final_0.0_nbartzis' # MSL 1 in Dampier 10000 yr
     28timestamp_dir = '20070613_061053_run_final_3.6_onslow_nbartzis' # HAT 1 in Onslow 10000 yr
     29#timestamp_dir = '20070613_061125_run_final_0.0_onslow_nbartzis' # MSL 1 in Onslow 10000 yr
     30#timestamp_dir = '20070613_044204_run_final_3.6_exmouth_nbartzis' # HAT 1 in Exmouth 10000 yr
     31#timestamp_dir = '20070613_225854_run_final_0.0_exmouth_nbartzis' # MSL 1 in Exmouth 10000 yr
     32
     33file_loc = project_urs.output_dir + timestamp_dir + sep
    2834swwfile = file_loc + project_urs.scenario_name #+ '.sww'
    2935buildings_filename = project_urs.buildings_filename
Note: See TracChangeset for help on using the changeset viewer.