Changeset 5781


Ignore:
Timestamp:
Sep 24, 2008, 11:56:50 AM (15 years ago)
Author:
kristy
Message:
 
Location:
anuga_work/production/perth
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/perth/build_boundary_27255.py

    r5761 r5781  
    44
    55Input: order_filename from project.py
     6       event_number needs to be reflected in the project file
    67Output: creates a sts file and csv files stored in project.boundaries_dir_event.
    78The run_perth.py is reliant on the output of this script.
  • anuga_work/production/perth/build_boundary_27283.py

    r5761 r5781  
    44
    55Input: order_filename from project.py
     6       event_number needs to be reflected in the project file
    67Output: creates a sts file and csv files stored in project.boundaries_dir_event.
    78The run_perth.py is reliant on the output of this script.
  • anuga_work/production/perth/build_perth.py

    r5761 r5781  
    130130if access(project.topographies_dir,F_OK) == 0:
    131131    mkdir (project.topographies_dir)
    132 G_clipped.export_points_file(project.combined_dir_name + '.pts')
    133 #G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC
     132#G_clipped.export_points_file(project.combined_dir_name + '.pts')
     133G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC
    134134
    135135import sys
  • anuga_work/production/perth/export_results_all.py

    r5761 r5781  
    2222#time_dir = '20080526_104946_run_final_0.6_test_kvanputt'
    2323#time_dir = '20080530_170833_run_final_0.6_exmouth_kvanputt'
    24 time_dir1 = '20080815_103442_run_final_0.0_polyline_alpha0.1_kvanputt'
    25 time_dir2 = '20080909_151438_run_final_0.0_polyline_alpha0.1_kvanputt'
     24#time_dir1 = '20080815_103442_run_final_0.0_polyline_alpha0.1_kvanputt'
     25time_dir1 = '20080912_154439_run_final_0.6_27255_alpha0.1_kvanputt'
    2626
    2727#cellsize = 20
     
    2929#timestep = 0
    3030#area = ['Geordie', 'Sorrento', 'Fremantle', 'Rockingham']
    31 area = ['All']
    32 #var = [1,2] # Absolute momentum and depth
     31area = 'All'
     32which_area = area
     33var = [1,2] # Absolute momentum and depth
    3334#var = [2] # depth
    34 var = [0,4] #stage and elevation
     35#var = [0,4] #stage and elevation
    3536
    36 time_dir = [time_dir1, time_dir2]
     37time_dirs = [time_dir1]
    3738for time_dir in time_dirs:
    3839
    3940    name1 = directory+time_dir+sep+project.scenario_name
    40     name2 = directory+time_dir+sep+'sww2'+sep+project.scenario_name+'_time_39900_0' #need to get assistance on how to make this into anything
     41    #name2 = directory+time_dir+sep+'sww2'+sep+project.scenario_name+'_time_39900_0' #need to get assistance on how to make this into anything
    4142
    42     names = [name1, name2]
     43    names = [name1]
    4344    for name in names:
    4445
    45         for which_area in area:
    46             if which_area == 'All':
    47            
    48             if which_area == 'Geordie':
    49                 easting_min = project.xminGeordie
    50                 easting_max = project.xmaxGeordie
    51                 northing_min = project.yminGeordie
    52                 northing_max = project.ymaxGeordie
    53 
    54             if which_area == 'Sorrento':
    55                 easting_min = project.xminSorrento
    56                 easting_max = project.xmaxSorrento
    57                 northing_min = project.yminSorrento
    58                 northing_max = project.ymaxSorrento
    59 
    60             if which_area == 'Fremantle':
    61                 easting_min = project.xminFremantle
    62                 easting_max = project.xmaxFremantle
    63                 northing_min = project.yminFremantle
    64                 northing_max = project.ymaxFremantle
    65 
    66             if which_area == 'Rockingham':
    67                 easting_min = project.xminRockingham
    68                 easting_max = project.xmaxRockingham
    69                 northing_min = project.yminRockingham
    70                 northing_max = project.ymaxRockingham
     46##        for which_area in area:
     47##            if which_area == 'All':
     48##           
     49##            if which_area == 'Geordie':
     50##                easting_min = project.xminGeordie
     51##                easting_max = project.xmaxGeordie
     52##                northing_min = project.yminGeordie
     53##                northing_max = project.ymaxGeordie
     54##
     55##            if which_area == 'Sorrento':
     56##                easting_min = project.xminSorrento
     57##                easting_max = project.xmaxSorrento
     58##                northing_min = project.yminSorrento
     59##                northing_max = project.ymaxSorrento
     60##
     61##            if which_area == 'Fremantle':
     62##                easting_min = project.xminFremantle
     63##                easting_max = project.xmaxFremantle
     64##                northing_min = project.yminFremantle
     65##                northing_max = project.ymaxFremantle
     66##
     67##            if which_area == 'Rockingham':
     68##                easting_min = project.xminRockingham
     69##                easting_max = project.xmaxRockingham
     70##                northing_min = project.yminRockingham
     71##                northing_max = project.ymaxRockingham
    7172
    7273             
     
    9394                    quantityname = 'elevation'  #Elevation
    9495
    95                 if which_area == 'All'
    96                     print 'start sww2dem',which_area
     96                else:
     97                    print 'start sww2dem',area
    9798                    sww2dem(name, basename_out = outname,
    9899                                quantity = quantityname,
     
    102103                                verbose = True,
    103104                                format = 'asc')
    104 
    105                 else:
    106                     print 'start sww2dem',which_area, easting_min
    107                     sww2dem(name, basename_out = outname,
    108                                 quantity = quantityname,
    109                                 #timestep = timestep,
    110                                 cellsize = cellsize,     
    111                                 easting_min = easting_min,
    112                                 easting_max = easting_max,
    113                                 northing_min = northing_min,
    114                                 northing_max = northing_max,       
    115                                 reduction = max,
    116                                 verbose = True,
    117                                 format = 'asc')
    118 
     105##
     106##                else:
     107##                    print 'start sww2dem',which_area, easting_min
     108##                    sww2dem(name, basename_out = outname,
     109##                                quantity = quantityname,
     110##                                #timestep = timestep,
     111##                                cellsize = cellsize,     
     112##                                easting_min = easting_min,
     113##                                easting_max = easting_max,
     114##                                northing_min = northing_min,
     115##                                northing_max = northing_max,       
     116##                                reduction = max,
     117##                                verbose = True,
     118##                                format = 'asc')
     119##
  • anuga_work/production/perth/project.py

    r5761 r5781  
    4141scenario = 'perth_tsunami_scenario'
    4242##One or all can be changed each time the run_scenario script is excuted
    43 tide = 0.6
    44 event_number = 27255
    45 #event_number = 27283
     43tide = 0 #0.6
     44#event_number = 27255
     45event_number = 27283
    4646alpha = 0.1
    4747friction=0.01
  • anuga_work/production/perth/run_perth.py

    r5761 r5781  
    130130    # Setup initial conditions
    131131    #-------------------------------------------------------------------------
    132    print 'Setup initial conditions'
     132    print 'Setup initial conditions'
    133133
    134134    #following sets the stage/water to be offcoast only
Note: See TracChangeset for help on using the changeset viewer.