Changeset 7327 for anuga_work


Ignore:
Timestamp:
Aug 6, 2009, 9:16:24 AM (16 years ago)
Author:
jgriffin
Message:
 
Location:
anuga_work/production/new_south_wales
Files:
8 edited

Legend:

Unmodified
Added
Removed

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

    r7225 r7327  
    2525output_dir = "anuga\\outputs\\"
    2626
    27 time_dir5 = '20090529_143332_run_final_0.0_51378_jgriffin'
    28 time_dir1 = '20090529_143419_run_final_0.0_51077_jgriffin'
    29 time_dir2 = '20090529_143442_run_final_0.0_51347_jgriffin'
    30 time_dir3 = '20090529_143458_run_final_0.0_51292_jgriffin'
    31 time_dir4 = '20090529_143527_run_final_0.0_51424_jgriffin'
     27time_dir1 = '20090529_143527_run_final_0.0_51424_jgriffin'
    3228
    33 
    34 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5]#, time_dir6, time_dir7, time_dir8, time_dir9]
     29time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5]#, time_dir6, time_dir7, time_dir8, time_dir9]
    3530
    3631for time_dir in time_dirs:
    … …  
    4338    #ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_conditions_ocean1"
    4439   
    45     print 'Process: Create File GDB'
    46     gp.CreateFileGDB_management(folder, "raster")
     40##    print 'Process: Create File GDB'
     41##    gp.CreateFileGDB_management(folder, "raster")
    4742
    4843    gp.Workspace = raster_gbd
    … …  
    5752
    5853    generate_filename = []
    59     input_ascii = glob.glob(folder + '*max.asc')
     54    input_ascii = glob.glob(folder + '*elevation_max.asc')
    6055    print time_dir
    6156

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

    r7252 r7327  
    11"""
    22Generates ascii grids of nominated areas -
    3 Input: sww file from run_perth.py
     3Input: sww file from run_model.py
    44       boundaries for grids from project.py
    55Outputs: ascii grids of specified variables
    6 Stored in the 'outputs_dir' folder for respective .sww file
     6Stored in the 'outputs_folder' folder for respective .sww file
    77
    88Note:
    … …  
    2323
    2424#Specify output directories
    25 time_dir1 = '20090601_172248_run_final_0.0_51204_jgriffin'
    26 time_dir2 = '20090529_143458_run_final_0.0_58346_jgriffin'
     25time_dir1 = '20090529_143527_run_final_0.0_51424_jgriffin'
     26#time_dir2 = '20090529_143458_run_final_0.0_58346_jgriffin'
    2727
    28 time_dirs = [time_dir1, time_dir2]   
     28time_dirs = [time_dir1]#, time_dir2]   
    2929
    3030cellsize = 20   #dependent on data resolution in area of interest.
    … …  
    5151
    5252# one or more key strings from var_equations above
    53 var = ['depth', 'speed','stage']
    54 
     53#var = ['depth', 'speed','stage']
     54var = ['elevation']
    5555######
    5656# Start running the various conversions we require.

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

    r7252 r7327  
    2626# Model specific parameters.
    2727# One or all can be changed each time the run_model script is executed
    28 tide = 0.0                # difference between MSL and HAT (1.0)
     28tide = 1.0                # difference between MSL and HAT (1.0)
    2929
    3030# the event number or the mux file name
    … …  
    117117# Used in run_building_inundation.py
    118118# Format latitude,longitude etc (geographic)
    119 ##building_exposure_filename = '' # from NEXIS
     119building_exposure_filename = '' # from NEXIS
    120120
    121121# AREA OF IMAGES - Extent of each image to find out highest runup
    … …  
    152152#-------------------------------------------------------------------------------
    153153
    154 # ASCII export grid for Busselton
    155 xminBusselton = 340000
    156 xmaxBusselton = 352000
    157 yminBusselton = 6271500
    158 ymaxBusselton = 6280000
    159 
    160 # ASCII export grid for Bunbury
    161 xminBunbury = 369000
    162 xmaxBunbury = 381000
    163 yminBunbury = 6308000
    164 ymaxBunbury = 6316500
     154# ASCII export grid for Batemans Bay
     155##xminBatemans_Bay =
     156##xmaxBatemans_Bay =
     157##yminBatemans_Bay =
     158##ymaxBatemans_Bay =
     159
     160
    165161
    166162################################################################################
    … …  
    257253# The absolute pathname for the building file
    258254# Used for run_building_inundation.py
    259 ##building_exposure = join(gauges_folder, building_exposure_filename)
     255building_exposure = join(gauges_folder, building_exposure_filename)
    260256
    261257# The absolute pathname for the image file

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

    r7059 r7327  
    1 """Run a tsunami inundation scenario for Busselton, WA, Australia.
     1"""Run a tsunami inundation scenario for Batemans Bay, NSW, Australia.
    22
    33The scenario is defined by a triangular mesh created from project.polygon, the
    … …  
    165165##
    166166##for t in domain.evolve(yieldstep=5,
    167 ##                       finaltime=6000,
    168 ##                       skip_initial_step=True):
     167##                       finaltime=project.finaltime,
     168##                       skip_initial_step=False):
    169169##    print domain.timestepping_statistics()
    170170##    print domain.boundary_statistics(tags='ocean')

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

    r7225 r7327  
    2929##time_dirs = [time_dir1, time_dir2]
    3030
    31 time_dir1 = '20090615_160527_run_final_1.0_51436_jgriffin'
    32 time_dir2 = '20090615_160610_run_final_0.0_51436_jgriffin'
    33 time_dir3 = '20090615_160900_run_final_0.0_58242_jgriffin'
    34 time_dir4 = '20090615_160939_run_final_1.0_58242_jgriffin'
    35 time_dir5 = '20090615_161039_run_final_0.0_58284_jgriffin'
    36 time_dir6 = '20090615_161145_run_final_0.0_58349_jgriffin'
    37 time_dir7 = '20090615_161318_run_final_0.0_58187_jgriffin'
    38 time_dir8 = '20090615_161339_run_final_0.0_58113_jgriffin'
    39 time_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'
     31time_dir1 = '20090615_161420_run_final_0.0_58025_jgriffin'
    4432
    45 time_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]
     33
     34time_dirs = [time_dir1]
    4635
    4736for time_dir in time_dirs:
    … …  
    5443    #ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_conditions_ocean1"
    5544   
    56     print 'Process: Create File GDB'
    57     gp.CreateFileGDB_management(folder, "raster")
     45##    print 'Process: Create File GDB'
     46##    gp.CreateFileGDB_management(folder, "raster")
    5847
    5948    gp.Workspace = raster_gbd

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

    r7252 r7327  
    11"""
    22Generates ascii grids of nominated areas -
    3 Input: sww file from run_perth.py
     3Input: sww file from run_model.py
    44       boundaries for grids from project.py
    55Outputs: ascii grids of specified variables
    … …  
    2323
    2424#Specify output directories
    25 time_dir1 = '20090615_160610_run_final_0.0_51436_jgriffin'
    26 time_dir2 = '20090615_161318_run_final_0.0_58187_jgriffin'
     25time_dir1 = '20090615_161420_run_final_0.0_58025_jgriffin'
     26#time_dir2 = '20090615_161318_run_final_0.0_58187_jgriffin'
    2727
    28 time_dirs = [time_dir1, time_dir2]
     28time_dirs = [time_dir1]#, time_dir2]
    2929 
    3030
    … …  
    5252
    5353# one or more key strings from var_equations above
    54 #var = ['depth', 'speed','stage']
    55 var = ['depth', 'speed', 'stage']
     54var = ['elevation']
     55#var = ['depth', 'speed', 'stage']
    5656######
    5757# Start running the various conversions we require.

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

    r7252 r7327  
    2626#-------------------------------------------------------------------------------
    2727
    28 tide = 0.0   
     28tide = 1.0   
    2929             # difference between MSL and HAT in metres
    3030zone = 56               # specify zone of model
    31 #event_number = 58242    # Puysegur 1 in 10 000 # the event number or the mux file name
    32 #event_number = 51436    # New Hebrides 1 in 10 000
     31#event_number = 58242    # Puysegur 1 in 10 000 (Event 1)# the event number or the mux file name
     32event_number = 51436    # New Hebrides 1 in 10 000 (Event 2)
    3333#event_number = 58349    # Puysegur 1 in 5000
    34 #event_number = 58284    # Puysegur 1 in 2000
     34#event_number = 58284    # Puysegur 1 in 2000 (Event 3)
    3535#event_number = 58187    # Puysegur 1 in 1000
    3636#event_number = 58113    # Puysegur 1 in 500
    37 event_number = 58025    # Puysegur 1 in 200
     37#event_number = 58025    # Puysegur 1 in 200
    3838
    3939alpha = 0.1             # smoothing parameter for mesh
    … …  
    102102# Used in run_building_inundation.py
    103103# Format latitude,longitude etc (geographic)
    104 ##building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
     104building_exposure_filename = '' # from NEXIS
    105105
    106106# AREA OF IMAGES - Extent of each image to find out highest runup
    … …  
    137137#-------------------------------------------------------------------------------
    138138
    139 # ASCII export grid for Busselton
    140 xminBusselton = 340000
    141 xmaxBusselton = 352000
    142 yminBusselton = 6271500
    143 ymaxBusselton = 6280000
    144 
    145 # ASCII export grid for Bunbury
    146 xminBunbury = 369000
    147 xmaxBunbury = 381000
    148 yminBunbury = 6308000
    149 ymaxBunbury = 6316500
     139# ASCII export grid for Gosford
     140##xminGosford =
     141##xmaxGosford =
     142##yminGosford =
     143##ymaxGosford =
     144
     145
    150146
    151147################################################################################
    … …  
    242238# The absolute pathname for the building file
    243239# Used for run_building_inundation.py
    244 ##building_exposure = join(gauges_folder, building_exposure_filename)
     240building_exposure = join(gauges_folder, building_exposure_filename)
    245241
    246242# The absolute pathname for the image file

  • anuga_work/production/new_south_wales/gosford/run_model.py

    r7252 r7327  
    1 """Run a tsunami inundation scenario for Busselton, WA, Australia.
     1"""Run a tsunami inundation scenario for Gosford, NSW, Australia.
    22
    33The scenario is defined by a triangular mesh created from project.polygon, the
    … …  
    159159t0 = time.time()
    160160
     161##for t in domain.evolve(yieldstep=5,
     162##                       finaltime=project.finaltime,
     163##                       skip_initial_step=False):
     164##    print domain.timestepping_statistics()
     165##    print domain.boundary_statistics(tags='ocean')
     166
    161167for t in domain.evolve(yieldstep=project.yieldstep,
    162168                       finaltime=project.finaltime,
    … …  
    165171    print domain.boundary_statistics(tags='ocean')
    166172    print domain.volumetric_balance_statistics()
    167 
     173   
    168174print 'Simulation took %.2f seconds' % (time.time()-t0)
Note: See TracChangeset for help on using the changeset viewer.